基于MATLAB的声音的延迟与混响(4)

z=[zeros(3000,1);y]; y=[y;zeros(3000,1)]; a=0.5; N=5;

R=9000;%滤波器阶数设置，其值越高，回声越明显

Bz1=[1,zeros(1,N*R-1),-a^N];%多重回声滤波器系统函数分子 Az1=[1,zeros(1,R-1),-a];%多重回声滤波器系统函数分母 y2=filter(Bz1,Az1,y); Y2=fft(y2,6001); [h,w]=freqz(Bz1,Az1); subplot(3,2,1);plot(abs(h)); title('多重回声滤波器幅频响应'); subplot(3,2,2); plot(angle(h)); title('多重回声滤波器相频响应');

subplot(3,2,3:4);plot(y2);title('多重回声滤波器时域图 '); subplot(3,2,5);plot(abs(Y2));title('多重回声滤波器幅频 '); subplot(3,2,6);plot(angle(Y2));title('多重回声滤波器相频 '); sound(y2,fs,bits);

3.7 无限个回声滤波器程序设计

[y,fs,bits]=wavread('C:\\Users\\Administrator\\Desktop\\音号.wav',[50000 150000]); y=y(:,1);

z=[zeros(3000,1);y]; y=[y;zeros(3000,1)]; a=0.5;

R=5000;%滤波器阶数设置，其值越高，回声越明显

16

频信

Bz=[0,zeros(1,R-1),1];%无限个回声滤波器系统函数分子 Az=[1,zeros(1,R-1),-a];%无限个回声滤波器系统函数分母 y2=filter(Bz,Az,y); Y2=fft(y2,6001); [h,w]=freqz(Bz,Az);

subplot(3,2,1);plot(abs(h)); title('无限个回声滤波器幅频响应'); subplot(3,2,2); plot(angle(h)); title('无限个回声滤波器相频响应');

subplot(3,2,3:4);plot(y2);title('无限个回声滤波器时域图 '); subplot(3,2,5);plot(abs(Y2));title('无限个回声滤波器信号幅频 '); subplot(3,2,6);plot(angle(Y2));title('无限个回声滤波器信号相频 '); sound(y2,fs,bits);

3.8 全通结构滤波器程序设计

[y,fs,bits]=wavread('C:\\Users\\Administrator\\Desktop\\音频号.wav',[50000 150000]); y=y(:,1);

z=[zeros(3000,1);y]; y=[y;zeros(3000,1)]; a=0.5;

R=5000;%滤波器阶数设置，其值越高，回声越明显 Bz1=[a,zeros(1,R-1),1];%全通滤波器系统函数分子 Az1=[1,zeros(1,R-1),a];%全通滤波器系统函数分母 [h,w]=freqz(Bz1,Az1); yy2=filter(Bz1,Az1,y); YY2=fft(yy2,6001);

subplot(3,2,1);plot(abs(h));

17

信

title('全通滤波器幅频响应'); subplot(3,2,2); plot(angle(h)); title('全通滤波器相频响应');

subplot(3,2,3:4);plot(yy2);title('全通结构时域图 '); subplot(3,2,5);plot(abs(YY2));title('全通结构幅频 '); subplot(3,2,6);plot(angle(YY2));title('全通结构相频 '); sound(yy2,fs,bits);

4 结果分析与体会

运行结果：

18

对信号采样后，信号时域波形出现了离散化的分布，不像原始图像那样密集分布，频谱图也是同样的变换，幅度和频率都是根据采样区间变化的

19

从时域图出延时是让原来的波形向右移动3000，频域的幅度发生了变化。

20

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库基于MATLAB的声音的延迟与混响(4)在线全文阅读。