matlab三维图得几种画法,第三章MATLAB分析验证

对于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的图得，笔算多得结果是正确的