【文章內(nèi)容簡介】 : RRazzzH???? 1)( , a 1 （ 4）全通結(jié)構(gòu)的混響器的系統(tǒng)函數(shù) : 全通濾波器的傳遞函數(shù)公式為 H(Z)=y(Z)/X(Z)=(K+Z^(m))/(1K*Z^(m)) 其中 m為回聲延時取樣， k為反饋系數(shù)。 用直接 1型表示這個傳遞函數(shù)則為： y(n)=k*x(n)+x(nm)+k*y(nm) 可見其實際上是一個簡單的 IIR 濾波器，時間 n的輸出有時間 n 的輸入和 m點之前的輸入與輸出計算而得。由于這個 IIR 濾波器的頻率響應為水平直線，所以被稱為全通濾波器。 13 3 程序設計的步驟和過程 原始語音信號時域波形和頻譜 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\音 頻信號 .wav39。)。%用于讀取語音，采樣值放在向量 y 中， fs 表示采樣頻率 (Hz)， bits 表示采樣位數(shù)。 [50000 150000]表示讀取從 50000 點到 150000 點的值（若只有一個 N的點則表示讀取前 N點的采樣值）。 y=y(:,1)。%只取單信道 sound(y,fs,bits)。%播放原始信號 Y=fft(y)。%對原始信號取傅里葉變換 subplot(2,2,1:2)。plot(y)。title(39。原始信號時域波形 39。)。 subplot(2,2,3)。plot(abs(Y))。title(39。原始信 號幅頻 39。)。 subplot(2,2,4)。plot(angle(Y))。title(39。原始信號相頻 39。)。 采樣后語音信號的時域波形和頻譜 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。%只取原信號的一部分，即時間區(qū)間 [50000 150000] y=y(:,1)。 Y=fft(y,6001)。%抽取 6001點 subplot(2,2,1:2)。plot(y)。title(39。采樣信號時域 波形 39。)。 subplot(2,2,3)。plot(abs(Y))。title(39。采樣信號幅頻 39。)。 subplot(2,2,4)。plot(angle(Y))。title(39。采樣信號相頻 39。)。 sound(y,fs,bits)。%回放采樣信號 14 對采樣后的信號延時 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。%延遲 3000 Z=fft(z,6001)。 subplot(2,2,1:2)。plot(z)。title(39。延時后時域波形 39。)。 subplot(2,2,3)。plot(abs(Z))。title(39。延時后幅頻 39。)。 subplot(2,2,4)。plot(angle(Z))。title(39。延時后相頻 39。)。 sound(y,fs,bits)。%回放延遲信號 對采樣后的信號混響 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。%補 3000 y1=y+z。%將原始信號與延遲信號相加，產(chǎn)生混響信號 Y1=fft(y1,6001)。 subplot(2,2,1:2)。plot(y1)。title(39?；祉懙臅r域波形 39。)。 subplot(2,2,3)。plot(abs(Y1))。title(39。混響的幅頻 39。)。 subplot(2,2,4)。plot(angle(Y1))。title(39?；祉懙南囝l 39。)。 sound(y,fs,bits)。%回放混響信號 15 單回聲濾波器程序設計 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設置，其值越高，回聲越明顯 Bz1=[1,zeros(1,R1),a]。%單回聲濾波器系統(tǒng)函數(shù)分子 Az1=(1)。%單回聲濾波器系統(tǒng)函數(shù) 分母 y2=filter(Bz1,Az1,y)。%單回聲濾波器系統(tǒng)函數(shù) Y2=fft(y2,6001)。 [h,w]=freqz(Bz1,Az1)。%求設計的濾波器頻譜 subplot(3,2,1)。plot(abs(h))。 title(39。單回聲濾波器幅頻響應 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。單回聲濾波器相頻響應 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。單回聲濾波器時域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。單回聲濾波器幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。單回聲濾波器相頻 39。)。 sound(y2,fs,bits)。 多重回聲濾波器程序設計 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。 y=y(:,1)。 16 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 N=5。 R=9000。%濾波器階數(shù)設置，其值越高，回聲越明顯 Bz1=[1,zeros(1,N*R1),a^N]。%多重回聲濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。%多重回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz1,Az1,y)。 Y2=fft(y2,6001)。 [h,w]=freqz(Bz1,Az1)。 subplot(3,2,1)。plot(abs(h))。 title(39。多重回聲濾波器幅頻響應 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。多重回聲濾波器相頻響應 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。多重回聲濾波器時域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。多重回聲濾波器幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。多重回聲濾波器相頻 39。)。 sound(y2,fs,bits)。 無限個回聲濾波器程序設計 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設置，其值越高，回聲越明顯 17 Bz=[0,zeros(1,R1),1]。%無限個回聲濾波器系統(tǒng)函數(shù)分子 Az=[1,zeros(1,R1),a]。%無限個回聲濾波器系統(tǒng)函數(shù)分母 y2=filter(Bz,Az,y)。 Y2=fft(y2,6001)。 [h,w]=freqz(Bz,Az)。 subplot(3,2,1)。plot(abs(h))。 title(39。無限個回聲濾波器幅頻響應 39。)。 subplot(3,2,2)。 plot(angle(h))。 title(39。無限個回聲濾波器相頻響應 39。)。 subplot(3,2,3:4)。plot(y2)。title(39。無限個回聲濾波器時域圖 39。)。 subplot(3,2,5)。plot(abs(Y2))。title(39。無限個回聲濾波器信號幅頻 39。)。 subplot(3,2,6)。plot(angle(Y2))。title(39。無限個回聲濾波器信 號相頻 39。)。 sound(y2,fs,bits)。 全通結(jié)構(gòu)濾波器程序設計 [y,fs,bits]=wavread(39。C:\Users\Administrator\Desktop\ 音頻信號 .wav39。,[50000 150000])。 y=y(:,1)。 z=[zeros(3000,1)。y]。 y=[y。zeros(3000,1)]。 a=。 R=5000。%濾波器階數(shù)設置，其值越高，回聲越明顯 Bz1=[a,zeros(1,R1),1]。%全通濾波器系統(tǒng)函數(shù)分子 Az1=[1,zeros(1,R1),a]。%全通濾波器系統(tǒng)函數(shù)分母 [h,w]=freqz(Bz1,Az1)。 yy2=filter(Bz1,Az1,y)。 YY2=fft(yy2,6001)。 subplot(3,2,1)。plot(abs(h))。 18 title(39。全通濾波器幅頻響應 39。)。 subplot(3,2,2)。 plot(angle(h))。 title