gauss.cpp

高斯消去法解线性方程组,适用于工程计算程序里的部分数值计算

#include "iostream.h"
#include"stdlib.h"
#include "math.h"
#include "stdio.h"
#define L 4
void form_factor_table(int n,float M[L][L],float FT[L][L])//求系数矩阵的因子表子函数
{
	int i,j,k;
    //对M做n次Gauss消去
	for(j=0;j<n;j++)//将M的对角元化为1，然后下三角部分化为0
	{
		FT[j][j]=M[j][j];//存归格化因子到因子表的对角上
		for(k=n-1;k>=j;k--)
			M[j][k]=M[j][k]/M[j][j];//对角元化为1
		for(i=j+1;i<n;i++)//消去第j列
		{
            FT[i][j]=M[i][j];//存消去因子到因子表的下三角部分
			for(k=n-1;k>=j;k--)
				M[i][k]=M[i][k]-M[i][j]*M[j][k];//第j列第i个元素化为0
		}
	}
	for(i=0;i<n;i++)
		for(j=i+1;j<n;j++)
			FT[i][j]=M[i][j];//M的上三角部分存到因子表的上三角部分
}
void main()
{
	
	int i,j,n;
	float M[L][L]={{1,2,1,1},{2,1,0,0},{1,0,1,0},{1,0,0,1}};//系数矩阵
	float FT[L][L]={0};//因子表矩阵
	n=4;
    form_factor_table(n,M,FT);
	for(i=0;i<n;i++)
	{
		for(j=0;j<n;j++)
			printf("%f\t",FT[i][j]);
		printf("\n");
	}
	 system("pause"); 
}