迭代法解线性方程组MATLAB编程实验报告

迭代法解线性方程组MATLAB编程实验报告一、问题描述：编写函数来完成迭代法解线性方程组；总结matlab解非线性方程的命令，用实例来展示。二、实验步骤（过程）：（1）迭代法解线性方程组——以书本例6.1为例clc,clear;A=[1031;2-103;1310];b=[14-514]';[m,n]=size(A);k=4;U=zeros(n,n);L=zeros(n,n);D=zeros(n,n);fori=1:nD(i,i)=A(i,i);endDT=A-D;fori=1:n-1forj=2:mU(i,j)=-T(i,j);endendUL=-T-UP=D-L;P=inv(P);GS=P*Ub1=P*bD=inv(D);B=D*(L+U);b2=D*b;Z=[B,b2]G=[GS,b1]W=zeros(k+1,n);Q=zeros(k+1,n);fori=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(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,:)（二）总结matlab解非线性方程的命令（1）利用solve()例：求方程17433xxx（2）利用fsolve()求解方程组在212221*21xxxx在x1=x2=1附近的近似解。三、结论：