基于DSP的FCM图像分割

int  m=50,k,ran[3];
double center[3],u[3][84*64],sum=0,temp=0;//中心，隶属度 

for(i=0;i<3;i++)//初始化聚类中心 
	{
		ran[i]=rand()%(84*64);
		center[i]=image[ran[i]]; 
	} 
	
	while(m--)//迭代 
	{
	
		for(i=0;i<3;i++)//计算隶属度 
		{
			for(j=0;j<84*64;j++)
			{
				sum=0;
				for(k=0;k<3;k++)
				{
					sum=sum+pow(fabs(image[j]-center[i])/fabs(image[j]-center[k]),2);
				}
				u[i][j]=1/sum;
			}
		} 
				
		for(i=0;i<3;i++)//更新中心点 
		{
			sum=0;
			temp=0;
			for(j=0;j<84*64;j++)
			{
				sum=sum+pow(u[i][j],2)*image[j];
				temp=temp+pow(u[i][j],2);
			}
			center[i]=sum/temp;
		}
	}		
	for(j=0;j<N;j++)//分类结果 
	{
		if(u[0][j]<u[1][j])
		{
			if(u[1][j]<u[2][j])
				image[j]=100;
			else
				image[j]=225;
		}
		else
		{
			if(u[0][j]<u[2][j])
				image[j]=100;
			else		
			    image[j]=0;
	    }
	}

}