数字信号处理课设报告-基于matlab的数字音效处理器

数字信号处理课设报告-基于matlab的数字音效处理器内容摘要：

y3t=39。 y_3(t)39。 figure。 subplot(2,1,1)。 t=n*T。 plot(t,y3)。 xlabel(39。 t/s39。 )。 ylabel(y3t)。 axis([0,t(end),min(y3),*max(y3)]) subplot(2,1,2)。 plot(f,abs(fft(y3)))。 %播放声音 % sound(x,fs)。 sound(y1,fs)。 %低通 sound(y2,fs)。 %高通 sound(y3,fs)。 %带通 %变速 % w= M=w*fs。 %相乘的数为大于一的数为快放，小于一的为慢放 sound(x,M)。 提高部分 : %对声音信号放大和衰减 % %放大 % figure o=10 for xa=1:16384。 p(xa,1)=x(xa,1)*xa*o。 end subplot(2,1,1)。 plot(x)。 xlabel(39。 原始波形 39。 ) subplot(2,1,2)。 plot(n/Fs,p)。 xlabel(39。 变换波形 39。 ) sound(20*p,fs)。 %衰减 % figure o1= for xa=1:16384。 p(xa,1)=x(xa,1)/xa*o1。 end subplot(2,1,1)。 plot(x)。 5 xlabel(39。 原始波形 39。 ) subplot(2,1,2)。 plot(n/Fs,p)。 xlabel(39。 变换波形 39。 ) sound(20*p,fs) %%%倒放 figure w1=x(1:2350)。 w2=x(2351:3850)。 w3=x(3851:4850)。 w4=x(4851:5550)。 w5=x(5551:8550)。 w6=x(10000:10150)。 w7=x(10151:12500)。 w8=x(12501:16384)。 y=[w8。 w7。 w6。 w5。 w4。 w3。 w2。 w1]。 sound(30*y,fs)。 %%混频 [x5,fs5,nbits]=wavread(39。 man39。 ,16384)。 %读声音文件 c=10*x+*x5。 plot(c)； sound(c)。 %%%%%%%回声 figure z=[zeros(1000,1)。 x]。 %声音延时 x1=[x。 zeros(1000,1)]。 %使原声音长度与延时后相等 y1=x1+z。 plot(y1)。 sound(10*y1,fs)。 %%%%%%%%%%男女声变换 %%%%%%%%%%%%%%%%%%%%%% %%%%%%男声变女声 [y,fs]=wavread(39。 man39。 )。 % save Y1 %luyin p=fft(y)。 xaa=p39。 N=1000。 pa=[zeros(1,N),xaa(1:22020),zeros(1,N)]。 pu=pa39。 Y1=3*real(ifft(pu))。 plot(Y1) sound(20*Y1,fs)。 %%%%%%%%女声变男声 [y,fs]=wavread(39。 bb39。 )。 % save Y1 %luyin p=fft(y)。 xaa=p39。 6 xab(1:17797)=xaa(1:17797)。 %取一半 N=3500。 pa=[xab(N:17797),。