6.3符号微积分.pptVIP

例6-36 求 */17 例6-37求 在区间[-1，1]*[-1，1]上的二重积分。 syms x y f=int(int(exp(-x^2/3)*sin(x^2+2*y),x,-1,1),y,-1,1) f = ((-pi)^(1/2)*sin(2)*erf((i + 1/3)^(1/2)))/(2*(i + 1/3)^(1/2)) - ((-pi)^(1/2)*sin(2)*erf((1/3 - i)^(1/2)))/(2*(1/3 - i)^(1/2)) double(f) ans = 0.4658 syms x y z; F2=int (int (int(x^2+y^2+z^2,z,sqrt(x*y),x^2*y),y,sqrt(x),x^2), x,1,2) VF2=vpa(F2) 执行结果如下： VF2 = 224.921535733312 */17 补充 求积分 */17 例6-38 求 syms x y int(int(sin(pi*(x^2+y^2)),y,-sqrt(1-x^2),sqrt(1-x^2)),x,-1,1) 执行结果如下： Warning: Explicit integral could not be found. In sym.int at 58 ans = int(2^(1/2)*cos(pi*x^2)*FresnelS((1-x^2)^(1/2)*2^(1/2))+2^(1/2)*sin(pi*x^2)*FresnelC((1-x^2)^(1/2)*2^(1/2)),x = -1 .. 1) 由执行结果可以看出，结果中仍带有int，表明Matlab没有求出这一积分的值，因此采用极坐标可化为二次积分，程序如下： syms a r int(int(r*sin(pi*r^2),r,0,1),a,0,2*pi) 执行结果如下： ans =2 */17 例6-39 求曲线积分 ， 其中 为曲线 在第一象限内的一段。 syms t; int(cos(t)*sin(t),0,pi/2) 执行结果如下： ans = 1/2 * * 第6章 符号计算 (3) 6.3 符号微积分 6.3.1 符号极限 求函数极限的函数：limit limit(f) 求符号函数f,默认自变量-0的极限值。 limit(f,a) 求符号函数f在默认自变量-a的极限值。 limit(f,x,a) 计算当自变量x趋近于常数a时，f（x）函数的极限值。 */17 limit(f,x,a,right) 功能：求符号函数f的极限值。 right表示变量x从右边趋近于a。 limit(f,x,a,left) 功能：求符号函数f的极限值。 left表示变量x从左边趋近于a。 */17 例6-26 求两个极限： syms x limit(sin(x)/x,x,0) ans = 1 limit((1-1/x)^x,x,inf) ans = 1/exp(1) 例6-27定义法求函数f(x)=sin(x)的导数。 操作步骤： syms t x; limit((sin(x+t)-sin(x))/t,t,0) ans = cos(x) */17 和 6.3.2 符号求和 级数求和运算函数：symsum symsum(s,x,a,b) 功能：计算符号表达式s的级数和。 计算符号表达式s的对于变量x从a到b的和。省略x表示对默认自变量求和，省略a，b时，如果变量为k表示s对变量k求从0到k-1的和。 */17 */17 例6-28 计算 syms x k t symsum(k) ans = k^2/2 - k/2 symsum(x^k/sym(k!),k,0,inf) ans = exp(x) symsum([t,k^2]) %数组的自变量为t ans = [ (t*(t - 1))/2, k^2*t] , , 6.3.3 符号微分 求符号表达式的微分的函数：diff diff(f) 功能：求f对自由变量的一阶微分 diff(f,t) 功能：求f对符号变量t的一阶微分 diff(f,n) 功能：求f对自由变量的n阶微分 diff(f,t,n) 功能：求f对符号变量t的n阶微分 */17 diff函数的应用: 求一元函数的导数 求参数方程所确定的函数的导数 求多元函数的偏导数 */17 */17 例6-29 以验证罗必塔法则。 分析：罗必塔法则是指在求极限时，可