给定的技术指标是:
已知通带截止频率fp=5kHz,通带最大衰减Wp=2dB,阻带截止频率fs=12kHz,阻带最小衰减Ws=30dB,按照以上技术指标设计巴特沃斯低通滤波器。
程序如下:
>> Wp=2*pi*5000; %通带截止角频率
>> Ws=2*pi*12000; %阻带截止角频率
>> Rp=2; %通带最大衰减
>> Rs=30; %阻带最小衰减
>> [n,Wn] = buttord(Wp,Ws,Rp,Rs,'s'); %求巴特沃斯阶数和3db截止角频率
>> [b,a] = butter(n,Wn,'s'); %求传递函数
>> [z,p,k] = butter(n,Wn,'s'); %求零极点及增益
>> w=linspace(1,15000)*2*pi;
>> H =freqs(b,a,w); %频率响应
>> magH=abs(H); %频率响应的幅度
>> phaH=unwrap(angle(H)); %频率响应的相位(平滑处理)
>> plot(w/(2*pi),20*log10(magH)); %频率响应的幅度的曲线图
>> title('巴特沃斯低通滤波器的幅频特性');
>> xlabel('频率/Hz');
>> ylabel('幅度/db')