第3章 MATLAB符号计算-fanny.pptVIP

3.4.1 符号极限 假定符号表达式的极限存在，Symbolic Math Toolbox提供了直接求表达式极限的函数limit，函数limit的基本用法如表3.3所示。 表 达 式 函 数 格 式 说 明 limt(f) 对x求趋近于0的极限 limt(f,x,a) 对x求趋近于a的极限，当左、右极限不相等时极限不存在 limt(f,x,a, left) 对x求左趋近于a的极限 limt(f,x,a, right) 对x求右趋近于a的极限 【例3.15】 分别求1/x在0处从两边趋近、从左边趋近和从右边趋近的3个极限值。 f=sym(1/x) f = 1/x limit(f) %对x求趋近于0的极限 ans = NaN limit(f,x,0) %对x求趋近于0的极限 ans = NaN limit(f,x,0,left) %左趋近于0 ans = ?inf limit(f,x,0,right) %右趋近于0 ans = inf 程序分析：当左、右极限不相等时，表达式的极限不存在，为NaN。 采用极限方法也可以求函数的导数： 。 【例3.15续】 求函数cos(x)的导数。 syms t x limit((cos(x+t) ?cos(x))/t,t,0) ans = ?sin(x) 函数diff是用来求符号表达式的微分。 语法： diff(f) %求f对自由变量的一阶微分 diff(f,t) %求f对符号变量t的一阶微分 diff(f,n) %求f对自由变量的n阶微分 diff(f,t,n) %求f对符号变量t的n阶微分 【例3.16】 已知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 = 2*a*x+b diff(f,a) %对符号变量a求一阶微分 ans = x^2 diff(f,x,2) %对符号变量x求二阶微分 ans = 2*a diff(f,3) %对默认自由变量x求三阶微分 ans = 0 微分函数diff也可以用于符号矩阵，其结果是对矩阵的每一个元素进行微分运算。 【例3.16续】 对符号矩阵 求微分。 syms t x g=[2*x t^2;t*sin(x) exp(x)] %创建符号矩阵 g = [ 2*x, t^2] [ t*sin(x), exp(x)] diff(g) %对默认自由变量x求一阶微分 ans = [ 2, 0] [ t*cos(x), exp(x)] diff(g,t) %对符号变量t求一阶微分 ans = [ 0, 2*t] [ sin(x), 0] diff(g,2) %对默认自由变量x求二阶微分 ans = [ 0, 0] [?t*sin(x), exp(x)] 【例3.16续】 可以使用diff计算向量间元素的差值。 x1=0:0.5:2; y1=sin(x1) y1 = 0 0.4794 0.8415 0.9975 0.9093 diff(y1) %计算元素差 ans = 0.4794 0.3620 0.1560 ?0.0882 程序分析：计算出的差值比原来的向量少1列。 ? 积分分为定积分和不定积分。运用函数int可以求得符号表达式的积分，即找出一个符号表达式F，使得diff(F)=f，也可以说是求微分的逆运算。 语法： int(f, t) %求符号变量t的不定积分 int(f, t,a,b) %求符号变量t的积分 int(f, t, m, n) %求符号变量t的积分 说明：t为符号变量，若t省略则为默认自由变量；a和b为数值，[a,b]为积分区间；m和n为符号对象，[m,n]为积分区间。与符号微分相比，符号积分复杂得多。函数的积分有时可能不存在，即使存在，也可能由于限于很多条件，MATLAB无法顺利得出。当MATLAB不能找到积分时，它将给出警告提示并返回该函数的原表达式。 【例3.17】 求 和 的积分。 f=sym(cos(x)); int(f) %求不定积分 ans = si