《2022年实验_求解线性方程组直接法参考 .pdf》由会员分享，可在线阅读，更多相关《2022年实验_求解线性方程组直接法参考 .pdf（8页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、1 数值分析实验报告二求解线性方程组的直接方法（2 学时）班级专业姓名学号日期一 实验目的1掌握求解线性方程组的高斯消元法及列主元素法；2. 掌握求解线性方程组的克劳特法；3. 掌握求解线性方程组的平方根法。二 实验内容1用高斯消元法求解方程组（精度要求为610） ：1231231233272212240 xxxxxxxxx2用克劳特法求解上述方程组（精度要求为610） 。3. 用平方根法求解上述方程组（精度要求为610） 。4. 用列主元素法求解方程组（精度要求为610） ：1231231233432222325xxxxxxxxx三 实验步骤（算法）与结果1.#include Stdio.h
2、 #include math.h 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 8 页 - - - - - - - - - 2 void 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; a31=2,a32=-2,a33=4; b1=7,b2=-1,b3=0; for(m=1;m=3;m+
3、) 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; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师
4、精心整理 - - - - - - - 第 3 页，共 8 页 - - - - - - - - - 4 for(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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 8 页 - - - - - - - - - 5 #include Stdio.h #i
5、nclude math.h void 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=0; for(m=1;m=3;m+) for(n=1;n=n) lmn=amn; for(i=1;i=n-1;i+) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - -
6、 - - - 第 5 页，共 8 页 - - - - - - - - - 6 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+) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页
7、，共 8 页 - - - - - - - - - 7 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(); 运行结果：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 8 页 - - - - - - - - - 8 3 4 四 实验收获与教师评语名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 8 页 - - - - - - - - -