求逆矩阵.h

二维点集合聚类的图像化演示程序 用它可以很直观的看到聚类的整个过程

#include <iostream>
#include <cassert>
#include <iomanip>
#include <new>
using namespace std;
void f_matrix(float **matrix,float **identity_matrix,int lenth)
{
	int i,j;
	float temp;
//原单位矩阵转上三角矩阵
	for(i=0;i<lenth-1;i++)
	{
		int tempi;
		for(tempi=i;tempi<lenth;tempi++)
		{
			if(matrix[tempi][i]!=0)
				temp=matrix[tempi][i];
			else
				temp=1;
			for(j=0;j<lenth;j++)
			{
				matrix[tempi][j]/=temp;
				identity_matrix[tempi][j]/=temp;
			}
		}
		//首元化一结果
        //	/*
		cout<<"首元化后原矩阵为："<<endl;
		for(tempi=0;tempi<lenth;tempi++)
		{
			for(int tempj=0;tempj<lenth;tempj++)
			{
				cout<<setw(11)<<matrix[tempi][tempj];
			}
			cout<<endl;
		}
		cout<<endl;
		cout<<"首元化后单位矩阵为："<<endl;
		for(tempi=0;tempi<lenth;tempi++)
		{
			for(int tempj=0;tempj<lenth;tempj++)
			{
				cout<<setw(11)<<identity_matrix[tempi][tempj];
			}
			cout<<endl;
		}
		cout<<endl;
        //	*/
		//化零
        cout<<i<<endl;
		for(tempi=i+1;tempi<lenth;tempi++)
		{
			
			for(j=0;j<lenth;j++)
			{
				//cout<<"化零"<<endl;
				matrix[tempi][j]-=matrix[i][j];
				identity_matrix[tempi][j]-=identity_matrix[i][j];
			}
		}
	}
//原单位矩阵转单位阵
    //首元化一
	temp=matrix[lenth-1][lenth-1];
	//cout<<temp;
	for(j=0;j<lenth;j++)
	{
		matrix[lenth-1][j]/=temp;
		identity_matrix[lenth-1][j]/=temp;
	}
    //首元化一结果
    //	/*
	cout<<"首元化后原矩阵为："<<endl;
	for(i=0;i<lenth;i++)
	{
		for(j=0;j<lenth;j++)
		{
			cout<<setw(8)<<matrix[i][j];
		}
		cout<<endl;
	}
	cout<<endl;
	cout<<"首元化后单位矩阵为："<<endl;
	for(i=0;i<lenth;i++)
	{
		for(j=0;j<lenth;j++)
		{
			cout<<setw(8)<<identity_matrix[i][j];
		}
		cout<<endl;
	}
	cout<<endl;
	// */
	
	for(i=lenth-1;i>=0;i--)
	{
		for(int tempi=i-1;tempi>=0;tempi--)
		{
			temp=matrix[tempi][i];
			for(j=0;j<lenth;j++)
			{
				matrix[tempi][j]-=matrix[i][j]*temp;
				identity_matrix[tempi][j]-=identity_matrix[i][j]*temp;
			}
		}
	}
	//最后结果
    //	/*
	cout<<"最后结果："<<endl;
	cout<<"原矩阵为："<<endl;
	for(i=0;i<lenth;i++)
	{
		for(j=0;j<lenth;j++)
		{
			cout<<setw(11)<<matrix[i][j];
		}
		cout<<endl;
	}
	cout<<endl;
	cout<<"单位矩阵为："<<endl;
	for(i=0;i<lenth;i++)
	{
		for(j=0;j<lenth;j++)
		{
			cout<<setw(11)<<identity_matrix[i][j];
		}
		cout<<endl;
	}
	cout<<endl;
	// */
}