第3章MATLAB在高等数学中的应用_Matlab及在电子信息类中的应用汇编.ppt

3.4.3函数的数值积分与微分 1．函数的数值积分 (1) 低阶数值积分函数quad 格式一：q=quad(‘fun’,a,b) 功能：采用自适应的Simpson积分方法，返回函数‘fun’ 在上限 a 和下限 b之间的数值积分。当给定一个输入值向量，‘fun’ 必须返回一个输出向量。函数‘fun’可以是函数名、函数句柄或字符串。 格式二：q=quad(‘fud’,a,b,tol) 功能：按指定绝对误差tol返回数值积分值, tol 缺省值为 1e-6。 (2) 高阶数值积分函数quadl ,具有与quad同样的调用方法，在处理软奇异函数时比quad更好。 例3-50：求y=e-x+x2在[0,1]区间上的定积分。 (3) 梯形面积法的积分函数trapz( ) 格式一：T=trapz(Y) 功能：以单位间隔，采用计算若干梯形面积的和来计算某函数的近似积分。如果Y为向量，计算Y的积分；如果Y是矩阵，得一个每列积分的行向量；如果Y为多维数组，则沿第一个非单元素维计算。 格式二：T=trapz(X,Y) 功能：用梯形积分法，依据X计算Y的积分。如果X为矢量，则Y必须是同大小的矢量；如果X是一列向量，并且数组Y第一非单元素维长度为length(X)，则在该维中计算。 例3-51: x1=-1:0.17:2;y1=humps(x1); Areal=trapz(x1,y1) x2=-1:0.07:2;y2=humps(x2);Aarea=trapz(x2,y2) (4) 双重积分函数dblquad MATLAB提供了一个求双重积分的函数dblquad，其基本调用格式为： 格式：Q=dblquad(fun,xmin,xmax,ymin,ymax,tol) 功能：按指定精度tol，对指定函数 f(x, y)在[xmin, xmax]范围和[ymin, ymax]范围进行双重积分。精度tol缺省时默认精度为1e-6。 例3-52:在0xπ, πy2π区域内计算ysinx+xcosy的二重积分 (5) 不定积分的计算 对于函数不定积分的计算，可以采用定积分函数来求不定积分的数值解。方法是： 固定积分下限，用for循环，把积分上限逐步增加即可实现。 例3-53：试求humps函数以x=0为下限的不定积分。 2．函数的数值微分 (1)通过计算数组中元素间的差分函数diff来粗略计算微分函数。 例3-54：试利用diff函数计算y=f(x)的微分，已知该函数的数据如下： X: 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Y： -0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2 (2)对数据拟合后再利用polyder函数微分 为了实现对数值较准确的微分，最好先用最小二乘曲线拟合这种数据，然后对所 得到的多项式进行微分。 例3-54：试利用polyder函数计算y=f(x)的微分，已知该函数的数据如下： X: 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Y： -0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2 3.4.4 常微分方程的数值求解 1．初始值的常微分方程求解 格式：[T,Y] = ode45(‘odefun’,tspan,y0) 功能：返回由文件‘odefun’所定义的具有初始条件为y0、时间t变化范围为[t0, tfinal] 的微分方程y’= f(t,y)的解，其中tspan = [t0,tfinal]。向量T中的每一列对应着矩阵Y的每一列。 例3-55：用数值积分法求解微分方程：y’’=-y+1-t2/2π。设初始时间t0=0,终止时间tf=3π，初始条件y(0)=0,y’(0)=0。 任何高阶常微分方程都可以变换成一阶微分方程，即表示成右函数形式，这是利用龙格-库塔法求解微分方程的前提。 2．边界条件的常微分方程求解 对于如下的微分方程： 用于边界条件的常微分方程求解问题: 函数bvp4c 格式：sol = bvp4c(‘odefun’,bcfun,solinit) 其中‘odefun’为常微分方程函数，bcfun为边界条件函数，solinit为求解的初始值。输出sol是一个结构，它