matlab第8章实例解析.pptVIP

实例解析【例8-1】复化求积示例（具体题目见书本例8-1）。解：编写如下语句：v=simple(dsolve(Dv=g-c_d/m*v^2,v(0)=0,t))%微分方程解析解f=eval([@(t,g,m,c_d),char(v)]);%将符号表达式写成匿名函数g=9.81;m=68.1;c_d=0.25;v_3=subs(v,{t,g,m,c_d},{3,g,m,c_d})%第3秒的下降速度x1=trape(f,0,3,[],g,m,c_d)%复化梯形公式求积x2=simpson(@(t)f(t,g,m,c_d),0,3)%复化辛普森公式求积x3=Cotes(@(t)f(t,g,m,c_d),0,3)%复化Cotes公式求积运行结果：

x1=41.9480x2=41.9481x3=41.9481【例8-2】自适应求积方法示例（具体题目见书本例8-9）。解：首先编写如下被积函数描述文件Dampedsinewave.m，程序如下：functionf=Dampedsinewave(t,xi)alpha=atan(-xi/sqrt(1-xi^2));%计算alphaf=exp(-xi*t).*cos(t*sqrt(1-xi^2)+alpha)/cos(alpha);%求函数值主程序：formatlongxi=0.1;I1=adapt_trape(@Dampedsinewave,0,20,[],xi)%自适应梯形求积I2=adapt_simp(@Dampedsinewave,0,20,[],xi)%自适应辛普森求积I3=adapt_Cotes(@Dampedsinewave,0,20,[],xi)%自适应Cotes求积I4=Romberg(@Dampedsinewave,0,20,[],xi)%Romberg求积运行结果：I1=0.302173980392894I2=0.302174200120611I3=0.302174464839673I4=0.302174214269305【例8-3】Gauss求积方法示例。解：编写如下程序求得Gauss求积数值解与解析解的误差：formatshortef=@(x)sin(x)./x;%定义被积函数form=1:4I1=Gauss_legendre(f,1,5,m);%高斯-勒让德公式求积分I2=Gauss_Lobatto(f,1,5,m);%Gauss_Lobatto公式求积err(1,m)=abs(I1-double(int(sin(x)/x,1,5)));%计算绝对误差err(2,m)=abs(I2-double(int(sin(x)/x,1,5)));%计算绝对误差enderr%显示误差结果输出结果：err=1.9400e-0064.8918e-0102.8755e-0143.0192e-0112.9093e-0066.9217e-0101.8661e-0092.3371e-009【例8-4】计算无穷积分无穷区间逼近f=@(x)exp(-x.^2);%定义被积函数I=quad_inf(f,0,inf)%求无穷积分变量替换f=@(x)exp(-(x./(1-x)).^2)./(1-x).^2;%定义被积函数I=simpson(f,eps,1-eps)%复化辛普森公式求积分运行结果：I=0.8862【例8-5】quadgk()函数示例。1、无穷积分：2、间断函数积分：3、震荡积分：4、复数积分：解：编写如下语句：formatlongf1=@(x)exp(-x.^2);%定义被积函数I1=quadgk(f1,0,inf)%quadgk函数求解无穷积分f2=@(x)x.^5.*exp(-x).*sin(x);%定义被积函数[I2,errbnd]=quadgk(f2,1,10,Waypoints,[25])%其中2,5为间断点f3=@(x)exp(x).*cos(1000*x);I3_quad=quad(f3,0,pi)%quad函数求解I3_quadgk=quadgk(f3,0,pi,MaxIntervalCount,1000)%quadgk函数求解i=sqrt(-1);f4=@(x)exp(-x.^2-i*x).*sin((7+2i)*x);%定义被积函数I4=qu