积分与微分.doc.doc

一 求极限、导数和极值, 级数求和与Taylor展开 学习目的 学习MATLAB语言有关极限、导数、极值、极限与Taylor展开等命令，加深对微分学的一些基本概念及理论的理解，掌握数值运算和符号运算的基本方法及MATLAB语言的机器实现。 学习步骤与内容 一、微分、数值微分 若f(x)在x=a可导，设h0且足够小，向前差商、向后差商、中心差商分别为 二、二阶导数公式 三、求极限、导数的MATLAB命令 limit 符号求极限； taylor 符号Taylor展开； symsum符号级数求和； polyder 多项式求导； diff 数值微分或符号微分； gradient 数值梯度； 数值求导 dy=diff(x) 向量x的差分； Fx=gradient(F,x) ; 即F( (x)，其中x与F是同维向量。 [Fx,Fy]=gradient(F,x,y);即二元函数数值梯度 （Fx, Fy），F为m(n矩阵时，x,y 分别是n 维和 m维向量。 例　设y = x3,求y’(x)在x＝1, 1.1 ,1.2 ,1.3处的值。 clear;x=[1 1.1 1.2 1.3];y=x.^3; dy=diff(y)./diff(x) dy = 3.3100 3.9700 4.6900 得到 y’(1),y’(1.1)和y’(1.2)的近似值（向前差商）。若用梯度(中心差商）求解 dy=gradient(y,x) dy = 3.3100 3.6400 4.3300 4.6900 得到 y’(1),y’(1.1),y’(1.2)和y’(1.3)的近似值.第一和最后一个数据分别用前、后差商代替。误差较小： 3*x.^2 ans = 3.0000 3.6300 4.3200 5.0700 解析运算 limit(s,x,a) 表达式 s当x????a时的极限。 symsum(s,n,a,b) 表达式 s的通项当自变量 n由a到b的和。 diff(s,x,n) 表达式 s对x的n阶导数。 taylor(s,n,a) 表达式 s在a点Taylor展开到n次式。 例１　　sinx在０点展开到６次。 syms x; taylor(sin(x),6,0) ans = x-1/6*x^3+1/120*x^5 例２　计算　 syms n x y;limit((1+x/n)^n,n,inf) ans = exp(x) symsum((-1)^n*x^n/n,n,1,inf) ans = -log(1+x) diff(sin(x*y),x,2) ans = -sin(x*y)*y^2 练习： 二 求积分 学习目的 学习数值积分方法MATLAB语言有关积分的命令，加深对积分学的一些基本概念及理论的理解，掌握数值积分和符号积分运算的基本方法及MATLAB语言的机器实现。了解积分在计算面积、体积等问题中的运用。 步骤与内容 数学原理： 矩形法 一、求积分的MATLAB命令 trapz 梯形法积分； dblquad 矩形区域积分； quad 辛普森数值积分； int 符号积分； quad8 　高精度数值积分。 梯形积分法 z=trapz(x,y) x是积分区间的离散向量；y是被积函数（与x同维）；z返回积分的近似值 例 1 clear;x=-1:0.1:1; y=exp(-x.^2); trapz(x,y) ans = 1.4924 2. 变步长数值积分 z=quad(‘Fun’,A,B,Tol): 用辛普森（2阶）公式计算以fun.m文件命名的函数在区间（a,b）上的积分，自动选择步长，相对误差为tol,并输出积分值。 Fun—被积函数的M函数名； A,B—下限，上限； Tol—精度，缺省值为1e-3. 例 先写M函数：eg2_1fun.m function y=fun(x) y=exp(-x.^2); 然后在命令窗口用 z=quad(eg2_1fun,-1,1) z = 1.4936 重积分 dblquad只能求矩形区域二重积分。 z=dblquad(‘Fun’,a,b,c,d) Fun—被积函数f的M函数名。 a,b—x的下上限；c,d—y的下上限。 符号积分 int(s) 符号表达式s的不定积分； int(s,v)符号表达式s关于变量v的不定积分; int(s,a,b)符号表达式s的定积分,a,b为下上限； int(s,v,a,b)符号表达式s关于变量v从