数学积分模型实验.pptVIP

第六讲 积分实验 问题：我国第一颗人造地球卫星近地点 距地球表面439km,远地点距地球 表面2384km,地球半径为6371km, 求该卫星的轨道长度。 分析与建模 人造卫星轨道可视为平面上的椭圆，参数方程为 a,b 分别是长、短半轴，根据计算参数方程的弧长的的公式，椭圆长度可表为如下椭圆积分 求解分析 椭圆积分无法用解析方法计算。数据 a=6371+2384=8755,b=6371+439=6810 可以用梯形公式或辛普森公式计算。 求积分的MATLAB命令 trapz 梯形法积分； dblquad 矩形区域积分； quad 变步长数值积分； int 符号积分； quad8 高精度数值积分。 梯形积分法 例 2. 变步长数值积分 z=quad8(‘Fun’,A,B,Tol) Fun—被积函数的M函数名； A,B—下限，上限； Tol—精度，缺省值为1e-3. (eg4_2fun.m) 3. 重积分 (1) dblquad只能求矩形区域二重积分。 3. 重积分 矩形区域积分 z=dblquad(‘Fun’,a,b,c,d) Fun—被积函数的M函数名； a,b—x的下上限；c,d—y的下上限。 （2） M文件dblquad2.m给出的二重计算方法 用法 I=dblquad2(‘f_name’,a,b,’c_lo’,’d_hi’.m,n) ‘f_name’—为被积函数f(x,y)的字符串，x为标量，y为向量；a,b为x的下上限，‘c_lo’和‘d_hi’是y的下上上限函数c(x),d(x),都是x的标量函数，m,n分别为 x、y方向的等分数（缺省值为100)。例1 4. 符号积分 int(s) 符号表达式s的不定积分； int(s,v)符号表达式s关于变量v的不定积分; int(s,a,b)符号表达式s的定积分,a,b为下上限； int(s,v,a,b)符号表达式s关于变量v从a到b的定积分。 梯形公式程序（x5) 辛普森公式程序（x6f) function y=fun(x) a=6755;b=6810; y=sqrt(a^2*sin(x).^2+b^2*cos(x).^2); 在command窗口执行 l2=4*quad(x6f,0,pi/2,1e-6) 练习 一台重5400kg的汽车在以速度30m/s行使时突然熄火，设滑行方程为 x为滑行距离，v为车速，计算要滑行多长距离后，速度可降至15m/s. * * 实验目的 实验内容 2、学会用Matlab求数值积分. 1、学会用Matlab求符号积分. 1、求符号积分. 4、实验作业. 2、求数值积分. z=trapz(x,y) x是积分区间的离散向量；y是被积函数（与x同维）；z返回积分的近似值 z=dblquad(‘Fun’,a,b,c,d) Fun—被积函数f的M函数名。 a,b—x的下上限；c,d—y的下上限。 a=6755;b=6810;t=0:pi/10:pi/2; y=sqrt(a^2*sin(t).^2+b^2*cos(t).^2); l1=4*trapz(t,y)