基于某matlab的数字音效处理器——数字信号处理课设报告材料

1、一：应用背景利用所学习的数字信号处理知识,自己动手制作一个有趣的音效处理系统,看看能不能完成声音的逐渐放大和逐渐衰减、看看能不能让自己的声音发生一些改变变得尖声尖气或粗声粗气、看看改变声音播放速度有什么方法等等,你还可以自己想想还有什么有趣的变化,可以通过我们已有的知识让它实现.作为课程设计,以下要求分为根本必做局部和提升必做局部,在提升局部你可以选择全部内容和局部内容,当然分数值是不一样.:、基于MATLAB数字音效处理器2,1:实现步骤根本要求描述40分1语音信号的采集2分要求利用Windows下的录音机,录制一段自己的话音,时间在5s内,存为*.WAV的文件.然后在Matlab软件平台下

2、,利用函数wavread对语音信号进行采样,记住采样频率和采样点2语音信号的频谱分析10分要求首先画出语音彳t号的时域波形；然后对语音号进行快速傅里叶变换,得到信号的频谱特性,分析基频.3设计数字滤波器和画出其频率响应10分给出各滤波器的性能指标:1低通滤波器性能指标fb=1000Hz,fc=1200Hz,As=100dB,Ap=1dB.2高通滤波器性能指标fc=4800Hz,fb=5000HzAs=100dB,Ap=1dB.3带通滤波器性能指标=100dB,Ap=1dB.fb1=1200Hz,fb2=3000Hz,fc1=1000Hz,fc2=3200Hz,As4用滤波器对信号进行滤波5分要

3、求学生用自己设计的各滤波器分别对采集的信号进行滤波,在Matlab中,FIR滤波器利用函数fftfilt对信号进行滤波,IIR滤波器利用函数filter对信号进行滤波.5比拟滤波前后语音信号的波形及频谱10分要求在一个窗口同时画出滤波前后的波形及频谱,做出分析.6回放语音信号1分在Matlab中,函数sound可以对声音进行回放.其调用格式：soundx,fs,b1s;可以感觉滤波前后的声音有变化.7实现慢录快放和快录慢放功能2分3提升局部要求选择60分内容8实现对声音信号放大和衰减功能10分9实现对录音内容倒播放30分10实现混音音效效果10分11实现回音音效效果10分12实现男女变声音效效

4、果30分13设计系统界面10分为了使编制的程序操作方便,要求有水平的学生,设计处理系统的用户界面.2.2 试验程序根底局部：clear;clc;closeall;%卖取声音信号%x,fs,nbits=wavread('bb',16384);烟声音文件N=length(x);n=0:N-1;X=fft(x);%傅里叶变换Fs=2*fs;%胎频T=1/Fs;f=n/N*Fs;喏巴点数转换成频率subplot(2,1,1);plot(n,x);湎出原声音信号ylabel('原声音信号)；xlabel('时间/s');subplot(2,1,2);plot(f,

5、abs(X);%画出原声音信号的幅度谱ylabel('语音的幅度谱)；xlabel('频率/Hz');%触波器设计%氐通滤波器%fp1=1000;fs1=1200;wp1=2*fp1/Fs;ws1=2*fs1/Fs;rp=1;as=100;N1,wp1=ellipord(wp1,ws1,rp,as);B,A=ellip(N1,rp,as,wp1);y1=filter(B,A,x);Y1=abs(fft(y1);%低通滤波器设计与实现绘图部%figure;freqz(B,A);y1t='y_1(t)'figure;subplot(2,1,1);t=n*T;

6、plot(t,y1);xlabel('t/s');ylabel(y1t);axis(0,t(end),min(y1),1.2*max(y1)subplot(2,1,2);plot(f,abs(fft(y1);%高通滤波器%fp2=4800;fs2=5000;wp2=2*fp2/Fs;ws2=2*fs2/Fs;rp=1;as=100;N2,wp2=ellipord(wp2,ws2,rp,as);B2,A2=ellip(N2,rp,as,wp2,'high');y2=filter(B2,A2,x);%高通滤波器设计与实现绘图局部figure;freqz(B2,A2)

7、;figure;y2t='y_2(t)'subplot(2,1,1);t=n*T;plot(t,y2);xlabel('t/s');ylabel(y2t);axis(0,t(end),min(y2),1.2*max(y2)subplot(2,1,2);plot(f,abs(fft(y2);%带通滤波器%fpl=1200;fpu=3000;fsl=1000;fsu=3200;烦定低通滤波器通带截止频率和阻带截止频率%十算椭圆低通模拟滤波器的阶数和通带边界频率%十算低通滤波器模拟滤波器系统函数系数%滤波器软件实现渔标范围%设定高通滤波器通带截止频率和阻带截止频率%计

8、算椭圆高通模拟滤波器的阶数和通带边界频率%计算高通滤波器模拟滤波器系统函数系数%滤波器软件实现wp3=2*fpl/Fs,2*fpu/Fs;ws3=2*fsl/Fs,2*fsu/Fs;rp=1;as=100;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),;figure;subplot(2,1,1)

9、;t=n*T;plot(t,y3);xlabel('t/s');ylabel(y3t);axis(0,t(end),min(y3),1.2*max(y3)subplot(2,1,2);plot(f,abs(fft(y3);%播放声音%sound(x,fs);sound(y1,fs);%低通sound(y2,fs);%高通sound(y3,fs);%带通.展速w=0.8M=w*fs;%相乘的数为大于一的数为快放,小于一的为慢放sound(x,M);提升局部：.寸声音信号放大和衰减%.做大figureo=10forxa=1:16384;p(xa,1)=x(xa,1)*xa*o;en

10、dsubplot(2,1,1);plot(x);xlabel('原始波形')subplot(2,1,2);plot(n/Fs,p);xlabel('变换波形')sound(20*p,fs);%t减%figureo1=0.9forxa=1:16384;p(xa,1)=x(xa,1)/xa*o1;endsubplot(2,1,1);plot(x);xlabel('原始波形')subplot(2,1,2);plot(n/Fs,p);xlabel('变换波形')sound(20*p,fs)%倒放figurew1=x(1:2350);w2=

11、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);%i频x5,fs5,nbits=wavread('man',16384);%读声音文件c=10*x+0.5*x5;plot(c);sound(c);%figurez=zeros(1000,1);x;%声音延时x1=x;zeros(1000,1);%使原声音长度与延时后相等

12、y1=x1+z;plot(y1);sound(10*y1,fs);%变换%变女声y,fs=wavread('man');%saveY1%luyinp=fft(y);xaa=p'N=1000;pa=zeros(1,N),xaa(1:22001),zeros(1,N);pu=pa'Y1=3*real(ifft(pu);plot(Y1)sound(20*Y1,fs);%!%男声y,fs=wavread('bb');%saveY1%luyinp=fft(y);xaa=p'xab(1:17797)=xaa(1:17797);%取一半N=3500;p

13、a=xab(N:17797),zeros(1,17797);pu=pa'Y2=3*real(ifft(pu);plot(Y2)sound(1000*Y1,fs);2.3 绘图局部：语音信号的采集号信音声原0.10.0500.05-0.1谱度幅的音语00.20.40.60.81时间/s201510500200040006000800010000120001400016000频率/Hz低通滤波器：1y00.10.20.30.40.50.60.70.80.91t/s高通滤波器:-3带通滤波器:NormalizedFrequeadysa(mple)0.)0.3y-0040200.10.20.30.40.50.60.70.80.91t/s放大:原始波形4变换波形衰减:原始波形-3x10变换波形倒放:混频:回声0.10.080.060.040.02-0.