6--数值微积分、微分方程课题.pptVIP

数值微积分 数值微(差)分 数值积分 数值定积分 数值重积分 DX=diff(A,n,dim)：计算矩阵A的n阶差分，dim=1时(缺省状态)，按列计算差分；dim=2，按行计算差分。 例如： 求下列定积分 数值积分 求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。 Matlab 计算定积分函数介绍 1、梯形法：trapz trapz 举例 抛物线法 quad 举例 建立函数法求定积分 例 求定积分： (1) 建立被积函数文件fesin.m。 function f=fesin(x) f=exp(-0.5*x).*sin(x+pi/6); (2) 调用数值积分函数quad求定积分。 [S,n]=quad(fesin,0,3*pi) 洛巴托数值积分 调用格式： Q = quadl (FUN,a,b,TOL) 例：用函数trapz、quad和quadl分别数值积分将计算结果精确值进行比较. 变限积分求导举例 例：求 综合应用举例 例：求下列函数所围成的平面区域D的面积S，并作出图形 ： 二重积分的计算 抛物线法计算二重积分： dblquad dblquad 举例 dblquad 举例 dblquad 举例 数值实验 非矩形区域重积分dblquad2.m function S=dblquad2(fun,a,b,clo,dhi,n) if nargin6,n=100;end x=linspace(a,b,n+1); S=0; for i=1:n S=S+dblquad(fun,x(i),x(i+1),feval(clo, … (x(i)+x(i+1))/2),feval(dhi,(x(i)+x(i+1))/2)); end 解答 调用函数dblquad2(@fun1,-1,1,@fun2,@fun3) 三重积分 求得三元函数Fun(x,y,z)的重积分, 格式类似dblquad。 常微分方程(组)求解 微分方程的解 代入微分方程后能使方程两端恒等的函数称为微分方程的解。 通解：若微分方程的解中所含任意常数的个数恰好等于微分方程的阶数，则称此解为微分方程的通解。 特解：通解中的任意常数得到确定后，称其为微分方程的特解。 数值解：求得的特解函数是用“列表法”表示的函数（非解析式）。 符号解(解析解)：解函数以解析式表示。 定解条件 定解条件：确定通解中的常数，使其成为特解所需的额外附加条件称为定解条件。 初始条件：定解条件由系统在某一瞬间所处的状态给出时（t0时刻的函数值及各阶导数值），称其为初始条件。 定解问题：带定解条件的微分方程问题。 初值问题：带初始条件的微分方程问题。 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Matlab函数数值求解 Matlab提供的ODE求解器 参数说明及举例 数值求解举例 数值求解举例 再说定义函数文件 再说定义函数文件 举例说明 举例说明（单个方程） 举例说明（方程组） 思考 说明 例 高阶常微分方程 思考 对于高阶微分方程组的数值求解。应用如何处理？ 综合举例 导弹追踪问题 设位于坐标原点的甲舰向位于x轴上的A(1,0)处的乙舰发射导弹，导弹始终对准乙舰。如果乙舰以最大的速度 沿平行于y轴的直线行驶，导弹的速度是5 ，求导弹运行的曲线。当乙舰行驶多远时，导弹将它击中。 综合举例 设任意时刻t，乙舰的坐标为 ， 导弹的坐标为 。导弹速度恒为 ，从原点射出，且速度是横向和纵向距离的一阶导数的矢量和，则有： (1) 由于导弹头始终对准乙舰，导弹头的速度向量平行于乙舰与导弹头位置的差向量，即： (2) 综合举例 故有（1）（2）式可得 （3） 解（3）得到 （4） 由题目知乙舰以最大速度 沿直线 运动，为便于