cheby2_lowpass.m

设计了巴特沃斯、车比雪夫I

%设计车比雪夫Ⅱ型 低通滤波器
fs=360;%采样频率，the Nyquist frequency (180 Hz)
fn=fs/2;
Wp=10/fn;Ws=100/fn;%分别为通带截止频率和阻带截止频率
Rp=3;As=60; %通带最大衰减Rp=3 dB;
                                                %阻带最小衰减As=60
[n,Wn]=cheb2ord(Wp,Ws, Rp,As)   
[b, a]=cheby2(n,Rp,Wn,'low');                    %
%[num,den]=bilinear(b,a,1000); 
%sztitle=sprintf('% d阶车比雪夫Ⅱ型 低通', n);
%freq(b,a,sztitle); 
freqz(b,a);
sztitle=sprintf('% d 阶车比雪夫Ⅱ型滤波器',n)
title(sztitle)



load st300.txt;
in=st300(1:21600,2);
%绘制输入信号
figure(2)
subplot(3, 2, [1 2]);
plot(in); 
grid on;
ylabel('幅度');
title('输入信号');
%设计巴特沃思低通滤波器
%使用设计的滤波器滤波,产生输出信号
out=filter(b, a, in);
%绘制输出信号
subplot(3, 2, [5 6]);
plot(out); 
grid on;
ylabel('幅度');
title('低通滤波后输出信号');
fprintf('车比雪夫Ⅱ型低通滤波器阶数=%d',n)

%[H,w]=freqz(b, a,FH, 2*FH);
[H,w]=freqz(b, a, 250, 'whole', 1);%求频率响应;
Hr=abs(H); %求幅频
Hphase=angle(H); %求相频
Hphase=unwrap(Hphase); %解卷绕
subplot(3, 2, 3)
plot(w,Hr); 
axis([0 0.5 0 1.2]);
grid on;
sztitle=sprintf('% d阶车比雪夫Ⅱ型 低通滤波特征', n);
title(sztitle);
ylabel('幅频响应')
xlabel('归一化频率')
subplot(3, 2, 4)
plot(w,Hphase); 
axis([0 0.5 -7 0]);
grid on;
ylabel('相频响应')
xlabel('归一化频率')