elliphfilt.m

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

function [ output_args ] = Untitled5( input_args )
%UNTITLED5 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;
Wph=2*5000/ft;
Wsh=2*4800/ft;
Rp = 1;
Rs = 100;
[Nh,Wnh] = ellipord(Wph,Wsh,Rp,Rs);
[bh,ah] = ellip(Nh,Rp,Rs,Wnh,'high');
vh=filter(bh,ah,y);
subplot(2,2,3);
plot(vh);
xlabel('Time Index t');
ylabel('Amplitude');
title('ELLIP高通滤波时域图象');grid;
subplot(2,2,4);
plot(abs(fft(vh)));
xlabel('Frequence Index w');
ylabel('Amplitude');
title('高通滤波后频域图象');grid;
sound(vh,ft,bits);