ellipbfilt.m

来自「使用Matlab实现了FIR与IIR滤波器」· M 代码 · 共 36 行

function [ output_args ] = Untitled6( input_args )
%UNTITLED6 Summary of this function goes here
%  Detailed explanation goes here
[y,ft,bits]=wavread('hello.wav');
subplot(2,2,1);
plot(y);
xlabel('Time Index t');
ylabel('Amplitude');
title('原始时域图象');grid;
subplot(2,2,2);
plot(abs(fft(y)));
xlabel('Frequence Index w');
ylabel('Amplitude');
title('原始频域图象');grid;
Wp1b =2*1200/ft;
Wp2b=2*3000/ft;
Ws1b= 2*1000/ft;
Ws2b=2*3200/ft;
Wpb=[Wp1b,Wp2b];
Wsb=[Ws1b,Ws2b];
Rp = 1;
Rs = 100;
[Nb,Wnb] = ellipord(Wpb,Wsb,Rp,Rs);
[bb,ab] = ellip(Nb,Rp,Rs,Wnb);
vb=filter(bb,ab,y);
subplot(2,2,3);
plot(vb);
xlabel('Time Index t');
ylabel('Amplitude');
title('ELLIP带通滤波时域图象');grid;
subplot(2,2,4);
plot(abs(fft(vb)));
xlabel('Frequence Index w');
ylabel('Amplitude');
title('带通滤波后频域图象');grid;
sound(vb,ft,bits);