数值分析实验王哲岩.doc

《数值分析实验王哲岩.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、《数值分析》实验报告姓名王哲岩学号1405H0012院（系）数计专业一计科年级14级指导教师李雅琴2017年3月12日1.设计说明书（目的要求、实验内容、步骤）【要求】利用matlab通过编程的方法，使用高斯消元法计算方程组的解。【内容】计算A二［64357;537812;34569;129101211;98875］B二［55,77,58,129,94］的解【步骤】首先新建脚本，将设计好的函数输入进去，在命令行输入A和b的值,调用函数计算2.程序流程图（NS盒图）消元过程K=i+1J=i+1A(kJ)=A(k,j)+A(ij)*(-A(k,i)

2、/A(i,i))b(k)=b(k)+b(i)*(-A(k,i)/A(U))A(k,i)=O直到k=n^直到i二nJ回代过程i=n-li=i-lSum=0J=i+1sum=sum+A(i,j)*x(j)直到j二nx(i)=(b(i)-sum)/A(i,i)直到日1.源程序functionx=gauss(A,b)n=length(b);x=zeros(n,i);fori=i:n-ifork=i+i:nforj=i+i:nA(k,j)=A(k,j)+A(i,j)*(-A(k,i)/A(i,i))；endb(k)=b(k)+b(i)*(-A(k,i)/

3、A(i,i));A(k,i)二o;endenddisp(A)disp(b)pausex(n)=b(n)/A(n,n);fori=n-i:-i:isum=o；forj=i+i:nsum=sum+A(i,j)*x(j);endx(i)=(b(i)-sum)/A(i,i);end2.运行结果A二[64357;537812;34569;129101211;98875]1212345643575378123345694129101211598875B二[55,77,5&129,94]1234555775812994运算结果>>x=gauss(A,b)6.

4、00004.00003.00005.00007.00000-0.33334.50003.83336.16670030.500026.500042.5000000-1.7049-8.885200009.846255.000031.1667217.5000-12.29519.84623.00002.00004.00002.00001.00005.总结与感受通过本次实验,掌握了利用matlab工具实现高斯消元法的过程。通过绘制ns盒图又详细的了解了高斯消元法在实现下三角的时候的具体流程。