Maltab设计IIR滤波器

clear clc

[audio,Fs]=audioread('D:\yinpin~1.mp3');%声音读取

audio = audio(:,1); %双通道变单通道

n=length(audio);

T = 1/Fs;%采样间隔

t = (0:n-1)*T;%时间轴

f = (0:n-1)/n*Fs;%频率轴

%Fs=8000;%T=1/Fs;

%绘出原始音频时域波

figure(5)

subplot(2,2,1);

plot(t,audio);

xlabel('时间/s');

ylabel('幅度');

title('初始信号波形');

grid;

%绘出原始音频频域频谱

audio_fft=fft(audio,n)*T;

subplot(2,2,2);

audiof = abs(audio_fft);

plot(f(1:(n-1)/2),audiof(1:(n-1)/2));

title('初始信号频谱');

xlabel('频率/Hz');ylabel('幅度');grid;

%=================================低通滤波器软件实现==&