资源描述：
《语音信号处理课程考核报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、语音信号处理课程考核报告题目1:绘制同一段语咅的宽带语谱图和牢带语谱图题EI2:绘制语音的倒谱图学号:201113136068姓名:陈凯班级:电信11032014年6月23日评语:绘制同一段语音的宽带语谱图和窄带语谱图加深信号处理理论在语音信号中的应用；理解语谱图与时频分辨率的关系；使用函数specgram_ex3pl9.画出所采集语咅的宽带语谱图和窄带语谱图，分析语谱图与时频分辨率的关系；语谱图的时间分辨率和频率分辨率是由窗函数的特性决定的。时间分辨率高，可以看岀时间波形的每个周期及共振峰随时间的变化，但频率分辨
2、率低，不足以分辨由于激励所形成的细微结构，称为宽带语谱图；而窄带语谱图正好与Z相反。宽带语谱图可以获得较高的时间分辨率，反映频谱的快速时变过程；窄带语谱图可以获得较高的频率分辨率，反映频谱的精细结构。两者相结合，可以提供带两与语音特性相关的信息。语谱图上因其不同的灰度，形成不同的纹路，称Z为“声纹二声纹因人而异，因此可以在司法、安全等场合得到应用。程序代码:function[spec]=specgram_ex3p19(x,winlgh,frmlgh,sampling_rate);%specgram_hw3p20(x
3、,winlgh,frmlgh,sampling_rate)%%functiontocomputeaspectrogram%%x=inputwaveform%winlgh=windowlengthinsamples%frmlgh=framelengthinsamples%sampling_rate=samples/sec%PARAMETERS%sampling_rate=10000;%samplingratelfft=1024;%FFTlengthlfft2=lfft/2;%winlgh=200;%(128)wind
4、owlength(insamples)%frmlgh=10;%frameinterval(insamples)noverlap=winlgh-frmlgh;%x=x(1:4500);x=2.0*x/max(abs(x));etime=length(x)/sampling_rate;spec=abs(specgram(x,lftt,sampling_rate,winlgh,noverlap));subplot(211);plot((l:length(x))/sampling_rate,x)xlabelfTime(s)
5、');titleCSPEECH*);axis([0etime-2.52.5]);grid；subplot(212)imagesc(O:.010:etime,0:1000:(sampling_rate/2),log10(abs(spec)));axis('xy*)xlabel('Time(ms),),ylabel(,Frequency(Hz)*);titleC'SPECTROGRAM');结果图：宽带：specgram_ex3p19(y,50,10,11025)SPEECH4000鳥1000。(ZH)xoUBnba匸
6、窄带：specgram_ex3p19(y,200,IO,11025)SPEECH(z£guanbauz绘制语音的倒谱图复倒谱的是x(n)的Z变换取对数示的逆Z变换，其表达式如下：A%=Z-1[lnZ[x(n)]](7)倒谱c(n)定义为x(n)取Z变换后的幅度对数的逆Z变换，即cS)=h[lnlX(z)l](为在吋域上，语音产牛模型实际上是一个激励信号与声道冲激响应的卷积。对于浊音，激励信号可以由周期脉冲序列表示；对于清音，激励信号可以由随机噪声序列表示。声道系统相当于参数缓慢变化的零极点线性滤波器。这样经过同态处
7、理后，语音信号的复倒谱，激励信号的复倒谱，声道系统的复倒谱之间满足下面的关系：AAAs(n)=e(n)-^-v(n)⑼由于倒谱对应于复倒谱的偶部，因此倒谱与复倒谱具有同样的特点，很容易知道语音信号的倒谱，激励信号的倒谱以及声道系统的倒谱之间满足下血关系：Cs(")=cW)+cW)(10)浊音信号的倒谱屮存在看峰值，它的岀现位證等于该语音段的基音周期，而清音的倒谱中则不存在峰值。利用这个特点我们可以进行清浊音的判断，并且可以估计浊音的基音周期。clc；clear；tic,[y,fs]=wavread(1ah02•wa
8、v!)；L=length(y);fw=y.*hamming(L)；r=real(log(fft(fwzL)))pfw=cceps(fw)；rpfw=rceps(fw)；z=rpfw(1:30);p=pfw(31:L)logz=real(exp(fft(z,L)))；logp=real(fft(p))；plot(y)/title(1原始波形1)o.so2000
