基于MATLAB的混响效果设计

1、信号混响效果的设计信号混响效果的设计一一. .设计要求设计要求 选择一段语音信号作为分析对象，对其进行频谱分析，然后在时域用数字信号处理方法将信号加入3种混响，再分析其频谱，并与原始信号进行比较。二二. .设计思路设计思路 1、利用Windows下的录音机或其他软件，录制一段语音信号，时间控制在3s左右，并对录制的信号进行采样 2、语音信号的频谱分析，画出采样后的时域波形和频谱图 3、将信号加入延时和混响，再分析其时域波形和频谱图，并与原始信号频谱进行比较 4、设计几种特殊类型的混响器：无限个回声混响器，多重回声混响器，全通结构的混响器，对采样后的信号进行混响处理，画出处理后的波形图和频谱图

2、5、分析处理后信号的波形图和频谱图变化三三. .设计步骤设计步骤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、三种混响器的设计，并对信号进行处理、三种混响器的设计，并对信号进行处理（1 1）无限回声混响器）无限回声混

8、响器a=0.05; %a取小于等于1Bz=0,0,0,0,0,0,0,0,0,0,1; %分子的系数Az=1,0,0,0,0,0,0,0,0,0,-a; %分母的系数yy1=filter(Bz,Az,z1); %滤波器进行滤波YY1=fft(yy1,32000); %经无限回声滤波器后的信号做32000点的FFT变换（2 2）多重回声混响器）多重回声混响器a=0.05; %a取小于等于1N=5Bz1=1,0,0,0,0,0,0,0,0,0,-0.5N %分子的系数Az1=1,0,0,0,0,0,0,0,0,0,-0.5; %分母的系数yy2=filter(Bz1,Az1,z1); %滤波器进行

9、滤波YY2=fft(yy2,32000); %经多重回声滤波器后的信号做32000点的FFT变换（3 3）全通结构的混响器）全通结构的混响器a=0.05; %a取小于等于1Bz1=a,0,0,0,0,0,0,0,0,0,1; %分子的系数Az1=1,0,0,0,0,0,0,0,0,0,a; %分母的系数yy3=filter(Bz1,Az1,z1); %滤波器进行滤波YY3=fft(yy3,32000); %经全通结构的混响器后的信号做32000点的FFT变sound(yy1,fs,bits);sound(yy2,fs,bits);sound(yy3,fs,bits);（4 4）画出经混响器处理

10、后信号的时域波形和频谱图）画出经混响器处理后信号的时域波形和频谱图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(全通结构滤波器的时域图) %全通结构的混响器时域波形subplot(2,1,2);plot(n1(1:1000),YY3(1:1000); title(全通结构滤波器的频谱图) %全通结构的混响器频谱图四四. .运行结果运行结果1.1.原始信号时域波形和频谱图原始信号时域波形和频谱图2 2、采样后信号的时域波形和频谱图、采样后信号的时域波形和频谱图3 3、对采样后信号延时的时域波形和频谱图、对采样后信号延时的时域波形和频谱图4 4、混响后的时域波形和频谱图、混响后的时域波形和频谱图5 5、原始信号的时域波形