%控制器Bode圖的顯示子程序 
Kfr=frsp(K,w);
figure(2);
subplot(2,1,1)
vplot('liv,lm',Kfr);
title('Bode plot of Controller')
xlabel('Frequency[rad/s]')
ylabel('Gain')
grid;
subplot(2,1,2)
vplot('liv,p',Kfr);%X加逗號
xlabel('Frequency[rad/s]')
ylabel('Phase(degrees)')
grid;