matrixmul.cpp

FCM算法的实现及测试

#include "data.h"
double **MatrixMul(double **A,int row1,int col1,double **B,int row2,int col2)	
{		
	//矩阵相乘（矩阵1，行数，列数，矩阵2，行数，列数），
	//函数内已开僻内存，直接返回给定义指针
	if(row2!=col1)
	{
		cout<<"error!两矩阵不满足相乘条件！"<<endl;
		return NULL;
	}

	int i,j,k;
	double **p;
	p=new double*[row1];
	for(i=0;i<row1;i++)
	{
		p[i]=new double[col2];
	}

	for(i=0;i<row1;i++)
	{
		for(j=0;j<col2;j++)
		{
			p[i][j]=0.0;
			for(k=0;k<col1;k++)
			{
				p[i][j]+=A[i][k]*B[k][j];
			}		
		}
	}
	
	return p;
	
	//内存释放
	for(i=0;i<row1;i++)
	{
		delete[]p[i];
	}
	delete[]p;
}