ch5MATLAB符号计算选读.ppt

* ? 采用极限方法也可以用来求函数的导数： 【例5.23】用求极限的方法求函数sin(x)的导数。 syms x t limit((sin(x+t)-sin(x))/t,t,0) ans = cos(x) * ? 2)多元函数的极限（这里主要指的是二元函数） 多元函数极限一般可表示为 在MATLAB中，同样是通过调用函数limit来实现，调用的语法格式如下： (1) limit(limit(f, x,x0 ), y, y0) (2) limit(limit(f, y,y0 ), x,x0 ) 这里f为函数表达式，一般为符号表达式，x,y是变量，x0,y0 是极限点，如果x0 或y0 不是确定的值，而是另外一个变量的函数，如x-g(y) ，则顺序不能改变。 * ? 【例5.24】求极限 syms x y a f=exp(-1/(x^2+y^2))*(sin(x))^2/x^2*(1+1/y^2)^(x+a^2*y^2) f = exp(-1/(x^2+y^2))*sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2) limit(limit(f,x,1/sqrt(y)),y,inf) ans = exp(a^2) * ? 5.4.2符号微分 函数diff是用来求符号表达式的微分。 1）一元函数的微分 语法： diff(f) %求f对自由变量的一阶微分 diff(f,t) %求f对符号变量t的一阶微分 diff(f,n) %求f对自由变量的n阶微分 diff(f,t,n) %求f对符号变量t的n阶微分 * ? 【例5.25】已知f(x)＝ax2+bx+c，求f(x)的微分。 f=sym(a*x^2+b*x+c) f = a*x^2 + b*x + c diff(f) %对默认自由变量x求一阶微分 ans = b + 2*a*x syms x a b c t limit((a*(x+t)^2+b*(x+t)+c-a*x^2-b*x-c)/t,t,0) ans = 2*a*x+b * ? diff(f,a) ans = x^2 diff(f,2) %对默认自由变量x求二阶微分 ans = 2*a diff(f,a,2) %对符号变量a求二阶微分 ans = 0 diff(f,3) %对默认自由变量x求三阶微分 ans = 0 * ? 微分函数diff也可以用于符号矩阵，其结果是对矩阵的每一个元素进行微分运算。 【例5.26】对符号矩阵 求微分。 clear syms x t A=[2*x t^2;t*cos(x) exp(x)] %创建符号矩阵 A = [ 2*x, t^2 ] [ t*cos(x), exp(x) ] * ? diff(A) %对默认自由变量x求一阶微分 ans = [ 2, 0 ] [ -t*sin(x), exp(x)] diff(A,t) %对符号变量t求一阶微分 ans = [ 0, 2*t] [ cos(x), 0] diff(A,2) %对默认自由变量x求二阶微分 ans = [ 0, 0 ] [ -t*cos(x), exp(x) ] * ? 2）多元函数的微分（这里主要指的是二元函数） 语法： diff(diff(f,x,m),y,n) %先对符号变量x求m阶微分, 再对符号变量y求n阶微分 diff(diff(f,y,n),x,m) %先对符号变量y求n阶微分, 再对符号变量x求m阶微分的值 syms x y f=(x^2-2*x)*exp(-x^2-y^2-x*y) f = -exp(- x^2 - x*y - y^2)*(- x^2 + 2*x) * ? fx=diff(f,x) fx = exp(- x^2 - x*y - y^2)*(2*x - 2) + exp(- x^2 - x*y - y^2)*(- x^2 + 2*x)*(