发布时间 : 星期五 文章数字信号课程设计报告更新完毕开始阅读51a8b28ff8c75fbfc67db213
多重回声滤波器幅频响应:
(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)时域图与频域图(延时与混响)
原始信号时域图:
延时时域图:
混响时域图: