chebyshev1gt.m

用MATLAB软件分析音频信号通过高通滤波前后的频谱

clear; close all     %清除工作空间中的变量                                
%读取wav文件
[y,fs,bits] =wavread('yinyue.wav');
sound(y,fs,bits); %对声音进行播放
Y=fft(y,4096); %4096点DFT，用FFT实现
%画出滤波前的音频信号的波形及频谱
subplot(221);plot(y);title('滤波前音频信号波形');
subplot(222);plot(abs(Y));title('滤波前音频信号频谱');
%高通chebyshev1型数字滤波器的设计
wp=80;ws=60;rp=1;rs=40;Fs=200;      %输入指标
% 求滤波器阶次和通带截止频率wn
[n,wn]=cheb1ord(wp/(Fs/2),ws/(Fs/2),rp,rs,'z'); 
[num,den]=cheby1(n,rp,wn,'high');    %设计滤波器,得出系数
[H,W]=freqz(num,den);%求数字滤波器的频率响应
x=filter(num,den,y);%对音频信号进行滤波
X=fft(x,4096);%4096点DFT，用FFT实现
%画出滤波后的音频信号的波形及频谱
subplot(223);plot(x);title('滤波后音频信号波形');
subplot(224);plot(abs(X));title('滤波后音频信号频谱');
sound(x,fs,bits);%回放音频信号