matrix.h

最优化作业c++源代码?钣呕饕礳++源代码最优化作业c++源代码

#include<math.h>
int const n=2;
int const m=3;
/*两个矩阵相乘*/
void matrix_mul(double (*x)[n],double (*y)[n],double (*temp)[n])
{
	int i,j,k;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
	        temp[i][j]=0;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			for(k=0;k<n;k++)
				temp[i][j]+=x[i][k]*y[k][j];
}
void matrix_mul(double *z,double(*x)[n],double *temp1)
{
	int i,j;
	for(i=0;i<n;i++)
		temp1[i]=0;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			temp1[i]+=x[i][j]*z[j];
}
void matrix_mul(double(*x)[n],double *z,double *temp1)
{
	int i,j;
    for(i=0;i<n;i++)
		temp1[i]=0;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			temp1[i]+=x[i][j]*z[j];
}
void matrix_mul(double (*x)[m],double (*y)[n],double (*temp)[n])
{
	int i,j,k;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
	        temp[i][j]=0;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			for(k=0;k<m;k++)
				temp[i][j]+=x[i][k]*y[k][j];
}
void matrix_mul(double (*x)[n],double (*y)[m],double (*temp)[n])
{
	int i,j,k;
	for(i=0;i<m;i++)
		for(j=0;j<m;j++)
	        temp[i][j]=0;
	for(i=0;i<m;i++)
		for(j=0;j<m;j++)
			for(k=0;k<n;k++)
				temp[i][j]+=y[i][k]*x[k][j];
}
//用Gass_Jordan法求矩阵的逆矩阵
void matrix_inv(double a[n][n],double c[n][n])
{
	double m[n];//辅助乘数
	double T;//存储换行临时变量
	double temp=0;
	double t;//最大列主元
	int tap1=0,tap2=0;//最大列主元下标
	double a1[n][n];
    for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			a1[i][j]=a[i][j];
			if(i==j)
			{
				c[i][j]=1;
			}
			else
			{
				c[i][j]=0;
			}
		}
	}
	for(int s=0;s<n;s++)
	{
		t=a[s][s];//赋初值
		for(int p=s;p<n;p++)//选取最大列主元
		{
			if(fabs(a1[p][s])>t)
			{
				t=a1[p][s];
				tap1=p;
				tap2=s;
			}
		}
		if(t==0)
		{cout<<"列主元为零，失败！"<<endl;break;}
		if(tap1!=tap2)//换行
		{
			for(int i=0;i<n;i++)
			{
				T=a1[tap1][i];
				a1[tap1][i]=a1[s][i];
				a1[s][i]=T;
				T=c[tap1][i];//逆矩阵换行
				c[tap1][i]=c[s][i];
				c[s][i]=T;
			}
		}
		tap1=tap2=0;//置零
		for(int j=0;j<n;j++)//单位化,(注意：逆矩阵每列都单位化，而原矩阵则不受影响)
		{
			a1[s][j]=a1[s][j]/t;
			c[s][j]=c[s][j]/t;
		}		
		for(int i=0;i<n;i++)//消元
		{
			if(i!=s)
			{
				m[i]=a1[i][s]/a1[s][s];
			    for(int j=0;j<n;j++)//注意：逆矩阵每列都消元，而原矩阵则不受影响
				{
					a1[i][j]=a1[i][j]-m[i]*a1[s][j];
					c[i][j]=c[i][j]-m[i]*c[s][j];
				}
			}
		}
	}
}
//两个矩阵相加
void matrix_plus(double (*x)[n],double (*y)[n],double (*temp)[n])
{
	int i,j;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			temp[i][j]=x[i][j]+y[i][j];
}
//两个矩阵相减
void matrix_minus(double (*x)[n],double (*y)[n],double (*temp)[n])
{
	int i,j;
	for(i=0;i<n;i++)
		for(j=0;j<n;j++)
			temp[i][j]=x[i][j]-y[i][j];
}