matlab上机实践五(计算物理班).docVIP

实验五 MATLAB解方程实验要求：为达到理想的实验效果，同学们务必做到：实验前认真准备，要根据实验目的和实验内容，复习好实验中可能要用到的命令，想好编程的思路，做到胸有成竹，提高上机效率。实验过程中积极思考，要深入分析命令、程序的执行结果以及各种屏幕信息的含义、出现的原因并提出解决办法。实验后认真总结，要总结本次实验有哪些收获，还存在哪些问题，并写出实验报告。实验报告应包括实验目的、实验内容、流程图（较大程序）、程序（命令）清单、运行结果以及实验的收获与体会等内容。同学们在上机过程中会碰到各种各样的问题，分析问题和解决问题的过程就是积累经验的过程。只要同学们按照上面3点要求去做，在学完本课程后就一定会有很大的收获。实验仪器：计算机实验时间：2018.5.实验原理：Matlab线性方程组求解：直接解法：左除，矩阵求逆；矩阵分解求解：LU分解迭代解法：jacobi迭代；gauss-serdel迭代2． Matlab非线性方程求解： 单变量非线性方程求解： z=fzero(@fname,x0,tol,trace) z=fzero(fname,x0,tol,trace)非线性方程组求解： X=fsolve(fun,X0,option) X=fsolve(@fun,X0,option)3.常微分方程（组）初值问题的数值解法龙格－库塔法 [t,y]=solver(fname,tspan,y0) [t,y]= solver(@fname,tspan,y0) 实验目的： 1.掌握matlab求解线性方程组2.掌握matlab求解非线性方程（组）3. 掌握常微分方程（组）初值问题的数值解法实验内容：求解线性方程组利用矩阵除法，矩阵求逆，矩阵分解，jacobi迭代法和和Gauss-serdel迭代法（初值x0=1，y0=1，z0=0.5）求解方程组。程序1：A=[2,3,4;3,7,4;1,-6,1];b=[8,3,5]’;x=A\b结果1：x = -2.4857 -0.6286 3.7143程序2：x=inv(A)*b结果2：x = -2.4857 -0.62863.7143程序3：A=[2,3,4;3,7,4;1,-6,1];b=[8,3,5]’;[L,U]=lu(A);x=U\(L\b)结果3：x = -2.4857 -0.6286 3.7143程序4：function [y,n]=jacobi(A,b,x0,eps)if nargin==3 eps=1.0e-6;elseif nargin3 error returnend D=diag(diag(A)); %求A的对角矩阵L=-tril(A,-1); %求A的下三角阵U=-triu(A,1); %求A的上三角阵B=D\(L+U);f=D\b;y=B*x0+f;n=1; %迭代次数while norm(y-x0)=eps x0=y; y=B*x0+f; n=n+1;endA=[2,3,4;3,7,4;1,-6,1];b=[8,3,5]’;[x,n]=jacobi(A,b,[1,1,0.5]’,1.0e-6)结果：x = NaN NaN NaNn = 1320程序5：function [y,n]=gauseidel(A,b,x0,eps)if nargin==3 eps=1.0e-6;elseif nargin3 error returnend D=diag(diag(A)); %求A的对角矩阵L=-tril(A,-1); %求A的下三角阵U=-triu(A,1); %求A的上三角阵G=(D-L)\U;f=(D-L)\b;y=G*x0+f;n=1; %迭代次数while norm(y-x0)=eps x0=y; y=G*x0+f; n=n+1;endA=[2,3,4;3,7,4;1,-6,1];b=[8,3,5]’;[x,n]= gauseidel(A,b,[1,1,0.5])结果：x = NaN NaN NaNn = 501求解非线性方程1）在附近的根。程序：function fx=funx(x)fx=x.^2+sin(x)-exp(x);z=fzero(@funx,1)结果：z = -1.10672) 在给定的初值，，下，求方程组的数值解。程序：function q=myfun(p)x=p(1);y=p(2);z=p(3);q(1)=sin(x)+y