差分方法实验报告解读

1、实验报告课程名称：计算方法院 系：数学科学系专业班级： 数应1001学 号：1031110139学生姓名： 姚海保指导教师：沈林开课时间：2012至2013学年第一学期一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后， 每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告， 不得 抄袭，不得缺交。学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。 字迹工 整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求1. 实验报告批改要深入细致， 批改过程中要发现和纠正学生实验报告中的问 题，给出评语和实验报告

2、成绩，签名并注明批改日期。实验报告批改完成后，应 采用适当的形式将学生实验报告中存在的问题及时反馈给学生。2. 实验报告成绩用百分制评定， 并给出成绩评定的依据或评分标准 （附于实 验报告成绩登记表后） 。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验 报告的学生应及时批评教育， 并对该次实验报告的分数以零分处理。 对单独设课 的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之 一以上，不得同意其参加本课程的考核。3. 各实验项目的实验报告成绩登记在实验报告成绩登记表中。 本学期实验项 目全部完成后，给定实验报告综合成绩。4. 实验报告综合成绩应按课程教学大纲规定比例（一

3、般为 10-15%）计入实 验课总评成绩； 实验总评成绩原则上应包括考勤、 实验报告、 考核（操作、理论） 等多方面成绩；5. 实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订， 按如下顺序装订成册： 实验报告封面、 实验报告成绩登记表、 实验报告成绩评定 依据、实验报告（按教学进度表规定的实验项目顺序排序） 。装订时统一靠左侧 按“两钉三等分”原则装订。实验名称MATLAB基本操作实验时间2012年9月13日学生姓名姚海保实验地点9#405数学实验室1、实验所用软件xp操作系统、Matlab实验软件2、实验目的熟悉Matlab编程环境，理解和掌握 Matlab执行命令的方式，会

4、做一些简单函数 的图像。通过用Matlab编程解决数值分析问题，更深一步的体会计算方法这门课的重 要性。3、实验内容(一) 、Matlab操作界面1. 命令窗口( comma nd wi ndoW2. 命令历史窗口( comma nd history)3. 工作空间管理窗口( workspace)4. 当前路径窗口( current directory)(二) 、具体练习1471、简单矩阵A= 2 5 8的输入步骤。.3 69一tt2、画出衰减振荡曲线y = e3si n3t及其它的包络线y=e3。t的取值范围是0,4兀。( / 2 + 2 )3、画出z-. 所表示的二维曲面。x, y的取值范

5、围是-8,8。7x + y4、复数矩阵的生成及运算.4、实验方法、步骤1. 了解matlab的硬件和软件必备环境；2. 启动 matlab ;3. 熟悉标题栏，菜单栏，工具栏，元素选择窗口，状态栏，控制栏以及系统布局区;4. 学习一些简单函数图形的绘制及命令的编写5、实验数据记录与分析1471、简单矩阵A=258的输入步骤。369 一(1) 在键盘上输入下列内容A = 1,4,7; 2,5,8; 3,6,9(2) 按【Enter】键，指令被执行。(3) 在指令执行后，MATLA指令窗中将显示以下结果：A = 1472 583 69tt2、画出衰减振荡曲线y二e刁sin3t及其它的包络线y。二。

6、t的取值范围是0,4：t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*si n(3*t);plot(t,y,-r,t,y0,:b,t,-y0,:b)grid10.80.60.40.20-0.2-0.4-0.6y i-0.8-1012r123、画出z/M J：十马所表示的三维曲面。x,y的取值范围是-8,8Jx2+y2clear;x=-8:0.5:8;y=x;X=on es(size(y)*x;Y=y* on es(size(x);R=sqrt(X.A2+Y.A2)+eps;%Z=si n(R)./R;%surf(X,Y,Z);%colormap(cool)%x

7、label(x),ylabel(y),zlabel(z)-10 -104、复数矩阵的生成及运算A=1,3;2,4-5,8;6,9*iB=1+5i,2+6i;3+8*i,4+9*iC=A*BA = 1.0000 - 5.0000i 3.0000 - 8.0000i2.0000 - 6.0000i 4.0000 - 9.0000iB =1.0000 + 5.0000i 2.0000 + 6.0000i3.0000 + 8.0000i 4.0000 + 9.0000iC =1.0e+002 *0.99001.1600 - 0.0900i1.1600 + 0.0900i 1.37006、实验结论通过本

8、次实验，我熟悉了 Matlab编程,掌握了基本的matlab编程命令，并熟悉 了一些作图方法，让我对这款软件充满了兴趣。指导教师评语和成绩评定指导教师签字：年 月 日实验名称插值算法实验时间2012年9月27日学生姓名姚海保实验地点9#405数学实验室1实验所用软件WINxp操作系统、Matlab2、实验目的设计一个程序，通过键盘指令输入参数，可以实现 Taylor插值方法、线性插值方 法、抛物插值方法以及一般情形的Lagrange插值方法的程序。使学生掌握插值算法的 基本理论，并提高学生在程序设计中发现问题、分析问题和解决问题的能力。3、实验内容1) .已知函数f(x)在1,3上具有二阶连续

9、导数，f(x)兰5，且满足条件 f (1) -1, f (3) - 2.求线性插值多项式和函数值f (1.5)，并估计其误差.2) .求函数f(x)=e*在0,1上线性插值多项式，并估计其误差3) .给出节 点数据 f(-2.00) =17.00，f(0.00)=1.00 ，f(1.00)=2.00，f(2.00)-17.00，作三次拉格朗日插值多项式计算f(0.6)，并估计其误差.4、实验方法、步骤1).解输入程序X=1,3;Y=1,2;101=poly(X(2)/(X(1)-X(2),111 =poly(X(1)/( X(2)- X(1), l0=poly2sym (101),1仁poly

10、2sym (111), P=l01* Y(1)+ 111* Y(2),L=poly2sym (P),x=1.5; Y = polyval(P,x)运行后输出基函数lo和l 1及其插值多项式的系数向量 P (略八 插值多项式L和插值丫 为l0 =l1 = L = Y =-1/2*x+3/21/2*x-1/21/2*x+1/21.2500输入程序 M=5;R 1=M *abs(x-X(1)* (x-X(2)/2运行后输出误差限为R1 =1.87502).解输入程序 X=0,1; Y =exp(-X),l0 仁 poly(X (2)/( X(1)- X( 2),l11= poly(X(1)/( X(

11、2)- X(1), l0=poly2sym (l01),l1= poly2sym (l11), P = l01* Y(1)+ l11* Y(2), L=poly2sym (P),运行后输出基函数l0和l1及其插值多项式的系数向量 P和插值多项式L为l0 = l1 =P =-x+1 x-0.6321 1.0000L =-1423408956596761/2251799813685248*x+1输入程序 M=1;x=0:0.001:1; R1= M*max(abs(x-X(1).*(x-X(2)./2运行后输出误差限为R1 =0.1250.3).解输入程序 X=-2,0,1,2; 丫 =17,1,

12、2,17; p仁poly(X(1); p2=poly(X (2); p3=poly(X (3); p4=poly(X )；101 =conv ( conv (p2, p3), p4)/(X(1)- X(2)*(X(1)- X(3)* ( X(1)-X），111=conv (conv(p1,p3),p4)/(X(2)- X(1)*(X(2)- X(3)* ( X(2)-X），121 =conv (conv(p1,p2),p4)/(X(3)- X(1)*(X(3)- X(2)* ( X(3)-X），131 =conv (conv(p1,p2),p3)/(X(4)- X(1)*(X(4)- X(2)

13、* ( X(4)-X(3),I0=poly2sym (l01),l1= poly2sym (I11),l2=poly2sym (121), I3=poly2sym (l31),P = l01* Y(1)+ 111* Y(2) + 121* Y(3) + l31* 丫,运行后输出基函数丨0,丨1, l2和l 3及其插值多项式的系数向量P (略)为l0 =-1/24*xA3+1/8*xA2-1/12*x，l1 =1/4*xA3-1/4*xA2-x+1l2 =-1/3*xA3+4/3*x，l3 =1/8*xA3+1/8*xA2-1/4*x输入程序 L=poly2sym (P),x=0.6; 丫 =

14、polyval(P,x)运行后输出插值多项式和插值为L =丫 =xA3+4*xA2-4*x+10.2560.输入程序 syms M; x=0.6;R3=M*abs(x-X(1)*(x-X(2) *(x-X(3) *(x-X(4)/24运行后输出误差限为R3 =91/2500*M即R 賂 0.04f(4)(，(2.00, 2.00).5、实验数据记录与分析1) 运行后输出基函数丨。和丨1及其插值多项式的系数向量 P (略八 插值多项式L和插 值丫为l0 =l1 = L = 丫 =-1/2*x+3/21/2*x-1/21/2*x+1/21.2500运行后输出误差限为R1 =1.87502) 运行后

15、输出基函数丨。和丨1及其插值多项式的系数向量 P和插值多项式L为l0 = l1 =P =-x+1 x-0.6321 1.0000L =-1423408956596761/2251799813685248*x+1运行后输出误差限为R1 =0.1250.3) 运行后输出基函数I。，丨1,丨2和13及其插值多项式的系数向量P (略)为10 =-1/24*xA3+1/8*xA2-1/12*x，I1 =1/4*xA3-1/4*xA2-x+1I2 =-1/3*xA3+4/3*x，I3 =1/8*xA3+1/8*xA2-1/4*x输入程序 L=poly2sym (P),x=0.6; Y = polyval(

16、P,x)运行后输出插值多项式和插值为L =Y =xA3+4*xA2-4*x+10.2560.输入程序 syms M; x=0.6;R3=M*abs(x-X(1)*(x-X(2) *(x-X(3) *(x-X(4)/24运行后输出误差限为R3 =91/2500*M即R ： 0.04f ()，(-2.00, 2.00).6、实验结论通过本次实验，运用了 matlab进行Taylor插值、线性插值、抛物插值以及一般 情形的Lagrange插值的简单程序，掌握了插值算法的基本理论，我在本次实验中获得 了巨大的收获，对数值分析这门课程有了更深刻的理解。指导教师评语和成绩评定指导教师签字：年 月 日实验名

17、称数值积分与数值微分实验时间2012年10月25日学生姓名姚海保实验地点9#405数学实验室1实验所用软件WIN7操 作系统、Matlab2、实验目的1. 掌握复合梯形法求积分、复合辛普生法方法求积分公式的程序设计2. 掌握Romberg积分的用法。3. 掌握高斯公式的程序设计3、实验内容1 1X2（1）辛普森公式计算1 =亠e 2 dx，取n 20 001个等距节点，并将计算结果与 P2兀o精确值比较，然后再取n =13计算，观察n对误差的影响.3T（2估计用辛普森公式计算定积分1 = F esinxdx时的误差，取h 兀/40Lo（3编写复合辛普森（Simpson）数值积分的MATLA主程

18、序4、实验方法、步骤1. 了解matlab的硬件和软件必备环境；2. 启动 matlab ；3. 学习一些简单的matlab编程。4. 了解各种求积分的方法的原理。5.编写程序。5、实验数据记录与分析1 1x!（1）辛普森公式计算1 -e 2dx，取n- 20 001个等距节点，并将计算结果与2応0精确值比较，然后再取n =13计算，观察n对误差的影响.解 由n = 2m 仁20 001，得m = 10 000 .根据辛普森(Simpson)公式编写并输入 下面的程序 a=0;b=1;m=10000; h=(b-a)/(2*m); x=a:h:b; y=exp(-x42)./2)./(sqrt

19、(2*pi);z1=y(1)+y(2*m+1); z2=2*sum(y(2:2:2*m); z3=4*sum(y(3:2:2*m);z=(z1+z2+z3)*h/3, syms t,f=exp(-tA2)/2)/(sqrt(2*pi);intf=in t(f,t,a,b), Fs=double(i ntf); Juewucha=abs(z-Fs)运行后屏幕显示用辛普森公式(9.11)计算定积分I的近似值z和精确值intf及其绝对 误差Juewucha (取n = 20 001个等距节点).(2估计用辛普森公式计算定积分I = 2 esinxdx时的误差，取hhQ40根据估计误差公式，先输入求f

20、(x)的程序syms x,y=exp(si n( x); yx4=diff(y,x,4)运行后输出被积函数的四阶导函数然后在输入误差估计程序h=pi/40; x=0:0.00001:pi/2;yx4=si n(x).*exp(si n(x)-4*cos(x).A2.*exp(si n(x)+3*si n(x).A2.*exp(si n(x)-6*si n(x).*cos(x).A2.*exp(si n(x)+cos(x).A4.*exp(si n(x);juyx4= abs(yx4); RS=(hA4)*(pi/2)*max(juyx4)/180运行后屏幕显示误差估计值RS =3.610450

21、295892220e-006(3)写复合辛普森(Simpson)数值积分的MATLA主程序fun cti on y=comsimps on(fun, a,b ,n)z1=feval (f un, a)+ feval (fun ,b);m=n/2;h=(b-a)/(2*m); x=a;z2=0; z3=0; x2=0; x3=0;for k=2:2:2*mx2=x+k*h; z2= z2+2*feval (fu n, x2);endfor k=3:2:2*mx3=x+k*h; z3= z3+4*feval (fu n, x3);endy=(z1+z2+z3)*h/3;6、实验结论通过本次实验果，

22、了解各种求积分的方法的原理。但在试验中还存在一些问题如：对matlab 些操作命令执行不大完整，容易出错等。所以在以后的学习中要努力学习 并在学习生活中运用，以达到实践与理论相结合。指导教师评语和成绩评定指导教师签字：年 月日实验名称常微分方程差分方法实验时间2012年11月8日学生姓名姚海保实验地点9#405数学实验室1实验所用软件WIN7操 作系统、Matlab2、实验目的1. 编写程序实现Euler求解方法及改进Euler求解方法2. 掌握龙格-库塔方法的用法3. 掌握方程组和高阶方程的程序实现3、实验内容（1）用欧拉方法求初值问题浮=x-y , OxEI,1 dx yL=o 二1的数值

23、解，分别取h =0.075 0,0.007 5，并计算误差，画出精确解和数值解的图形（2）用改进的欧拉方法梯形格式求解区间0,2上的初值问题叟= x0=0;y0=1;b=1;h=0.0750;P=Eulerli1(xO,yO,b,h)在MATLA工作窗口输入下面的程序h1=0.0075; P1=Eulerli1(x0,y0,b,h1)legend(h1=0.0075时，dy/dx=x-y , y(0)=1 在0,1上的数值解,精确解y=x-1+2 exp(-x)(2)方程为dy = 3y8x7, y(0) =1求解区间0,2取步长h =0.05,精度为10J dx输入程序 x0=0;y0=1;

24、 b=2; tol=0.1; h=0.05;X,Yt, n,Pt=odtixi ng1(fu nfcn,x0,b,y0,h,tol)hold onS1=8/3*X-29/9+38/9*exp(-3*X);plot(X,S1,b-), hold offlegend(h=0.05,用梯形公式计算 dy/dx=8x-3y-7，y(0)=1 在0,2上的数值 解,dy/dx=8x-3y-7，y(0)=1 在0,2上的精确解)juwYt=S1-Yt; xiwYt=juwYt./Yt; Lt=Pt,S1,juwYt,xiwYt运行后屏幕显示取精度为10d，分别用梯形公式和向前欧拉公式求解此初值问题在区 间

25、0,2上的自变量X处数值解Y (i=t,q)和精确解S,步长H Y的相对误差xiwYi和绝对误差juwYi (略)及其数值解和精确解的图形(3)微分方程组为dy _ 2xy1 dx1 + x2J y | x 鱼=0取c1 =1/4, c2 =3/4, a2 = b21 = 2 /3， h = 1/4二阶龙格一库塔方法的MATLA主程序：fun ctio n k,X,Y,fxy,wch,wucha,P=RK2(fu nfcn ,fu n,x0,b,C,y0,h) x=x0; y=y0;p=128;n=fix(b-x0)/h);fxy=zeros(p,1);wucha=zeros(p,1);wch

26、=zeros(p,1);X=zeros(p,1);Y=zeros(p,le ngth(y);k=1; X(k)=x; Y( k,:)=y;%绘图clc,x,h,y%计算%fxy=fxy(:);for k=2:n+1x=x+h;a2=C(3);b2 仁C(4);c仁 C(1);c2=C (2);x1=x+a2*h;k仁feval(f unfcn, x,y);y仁 y+b21*h*k1;k2=feval(fu nfcn,x1,y1);fxy(k)=feval(fu n, x);y=y+h*(c1*k1+c2*k2);X(k)=x;Y( k,:)=y;k=k+1;plot(X,Y,mh,X,fxy,

27、bo)y(x0)=y0 在x0,b上grid,xlabel( 自变量 X),ylabel( 因变量 Y) legend(用二阶龙格-库塔方法计算dy/dx=f(x,y), 的数值解,y/dx=f(x,y), y(xO)=yO 的精确解 y=f(x)endfor k=2:n+1wucha(k)=norm(Y(k-1)-Y(k); wch(k)=n orm(fxy(k)-Y(k); endX=X(1:k);Y=Y(1:k,:); fxy=fxy(1:k,:); n=1:k;wucha=wucha(1:k,:); wch=wch(1:k,:); P= n,X,Y,fxy,wch,wucha;再在MA

28、TLA工作窗口输入下面的程序 x0=0;b=2;C=1/4,3/4,2/3,2/3;y0=0;h=1/4;k,X,Y,fxy,wch,wucha,P=RK2(fu nfcn ,fu n,x0,b,C,y0,h)6、实验结论Euler求解方法,龙通过本次试验，基本掌握了编程实现Euler求解方法及改进格-库塔方法的用法和方程组，高阶方程的程序实现。指导教师评语和成绩评定指导教师签字：年 月 日实验名称方程求根的迭代法实验时间2012年11月29日学生姓名姚海保实验地点9#405数学实验室1实验所用软件WIN7操 作系统、Matlab2、实验目的1. 掌握二分法与牛顿法的基本原理及应用2. 掌握牛

29、顿法、弦截法的程序设计并用具体例子加以实现3、实验内容(1) 编与二分法的主程序。(2) 编写牛顿切线法的MATLA主程序(3) 用牛顿切线法求方程2x3 - 3x2 +1 =0在x0 = -0.4和0.9的近似根，要求精度 呂=10.4、实验方法、步骤1. 了解matlab的硬件和软件必备环境；2. 启动 matlab ;3. 学习一些简单的matlab编程。4. 了解二分法与牛顿法的的基本原理及应用5. 了解牛顿法，弦截法的程序设计。6. 编写程序5、实验数据记录与分析(1) 编写程序如下：function k,x,wuca,yx=erfen(a,b,abtol)a(1)=a; b(1)=

30、b;ya=fun(a(1); yb=fun(b(1); %程序中调用的 fun.m 为函数if ya* yb0,disp(注意：ya*yb0,请重新调整区间端点a和b.), return endmax 1=-1+ceil(log(b-a)- log(abtol)/ log(2); % ceil向取整for k=1: max1+1a;ya=fu n( a); b;yb=fu n( b); x=(a+b)/2; yx=fun(x); wuca=abs(b-a)/2; k=k-1; k,a,b,x,wuca,ya,yb,yxif yx=0a=x; b=x;elseif yb*yx0b=x;yb=yx

31、;elsea=x; ya=yx;endif b-a abtol , retur n, endendk=max1; x; wuca; yx=fu n( x);(2) 编程如下：function k,xk,yk,piancha,xdpiancha=newtonqx(xO,tol,ftol,gxmax) x(1)=x0;for i=1: gxmax x(i+1)=x(i)-fnq(x(i)/(dfnq(x(i)+eps); pia ncha=abs(x(i+1)-x(i);xdpiancha= piancha/( abs(x(i+1)+eps); i=i+1;xk=x(i);yk=f nq (x(i

32、); (i-1) xk yk pia ncha xdpia ncha if (abs(yk)ftol)&( (pia nchagxmaxdisp(请注意：迭代次数超过给定的最大值gxmax )k=i-1; xk=x(i);(i-1) xk yk pia ncha xdpia ncha return;end(i-1),xk,yk,pia ncha,xdpia ncha;(3)方程：2x3 -3x2 1-0在MATLA工作窗口输入程序 k,xk,yk,pia ncha,xdpia ncha=newto nqx(-0.4,0.001,0.001,100) 运行后输出初始值X。二(.4的迭代结果.同理

33、在MATLAB：作窗口输入程序 k,xk,yk,pia ncha,xdpia ncha=newto nqx(0.9,0.001,0.001,100)运行后输出初始值x=0.9的迭代结果.6、实验结论通过本次试验 掌握了二分法与牛顿法的基本原理及应用和牛顿法弦截法的程序设计并能运用熟练举出例子等。指导教师评语和成绩评定指导教师签字：年 月 日实验名称线性方程组的迭代法实验时间2012年12月13日学生姓名姚海保实验地点9#405数学实验室1实验所用软件WIN7操 作系统、Matlab2、实验目的设计程序实现雅可比迭代法，高斯-赛德尔迭代法等算法。掌握方程求根的基本理 论。并用实例在计算机上计算。

34、3(、实验内容1)编写用雅可比迭代解线性方程组 AX = b的MATLAB程序(2) 编写用高斯-塞德尔迭代定义解线性方程组 AX= b的MATLA主程序(3) 用高斯-塞德尔迭代定义的MATLAB程序解下列线性方程组,取初始值Xi(0),x20),x30) =(0,0,0)，要求当 x(k 也-x(k)L10 时，迭代终止3x1 +4x2 5x3 +7x4 =5,2Xi 8x2 *3x3 2x4 =2,4x=0disp(请注意：系数矩阵A不是严格对角占优的，此雅可比迭代不一 定收敛)returnendendif a(i)wucha )&( xdwcXwucha)disp(请注意：雅可比迭代次

35、数已经超过最大迭代次数max1 )enda,X=X;jX=X1,(2) 编写程序如下fun ction X=gsdddy(A,b,X0,P,wucha,max1)D=diag(diag(A);U=-triu(A,1);L=-tril(A,-1); dD=det(D);if dD=0disp(请注意：因为对角矩阵D奇异，所以此方程组无解.)elsedisp(请注意：因为对角矩阵D非奇异，所以此方程组有解.)iD=inv(D-L); B2=iD*U;f2=iD*b;jX=Ab; X=X0;n m=size(A);for k=1:max1X1= B2*X+f2; djwcX=norm(X1-X,P)

36、;xdwcX=djwcX/( norm(X,P)+eps);if (djwcXwucha)|(xdwcXwucha)returnelsek,X1,k=k+1;X=X1;endendif (djwcXwucha)|(xdwcX A=3 4 -5 7;2 -8 3 -2;4 51-13 16;7 -2 21 3;b=5;2;-1;21;X0=0 0 0 0;X=gsdddy(A,b,X0,i nf,0.001,100)运行后输出结果请注意：因为对角矩阵D非奇异，所以此方程组有解.请注意：高斯-塞德尔迭代的记过没有达到给定的精度，并且迭代次数已经 超过最大迭代次数max1,方程组的精确解jX和迭代向

37、量X如下：jX =0.1821 -0.25710.72861.3036X = 1.0e+142 *0.28830.10620.3622 -3.13746、实验结论通过本次试验，设计程序实现了雅可比迭代法，高斯-赛德尔迭代法等算法。基本掌握 方程求根的基本理论。并用实例在计算机上计算。实验的结果与预期中的基本一样， 基本达到了实验要求。指导教师评语和成绩评定指导教师签字：年 月 日实验名称线性方程组的直接法实验时间2012年12月20日学生姓名姚海保实验地点9#405数学实验室1实验所用软件WIN7操 作系统、Matlab2、实验目的1. 熟悉Gauss列主元消去法，编出实用程序。2. 认识选主

38、元技术的重要性。4. 编制程序，用Gauss列主元消去法求解线性方程组 Ax = b。5. 编制程序，用平方根法求解线性方程组三角方程组3、实验内容（1）编写用列主元消元法解线性（2）用列主元消元法解线性方程:1方程组 AX = b的MATLA程序组的 MATLAB?序解方程组-X2 X3 +X4 =0,Xi - X2 + X3 -3X4 = 1,.2x1 2x2 4x3 +6x4 =T,% -2x2 -4x3 + x4 =-1.(3) 编写解上三角形线性方程组 AX二b的MATLABS序(4) 用解上三角形线性方程组的 MATLAB?序解方程组5xi x2 十2 x3 +3x4 =20, -2x2 +7x3 4x4 = -7,16x3 +5x4 = 4,3x4 = 6.4、实验方法、步骤1. 了解matlab的硬件和软件必备环境;2. 启动 matlab ；3. 学习一些简单的matlab编程。4. 了解Gauss列主元消去法的原理。5.编写程序。5、实验数据记录与分析(1)编写程序如下fun ctio n RA,RB ,n ,X=liezhu(A,b)B=A b; n=length(b); RA=rank(A);RB=ra nk(B);zhica=R