求救于各位大神！！感谢！！

jm258

wn=6;
kosi=[0.1:0.1:1.0,2];
w=logspace(-1,1,100);
figure(1)
for kos=kosi
    num=[2*kos*wn,wn .^2]
    den=[1,2*kos*wn,wn .^2]
    [mag,pha,w]=bode(num,den,w)
    subplot(2,1,1);
    hold on
    semilogx(w1,mag);
    subplot(2,1,2);
    hold on
    semilogx(w1,pha);
    subplot(2,1,1);
    grid on
    title('Bode Plot');
    xlabel('Frequency (rad/sec)');
    ylabel('Gain dB')
    subplot(2,1,2);
    grid on
    xlabel('Frequency (rad/sec)');
    ylabel('Phase deg')
    hold off

究竟这个程序出现了什么错误呢？无法得到振动系统bode图