任意N阶复数矩阵求逆任意N阶复数矩阵求逆
#include
#include
#include"stdlib.h"
#include "stdio.h"
using namespace std;
template
void reciprocal_complex_matrix(double c0_real[n][n],double c0_imag[n][n],double
cr_real[n][n],double cr_imag[n][n])
{
int i,j,k,m;
double t;
double M...
任意N阶复数矩阵求逆
#include
#include
#include"stdlib.h"
#include "stdio.h"
using namespace std;
template
void reciprocal_complex_matrix(double c0_real[n][n],double c0_imag[n][n],double
cr_real[n][n],double cr_imag[n][n])
{
int i,j,k,m;
double t;
double M[3*n][3*n]={0};
for(i=0;ifabs(M[j][j]))
m=i;
else
m=j;
}
for(k=0;k<3*n;k++)
{
t=M[j][k];
M[j][k]=M[m][k];
M[m][k]=t;
}
for(i=j+1;i<2*n;i++)
for(k=3*n-1;k>=j;k--)
M[i][k]=M[i][k]-M[j][k]*M[i][j]/M[j][j];
}
for(j=2*n-1;j>0;j--)
for(i=0;i=0;k--)
M[i][k]=M[i][k]-M[j][k]*M[i][j]/M[j][j];
for(i=0;i<2*n;i++)
for(k=3*n-1;k>=0;k--)
M[i][k]=M[i][k]/M[i][i];
for(i=0;i(c0_real,c0_imag,cr_real,cr_imag);
for(int i=0;i
本文档为【任意N阶复数矩阵求逆】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑,
图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。