MATLAB应用第4章第4版 微积分问题的计算机求解.docVIP

第章 微积分问题的计算机求解 4.1　极限问题的解析解4.1.1 单变量函数的极限 其中， x0可以是一个确定的值，也可以是无穷大。 对某些函数来说，还可以定义单边极限（或称左右极限）， 或 前者表示x从左侧趋近于x0点，称为左极限，后者相应地称为右极限。 极限问题在MATLAB符号运算工具箱中可以使用1imit函数直接求出，该函数的调用格式为 1．limit(f,x,a)：求 2．limit(f,a)：求f(x)的极限值。符号函数f(x)的变量为函数findsym(f)确定的默认自变量a，即x→a 3．limit(f)：求f(x)的极限值。符号函数f(x)的变量为函数findsym(f)确定的默认自变量。没有指定变量的目标值，系统默认自变量趋近于0，即a＝0。 4．limit(f,x,a,’right’)：求。’right’表示变量x从右边趋近于a 5．limit(f,x,a,’left’)：求。’left’表示变量x从左边趋近于a 例1 求解极限问题。 syms x a b; f = x*(1+a/x) ^ x*sin(b/x); L = limit(f,x,inf) 例2．求 syms a m x f = (x ^ (1/m) - a ^ (1/m)) / (x-a) limit (f,x,a) 例3．求 syms x f = x * (sqrt(x ^ 2 + 1) - x) limit(f, x, inf, left) %为什么不right？ 例4．求 syms a m x f = (sqrt(x) - sqrt(a) + sqrt(x - a)) / sqrt(x * x - a * a) limit(f, x, a, right) 例5 求解单边极限问题 syms x; limit((exp(x ^ 3) - 1)/(1 - cos(sqrt(x - sin(x)))), x, 0, right) 例6 分别求出tanu函数关于π / 2点处的左右极限。 syms u; f = tan(u); L1 = limit(f, u, pi/2, left) L2 = limit(f, u, pi/2, right ) 4. 1. 2 多变量函数的极限 limit函数直接求解。假设有二元函数f (x, y)，若要求出二元函数的极限 则可以嵌套使用limit函数。例如， L1 = limit (limit (f, x, x0), y, y0) 或 L2—limit (limit (f, y, y0), x, x0) 如果x0或y0不是确定的值，而是另一个变量的函数，例如 x → g(y)，则上述的极限求取顺序不能交换。 例7 求出二元函数极限值 syms x y a; f = exp( -1 / (y ^ 2 + x ^ 2)) * sin(x) ^ 2 / x ^ 2 * (1 + 1 / y ^ 2 ) ^ (x + a ^ 2 * y ^ 2); L = limit(limit(f, x, 1 / sqrt(y)), y, inf) 4. 2 函数导数的解析解 4. 2. 1 函数的导数和高阶导数 diff(s)：没有指定变量和导数阶数，则系统按findsym函数指示的默认变量对符号表达式s求一阶导数。 diff(s,v)：以v为自变量，对符号表达式s求一阶导数。 diff(s,n)：按findsym函数指示的默认变量对符号表达式s求n阶导数，n为正整数。 diff(s,v,n)：以v为自变量，对符号表达式s求n阶导数。 例1．，求 syms x y=sqrt(1+exp(x)) diff(y) %求1。未指定求导变量和阶数，按默认规则处理 例2．，求y’’、 y’’’。 syms x y=x*cos(x) diff(y,x,2) %求2。求y对x的2阶导数 diff(y,x,3) %求2。求y对x的3阶导数 例3 给定函数，试求出 syms x; f = sin(x) / (x ^ 2+4 * x + 3); f4 = diff(f, x, 4) 4. 2. 2 隐函数的偏导数 ，则可以通过隐函数对它们的偏导数求出自变量之间的偏导数。具体可以用下面的公式求出： 例4．由定义，求， syms a x y z f=x^2+y^2+z^2-a^2 zx=-diff(f,x)/ diff(f,z) %求z对x的偏导数 zy=-diff(f,y)/ diff(f,z) %求z对y的偏导数 例5二元函数，求 syms x y; f = (x ^ 2 - 2 * x) * exp(- x ^ 2 -