实验09 数值微积分与方程数值解(第6章).docVIP

PAGE PAGE 1 《 《数学软件》课内实验 王平 实验09 数值微积分与方程数值求解 （第6章 MATLAB数值计算） 一、实验目的 1. 掌握求数值导数和数值积分的方法。 2. 掌握代数方程数值求解的方法。 3. 掌握常微分方程数值求解的方法。 二、实验内容 1. 求函数在指定点的数值导数 程序及运行结果： 2. 用数值方法求定积分 (1) 的近似值。 程序及运行结果： (2) 程序及运行结果： 3. 分别用3种不同的数值方法解线性方程组 程序及运行结果： 4. 求非齐次线性方程组的通解 程序及运行结果（提示：要用教材中的函数程序line_solution）： 5. 求代数方程的数值解 (1) 3x+sinx-ex=0在x0=1.5附近的根。 程序及运行结果（提示：要用教材中的函数程序line_solution）： (2) 在给定的初值x0=1，y0=1，z0=1下，求方程组的数值解。 程序及运行结果： 6. 求函数在指定区间的极值 (1) 在(0,1)内的最小值。 程序及运行结果： (2) 在[0,0]附近的最小值点和最小值。 程序及运行结果： 7. 求微分方程的数值解，并绘制解的曲线 程序及运行结果（注意：参数中不能取0，用足够小的正数代替）： 令y2=y,y1=y ，将二阶方程转化为一阶方程组： 8. 求微分方程组的数值解，并绘制解的曲线 程序及运行结果： 三、实验提示 四、教程：第6章 MATLAB数值计算（2/2） 6.2 数值微积分 p155 6.2.1 数值微分 1. 数值差分与差商 对任意函数f(x)，假设h0。 向前差分： 向后差分： 中心差分： 当步长h充分小时，有 向前差商： 向后差商： 中心差商： 2. 数值微分的实现 MATLAB没有提供求数值导数的函数，只有计算向前差分的函数diff。 DX=diff(X)：计算向量X的向前差分，DX(i)=X(i+1)-X(i)，i=1,2,…,n-1。 DX=diff(X,n)：计算X的n阶向前差分。 例如，diff(X,2)=diff(diff(X))。 DX=diff(A,n,dim)：计算矩阵A的n阶差分，dim=1时(缺省状态)，按列；dim=2，按行。 例6.18 （向前差分）求1~3阶差分p156 设x由[0,2π]间均匀分布的10个点组成，求sinx的1~3阶差分。 X=linspace(0,2*pi,10); Y=sin(X); DY=diff(Y) D2Y=diff(Y,2) D3Y=diff(Y,3) DY = 0.6428 0.3420 -0.1188 -0.5240 -0.6840 -0.5240 -0.1188 0.3420 0.6428 D2Y = -0.3008 -0.4608 -0.4052 -0.1600 0.1600 0.4052 0.4608 0.3008 D3Y = -0.1600 0.0556 0.2452 0.3201 0.2452 0.0556 -0.1600 例6.19 （数值导数）3种方法求导p157 用不同的方法求f(x)的数值导数，并在同一个坐标系中做出f (x)的图像。 f=inline(sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+5*x+2); g=inline((3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-x+12)/2+1/6./(x+5).^(5/6)+5);%导函数 x=-3:0.01:3; p=polyfit(x,f(x),5);%用5次多项式p拟合f(x) dp=polyder(p); %对拟合多项式p求导数dp dpx=polyval(dp,x);%求dp在假设点的函数值 dx=diff(f([x,3.01]))/0.01;%直接对f(x)求数值导数 gx=g(x);%求函数f的导函数g在假设点的导数 plot(x,dpx,x,dx,.,x,gx,-); 6.2.2 数值积分 p157 1. 数值积分基本原理 求解定积分的数值方法： 梯形法 辛普生(Simpson)法 牛顿－柯特斯(Newton-Cotes)法 基本思想： 将整个积分区间[a,b]分成n个子区间 [xi, xi+1]，i=1,2,…,n，其中x1=a，xn+1=b 这样求定积分问题就分解为求和问题。 2. 数值积分的实现 (1) 被积函数是一个解析式 调用格式： quad(fname,a,b,tol,trace) quadl(fname,a,b,tol,trace) fname是被积函数名。 a和b分别是定