example8_12a.m

来自「<MATLAB 7辅助控制系统设计与仿真>书籍内的源代码.」· M 代码 · 共 13 行

num=400;
den=[1 2 400];
G=tf(num,den);  %系统传递函数
nW1=100*conv([0.005,1],[0.005,1]);
dW1=conv([0.2,1],[0.2,1]);
nW3=[1,0,0];
dW3=[40000];
[m1,p1,w1]=bode(tf(nW1,dW1));   %加权函数W1
semilogx(w1,20*log10(m1(:)))
figure;
[m3,p3,w3]=bode(tf(nW3,dW3));  %加权函数W3
semilogx(w3,20*log10(m3(:)))