《MATLAB 8.X程序设计及典型应用》课件第4章.pptVIP

在指令窗中执行exm4_31，结果为4.5.2数值积分

1.一重积分

常用的一重积分函数有两个，即quad()和quadl()。

quad()函数是基于变步长辛普生法的求定积分函数，调用格式如下：

[I，n]=quad(fun，a，b，tol)quadl()是基于洛巴托法求定积分函数，调用格式如下：

[I，n]=quadl(fun，a，b，tol)

功能同函数quad()。该函数可以更精确地求出定积分的值，且一般情况下函数调用的步数明显小于quad()函数，从而保证能以更高效率求出所需的定积分值。tol缺省时为。

对于以表格形式定义的函数，MATLAB计算定积分函数为trapz()，格式为

St=trapz(x，y)

采用梯形法求函数y关于自变量x的积分。2.二重积分

采用递推自适应辛普森法计算二重dblquad()和三重triplequad()定积分的调用格式分别如下：

S=dblquad(fun，xmin，xmax，ymin，ymax，tol)

S=triplequad(fun，xmin，xmax，ymin，ymax，zmin，zmax，tol)

tol缺省时，积分的绝对精度为。【例4-33】求的结果。

(1)编写文件名为fin的函数文件。(2)调用dblquad函数求二重定积分。编写文件名为exm4_35的脚本文件：4.6初值问题的常微分方程数值解

凡含有参数、未知函数和未知函数导数(或微分)的方程，称为微分方程。未知函数是一元函数的微分方程称做常微分方程(OrdinaryDifferentialEquation，ODE)。许多物理模型都可以表示为常微分方程的形式。常微分方程的求解问题可以分为初值问题和边值问题。4.6.1ode()指令的执行机理

ode()指令是MATLAB中针对一阶微分方程或方程组的初值问题编写的指令。

对于微分方程和初值条件为：

的n阶常微分方程的初值问题，可将其转化为一阶常微分方程的初值问题，再应用ode()指令求解。具体步骤为：

(1)将n阶常微分方程变换为一阶常微分方程组，令则原n阶常微分方程以及初值问题变换为

和(2)将一阶常微分方程组编写成函数文件，并执行ode()

指令。

ode()指令主要采用龙格-库塔(Runge-Kutta)数值积分法。常用的ode()指令见表4-2。4.6.2ode()指令求解微分方程示例

以ode45()为例，介绍ode()指令的调用格式：

[X，Y]=ode45(‘DY’，xspan，Y0)

计算常微分方程DY的数值解，其中xspan为一维行向量，可以取两种形式：

(1)当时，表示微分方程的积分上、下限分别为x0和xf；

(2)当时，计算出这些离散点处的微分方程的数值解。Y0为微分方程的初值向量。

【例4-34】求微分方程组在区间[012]内

的数值解，初始条件为。

(1)编写文件名为ydot1关于微分方程组的函数文件：

functiondy=ydot1(x，y)

dy=[y(2)*y(3)；-y(1)*y(3)；2*y(1)*y(2)]；

(2)执行ode()指令求解该方程组的数值解，并图示。编写文件名为exm4_34的脚本文件：图4.4图示微分方程组数值解【例4-35】求解描述振荡器工作的范德波(VandePol)微分方程在初始条件

下在范围内的数值解，并分别绘制其解的曲线。参数。(1)将二阶微分方程减阶为一阶微分方程组。令

，则原方程演变为一阶微分方程组：图4.5图示范德波(VandePol)微分方程的数值解5.多项式部分分式展开

执行部分分式展开的函数是residue()，格式为

[r，p，k]=residue(b，a)

[b，a]=residue(r，p，k)

对多项式进行部分分式的展开和其逆运算。a为分式的分母系数向量，b为分式的分子系数向量。p为极点(Pole)，r为留数(Residue)，k(x)为直项(Directterm)。对于多项式b(x)和不含重根的n阶多项式a(x)之比，满足：

如果a(x)有m阶重根，则对应的部分分式可以写成：

【例4-24】对有理多项式

进行部分分式展开。4.3多项式插值和拟合

4.3