MATLAB IIR 滤波器设计流程图及代码示例

本教程将详细介绍使用 MATLAB 设计 IIR 滤波器的流程，包括代码示例和流程图。

流程图:

开始节点 -> 设置通带截至频率/Hz -> 设置阻带截至频率/Hz -> 设置阻带波纹/dB -> 设置通带波纹/dB -> 求通带截止角频率 -> 求阻带截止角频率 -> 求低通滤波器阶数和截止频率 -> 求S域频率响应参数 -> 双线性变换实现S域到Z域变换 -> 根据参数求频率响应 -> 绘制滤波器 -> 利用滤波器对加噪语音信号滤波 -> 对滤波后的语音信号进行傅里叶变换 -> 绘制加噪语音信号时域波形 -> 绘制滤波后语音信号时域波形 -> 绘制加噪语音信号频谱 -> 绘制滤波后语音信号频谱 -> 发声滤波后的语音信号 -> 导出去噪语音信号(单轨) -> 结束节点

代码示例:

Filter_p=0.5*1e4;              % 设置通带截至频率/Hz
Filter_s=0.7*1e4;              % 设置阻带截至频率/Hz
As=20;                         % 阻带波纹/dB
Ap=1;                          % 通带波纹/dB
wp=2*pi*Filter_p/fs;
ws=2*pi*Filter_s/fs; 
filter_p=2*fs*tan(wp/2);
filter_s=2*fs*tan(ws/2);
[n,wn]=buttord(wp,ws,Ap,As,'s');% 求低通滤波器阶数和截止频率
[b,a]=butter(n,wn,'s');         % 求S域频率响应参数 
[B,A]=bilinear(b,a,1);          % 双线性变换实现S域到Z域变换 
[h,w]=freqz(B,A);               % 根据参数求频率响应
figure
plot(w*fs*0.5/pi,abs(h));       % 绘制滤波器
title('IIR滤波器');
xlabel('频率/Hz');ylabel('幅度');grid on;axis tight;
% 滤波
x_denoise=filter(B,A,x_noise);  % 利用滤波器对加噪语音信号滤波 
X_denoise=fft(x_denoise); 
figure  
subplot(2,2,1);plot(t,x_noise);
title('加噪语音信号时域波形');
xlabel('时间/s');ylabel('幅度');grid on;axis tight;
subplot(2,2,2);plot(t,x_denoise);
title('滤波后语音信号时域波形');
xlabel('时间/s');ylabel('幅度');grid on;axis tight;
subplot(2,2,3);plot(f,fftshift(abs(X_noise)));
title('加噪语音信号频谱');
xlabel('频率/Hz');ylabel('幅度');grid on;axis tight; 
subplot(2,2,4);plot(f,fftshift(abs(X_denoise)));
title('滤波后语音信号频谱');
xlabel('频率/Hz');ylabel('幅度');grid on;axis tight;
% 发声，导出去噪音频
sound(x_denoise,fs);                           % 发声滤波后的语音信号         
audiowrite('yinpin_denoise.mp4',x_denoise,fs); % 导出去噪语音信号(单轨)

代码说明:

代码中使用了 buttord 函数计算低通滤波器的阶数和截止频率，butter 函数求解 S 域频率响应参数，bilinear 函数实现 S 域到 Z 域的双线性变换。
freqz 函数根据参数计算频率响应，并绘制滤波器的频率响应曲线。
使用 filter 函数对加噪语音信号进行滤波，并利用 fft 函数对滤波后的信号进行傅里叶变换。
最后使用 sound 函数播放滤波后的语音信号，并使用 audiowrite 函数将去噪后的语音信号导出为音频文件。

注意事项:

代码中需要根据实际情况修改 fs、x_noise 和 t 等参数。
滤波器的参数需要根据实际需求进行调整。
滤波器设计是一个复杂的工程问题，需要根据具体的应用场景进行优化。

更多信息:

MATLAB 官方文档: https://www.mathworks.com/help/matlab/
数字信号处理相关书籍和网站

希望本教程对您有所帮助。如果您有任何问题，请随时提出。