如何画一阶butterworth的响应函数？

re-us

如何画下面的butterworth的响应函数。一阶的，要求中心频率的响应为1，两侧的响应为0.5和0.5.即图中的虚线。知道有butter这个函数，但是没有找到如何画响应函数的命令？

谢谢

songzy41

采样频率为1（1/天），通带是1/60~1/30（1/天）。所以滤波器设计程序可以有：
fs=1;
fl=1/60; fh=1/30;
[b,a]=butter(1,[fl fh]/(fs/2));
fprintf('B=%5.6f  %5.6f  %5.6f\n',b);
fprintf('A=%5.6f  %5.6f  %5.6f\n',a);
[H,w]=freqz(b,a);
plot((w/2/pi),H);
axis([0 0.1 0 1]);
set(gca, 'XTickMode', 'manual', 'XTick', [0, fl, fh, 0.1]);
grid;
求出的滤波器系数为
B=0.049798  0.000000  -0.049798
A=1.000000  -1.879583  0.900404
其频率响应图为：