3.6符号微积分.pptVIP

* 3.6 符号微积分 微积分是高等数学中最重要的基础内容之一，它被广泛地应用于许多的工程学科中。MATLAB的符号数学工具箱为我们提供了快速、简便地计算微积分的工具。MATLAB符号数学工具箱中的符号微积分包括符号极限、符号微分、符号积分等。 1 符号极限 极限是微积分学的基础。在MATLAB中，符号极限由函数“limit”来实现的。Limit函数的调用格式如下： 【例3.6.1】 syms x; limit((x-2)/(x^2-4)) ans = 1/2 Limit(F) Limit(F,a) 返回符号表达式F的默认变量趋向于a 时的极限 【例3.6.2】 syms x; limit((x-2)/(x^2-4),2) ans = 1/4 【例3.6.3】 syms x h; limit((sin(x+h)-sin(x))/h,h,0) ans = cos(x) Limit(F,x,a) Limit(F,x,a,’right’) Limit(F,x,a,’left’) 【例3.6.4】 syms x; limit(1/x,x,0,left) ans = -inf Limit(F,x,a,’right’) Limit(F,x,a,’left’) 【例3.6.4】 syms x; limit(2/(1+exp(1/x))+sin(x)/(x),x,0,left) ans = 3 %求 的左右极限 limit(2/(1+exp(1/x))+sin(x)/(x),x,0,right) ans = 1 【例3.6.5】 syms x a v; v=[(1+a/x)^x,exp(-x)]; limit(v,x,+inf) ans = [ exp(a), 0] 2 符号微分 微分是高等数学中最基础的内容之一。在MATLAB中，符号微分由函数diff(differentiation)来实现的。实际上，在数值运算一章中已介绍过diff函数了。Diff函数可同时计算数值微分与符号微分。当输入的参数是数值时，MATLAB能非常巧妙地对其进行数值微分；当输入的参数是符号字符串时，MATLAB同样能非常巧妙地对其进行符号微分。 Diff(f) 对findsym函数返回的独立变量求导数， f为符号表达式 Diff(f,’a’) 对a变量求导数，f为符号表达式 Diff(f,n) 对findsym函数返回的独立变量求n次导 数，f为符号表达式 Diff(f,’a’,n) 或Diff(f,n,’a’) 对变量a求n次导数， f为符号表达式 1）一元函数的导数 一元函数导数的调用格式如下： 【例3.6.6】 syms x; f=sym((x-1)^3/(x-1)); B=diff(f) B = 2*x-2 【例3.6.7】 syms x; f=sym((x-1)^3/(x-1)); B=diff(f,2) B = 2 同样地，函数Diff也可对符号矩阵进行运算。此时，它是对符号矩阵中的每个元素进行微分。 【例3.6.8】 A=sym([cos(x),sin(x);x^ 2+x+1 tan(x)]); B=diff(A) B = [ -sin(x), cos(x)] [ 2*x+1, 1+tan(x)^2] 2）多元函数的偏导数 在多元函数中，仿照单元函数的极限、可微的概念引入了偏导数。多元函数偏导数在非线性方程的求解和变分原理中有极其重要的应用，在MATLAB中，此问题的实现由函数jacobian完成。 jacobian(f,v) 计算数量或向量f对向量v的Jacobi 矩阵。注意当f为数量时，函数返 回f的梯度。 设 则 【例3.6.9】求各函数的Jacobi矩阵。 (1) 解：(1) x=sym([x]); y=sym([y]); z=sym([z]); jacobian([x^2+y^2;x^2-y^2],[x,y]) a