第2章MATLAB数值计算 数据处理与多项式计算 数值微积分.pptVIP

例 6.18 设 x 由 [0,2π] 间均匀分布的 10 个点组成，求 sin x 的 1~3 阶差分。 命令如下： X=linspace(0,2*pi,10); Y=sin(X); DY=diff(Y); % 计算 Y 的一阶差分 D2Y=diff(Y,2); % 计算 Y 的二阶差分，也可用命令 diff(DY) 计算 D3Y=diff(Y,3); % 计算 Y 的三阶差分，也可用 diff(D2Y) 或 diff(DY,2) 例 6.19 用不同的方法求函数 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 数值积分 1. 数值积分基本原理 求解定积分的数值方法多种多样，如简单 的梯形法、辛普生 (Simpson) ? 法、牛顿－柯 特斯 (Newton-Cotes) 法等都是经常采用的方 法。它们的基本思想都是将整个积分区间 [a,b] 分成 n 个子区间 [x i ,x i+1 ] ， i=1,2,…,n ，其 中 x 1 =a ， x n+1 =b 。这样求定积分问题就分解 为求和问题。 2. 数值积分的实现 (1) 被积函数是一个解析式 MATLAB 提供了 quad 函数和 quadl 函数来求 定积分。它们的调用格式为： quad(filename,a,b,tol,trace) quadl(filename,a,b,tol,trace) 例 6.20 用两种不同的方法求定积分。 先建立一个函数文件 ex.m ： function ex=ex(x) ex=exp(-x.^2); 然后在 MATLAB 命令窗口，输入命令： format long I=quad(ex,0,1) % 注意函数名应加字符引号 I = 0.74682418072642 I=quadl(ex,0,1) I = 0.74682413398845 也可不建立关于被积函数的函数文件，而使用语句函数 ( 内联函数 ) 求解，命令如 下： g=inline(exp(-x.^2)); % 定义一个语句函数 g(x)=exp(-x^2) I=quadl(g,0,1) % 注意函数名不加 号 I = 0.74682413398845 format short (2) 被积函数由一个表格定义 在科学实验和工程应用中，函数关系往往是不知道 的，只有实验测定的一组样本点和样本值，这时， 就无法使用 quad 函数计算其定积分。在 MATLAB 中，对由表格形式定义的函数关系的求定积分问 题用 trapz(X,Y) 函数。其中向量 X 、 Y 定义函数关 系 Y=f(X) 。 X 、 Y 是两个等长的向量： X=(x1,x2,…,xn) ， Y=(y1,y2,…,yn) ，并且 x1x2…xn ，积分区间是 [x1 ， xn] 。 例 6.21 用 trapz 函数计算定积分。 在 MATLAB 命令窗口，输入命令： X=0:0.01:1; Y=exp(-X.^2); trapz(X,Y) ans = 0.7468 (3) 二重积分数值求解 使用 MATLAB 提供的 dblquad 函数就可以直 接求出上述二重定积分的数值解。该函数 的调用格式为： I=dblquad(f,a,b,c,d,tol,trace) 该函数求 f(x,y) 在 [a,b] × [c,d] 区域上的二重定 积分