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

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); 3. 语音信号数据转换clear; close all; clc; [x,FS] = audioread(guxiang1.wav); [x,NBITS] = audioread(guxiang1.wav); % 将 WAV 文件转换成变量 N=length(x); % 计算数据点数 % 不是偶数点化成偶数点 if mod(N,2)==0; N=N; else x(N)=[]; N=N-1; end; % 原信号波形频谱分析 tx=(0:N-1)/FS; % 计算原信号数据点时刻 subplot(3,2,1);plot(tx,x); % 绘制原信号波形 title( 原信号波形图 ); % 加标题 xf=fft(x); % 求原信号频谱 fx=(0:N/2)*FS/N; % 确定频谱图频率刻度 subplot(3,2,2);plot(fx,abs(xf(1:N/2+1))); % 绘制原信号频谱 title( 原信号频谱图 ); % 加标题 % 实现数据抽取 k=[1:N/2]; % 确定抽取位置 y=x(2*k); % 实现抽取后的数据 M=length(y); % 计算抽取后数据点数 % 抽取数据在原采样频率 FS 下的波形频谱分析 ty=(