数字信号在语音处理上的应用程序.doc

 PAGE \* MERGEFORMAT 9 2016年（秋季）学期 现代数字信号处理实践仿真 分析报告 学院： 电气工程与信息工程学院 专业： 电力电子与电力传动 姓名： 房理想 学号： 162080804005 语音信号时域和频谱 [x,fs] = audioread(guxiang1.wav); [x,bits] = audioread(guxiang1.wav); fs=25600; %语音信号采样频率为25600 y1=fft(x,4096);%对信号做2048点FFT变换? f=fs*(0:2047)/4096; figure(1) magy1=abs(y1); angy1=angle(y1); subplot(3,1,1),plot(x);title(原始信号波形) subplot(3,1,2),plot(magy1);title(原始信号幅值) subplot(3,1,3),plot(angy1);title(原始信号相位) figure(2) freqz(x)%绘制原始语音信号的频率响应图 title(频率响应图) figure(3) plot(f,abs(y1(1:2048)));title(原始语音信号频谱) xlabel(Hz);ylabel(fudu); axis([0 4500 0 100]) 1.语音信号倒谱分析 x=audioread(guxiang1.wav) plot(x); [x,fs] = audioread(guxiang1.wav); p = audioplayer(x,fs); play(p); [x,nb]=audioread(guxiang1.wav); ms2=floor(fs*0.002); ms10=floor(fs*.01); ms20=floor(fs*0.02); ms30=floor(fs*0.03); % 画波形图 t=(0:length(x)- 1)/ fs; subplot(2,1,1); plot(t,x); legend( Waveform ); xlabel( Time (s) ); ylabel( Amplitude ); % 加窗 w=hamming(ms30); pos=1; fx=[]; while(pos+ms30)=length(x) y=x(pos:pos+ms30- 1);y=y- mean(y); Y=fft(y.*w); %作加窗信号的傅利叶变换 C=fft(log(abs(Y)+eps)); %倒谱分析 [c,fxval]=max(abs(C(ms2:ms20))); %找出 2ms(=500Hz)和 20m(=50Hz)之间的最大值 fx=[fx fs/ (ms2+fxval- 1)]; pos=pos+ms10; end; %画出倒谱图 t=(0:length(fx)- 1)*0.01 subplot(2,1,2); plot(t,fx); legend( FX Trace ); xlabel( Time(s) ); ylabel( Frequency(Hz) ); 2. 语音信号数字滤波 clear;clc; [X,FS] = audioread(guxiang1.wav); [X,NBITS] = audioread(guxiang1.wav); %利用 FDATool 设计一个 LowpassButterworth 滤波器 %指标 FS=22050Hz Fp=1000Hz Ap=1dB Fs=3000Hz As=20dB B =[0.0062,0.0187,0.0187,0.0062]; % 分子系数 A =[1,-2.1706,1.6517,-0.4312]; % 分母系数 Y=filter(B,A,X); % 实现数字滤波 t=(0:length(X)-1)/FS; % 计算数据时刻 subplot(2,2,1);plot(t,X); % 绘制原波形图 title( 原信号波形图 ); % 加标题 subplot(2,2,3);plot(t,Y); % 绘制滤波波形图 title( 滤波后波形图 ); % 加标题 xf=fft(X); % 作傅里叶变换求原频谱 yf=fft(Y); % 作傅里叶变换求滤波后频谱 fm=3000*length(xf)/FS; % 确定绘频谱图的上限频率 f=(0:fm)*FS/length(xf); % 确定绘频谱图的频率刻度 subplot(2,2,2);plot(