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

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

2、5、分析處理后信號(hào)的波形圖和頻譜圖變化三三. .設(shè)計(jì)步驟設(shè)計(jì)步驟1 1、讀取、讀取3s3s的語音信號(hào)并畫出時(shí)域波形和頻譜圖的語音信號(hào)并畫出時(shí)域波形和頻譜圖x1=wavread(yuying111.wav); %讀取信號(hào)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(原始信號(hào)波形);subplot(2,1,2);plot(X); tit

3、le(原始信號(hào)頻譜);2 2、對語音信號(hào)進(jìn)行采樣并畫出采樣后信號(hào)的時(shí)域波、對語音信號(hào)進(jìn)行采樣并畫出采樣后信號(hào)的時(shí)域波和頻譜圖和頻譜圖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); %對信號(hào)做32000點(diǎn)FFT變換figure(2);subplot(2,1,1);plot(x);title(采樣后信號(hào)的時(shí)域波形);subplot(2,1,2);plot(n1(1:1000),Y(1:1000); titl

4、e(采樣后信號(hào)的頻譜圖); %信號(hào)幅度3 3、對采樣后的信號(hào)延時(shí)，并畫出延時(shí)后的時(shí)域波形、對采樣后的信號(hào)延時(shí)，并畫出延時(shí)后的時(shí)域波形和頻譜圖和頻譜圖z1=zeros(1000,1);x; %對信號(hào)進(jìn)行延時(shí)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(延時(shí)后的時(shí)域圖1); %畫出延時(shí)后的信號(hào)時(shí)域圖 subplot(3,1,2);plot(z2); title(延時(shí)后的時(shí)域圖2);su

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

6、00,1); %使語音信號(hào)與延時(shí)后信號(hào)同等長度x2=x;zeros(2000,1);x3=x;zeros(3000,1);y1=x1+z1; %信號(hào)的混響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(混響的時(shí)域圖1); %混響時(shí)域圖subplot(3,1,2);plot(y2); title(混響的時(shí)域圖2);subplot(3,1,3);plot(y3); title(混響的時(shí)域圖3);Y1=fft(y1,3200

7、0); %對混響信號(hào)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è)計(jì)，并對信號(hào)進(jìn)行處理、三種混響器的設(shè)計(jì)，并對信號(hào)進(jìn)行處理（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); %濾波器進(jìn)行濾波YY1=fft(yy1,32000); %經(jīng)無限回聲濾波器后的信號(hào)做32000點(diǎn)的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); %濾波器進(jìn)行

9、濾波YY2=fft(yy2,32000); %經(jīng)多重回聲濾波器后的信號(hào)做32000點(diǎn)的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); %濾波器進(jìn)行濾波YY3=fft(yy3,32000); %經(jīng)全通結(jié)構(gòu)的混響器后的信號(hào)做32000點(diǎn)的FFT變sound(yy1,fs,bits);sound(yy2,fs,bits);sound(yy3,fs,bits);（4 4）畫出經(jīng)混響器處理

10、后信號(hào)的時(shí)域波形和頻譜圖）畫出經(jīng)混響器處理后信號(hào)的時(shí)域波形和頻譜圖figure(8); subplot(2,1,1);plot(yy1); title(無限個(gè)回聲濾波器時(shí)域圖); %無限回聲濾波器時(shí)域波形subplot(2,1,2);plot(n1(1:1000),YY1(1:1000); title(無限個(gè)回聲濾波器頻譜圖 ); %無限回聲濾波器頻譜圖figure(9)subplot(2,1,1);plot(yy2); title(多重回聲濾波器的時(shí)域圖) %多重回聲濾波器的混響器時(shí)域波形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)濾波器的時(shí)域圖) %全通結(jié)構(gòu)的混響器時(shí)域波形subplot(2,1,2);plot(n1(1:1000),YY3(1:1000); title(全通結(jié)構(gòu)濾波器的頻譜圖) %全通結(jié)構(gòu)的混響器頻譜圖四四. .運(yùn)行結(jié)果運(yùn)行結(jié)果1.1.原始信號(hào)時(shí)域波形和頻譜圖原始信號(hào)時(shí)域波形和頻譜圖2 2、采樣后信號(hào)的時(shí)域波形和頻譜圖、采樣后信號(hào)的時(shí)域波形和頻譜圖3 3、對采樣后信號(hào)延時(shí)的時(shí)域波形和頻譜圖、對采樣后信號(hào)延時(shí)的時(shí)域波形和頻譜圖4 4、混響后的時(shí)域波形和頻譜圖、混響后的時(shí)域波形和頻譜圖5 5、原始信號(hào)的時(shí)域波形