MATLAB第3章第3节I.pptVIP

几种求积指令的比较 计算积分 （很显然这个积分的值是2） 编写M文件如下： x=linspace(0,pi,50); format long y=sin(x); T1=cumsum(y(1:49))*pi/(49); %左矩形法 T1=T1(49) T2=trapz(x,y); %梯形法 3.3.4 用Matlab作数值积分 3.3.4 用Matlab作数值积分 T3=quad(sin,0,pi); %辛普森法 % r1 r2 r3是几种方法的绝对误差 r1=abs(T1-2); r2=abs(T2-2); r3=abs(T3-2); 计算结果： r1 = 6.851506759375514e-004 r2= 6.851506759371073e-004 r3= 3.601569042999131e-009 3.3.5 实验 利用梯形公式和辛普森公式计算3.3.1中 卫星轨道长度问题，编写M文件如下： function y=x5(t) a=8755;b=6810; y=sqrt(a^2*sin(t).^2+b^2*cos(t).^2); 在命令窗口键入以下程序： t=0:pi/10:pi/2; y1=x5(t); l1=4*trapz(t,y1); 定义一个函数 调用这个函数 3.3.5 实验 l2=4*quad(x5,0,pi/2,1e-6); format long 用长格式输出结果： l1=4.908996526785276e+004 l2= 4.908996526863898e+004  Harbin University of Science and Technology Experiments in Mathematics  Harbin University of Science and Technology Experiments in Mathematics  Harbin University of Science and Technology Experiments in Mathematics 3.3 数值积分 3.3.1 卫星轨道长度问题 3.3.2 公式的导出 3.3.3 误差估计和收敛性 3.3.4 用Matlab作数值积分 3.3.5 实验 3.3.1 卫星轨道长度问题 问题： 人造地球卫星轨道可视为平面上的椭圆。我国第一颗人造地球卫星近地点距地球表面439km ，远地点距地球表面2384km，地球半径为 6371km。 求: 该卫星的轨道长度。 3.3.1 卫星轨道长度问题 模型 a,b分别是长半轴和短半轴; 焦距为c， 地球半径为 r ， 近地点和远地点与地球表面的距离分别是 和 。 图3.4 卫星轨道的示意图 3.3.1 卫星轨道长度问题 椭圆的参数方程为 弧长的公式 椭圆长度 椭圆积分 无法用解析方法计算, 讨论用数值方法来计算。 3.3.2 公式的导出 给定函数f(x)，关于积分，有Newton-Leibniz 公式 但是，在下列情况下， 函数在离散点处给出； 被积函数的原函数无法用初等函数表示； 被积函数的原函数虽有初等函数表达式，但过于复杂； 就必须借助数值方法来求函数的积分 3.3.2 公式的导出 用数值方法近似地求一个函数 在区间(a,b)上的 定积分的基本思路，可以归结到定积分的定义 (3-4) 3.3.2 公式的导出 取等距步长，当n充分大时， 就是 的数值积 分， 是第k小区间中x的取值，显然， 取值不同 ，数值积分 的结果就不同。 这种做法相当于用相对简单的阶梯函数 (k=1,…,n) 代替 作积分。 实际上各种不同的数值积分方法就在于，研究用什么样的简单函数代替 ，使得既能保证结果有一定的精度，计算量又小。 3.3.2 公式的导出 牛顿-柯特斯(Newton-Cote’