上机习题3().docVIP

实验三（2） MATLAB7.0符号运算（ = 2 \* ROMAN II） 实验目的： ① 掌握MATLAB7.0subs、simple、simplify等函数的应用； ② 掌握利用MATLAB7.0计算极限和级数的方法，计算复合函数的方法； = 3 \* GB3 ③了解利用MATLAB7.0计算taylor级数、Fourier变换和Laplace变换的方法。 实验要求:给出程序和实验结果。 实验内容： 一、指出下面各条语句的输出结果 f=sym(2*u); subs(f,u,2) f2=subs(f,u,u+2) a=3; subs(f2,u,a+2) subs(f2,u,a+2) syms x y f3=subs(f,u,x+y) subs(f3,[x,y],[1,2]) subs(f3,[x,y],[x+y,x+y]) 二、简化 三、计算 和 四、1 计算级数 及其前100项的部分和； 2 计算函数级数 五、(1)设求复合函数 f(g(z)); (2)设求复合函数 f(g(z))。 六、设a、b定义如下，试上机输出factor(a)和factor(b)的结果，并指出那个结果才是12345678901234567890的因式分解，为什么？ a=sym(‘12345678901234567890’); b=sym(12345678901234567890); 七、替换函数subs的使用。 (1)试指出下面程序中的f1、f2、f3、f4、f5的值。 syms a x; f=a*sin(x)+5; f1=subs(f,sin(x),sym(y)) f2=subs(f,{a,x},{2,sym(pi/3)}) f3=subs(f,{a,x},{2,pi/3}) f4=subs(subs(f,a,2),x,0:pi/6:pi) f5=subs(f,{a,x},{0:6,0:pi/6:pi}) (2)试指出下面替换的结果，并说明理由。 syms x y; f=2*x+y; subs(f,{x,y},{x+y,x-y}) 八、计算sin(t)在t=0处的taylor级数、Fourier变换和Laplace变换。 九、计算的和。 实验结果： 1、f=sym(2*u) f = 2*u subs(f,u,2) ans = 4 f2=subs(f,u,u+2) f2 = 2*(u+2) a=3 a = 3 subs(f2,u,a+2) ans = 14 subs(f2,u,a+2) ans = 2*((a+2)+2) syms x y f3=subs(f,u,x+y) f3 = 2*x+2*y subs(f3,[x,y],[1,2]) ans = 6 subs(f3,[x,y],[x+y,x+y]) ans = 6*x+4*y 2、syms x f=[x^(-3)+6*x^(-2)+12*x^(-1)+8]^(1/3); simplify(f) ans = ((2*x+1)^3/x^3)^(1/3) 3、syms x h L=[log(x+h)-log(x)]/h; limit(L,h,0) ans = 1/x syms x n M=(1-x/n)^n; limit(M,n,Inf) ans = exp(-x) 4、（1）、syms n S=1/(n^2); symsum(S,1,inf) ans = 1/6*pi^2 symsum(S,1,100) ans = 1589508694133037873112297928517553859702383498543709859889432834803818131090369901/972186144434381030589657976672623144161975583995746241782720354705517986165248000 （2）、syms x n S=x/(n^2); symsum(S,n,1,inf) ans = 1/6*x*pi^2 5、（1）syms x u y t f=(x^2)/(1+u) f = x^2/(1+u) g=cos(y+2*t) g = cos(y+2*t) compose(f,g,x,y) ans = cos(y+2*t)^2/(1+u) （2）、 syms x u y t f=(x^2)/(1+u); g