语音信号处理及MATLAB实现-数字信号处理 课程设计
3.6对加高频噪音的程序进行滤波处理
用椭圆滤波器对信号的高频噪音进行滤波,程序如下: %加高通滤波 wp1=100*pi*2/fs; ws1=60*pi*2/fs; Rp1=1; As1=100;
[N1,wso]=ellipord(wp1,ws1,Rp1,As1); [B1,A1]=ellip(N1,Rp1,As1,wso,'high'); figure(9) freqz(B1,A1); grid on; ylim([0,1]); y4=filter(B1,A1,x4); X1=fft(y4,20480); %滤波后时域图 subplot(2,1,1); plot(y4); grid on; %axis tight; ylim([-0.03,0.03]);
title('滤波后信号时域波形'); xlabel('time(s)'); ylabel('幅度'); %滤波后频谱 subplot(2,1,2);
plot(abs(X1(1:10240))); title('滤波后信号频谱'); grid on;