MATLAB第三次上机实验报告.docVIP

电子科技大学 电子工程 学院 标 准 实 验 报 告 （实验）课程名称 MATLAB与数值分析 学生姓名： 学 号： 指导教师： 实验名称 实验 通过上机实验，使学生对数值积分、微分方程求解方法有一个初步的理解。 实验涉及的核心知识点:数值积分和积分方程的数值求解、常微分方程数值求解。 实验重点与难点:算法设计和MATLAB编程 实验内容 1. 选用复合Simpson公式，计算 并用Matlab的符号运算工具箱计算其精确值。比较结果，找出问题原理，提出解决问题的方法。 2. 求积分方程 的数值解和精确解，分析二者的差异。 3. 利用Euler法对不同的步长求下面初值问题的数值解： 并通过绘图，与方程的解析解进行比较。 实验数据及结果分析 1. 选用复合Simpson公式，计算 并用Matlab的符号运算工具箱计算其精确值。比较结果，找出问题原理，提出解决问题的方法。 function s=simpr1(f,a,b,n) h=(b-a)/(2*n); s1=0; s2=0; for k=1:n x=a+h*(2*k-1); s1=s1+feval(f,x); end for k=1:(n-1) x=a+h*k*2; s2=s2+feval(f,x); end s=h*(feval(f,a)+feval(f,b)+4*s1+2*s2)/3; function y=f(x) y=(log(1+x))*(log(1-x)); (1)在主界面中输入 simpr1(f,-1,1,100) 所得结果为： ans = -Inf (2)在主界面中输入 simpr1(f,-0.9999,0.9999,100) 所得结果为 ans = -1.1104 2. 求积分方程 的数值解和精确解，分析二者的差异。 对原方程进行化简和变形，得： 原问题转化为一个常微分初值问题。 function [x,y]=euler(fun,x0,y0,xN,N) y(1)=y0; h=(xN-x0)/N; x=x0:h:xN; for n=1:N x(n+1)=x(n)+h; y(n+1)=y(n)+h*feval(fun,x(n),y(n)); end plot(x,y,b) f=inline(dsolve(Dy=y,y(0)=1,x)); F=inline(y,x,y) y=euler(F,0,1,5,500); hold on X=linspace(0,5); Y=f(X); plot(X,Y,r) legend(Euler近似,y=e^t,0) 3. 利用Euler法对不同的步长求下面初值问题的数值解： 并通过绘图，与方程的解析解进行比较。 function E=Euler_1(fun,x0,y0,xN,N) x=zeros(1,N+1);y=zeros(1,N+1); x(1)=x0;y(1)=y0; h=(xN-x0)/N; for n=1:N x(n+1)=x(n)+h; y(n+1)=y(n)+h*feval(fun,x(n),y(n)); end T=[x,y]; plot(0:N,0:N,x,y,o,x,y,-) function z=f(x,y) z=y-20; 在主界面中输入 Euler_1(f,-10,10,10,10) 在主界面中输入 Euler_1(f,-10,10,10,20) 在主界面中输入 Euler_1(f,-10,10,10,30) 该常微分方程应为，图像如下 总结及心得体会 由上面的实验结果可以看出有时候复合Simpson公式并不能求出所需解。本题的原因是在取1，-1时式子没有意义，无法进行运算。我运用的解决方法是分别用0.9999，-0.9999代替。精度越高，此方法所得结果越接近于真实值，但是和真实值还是有差距。 由图像可知，Euler公式能够较精确的求出常微分方程的解。适当的步长h，能比较精确的表示该常微分方程的解。 由三个图像比较可知，用Euler公式求解常微分方程初值问题时，步长越短，其数值解越接近其解析解。 此次MATLAB实验涉及的核心知识点是数值积分和积分方程的数值求解、常微分方程数值求解。实验重点与难点是算法设计和MATLAB编程 通过此次上机实验，我明白了MATLAB求积分、常微分方程初值问题的原理，明白了这些原理方法在实际应用中需要结合实际操作进行适当处理，而且要通过更好的算法，尽量减小误差，得到更加精确的解。