数值计算方法实验报告--矩阵直接LU分解法方程组和微分方程的数值解法上机实现.docx

《数值计算方法实验报告--矩阵直接LU分解法方程组和微分方程的数值解法上机实现.docx》由会员分享，可在线阅读，更多相关《数值计算方法实验报告--矩阵直接LU分解法方程组和微分方程的数值解法上机实现.docx（3页珍藏版）》请在第一文库网上搜索。

1、数值计算方法实验报告实验报告题目：矩阵直接1U分解法方程组和微分方程的数值解法上机实现实验要求矩阵直接1U分解法方程组和微分方程的数值解法二、 实验分析（包括数学原理，小组分析讨论后确定实验方案和实现思路）用matIab对矩阵三角分解法求解方程组三、 实脸步骤(过程)(包括程序及上机的实现的结果)functionX=1U1(a,b)b=b,;a=a,b,n=1ength(b);x=zeros(n,1);y=zeros(n,1);U=zeros(n);1=eye(n);fork=1:nU(1,k)=a(1,k)；1(k,1)=a(k,1)U(1,1);endfori=2:nfork=i:n1u=

2、0;1u1=0;forj=1:i-11u=1u+1(i,j)*U(j,k);1u1=1u1+1(k,j)*U(j.i);endU(i,k)=a(i,k)-1u;1(k,i)=(a(k,i)-Iu1)/U(i,i);endendfori=1:n1y=0;forj=1:i1y=1y+1(i,j)*y(j);endy(i)=B(i)-1y;endfori=n:-1:11y1=0;forj=i+1:n1y1=1y1+U(izj)*X(j);endX(i)=(y(i)-1y1)U(i,i);end输入实验内容:1-723-32(57输入：A=10,-7,0;-3,2,6;5,-1,5;B=7,4,6;X=1U1(A,B)A=10-707-32645-1561=1.000000-0.30001.00000.5000-25.00001.0000U=10.0000-7.000000-0.10006.000000155.0000X=-0.0000-1.00001.0000四、总结（包括实验过程遇到的情况等，组长总结组员在整个过程的参与情况）在学习中，通过小组讨论探讨以及网上查找资料等，学习了用1U分解法解方程组