《迭代法解线性方程组MATLAB编程实验报告》由会员分享,可在线阅读,更多相关《迭代法解线性方程组MATLAB编程实验报告(3页珍藏版)》请在人人文库网上搜索。
1、迭代法解线性方程组MATLAB编程实验报告一、 问题描述:编写函数来完成迭代法解线性方程组;总结matlab解非线性方程的命令,用实例来展示。二、 实验步骤(过程):(1)迭代法解线性方程组以书本例6.1为例clc,clear;A=10 3 1;2 -10 3;1 3 10;b=14 -5 14;m,n=size(A);k=4;U=zeros(n,n);L=zeros(n,n);D=zeros(n,n);for i=1:nD(i,i)=A(i,i);endDT=A-D;for i=1:n-1for j=2:mU(i,j)=-T(i,j);endendUL=-T-U P=D-L;P=inv(P)。
2、;GS=P*Ub1=P*bD=inv(D);B=D*(L+U);b2=D*b;Z=B,b2G=GS,b1W=zeros(k+1,n);Q=zeros(k+1,n); for i=2:k+1j=1;W(i,j)=Z(1,j)*W(i-1,j)+Z(1,j+1)*W(i-1,j+1)+Z(1,j+2)*W(i-1,j+2)+Z(1,j+3);W(i,j+1)=Z(2,j)*W(i-1,j)+Z(2,j+1)*W(i-1,j+1)+Z(2,j+2)*W(i-1,j+2)+Z(2,j+3);W(i,j+2)=Z(3,j)*W(i-1,j)+Z(3,j+1)*W(i-1,j+1)+Z(3,j+2)*W(。
3、i-1,j+2)+Z(3,j+3);P(i,j)=G(1,j)*P(i-1,j)+G(1,j+1)*P(i-1,j+1)+G(1,j+2)*P(i-1,j+2)+G(1,j+3);P(i,j+1)=G(2,j)*P(i-1,j)+G(2,j+1)*P(i-1,j+1)+G(2,j+2)*P(i-1,j+2)+G(2,j+3);P(i,j+2)=G(3,j)*P(i-1,j)+G(3,j+1)*P(i-1,j+1)+G(3,j+2)*P(i-1,j+2)+G(3,j+3);endJ=W(k+1,:)s=P(k+1,:)(2) 总结matlab解非线性方程的命令(1) 利用solve()例:求方程(2) 利用fsolve()求解方程组在在x1=x2=1附近的近似解。三、 结论。