MATLAB次上机作业曹红兴.docVIP

第五次上机作业 时间：2013-6-25 E-207（周二）晚上6：00~10：00 班级：071111学号名：曹红兴 符号运算： 求符号矩阵的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 程序： A = sym ( [a11 a12 a13 ; a21 a22 a23 ; a31 a32 a33] ); B = det(A); C = inv(A); [ Bs,e] = subexpr( B,e) [ Cs,d] = subexpr( C,d ) 运行结果： Bs = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 e = [ empty sym ] Cs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 对函数，当为正实数，求。 程序： syms k a z f = a^k * z^(-k); S = symsum(f,k,0,inf) 运行结果： S = piecewise([a = z, Inf], [a 0 and z = 0, 0], [a z and z 0, limit((a^k*z)/(z^k*(a - z)), k = Inf) - z/(a - z)]) 对，求。 程序： syms k x f = 2*((x-1)/(x+1))^(2*k+1)/(2*k+1) S = symsum(f,k,0,inf) 运行结果： f = (2*((x - 1)/(x + 1))^(2*k + 1))/(2*k + 1) S = piecewise([abs(x + 1)*(abs(x - 1) - abs(x + 1)) 0, 2*atanh((x - 1)/(x + 1))]) 求的导数，， 程序： syms x f = abs(sin(x)) g = diff(f,x) g1 = subs(g,x,-eps) g2 = subs(g,x,pi/2) 运行结果： f = abs(sin(x)) g = sign(sin(x))*cos(x) g1 = -1 g2 = 6.1232e-017 计算二重积分 程序： syms x y f = int(int((x^2+y^2),y,1,x^2),x,1,2) vpa(f) 运行结果： f = 1006/105 ans = 9.5809523809523809523809523809524 求的Fourier变换，并画出时的频谱 程序； syms A t w tao; ft = A*(1-abs(t)/tao)*(heaviside(t+tao)-heaviside(t-tao)); Fw = fourier(ft,t,w) h = subs(Fw,A,2); h = subs(h,tao,2); x = linspace(-20,20,1000)； sn = subs(h,w,x); plot(x,sn) 运行结果： Fw = A*((sign(tao)/w^2 + 1/w^2 + pi*dirac(-w, 1)*i - (exp(tao*w*i)*sign(tao))/w^2 + (tao*exp(tao*w*i)*sign(tao)*i)/w)/tao - (1/w^2 - sign(tao)/w^2 + pi*dirac(-w, 1)*i + ((1/exp(tao*w*i))*sign(tao))/w^2 + (tao*(1/exp(tao*w*i))*sign(tao)*i)/w)/tao + (1/exp(tao*w*i))*(- pi*dirac(-w) + i/w) - exp(tao*w*i)*(- pi*dirac(-w) + i/w)) Warning: Imaginary parts of complex X and/or Y arguments ignored In workfive at 10