語音信號的采集與時頻域分析系統(tǒng)的設(shè)計-資料下載頁
2025-06-26 09:17本頁面
【正文】 N1)n Zm(k)=0。 for i=0:N1。 Zm(k)=Zm(k)+w(k+i)。 end j=j+1。 k=k+N/2。 %每次移動半個窗 end for w=1:j Q(w)=Zm(230*(w1)+1)/(2*N)。 %短時平均過零率 end subplot(2,1,1),plot(Q),xlabel(39。采樣點39。),ylabel(39。短時平均過零率39。)。grid。程序四：；；；；；；；到譜圖程序；；；；；；；；；；；；；[x,fs,nb]=wavread(39。F:\xinhaochuli\39。)。tt=floor(fs*)。tt1=floor(fs*.01)。tt2=floor(fs*)。tt3=floor(fs*)。%畫波形圖t=(0:length(x)1)/fs。subplot(2,1,1)。plot(t,x)。legend(39。Waveform39。 )。title(39。波形圖39。)。xlabel(39。 時間39。 )。ylabel(39。幅度39。 )。%加窗w=hamming(tt3)。pos=1。fx=[]。while(pos+tt3)=length(x)y=x(pos:pos+tt31)。y=ymean(y)。Y=fft(y.*w)。 %作加窗信號的傅利葉變換C=fft(log(abs(Y)+eps))。 %倒譜分析[c,fxval]=max(abs(C(tt:tt2)))。%找出 2ms(=500Hz)和 20m(=50Hz)之間的最大值fx=[fx fs/(tt+fxval1)]。pos=pos+tt1。end。%畫出倒譜圖t2=(0:length(fx)1)*。subplot(2,1,2)。plot(t2,fx)。legend(39。 FX Trace39。 )。title(39。倒譜圖39。)。ylabel(39。 頻率(Hz)39。 )。程序五：；；；；濾波器設(shè)計與濾波；；；；；；；；；；；；[z,fs,bits]=wavread(39。F:\xinhaochuli\39。)。y=awgn(z,50)。ws=*pi。wp=*pi。Rp=。As=50。[N,wn]=ellipord(wp/pi,ws/pi,Rp,As)。[b,a]=ellip(N,Rp,As,wn)。[Rows,Cols]=size(y)。y(:,1)=y(:,1) mean(y(:,1))。Y=abs(fft(y)+eps)。 %為防止出現(xiàn) 0 值Y(1)=Y(2)。L=Rows。w=0:(L1)。w=2*pi/ Rows*w。x=filter(b,a,y)。X=abs(fft(x)+eps)。 %為防止出現(xiàn) 0 值X(1)=X(2)。figure(1)。subplot(3,2,1)。plot(z)。title(39。 原始信號波形 39。 )。subplot (3,2,2)。plot (w/ pi,Y(1:L))。title(39。 原始信號頻譜 39。 )。 xlabel(39。w(單位 pi)39。 )。ylabel(39。 幅值 39。 )。subplot(3,2,3)。plot(y)。title(39。 濾波前信號波形 39。 )。subplot (3,2,4)。plot (w/ pi,Y(1:L))。title(39。 濾波前信號頻譜 39。 )。 xlabel(39。w(單位 pi)39。 )。ylabel(39。 幅值 39。 )。subplot(3, 2 ,5)。plot(x)。title(39。 濾波后信號波形 39。 )。subplot(3, 2 ,6)。plot(w/ pi,X(1:L))。title(39。 低通濾波后信號頻譜 39。 )。 xlabel(39。w(單位 pi)39。 )。ylabel(39。 幅值 39。 )。燕山大學(xué)課程設(shè)計評審意見表指導(dǎo)教師評語：成績： 指導(dǎo)教師： 年 月 日答辯小組評語：成績： 評閱人： 年 月 日課程設(shè)計總成績：答辯小組成員簽字： 年 月 日 第 頁 共 18 頁
研究報告相關(guān)推薦
鄂ICP備17016276號-1