[工学]第六讲 MATLAB数值计算.pptVIP

例5.41 用不同的方法求函数f(x)的数值导数，并在同一个坐标系中做出f(x)的图象。 程序如下： f=inline(sqrt(x.^3+2*x.^2-x+12)+(x+5).^(1/6)+5*x+2); g=inline((3*x.^2+4*x-1)./sqrt(x.^3+2*x.^2-x+12)/2+1/6./(x+5).^(5/6)+5); x=-3:0.01:3; p=polyfit(x,f(x),5); %用5次多项式p拟合f(x) dp=polyder(p); %对拟合多项式p求导数dp dpx=polyval(dp,x); %求dp在假设点的函数值 dx=diff(f([x,3.01]))/0.01; %直接对f(x)求数值导数 gx=g(x); %求函数f的导函数g在假设点的导数 plot(x,dpx,x,dx,g.,x,gx,r-); %作图 5.6.2 数值积分 (1)被积函数是一个解析式 quad, quadl, odblquad q1 = quad (fun,a,b,tol,trace,p1,p2…) q2 = quadl (fun,a,b,tol,trace,p1,p2…) q3 = dblquad (fun,xmin,xmax,ymin,ymax,tol,method) 求被积函数fun在[a,b]上的定积分，tol是计算精度，缺省值是1.0e-6。trace非0时，显示中间计算结果。注意，调用quad函数时，先要建立一个描述被积函数f(x)的函数文件或语句函数。当被积函数fun含有一个以上的变量时，quad函数的调用格式为：quad(f,a,b,tol,trace,p1,p2) q1,q2,q3分别为自适应递推辛普生(simpson)方法、自适应Lobatto方法、二重积分结果。 注：quad8（自适应递推牛顿-柯西（Newton-Cotes)法）已废弃，推荐用quadl 例5.42 求积分 Q1=quad(sqrt(4*cos(2*t).^2+sin(t).^2+1), 0, 3*pi) %Simpson方法计算结果 Q2=quadl(sqrt(4*cos(2*t).^2+sin(t).^2+1), 0, 3*pi) %Lobatto方法计算结果 例5.43 求积分 Q3=dblquad(y*sin(x)+x*cos(y), 3*pi, 4*pi, pi, 2*pi) 例5.44 用trapz函数计算积分。 在MATLAB命令窗口，输入命令： X=0:0.01:1;Y=exp(-X.^2); trapz(X,Y)；%梯形法求数值积分 (2)被积函数由一个表格定义 MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。其中向量X、Y定义函数关系Y=f(X)。 (3)二重积分 例5.44计算二重积分。 建立一个函数文件fixy.m： function f=f(x,y) f=exp(-x.^2-y.^2); return 建立一个命令文件ftxy1.m： for i=1:20 int2(i)=quad(fixy,0,1,[],[],x(i)); %在二维函数fixy中以x=x(i)代入并对y积分。 end 在MATLAB命令窗口，输入命令： x=linspace(0,1,20); ftxy1 trapz(x,int2) 实际上，MATLAB提供了计算二重积分的函数： dblquad(f,a,b,c,d,tol,trace) 该函数求f(x,y)在[a,b]×[c,d]区域上的二重积分。参数tol，trace的用法与函数quad完全相同。 如果直接使用这里介绍的二重积分函数dblquad来求解本例就非常简单，命令如下： g=inline(exp(-x.^2-y.^2)); dblquad(g,0,1,0,1) %直接调用二重积分函数求解 5.7 常微分方程的数值求解 基于龙格－库塔法，MATLAB提供了求常微分方程数值解的函数，一般调用格式为： [X,Y]=ode23(f,[x0,xn],y0) [X,Y]=ode45(f,[x0,xn],y0) 其中X、Y是两个向量，X对应自变量x在求解区间[x0，xn]的一组采样点，其采样密度是自适应的，无需指定；Y是与X对应的一组解，f是一个函数，[x0，xn]代表自变量的求解区间，y0=y(x0)，由方程的初值给定。函数在求解区间[x0，xn]内，自动设立采样点向量X，并求出解函数y在采样点X处的样本值。 例5.45 求微分方程