Question

matlab中画波特图的某程序的一些疑问

b=[0.4914];
a=[1,0.9884,1.2386,0.4914];

[h w]=freqs(b,a);
phaiw=angle(h);

subplot(211);
plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('幅度');
title('幅频响应');
subplot(212)
plot(w,unwrap(phaiw) , '-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('相位');
title('相频响应');
这个程序中的w应该是角频率了，那是不是传递函数也必须要H（w）?如果是H（s）需要除以个2π？

Answer 1

b=0.4914;
a=[1,0.9884,1.2386,0.4914];
[h, w]=freqs(b,a);
freqs(b,a)
phaiw=angle(h);
figure
subplot(211);
plot(w,abs(h) ,'-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('幅度');
title('幅频响应');
subplot(212)
plot(w,unwrap(phaiw) , '-k','linewidth',0.5,'markersize',3);
hold on;
xlabel('角频率(rad)');
ylabel('相位');
title('相频响应');
按上面语句试试，图2可参照第一个改

Answer 2

freqs 函数返回的是 角频率 rad/s, 输入的是H(s)的系数，
注意是系数，除什么都无关结果（分子分母消掉了）