MATLAB第3篇第3节.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); 在命令窗口键入以下程序：