数值分析实验报告.doc

数值分析数值分析试验报告东北高校20__年5月

一、试验要求)建立矩阵求逆的算法公式2)编制下三角部分消元，上三角部分消元和对角元单位化等子程序3)应用结构程序设计编程出求n阶矩阵的逆矩阵的通用程序

二、试验目的与意义)通过该课题的试验，把握求逆矩阵的程序设计方法2)学会利用矩阵求逆运算求解线性方程组的方法，提高数值方法的运用和编程的力量

三、试验环境编译环境：VS20__编译语言：C+

四、试验过程与分析XXX;classJuZhenprivate:;voidJuZhen:setSizevoidJuZhen:showintn;cout请?输?入?矩?阵的?阶:;cinn;size=n;doubledata020;intsize;voidqiuNi;voidsetSize;voidshow;voidchuShi;public:数值分析voidJuZhen:chuShiinti,j;voidJuZhen:qiuNiinti,j,k;intma_I=0;for(i=;isize;i+)if(ma_I!=0)doubletemp;for(j=0;j2_size;j+)if(fabs(datama_I0)fabs(datai0)ma_I=i;for(i=0;isize;i+)cout请?输?入?i+行D的?元a素?:;for(j=0;j2_size;j+)if(jsize)cindataij;dataij=

1.0;dataij=0.0;elseif(j=i+size)elsefor(i=0;isize;i+)cout_endl;coutendl;for(j=size;j2_size;j+)coutendl;XXX

(0);XXX(ios:right);coutdataij;数值分析voidmainwhile

(1)JuZhena;XXX;XXX;XXXNi;doubletemp2;for(i=0;isize;i+)cout逆?矩?阵为a：阰endl;show;if(dataii!=0)for(j=0;j2_size;j+)dataij_=temp2;for(j=0;jsize;j+)if(j!=i)doubletemp3=dataji;for(k=0;k2_size;k+)datajk-=temp3_dataik;cout此?矩?阵无T逆?!endl;return;temp2=

XXX;elsetemp=data0j;data0j=datama_Ij;datama_Ij=temp;

五、试验结果与分析数值分析如图，输入矩阵的阶，然后依次输入矩阵中元素，即可得到逆矩阵，若该矩阵无逆，则输出“此矩阵无逆”。

六、附录参考资料：)张X，闫XX数值分析北京，冶金工业出版社，20__