09实验三及答案.docVIP

实验三、 1、设x为符号变量，令y=x4+2x2+1，z=y2，w=sin(z)，试求w关于x的符号表达式，并求当x=5时，w的小数点之后具有1位、2位、10位、20位和50位有效数字的数值解。 解： syms x y=x^4+2*x^2+1 y = x^4+2*x^2+1 z=y^2 z = (x^4+2*x^2+1)^2 w=sin(z) w = sin((x^4+2*x^2+1)^2) u=subs(w,5) u = -0.0673 vpa(u,1) ans = -.7e-1 vpa(u,2) ans = -.67e-1 vpa(u,10) ans = -.6734017262e-1 vpa(u,20) ans = -.67340172619870686255e-1 vpa(u,50) ans = -.67340172619870686254728298081317916512489318847656e-1 2、求级数的级数和以及级数的级数和。 解： syms n k r1=symsum((-1)^(n-1)/(n*(n+1)),n,1,inf) r1 = -1+2*log(2) r2=symsum((-1)^(n-1)/((2*n+1)*(2*n+3)),n,0,inf) r2 = 1/2-1/4*pi 3、将二元函数展开为x的4阶、y的7阶麦克劳林级数，在x=3处展开为6阶泰勒级数。 syms x y z=sin(2*x+5*y^2) z = sin(2*x+5*y^2) taylor(z,5,x) ans = sin(5*y^2)+2*cos(5*y^2)*x-2*sin(5*y^2)*x^2-4/3*cos(5*y^2)*x^3+2/3*sin(5*y^2)*x^4 taylor(z,8,y) ans = sin(2*x)+5*cos(2*x)*y^2-25/2*sin(2*x)*y^4-125/6*cos(2*x)*y^6 taylor(z,7,x,3) ans = sin(6+5*y^2)+2*cos(6+5*y^2)*(x-3)-2*sin(6+5*y^2)*(x-3)^2-4/3*cos(6+5*y^2)*(x-3)^3+2/3*sin(6+5*y^2)*(x-3)^4+4/15*cos(6+5*y^2)*(x-3)^5-4/45*sin(6+5*y^2)*(x-3)^6 4、对符号表达式，分别进行如下变换。 （1）关于y的Fourier变换及其逆变换。 （2）关于x的Laplace变换及其逆变换。 （3）分别关于x和y的Z变换及其逆变换。 解: syms x y syms x y u v w z=(2*x+3*y)*exp(-2*x^2-3*y^2) z = (2*x+3*y)*exp(-2*x^2-3*y^2) fourier(z,y,u) ans = 1/6*exp(-2*x^2)*3^(1/2)*exp(-1/12*u^2)*pi^(1/2)*(4*x-i*u) ifourier(z,y,u) ans = 1/12*3^(1/2)/pi^(1/2)*(4*x+i*u)*exp(-2*x^2-1/12*u^2) laplace(z,x,u) ans = 1/2*exp(-3*y^2)-1/8*erfc(1/4*u*2^(1/2))*exp(-3*y^2+1/8*u^2)*2^(1/2)*pi^(1/2)*(u-6*y) ilaplace(z,x,u) ans = 2*ilaplace(exp(-2*x^2-3*y^2)*x,x,u)+3*y*ilaplace(exp(-2*x^2-3*y^2),x,u) ztrans(z,x,u) ans = -2*u*diff(ztrans(exp(-2*x^2-3*y^2),x,u),u)+3*y*ztrans(exp(-2*x^2-3*y^2),x,u) iztrans(z,x,u) ans = iztrans((2*x+3*y)*exp(-2*x^2-3*y^2),x,u) ztrans(z,y,u) ans = 2*x*ztrans(exp(-2*x^2-3*y^2),y,u)-3*u*diff(ztrans(exp(-2*x^2-3*y^2),y,u),u) iztrans(z,y,u) ans = iztrans((2*