在命令行中输入:
designfilt
然后就会弹出滤波器的种类,选定后会让你填参数。
参数填好之后点确定,会在命令行窗口生成这个滤波器,复制到编辑器里直接用就好。
%读取音频文件%
[x,Fs]=audioread('C:\Users\16000\Desktop\testfile.wav');
%x是一个序列
%Fs是这个音频的采样率
%复制过来的滤波器%
High_Filt=designfilt('highpassfir', 'StopbandFrequency', 500, 'PassbandFrequency', 600, 'StopbandAttenuation', 60, 'PassbandRipple', 1, 'SampleRate', 48000);%改成单声道%
x=x(:,1);
%滤波操作%
y_1 = filter(High_Filt,x);sigLength=length(y_1);
t=(0:length(y_1)-1)/Fs; %时间轴
subplot(2,1,1);
plot(t,y_1);xlabel('Time(s)');%在第一个窗口画波形
grid on;
使用designfilt方法的好处是自动生成,而且不用在filter命令里那里给滤波器填参数。
%部分重要参数说明%
Fs=Sample Rate;
Passband frequency是滤波器起始点
Stopband frequency是滤波器结束点
例如,Passband=400;Stopband=500的highpassfir就是说在400Hz到500Hz之间把波滤干净。
最后总结一下:
- 在命令行中输入designfilt;
- 设置好后将生成的代码复制进编辑器区域;
- 使用audioread(‘路径’)读入音频,拿到信号x和采样率Fs;
- y=Filter(‘滤波器名称’,x);
- 画图
一些其他看起来有点复杂的东西:
help filterDesigner %滤波器涉及图形化界面
help butter %butterworth滤波器