关于用matlab画bode图的问题

eric415

已知一个系统K=[3 -1 0 0;-1 3 -2 0;0 -2 4 -1;0 0 -1 2];M=[2 0 0 0;0 3 0 0;0 0 5 0;0 0 0 3];C=[0 0 0 0;0 2 0 0;0 0 3 0;0 0 0 0];
可以用bode命令画出频响函数图，但是是分开的。
A=[zeros(4,4) ones(4,4);-inv(M)*K -inv(M)*C];
B=[zeros(4,4);inv(M)];
C=[1 zeros(1,7)];
D=zeros(1,4);
sys=ss(A,B,C,D);
bode(sys)


有没有什么方法让这4个频率对应的尖峰画到一个图里，变成一条曲线。并标注横向坐标值。
类似这种：

                               
登录/注册后可看大图

ChaChing

[mag,phase,w] = bode(sys) or [mag,phase] = bode(sys,w)
然后各自画+hold on

eric415

画完了是这样的

ChaChing

LZ怎么画! 代码?

eric415

[mag,phase,w] = bode(sys);
for i=1:4
     for k=1:43
          y(i,k)=mag(1,i,k);
     end
end
magdb1= 20*log10(y(1,:));
magdb2= 20*log10(y(2,:));
magdb3= 20*log10(y(3,:));
magdb4= 20*log10(y(4,:));
hold on;
plot(w,20*log10(magdb1))
hold on;
plot(w,20*log10(magdb2))
hold on;
plot(w,20*log10(magdb3))
hold on;
plot(w,20*log10(magdb4))

ChaChing

刚才发现这个忘了回应~!
LZ给的代码好像根本有问题!? 无法画出3F的图!