第8章 matlab数值微分与积分.ppt

第8章 Matlab数值微分和积分;本章知识要点;数值微分、数值积分在化工计算中的作用;单变量函数的极限 格式1： L= limit( fun, x, x0) 格式2： L= limit( fun, x, x0, ‘left’ 或 ‘right’);例: 试求解极限问题 syms x a b; L=limit(‘x*(1+a/x)^x*sin(b/x)’,x,inf) L = exp(a)*b 例:求解单边极限问题 syms x; limit((exp(x^3)-1)/(1-cos(sqrt(x-sin(x)))),x,0,right) ans = 12;多变量函数的极限： 格式： L1=limit(limit(f,x,x0),y,y0) 或 L1=limit(limit(f,y,y0), x,x0) 如果x0 或y0不是确定的值，而是另一个变量的函数，如x-g(y)，则上述的极限求取顺序不能交换。;例：求出二元函数极限值 syms x y a; f=exp(-1/(y^2+x^2)) … *sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2); L=limit(limit(f,x,1/sqrt(y)),y,inf) L = exp(a^2);8.2 函数导数的解析解;syms x; f1=diff(sin(x)/x) 得结果: f1=cos(x)/x-sin(x)/x^2 ;多元函数的偏导： 格式： f=diff(diff(f,x,m),y,n) 或 f=diff(diff(f,y,n),x,m) 例： 求偏导数 syms x y; z=(x^2-2*x)*exp(-x^2-y^2-x*y); zx=diff(z,x) zx =(2*x-2)*exp(-x^2-y^2-x*y)+(x^2-2*x)*(-2*x-y)*exp(-x^2-y^2-x*y); zy=diff(z,y) zy = (x^2-2*x)*(-2*y-x)*exp(-x^2-y^2-x*y) 直接绘制三维曲面 [x,y]=meshgrid(-3:.2:3,-2:.2:2); z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y); surf(x,y,z), axis([-3 3 -2 2 -0.7 1.5]) ;隐函数的偏导数： 格式：F=diff(f,xj)/diff(f,xi) ;例： z=f(x,y)=-x2-y2-x*y; syms x y; f=x^2-y^2-x*y; diff(f,x)/diff(f,y) 或者： syms x y; diff(x^2-y^2-x*y,x)/diff(x^2-y^2-x*y,y) ans = (-2*x-y)/(-2*y-x) ;对于列表型函数往往需要用数值方法计算函数的微分 数值微分的基本方法 差分 利用插值（拟合）多项式求微分 利用三次样条插值（拟合）函数求微分 数值微分可以放大误差，应谨慎使用;函数diff;Matlab数值微分实现方法;例题;;8.2 数值积分;数值积分;梯形法数值积分：trapz(); 格式： S=trapz(x,y) 例： x1=[0:pi/30:pi]; y=[sin(x1) cos(x1) sin(x1/2)]; S1=trapz(x1,y) S1 = 1.9982 0.0000 1.9995;等同于： x=[0:pi/30:pi]; y1=sin(x); y2=cos(x) ; y3=sin(x1/2); S1=trapz(x,y1) ; S2=trapz(x,y2) ; S3=trapz(x,y3) ; S1, S2, S3 ;调用格式： q=quad(’fun‘,a,b) q=quad(’fun‘,a,b,tol) 输入参数： fun 被积函数。在定义fun时，被积函数表达 式必须是向量形式，即表达式必须使用 点运算符(.*、./和.^)以支持向量 a,b 即积分限[a,b] tol 绝对误差限，默认值为1.e-6 输出参数： q 积分结果 自适应Lobatto法数值积分：quadl() 调用格式同quad; 例：; 用 inline 函数定义被积函数： f=inline(2/sqrt(pi)*exp(-x.^2),x); y=quad(f,0,1.5) y = 0.9661 y=quad(2/sqrt(pi)*e