C常用矩阵子函数

相关主题

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

double scalar(double MA[R1][C1],double k)

{

int i,j,R1,C1;

double k,MA2[][];

for (i=0;i

{

for (j=0;j

MA2[i][j]=k*MA[i][j];

}

return MA2[][];

}

double Add(double MA1[R1][C1],double MA2[R2][C2]) {

inti,j;

double MA[3][3];

for(i=0;i

{

for(j=0;j

MA[i][j]=MA1[i][j]+MA2[i][j];

}

return MA[3][3];

}

double Sub(double MA1[R1][C1],double MA2[R2][C2]) {

inti,j;

double MA[3][3];

for(i=0;i

{

for(j=0;j

MA[i][j]=MA1[i][j]-MA2[i][j];

}

return MA[3][3];

}

double Mul(double MA1[R1][C1],double MA2[R2][C2]) {

inti,j,k;

double MA[3][3];

for(i=0;i

{

for(j=0;j

{

for(k=0;k

MA[i][j]+=MA1[i][k]*MA2[k][j];

}

return MA[3][3];

}

double Tra(double MA1[R1][C1])

{

inti,j;

double MA[R1][C1];

for(i=0;i

{

for(j=0;j

MA[j][i]=MA1[i][j];

}

return MA[R1][C1];

}

double Inv(a[M][M])

{

double a[M][2*M];

double b[N][2*N];

floatt,x;

inti,j,k,T=M;

for(i=0;i

for(j=0;j<(2*T);j++)

{ if (j

a[i][j]=b[i][j];

else if (j==T+i)

a[i][j]=1.0;

else

a[i][j]=0.0;

}

for(i=0;i

{

for(k=0;k

{

if(k!=i)

{

t=a[k][i]/a[i][i];

for(j=0;j<(2*T);j++)

{

x=a[i][j]*t;

a[k][j]=a[k][j]-x;

}

for(i=0;i

{

t=a[i][i];

for(j=0;j<(2*T);j++)

a[i][j]=a[i][j]/t;

}

for(i=0;i

y=y*a[i][i];

if(y==0)

printf("Error\n");

else

{

for(i=0;i

for(j=0;j

b[i][j]=a[i][j+T];

};

return b[i][j];

}