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

1、 2016年（秋季）学期现代数字信号处理实践仿真分析报告学院： 电气工程与信息工程学院 专业： 电力电子与电力传动 姓名： 房理想 学号： 162080804005 1. 语音信号时域和频谱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);

2、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');

3、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,

4、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(=5

5、0Hz)之间的最大值fx=fx fs/ (ms2+fxval- 1);pos=pos+ms10;end;%画出倒谱图t=(0:length(fx)- 1)*0.01subplot(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&#

6、39;);%利用 FDATool 设计一个 LowpassButterworth 滤波器%指标 FS=22050Hz Fp=1000Hz Ap=1dB Fs=3000Hz As=20dBB =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

7、);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(f,abs(xf(1:length(f); % 绘制原波形频谱图title(' 原信号频谱图 '); % 加标题subplot(2,2,4);plot(f,abs(yf(1:length(f);

8、 % 绘制滤波后频谱图title(' 滤波后信号频谱图 '); % 加标题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(

9、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=(0:M-1)/FS; % 计算数据点时刻subplot(3,2,3);plot(ty,y); % 绘制信号波形图title(' 原采样率下新波形图 '); % 加标题yf=fft(y); % 求频谱fy=(0:M/2)*FS/M; % 确定频谱图频率刻度subplot(3,2,4);plot(fy,abs(yf(1:M/2+1); % 绘制频谱图title(' 原采样率下新频谱图 '); % 加标题% 抽取数据在 FS/2 采样频率下的波形频谱分析tz=(0:M-1)/(FS/2); % 计算数据点时刻subplot(3,2,5);plot(tz,y); %