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