实验7微积分基本运算.docVIP

实验7 微积分基本运算 一、实验目的 学会函数的极值、微分、积分的方法. 二、实验内容与要求 1. 函数的极限 格式： limit(F,x,a) %计算符号表达式F=F(x)当x→a时的极限值. limit(F,x,a,right) %计算符号函数F当x→a+时的右极限. limit(F,x,a,left) %计算符号函数F当x→a－时的左极限. 【例】 syms x a t h n; L1=limit((cos(x)-1)/x) %缺省状态下，计算当x→0时的极限值 L2 = limit(1/x^3,x,0,right) L3 = limit(1/x,x,0,left) L4 = limit((log(x+h)-log(x))/h,h,0) v = [(1+a/x)^x, exp(-x)]; L5 = limit(v,x,inf,left) L6 = limit((1+2/n)^(3*n),n,inf) 计算结果为： L1 = 0 L2 = inf L3 = -inf L4 = 1/x L5 = [ exp(a), 0] L6 = exp(6) 2. 求单变量函数的极值 fminbnd(Scalar bounded nonlinear function minimization. ) 格式： x=fmin (F,a,b)%计算区间a-b上函数F取最小值时的x值. 说明：在5.3及以上版本命令fmin已改为fminbnd，常用格式如下. x=fminbnd(F,a,b) %计算区间a-b上函数F取最小值时的x值. [x,fval]= fminbnd(F,a,b) %计算在区间a-b上函数F的最小值fval和对应的x值. 【例】 求函数f(x)= 2x3－6x2－18x+7在区间(-2,4)的极小值，并作图. f=inline(2*x.^3-6*x.^2-18*x+7); %建立内联函数f(x). [x,fval] = fminbnd(f, -2, 4); %求函数f的最小值和对应的x值. fplot(f,[-2,4]) 结果为： x = 3.0000 fval = -47.0000 如图1.12所示. 图1.12 注意：用inline建立的函数f，在fminbnd和fplot命令中不用加单引号，而用M函数文件建立的函数则要加单引号. 问题：如何求函数f的最大值？ 3. 函数的微分 格式：diff(S,‘v’,n) %对符号表达式S中指定的符号变量v计算S的n阶导数，在缺省状态下，v=findsym(S)，n=1. 【例】 syms x y t D1 = diff(sin(x^2)*y^2,2) %计算 D2 = diff(D1,y) %计算 D3 = diff(t^6,6) 计算结果为： D1= -4*sin(x^2)*x^2*y^2+2*cos(x^2)*y^2 D2= -8*sin(x^2)*x^2*y+4*cos(x^2)*y D3= 720 4. 函数的积分 quad法数值积分 quad(adaptive Simpson quadrature.求[面]积) 格式：s = quad( fun,a,b) %近似地从a到b计算函数fun的数值积分，误差为10－6. s = quad( fun,a,b,tol) %用指定的绝对误差tol代替缺省误差. s =quad8( fun,a,b,…) %用高精度进行计算，效率可能比quad更好. 说明：quad8命令在6.x版本用quadl代替. quadl(Lobatto quadrature.洛巴托求积) s=quadl( fun,a,b,…) 【例】 fun = inline(3*x.^2./(x.^3-2*x.^2+3)); %构造一函数 S1 = quad(fun,0,2) S2 = quadl(fun,0,2) 计算结果为： S1 = 3.7224 S2 = 3.7224 注意：用inline构造函数比用function构造函数简单；命令quadl 最后是字母l，不是数字1. 符号函数的积分 int(integral) 格式：R = int(S,v) %对符号表达式S中指定的符号变量v计算不定积分. R = int(S,v,a,b) %对表达式S中指定的符号变量v计算从a到b的定积分. 【例】 syms x z t alpha I1 = int(-2*x/(1+x^3)^2) I2 = int(x/(1