MATLAB第 八讲 符号运算(续)学科知识.pptVIP

积分运算;2.计算积分; (1) x=sym(x); f=(3-x^2)^3; int(f) %求不定积分(1) (2) f=sin(x)^2; int(f) %求不定积分(2) g=simple(ans) %调用simple函数对结果化简;2。 符号函数的定积分：定积分在实际工作中有广泛的应用。在MATLAB中，定积分的计算使用函数：int(f,x,a,b) 例 求定积分。 (1) x=sym(x);t=sym(t); int(abs(1-x),1,2) %求定积分(1) (2) f=1/(1+x^2); int(f,-inf,inf) %求广义定积分(2) ; (3) f=x^3/(x-1)^10; I=int(f,2,3) %用符号积分的方法求定积分(3) double(I) %将上述符号结果转换为数值 (4) int(4*x/t,t,2,sin(x)) %求定积分(4) ;3. 求椭球的体积。 命令如下： syms a b c z; f=pi*a*b*(c^2-z^2)/c^2; V=int(f,z,-c,c) V = 4/3*pi*a*b*c; 例 求空间曲线c从点(0，0，0)到点(3，3，2)的长度。???曲线c的长度是曲线一型积分问题。按公式可转化为定积分问题。 命令如下： syms t; x=3*t;y=3*t^2;z=2*t^3; f=diff([x,y,z],t) %求x,y,z对参数t的导数 g=sqrt(f*f) %计算一型积分公式中的根式部分 l=int(g,t,0,1) %计算曲线c的长度;; 例1 求函数y= 的傅立叶变换及其逆变换。 syms x t; y=abs(x); Ft=fourier(y,x,t) %求y的傅立叶变换 fx=ifourier(Ft,t,x) %求Ft的傅立叶逆变换 Ft = -2/t^2 fx = -transform::fourier(1/t^2, t, x)/pi; 2. 拉普拉斯(Laplace)变换 逆变换： 在MATLAB中，进行拉普拉斯变换的函数是： laplace(fx,x,t) 求函数f(x)的拉普拉斯像函数F(t)。 ilaplace(Fw,t,x) 求拉普拉斯像函数F(t)的原函数f(x)。; 例 计算y=x2的拉普拉斯变换及其逆变换. 命令如下： x=sym(x); t=sym(t); y=x^2; Ft=laplace(y,x,t) %对函数y进行拉普拉斯变换 fx=ilaplace(Ft,t,x) %对函数Ft进行拉普拉斯逆变换 Ft = 2/t^3 fx = x^2; 3. Z变换 逆变换： 对数列f(n)进行z变换的MATLAB函数是： ztrans(fn,n,z) 求fn的Z变换像函数F(z) iztrans(Fz,z,n) 求Fz的z变换原函数f(n) ; 在前面讨论过求和函数sum，sum处理的级数是有穷级数。对于无穷级数求和，sum是无能无力的。求无穷级数的和需要使用符合表达式求和函数symsum。 1. 级数的符号求和 级数符号求和函数symsum，调用格式为： symsum(a,n,n0,nn) ; 例1 求级数之和。 (1) 常数项级数 n=sym(n); s1=symsum(1/n^2,n,1,inf) %求s1 s1 = pi^2/6 (2) s2=symsum((-1)^(n+1)/n,1,inf) %求s2。未指定求和变量，缺省为n s2 = log(2); (3) 函数项级数 s3=symsum(n*x^n,n,1,inf) %求s3。此处的求和变量n不能省略。 s3 = piecewise([abs(x) 1, x/(x - 1)^2]) (4) s4=symsum(n^2,1,100) %求s4。计算有限级数的和 s4 = 338350 (5) s5=symsum(sin(n*pi/6),1,inf)