数字信号处理课设报告-基于matlab的数字音效处理器Word文档格式.doc

1、sound（x,fs,bits）;可以感觉滤波前后的声音有变化。7）实现慢录快放和快录慢放功能（2分）3 提高部分要求（选择60分内容）8）实现对声音信号放大和衰减功能（10分）9）实现对录音内容倒播放（30分）10） 实现混音音效效果（10分）11）实现回音音效效果（10分）12）实现男女变声音效效果（30分）13）设计系统界面（10分）为了使编制的程序操作方便,要求有能力的学生,设计处理系统的用户界面。2.2试验程序基础部分：clear;clc;close all;%读取声音信号%x,fs,nbits=wavread（bb,16384）; %读声音文件 N=length（x）;n=0:N-

2、1;X= fft（x）; %傅里叶变换Fs=2*fs; %2倍频T=1/Fs;f=n/N*Fs; %把点数转换成频率subplot（2,1,1）;plot（n,x）; %画出原声音信号ylabel（原声音信号）;xlabel（时间/s subplot（2,1,2）;plot（f,abs（X）; %画出原声音信号的幅度谱语音的幅度谱频率/Hz% %滤波器设计% %低通滤波器% fp1=1000;fs1=1200; %设定低通滤波器通带截止频率和阻带截止频率 wp1=2*fp1/Fs; ws1=2*fs1/Fs;rp=1;as=100; N1,wp1=ellipord（wp1,ws1,rp,as）

3、; %计算椭圆低通模拟滤波器的阶数和通带边界频率 B,A=ellip（N1,rp,as,wp1）; %计算低通滤波器模拟滤波器系统函数系数y1=filter（B,A,x）; %滤波器软件实现Y1=abs（fft（y1）;% 低通滤波器设计与实现绘图部%figure;freqz（B,A）; y1t=y_1（t）;t=n*T;plot（t,y1）;t/sylabel（y1t）;axis（0,t（end）,min（y1）,1.2*max（y1）%坐标范围plot（f,abs（fft（y1）; %高通滤波器% fp2=4800;fs2=5000; %设定高通滤波器通带截止频率和阻带截止频率 wp2=2

4、*fp2/Fs; ws2=2*fs2/Fs; N2,wp2=ellipord（wp2,ws2,rp,as）; %计算椭圆高通模拟滤波器的阶数和通带边界频率 B2,A2=ellip（N2,rp,as,wp2,high %计算高通滤波器模拟滤波器系统函数系数 y2=filter（B2,A2,x）; %滤波器软件实现% 高通滤波器设计与实现绘图部分freqz（B2,A2）; y2t=y_2（t）plot（t,y2）;ylabel（y2t）;axis（0,t（end）,min（y2）,1.2*max（y2）plot（f,abs（fft（y2）;%带通滤波器% fpl=1200;fpu=3000;fsl

5、=1000;fsu=3200; wp3=2*fpl/Fs,2*fpu/Fs;ws3=2*fsl/Fs,2*fsu/Fs; N3,wp3=ellipord（wp3,ws3,rp,as）; %计算椭圆带通模拟滤波器的阶数和通带边界频率 B3,A3=ellip（N3,rp,as,wp3）; %计算带通滤波器模拟滤波器系统函数系数 y3=filter（B3,A3,x）; %滤波器软件实现% 带通滤波器设计与实现绘图部分 figure; freqz（B3,A3）;y3t=y_3（t）plot（t,y3）;ylabel（y3t）;axis（0,t（end）,min（y3）,1.2*max（y3）plot（

6、f,abs（fft（y3）; %播放声音%sound（x,fs）;sound（y1,fs）; %低通sound（y2,fs）; %高通sound（y3,fs）; %带通%变速%w=0.8 M=w*fs;%相乘的数为大于一的数为快放，小于一的为慢放sound（x,M）;提高部分:%对声音信号放大和衰减%放大%figureo=10for xa=1:16384; p（xa,1）=x（xa,1）*xa*o;end subplot（2,1,1）; plot（x）; xlabel（原始波形） subplot（2,1,2）; plot（n/Fs,p）;变换波形 sound（20*p,fs）; %衰减% fi

7、gureo1=0.9 p（xa,1）=x（xa,1）/xa*o1; sound（20*p,fs）%倒放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（manc=10*x+0.5*x5;plot（c）；sound（c）;%回声z=zeros（1000,1）;

8、x;%声音延时x1=x;zeros（1000,1）;%使原声音长度与延时后相等y1=x1+z;plot（y1）;sound（10*y1,fs）;%男女声变换% %男声变女声y,fs=wavread（% save Y1 %luyinp=fft（y）;xaa=pN=1000;pa=zeros（1,N）,xaa（1:22001）,zeros（1,N）;pu=paY1=3*real（ifft（pu）;plot（Y1）sound（20*Y1,fs）;%女声变男声xab（1:17797）=xaa（1:17797）;%取一半N=3500;pa=xab（N:17797）,zeros（1,17797）;Y2=3*real（ifft（pu）;plot（Y2）sound（1000*Y1,fs）;2.3 绘图部分：语音信号的采集 低通滤波器：高通滤波器：带通滤波器：放大：衰减：倒放：混频：回声：男声变女声：女声变男声：三 课设总结：参考文献美数字信号处理使用MATLABM.西安：西安交通大学出版社，2002. 11