基于MATLAB的混響效果設(shè)計

1、信號混響效果的設(shè)計信號混響效果的設(shè)計一一. .設(shè)計要求設(shè)計要求 選擇一段語音信號作為分析對象，對其進行頻譜分析，然后在時域用數(shù)字信號處理方法將信號加入3種混響，再分析其頻譜，并與原始信號進行比較。二二. .設(shè)計思路設(shè)計思路 1、利用Windows下的錄音機或其他軟件，錄制一段語音信號，時間控制在3s左右，并對錄制的信號進行采樣 2、語音信號的頻譜分析，畫出采樣后的時域波形和頻譜圖 3、將信號加入延時和混響，再分析其時域波形和頻譜圖，并與原始信號頻譜進行比較 4、設(shè)計幾種特殊類型的混響器：無限個回聲混響器，多重回聲混響器，全通結(jié)構(gòu)的混響器，對采樣后的信號進行混響處理，畫出處理后的波形圖和頻譜圖

2、5、分析處理后信號的波形圖和頻譜圖變化三三. .設(shè)計步驟設(shè)計步驟1 1、讀取、讀取3s3s的語音信號并畫出時域波形和頻譜圖的語音信號并畫出時域波形和頻譜圖x1=wavread(yuying111.wav); %讀取信號x,fs,bits=wavread(yuying111.wav,1024 5120);x=x(:,1); %只取單聲道sound(x,fs,bits);X=fft(x,32000);magX=abs(X);angX=angle(X);figure(1);subplot(2,1,1);plot(x);title(原始信號波形);subplot(2,1,2);plot(X); tit

3、le(原始信號頻譜);2 2、對語音信號進行采樣并畫出采樣后信號的時域波、對語音信號進行采樣并畫出采樣后信號的時域波和頻譜圖和頻譜圖x,fs,bits=wavread(yuying111.wav,10000 40000);x=x(:,1); %只取單聲道sound(5*x,fs); %對聲音的回放n1=0:2000;N=size(x,1);Y=fft(x,32000); %對信號做32000點FFT變換figure(2);subplot(2,1,1);plot(x);title(采樣后信號的時域波形);subplot(2,1,2);plot(n1(1:1000),Y(1:1000); titl

4、e(采樣后信號的頻譜圖); %信號幅度3 3、對采樣后的信號延時，并畫出延時后的時域波形、對采樣后的信號延時，并畫出延時后的時域波形和頻譜圖和頻譜圖z1=zeros(1000,1);x; %對信號進行延時z2=zeros(2000,1);x;z3=zeros(3000,1);x;Z1=fft(z1,32000);Z2=fft(z2,32000);Z3=fft(z3,32000);figure(3); subplot(3,1,1);plot(z1); title(延時后的時域圖1); %畫出延時后的信號時域圖 subplot(3,1,2);plot(z2); title(延時后的時域圖2);su

5、bplot(3,1,3);plot(z3); title(延時后的時域圖3); figure(4)subplot(3,1,1);plot(n1(1:1000),Z1(1:1000);title(延時后的頻譜圖1); %延時后的信號頻譜圖subplot(3,1,2);plot(n1(1:1000),Z2(1:1000);title(延時后的頻譜圖2); subplot(3,1,3);plot(n1(1:1000),Z3(1:1000);title(延時后的頻譜圖3); 4 4、對信號進行混響，并畫出混響后的時域波形和頻譜圖、對信號進行混響，并畫出混響后的時域波形和頻譜圖x1=x;zeros(10

6、00,1); %使語音信號與延時后信號同等長度x2=x;zeros(2000,1);x3=x;zeros(3000,1);y1=x1+z1; %信號的混響y2=x2+z2;y3=x3+z3;sound(y1,fs,bits);sound(y2,fs,bits);sound(y3,fs,bits);figure(5);subplot(3,1,1);plot(y1); title(混響的時域圖1); %混響時域圖subplot(3,1,2);plot(y2); title(混響的時域圖2);subplot(3,1,3);plot(y3); title(混響的時域圖3);Y1=fft(y1,3200

7、0); %對混響信號FFT變換Y2=fft(y2,32000); Y3=fft(y3,32000);figure(6);subplot(3,1,1);plot(n1(1:1000),Y1(1:1000); title(混響的頻譜圖1); %混響頻譜圖subplot(3,1,2);plot(n1(1:1000),Y2(1:1000); title(混響的頻譜圖2); subplot(3,1,3);plot(n1(1:1000),Y3(1:1000); title(混響的頻譜圖3); 5 5、三種混響器的設(shè)計，并對信號進行處理、三種混響器的設(shè)計，并對信號進行處理（1 1）無限回聲混響器）無限回聲混

8、響器a=0.05; %a取小于等于1Bz=0,0,0,0,0,0,0,0,0,0,1; %分子的系數(shù)Az=1,0,0,0,0,0,0,0,0,0,-a; %分母的系數(shù)yy1=filter(Bz,Az,z1); %濾波器進行濾波YY1=fft(yy1,32000); %經(jīng)無限回聲濾波器后的信號做32000點的FFT變換（2 2）多重回聲混響器）多重回聲混響器a=0.05; %a取小于等于1N=5Bz1=1,0,0,0,0,0,0,0,0,0,-0.5N %分子的系數(shù)Az1=1,0,0,0,0,0,0,0,0,0,-0.5; %分母的系數(shù)yy2=filter(Bz1,Az1,z1); %濾波器進行

9、濾波YY2=fft(yy2,32000); %經(jīng)多重回聲濾波器后的信號做32000點的FFT變換（3 3）全通結(jié)構(gòu)的混響器）全通結(jié)構(gòu)的混響器a=0.05; %a取小于等于1Bz1=a,0,0,0,0,0,0,0,0,0,1; %分子的系數(shù)Az1=1,0,0,0,0,0,0,0,0,0,a; %分母的系數(shù)yy3=filter(Bz1,Az1,z1); %濾波器進行濾波YY3=fft(yy3,32000); %經(jīng)全通結(jié)構(gòu)的混響器后的信號做32000點的FFT變sound(yy1,fs,bits);sound(yy2,fs,bits);sound(yy3,fs,bits);（4 4）畫出經(jīng)混響器處理

10、后信號的時域波形和頻譜圖）畫出經(jīng)混響器處理后信號的時域波形和頻譜圖figure(8); subplot(2,1,1);plot(yy1); title(無限個回聲濾波器時域圖); %無限回聲濾波器時域波形subplot(2,1,2);plot(n1(1:1000),YY1(1:1000); title(無限個回聲濾波器頻譜圖 ); %無限回聲濾波器頻譜圖figure(9)subplot(2,1,1);plot(yy2); title(多重回聲濾波器的時域圖) %多重回聲濾波器的混響器時域波形subplot(2,1,2);plot(n1(1:1000),YY2(1:1000); title(多重

11、回聲濾波器的頻譜圖) %多重回聲濾波器的頻譜圖figure(10)subplot(2,1,1);plot(yy3); title(全通結(jié)構(gòu)濾波器的時域圖) %全通結(jié)構(gòu)的混響器時域波形subplot(2,1,2);plot(n1(1:1000),YY3(1:1000); title(全通結(jié)構(gòu)濾波器的頻譜圖) %全通結(jié)構(gòu)的混響器頻譜圖四四. .運行結(jié)果運行結(jié)果1.1.原始信號時域波形和頻譜圖原始信號時域波形和頻譜圖2 2、采樣后信號的時域波形和頻譜圖、采樣后信號的時域波形和頻譜圖3 3、對采樣后信號延時的時域波形和頻譜圖、對采樣后信號延時的時域波形和頻譜圖4 4、混響后的時域波形和頻譜圖、混響后的時域波形和頻譜圖5 5、原始信號的時域波形