[理学]3、MATLAB符号计算.ppt

limit函数的调用格式为： limit(f,x,a) limit函数的另一种功能是求单边极限，其调用格式为： limit(f,x,a,right) 或 limit(f,x,a,left) 例: 求极限。 在MATLAB命令窗口，输入命令： syms a m x; f=(x^(1/m)-a^(1/m))/(x-a); limit(f,x,a) %求极限(1) f=(sin(a+x)-sin(a-x))/x; limit(f,x,0) %求极限(2) limit(f,inf) %求f函数在x→∞(包括+∞和-∞)处的极限 limit(f,x,inf,left) %求极限(3) g=(sqrt(x)-sqrt(a)-sqrt(x-a))/sqrt(x*x-a*a); limit(g,x,a, right) %求极限(4) MATLAB中的求导的函数为：diff(f,x,n) diff函数求函数f对变量x的n阶导数。参数x的用法同求极限函数limit，可以缺省，缺省值与limit相同，n的缺省值是1。 例：求函数的导数。 命令如下： syms a b t x y z; f=sqrt(1+exp(x)); % diff(f) %求(1) 。未指定求导变量和阶数，按缺省规则处理 f=x*cos(x); diff(f,x,2) %求(2) 。求f对x的二阶导数 diff(f,x,3) %求(3) 。求f对x的三阶导数 f1=a*cos(t);f2=b*sin(t); diff(f2)/diff(f1) %求(4)。按参数方程求导公式求y对x的导数 (diff(f1)*diff(f2,2)-diff(f1,2)*diff(f2))/(diff(f1))^3 %求(3)。求y对x的二阶导数 f=x*exp(y)/y^2; diff(f,x) %求(5)。z对x的偏导数 diff(f,y) %求(6)。z对y的偏导数 f=x^2+y^2+z^2-a^2; zx=-diff(f,x)/diff(f,z) %求(7)。按隐函数求导公式求z对x的偏导数 zy=-diff(f,y)/diff(f,z) %求(8)。按隐函数求导公式求z对y的偏导数 3.2 符号函数的定积分 定积分在实际工作中有广泛的应用。在MATLAB中，定积分的计算使用函数： int(f,x,a,b) 例：求定积分。 命令如下： x=sym(x);t=sym(t); int(abs(1-x),1,2) %求定积分(1) f=1/(1+x^2); int(f,-inf,inf) %求定积分(2) int(4*t*x,x,2,sin(t)) %求定积分(3) f=x^3/(x-1)^100; I=int(f,2,3) %用符号积分的方法求定积分(4) double(I) %将上述符号结果转换为数值 例：求级数之和。 命令如下： n=sym(n); s1=symsum(1/n^2,n,1,inf) %求s1 s2=symsum((-1)^(n+1)/n,1,inf) %求s2。未指定求和变量，缺省为n s3=symsum(n*x^n,n,1,inf) %求s3。此处的求和变量n不能省略。 s4=symsum(n^2,1,100) %求s4。计算有限级数的和 MATLAB中提供了将函数展开为幂级数的函数taylor，其调用格式为： taylor(f,v,n,a) 5.1. 傅立叶(Fo