数字信号课程设计报告

多重回声滤波器幅频响应：

（4）具体代码实现（多回声）

%原声音信号

[x,fs,bits]=wavread('E:\\mcpass.wav'); N=size(x,1);

t=(0:length(x)-1)/fs; f=fs*(0:(N+1)/2-1)/N+1; %多重回声 a=0.5; R=fs*a;

B=[1,zeros(1,R-2)]; A=[1,zeros(1,R-1),a];

yd = filter(B,A,x); %滤波器函数 %频率响应

[Hb,Wc]=freqz(B,A); %求多重回声滤波器幅频响应函数 %绘图部分 figure(2);

plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('多重滤波器频域响应图')

xlabel('频率(Hz)'); ylabel('幅度'); figure(1);

subplot(4,1,1);

plot(t,x);title('原信号时域'); xlabel('时间(s)'); ylabel('幅度'); X=fft(x,N+1); subplot(4,1,3);

plot(f,abs(X(1:(N+1)/2)));title('原信号频域'); xlabel('频率(Hz)'); ylabel('幅度'); figure(1);

subplot(4,1,2);

plot(t,yd);title('多重信号时域图') xlabel('时间(s)'); ylabel('幅度');

Yd=fft(yd,N+1); %多重回声滤波后信号FFT变换 subplot(4,1,4);

plot(f,abs(Yd(1:(N+1)/2)));title('多重信号频域图'); xlabel('频率(Hz)'); ylabel('幅度'); sound(x); sound(yd);

（5）时域图与频域图（全通滤波器）

原始信号时域图与全通滤波器时域图：

原始信号时域图与全通滤波器频域图：

全通滤波器的频域响应图：

（6）具体代码实现（全通滤波器）：

%原声音信号

[x,fs,bits]=wavread('E:\\mcpass.wav'); N=size(x,1);

t=(0:length(x)-1)/fs; f=fs*(0:(N+1)/2-1)/N+1; %全通结构混响 a=0.5; R=fs*a;

B=[a,zeros(1,R-2),1]; A=[1,zeros(1,R-2),a]; yd = filter(B,A,x); %频率响应

[Hb,Wc]=freqz(B,A); %求全通滤波器幅频响应函数 figure(2);

plot(fs*Wc/(2*pi),20*log10(abs(Hb)));title('全通结构频域响应图') xlabel('频率(Hz)');

ylabel('幅度'); figure(1);

subplot(4,1,1);

plot(t,x);title('原信号时域'); xlabel('时间(s)'); ylabel('幅度'); X=fft(x,N+1); subplot(4,1,3);

plot(f,abs(X(1:(N+1)/2)));title('原信号频域'); xlabel('频率(Hz)'); ylabel('幅度'); figure(1);

subplot(4,1,2);

plot(t,yd);title('全通滤波器时域图'); xlabel('时间(s)'); ylabel('幅度'); Yd=fft(yd,N+1); subplot(4,1,4);

plot(f,abs(Yd(1:(N+1)/2)));title('全通滤波器频域图'); xlabel('频率(Hz)'); ylabel('幅度'); sound(x); sound(yd);

（7）时域图与频域图（延时与混响）

原始信号时域图：

延时时域图：

混响时域图：