电子科技大学matlab-上机实验五.pdfVIP

Matlab 与数值分析实验报告（五） 沈汉男 一、 实验的性质和任务 通过上机实验，使学生对数值积分、微分方程求解方法有一个初步的理解。 二、 实验内容 1. 选用复合 Simpson 公式，计算 1 ln(1 x ) ln( 1 x)dx 1 并用 Matlab 的符号运算工具箱计算其精确值。比较结果，找出问题原理，提出解决问题的 方法。 2. 求积分方程 2 1 t t y(t ) e y(s)ds e e 1 0 的数值解和精确解，分析二者的差异。 3. 利用 Euler 法对不同的步长求下面初值问题的数值解： y (t ) y (t ) 20 y(0) 10 并通过绘图，与方程的解析解进行比较。 三、 实验过程 实验 1 a. 实验过程： （1）、创建文件 simpr1.m 如下： 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*2*k; s2=s2+feval( f ,x); end s=h*(feval( f ,a)+feval( f ,b)+4*s1+2*s2)/3 创建文件 f.m 如下： function y=f(x) y=(log(1+x))*(log(1-x)); 在主界面中输入 simpr1(f,-1,1,10) 所得结果为： s = -Inf ans = -Inf （2 ）在主界面中输入程序 clear all syms x y=int(log(1+x)*log(1-x),x,-1,1) eval(y) 所得结果为 y = 2*log(2)^2 - pi^2/3 - log(16) + 4 ans = -1.1016 b. 实验结果分析： 由上面的实验结果可以看出有时候复合 Simpson 公式并不能求出所需解 实验 2： a. 实验过程： 原积分方程为 2 1 t t y (t ) e y(s)ds e e 1 0 对其进行化简和变形，得： y y e 1 y (0)