ans = 0.9307
39．试用roots命令算出x?2x2?x?1?0的根。
>> d=[1 2 1 1] d =
1 2 1 1 >> r=roots(d) r =
-1.7549 -0.1226 + 0.7449i -0.1226 - 0.7449i
10．试用residue命令来计算下列表达式的部分分式展开：同时再利用residue命令，验算所得答案是否正确 >> collect((s+1)^3*(s^2+1)) ans =
1+s^5+3*s^4+4*s^3+4*s^2+3*s >> b=[0 0 0 3 5 2]; >> a=[1 3 4 4 3 1]; >> [r,p,k]=residue(b,a) r =
1.0000 -0.5000 0.0000 -0.5000 - 0.7500i -0.5000 + 0.7500i p =
-1.0000 -1.0000 -1.0000 0.0000 + 1.0000i 0.0000 - 1.0000i k =
3s2?5s?2(s?1)3(s2?1)
[] 验证：
>> [b,a]=residue(r,p,k) b =
-0.0000 0.0000 3.0000 5.0000 2.0000 a =
1.0000 3.0000 4.0000 4.0000 3.0000 1.0000 11．已知一个二阶的离散系统的脉冲传递函数为：
G(z)?0.9Z?0.05Z2?0.6z?0.8
试绘制离散系统当Ts=0.2s时的Bode图。 >> g=tf([0.9 0.05],[1 -0.6 0.8],0.2); >> bode(g)
12．
f(t)?e?2tcos(10t?15)?1,求L[d5f(t)/dt5]>> syms t
>> F=exp(-2*t)*cos(10*t+15)+1; >> F1=diff(F,t,5) F1 =
-92032*exp(-2*t)*cos(10*t+15)-60800*exp(-2*t)*sin(10*t+15) >> F2=simple(laplace(F1)) F2 =
-128*(719*cos(15)*s+6188*cos(15)-6240*sin(15)+475*sin(15)*s)/(s^2+4*s+104)