对于3—9题,特征方程为:3*s^4+10*s^3+5*s^2+s+2=0。所得结果为:该系统不稳定。
用MATLAB所写如下:
>
bac=[3 10 5 1 2];
>> p=roots(bac)
p =
-2.7362 -0.8767 0.1398 + 0.5083i
0.1398 - 0.5083i
由结果可得,并非所有特征方程的跟的实数部分都为负数,所以与线性系统稳定的充分必要条件想矛盾,所以肯定了笔算结果的正确性。该系统不稳定
对于3-11,开环传递函数为:k*(0.5*s+1)/s(s+1)(0.5*s^2+s+1)。笔算所得结果为:0<k<=1.07用MATLAB所写如下:for
k=0:1.07;
num=[0.5*k 1];
den=[0.5 1.5 2 (1+0.5*k) k];
num=[0.5*k 1];Y=tf(num,den);
pzmap(Y)
hold on
end
k=1.07;num=[(0.5*k) 1]
num =
0.5350 1.0000
>> den=[0.5 1.5 2 (1+0.5*k) k]
den =
0.5000 1.5000 2.0000 1.5350 1.0700
H=tf(num,den)
Transfer function:
0.535 s + 1
------------------------------------------
0.5 s^4 + 1.5 s^3 + 2 s^2 + 1.535 s + 1.07
>> sgrid
>> pzmap(H)
所得图如下:
for k=0:10
num=[(0.5*k) 1]
n=[0.5 1.5 2 (1+0.5*k) k]
tf(num,den)
sgrid
pzmap(H)
可以看到所得点均位于坐标轴左半边,所以当k=1.07时。系统是稳定的
有MATLAB的图得,笔算多得结果是正确的