求解线性方程组直接法.doc

《求解线性方程组直接法.doc》由会员分享，可在线阅读，更多相关《求解线性方程组直接法.doc（9页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、数值分析实验报告二求解线性方程组的直接方法 一 实验目的1掌握求解线性方程组的高斯消元法及列主元素法；2. 掌握求解线性方程组的克劳特法；3. 掌握求解线性方程组的平方根法。二 实验内容1用高斯消元法求解方程组（精度要求为）：2. 用列主元素法求解方程组（精度要求为）：三 实验步骤（算法）与结果1. #include Stdio.h#include math.hvoid main() float l44,u44; int a44,b4; float x5,z4; int n,m,i; l11=1,l22=1,l33=1; a11=3,a12=-1,a13=2,a21=-1,a22=2,a23=

2、-2; a31=2,a32=-2,a33=4; b1=7,b2=-1,b3=0; for(m=1;m=3;m+) for(n=1;n=3;n+) if(m=n) umn=amn; for(i=1;in) lmn=amn; for(i=1;i=n-1;i+) lmn-=lmi*uin; lmn/=unn; for(m=1;m=3;m+) for(n=1;n=3&n=m;n+) printf(l%d%d=%ft,m,n,lmn); printf(u%d%d=%ft,n,m,unm); for(m=1;m=3;m+) zm=bm; for(i=1;i=m-1;i+) zm-=lmi*zi; for(

3、m=1;m0;m-) xm=zm; for(i=m+1;i0&i=3;i+) xm-=umi*xi; xm=xm/umm; for(m=1;m=3;m+) printf(x%d=%ft,m,xm); getch();运行结果：2#include Stdio.h#include math.hvoid main() float l44,u44; int a44,b4; float x5,z4; int n,m,i; u11=1,u22=1,u33=1; a11=3,a12=-1,a13=2,a21=-1,a22=2,a23=-2; a31=2,a32=-2,a33=4; b1=7,b2=-1,b3

4、=0; for(m=1;m=3;m+) for(n=1;n=n) lmn=amn; for(i=1;i=n-1;i+) lmn-=lmi*uin; if(mn) umn=amn; for(i=1;i=m-1;i+) umn-=lmi*uin; umn/=lmm; for(m=1;m=3;m+) for(n=1;n=3&n=m;n+) printf(l%d%d=%ft,m,n,lmn); printf(u%d%d=%ft,n,m,unm); for(m=1;m=3;m+) zm=bm; for(i=1;i=m-1;i+) zm-=lmi*zi; zm/=lmm; for(m=1;m0;m-) xm=zm; for(i=m+1;i0&i=3;i+) xm-=umi*xi; for(m=1;m=3;m+) printf(x%d=%ft,m,xm); getch();运行结果：四 实验收获 通过做实验，增加了对程序的了解程度，对方程组的解法有了更深的了解，能够更熟练的应用LU分解法解决实际问题了。