基于倒譜的語(yǔ)音特性提取算法設(shè)計(jì)及其實(shí)現(xiàn)(參考版)
2025-06-22 12:45本頁(yè)面
【正文】 這期間身邊的同學(xué)也幫我一起解決了很多相關(guān)知識(shí)和疑問(wèn)。 致謝： 本論文的指導(dǎo)老師是西南大學(xué)電子信息工程學(xué)院的李太華老師。Fundamental Frequency F0= %\n39。F0=1/T0。[Cmax Cloc]=max(abs(z(41:320)))。r39。subplot(2,1,2)。)。title(39。%fft變換subplot(2,1,2)z2=fft(z11)。 zeros(1,10002*mcep) z1(mcep:1:2)39。Envelope Information in Quefrency Space39。plot(t(nn),abs(z1(nn)))。%加窗subplot(2,1,1)nn=1:mcep。)。 hold ontitle(39。)Y=fft(y)。xlabel(39。Amplitude39。z=rceps(y)39。 %MATLAB提供的倒頻譜函數(shù)plot(t(nn),abs(z(nn)))。dt=1/sf。t=0:1/sf:(nfft1)/sf。index2=index1+nfft1。index1=4000。39。end2 提取基因和共振峰的程序%% dpclear all。endif nargout 1 x=sparse(r,c,v)。m39。) v=1cos(v*pi/2)。endif any(w==39。 mn=b1+1。 c=[k2:k4 1:k3]。 k4=b4b1+1。 k2=b2b1+1。 fp=floor(pf)。 b4=min(fn2,ceil(bl(4)))1。 mx=fn2+1。 v=2*[ ones(1,b21) 1pf+fp pffp ones(1,fn2b31) ]。 r=[ones(1,b2) fp fp+1 p*ones(1,fn2b3)]。) pf=log((f0+(b2:b3)/n)/(f0+fl))/lr。if any(w==39。b2=ceil(bl(2))。lr=log((f0+fh)/(f0+fl))/(p+1)。 end endendf0=700/fs。 if nargin 5 fh=。 without even the implied warranty of% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the% GNU General Public License for more details.%% You can obtain a copy of the GNU General Public License from% ftp://% Free Software Foundation, Inc.,675 Mass Ave, Cambridge, MA 02139, USA.%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if nargin 6 w=39。 you can redistribute it and/or modify% it under the terms of the GNU General Public License as published by% the Free Software Foundation。%% To plot filterbanks . plot(melbankm(20,256,8000)39。% c=dct(z)。 z=log(x*(f(na:nb)).*conj(f(na:nb)))。 [x,na,nb]=melbankm(p,n,fs)。 f=fft(s)。ny39。ty39。y39。z39。m39。n39。t39。 f = f .* w(ones(nf,1),:)。f(:) = x(indf(:,ones(1,len))+inds(ones(nf,1),:))。indf= inc*(0:(nf1)).39。endnf = fix((nxlen+inc)/inc)。else len = nwin。nwin=length(win)。 either version 2 of the License, or% (at your option) any later version.%% This program is distributed in the hope that it will be useful,% but WITHOUT ANY WARRANTY。) function f=enframe(x,win,inc)%ENFRAME split signal up into (overlapping) frames: one per row. F=(X,WIN,INC)%% F = ENFRAME(X,LEN) splits the vector X(:) up into% frames. Each frame is of length LEN and occupies% one row of the output matrix. The last few frames of X% will be ignored if its length is not divisible by LEN.% It is an error if X is shorter than LEN.%% F = ENFRAME(X,LEN,INC) has frames beginning at increments of INC% The centre of frame I is X((I1)*INC+(LEN+1)/2) for I=1,2,...% The number of frames is fix((length(X)LEN+INC)/INC)%% F = ENFRAME(X,WINDOW) or ENFRAME(X,WINDOW,INC) multiplies% each frame by WINDOW(:)% Copyright (C) Mike Brookes 1997% Version: $Id: ,v 2006/06/22 19:07:50 dmb Exp $%% VOICEBOX is a MATLAB toolbox for speech processing.% Home page: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% This program is free software。title(39。幅值39。)。xlabel(39。A=size(ccc)。一維數(shù)組及其幅值的關(guān)系39。)。ylabel(39。MFCC39。plot(ccc_1)。%去除首尾兩幀，因?yàn)檫@兩幀的一階差分參數(shù)為0ccc = ccc(3:size(m,1)2,:)。enddtm = dtm / 3。end%求取差分系數(shù)dtm = zeros(size(m))。 c2 = c1.*w39。%fft快速傅立葉變換 t = t.^2。 .* hamming(256)。%對(duì)x 256點(diǎn)分為一幀% 計(jì)算每幀的MFCC參數(shù)for i=1:size(xx,1) y = xx(i,:)。xx=filter([1 ],1,xx)。w = w/max(w)。 dctcoef(k,:)=cos((2*n+1)*k*pi/(2*24))。bank=bank/max(bank(:))。)。bank=mel(24,256,fs,0,39。39。47(2)。參考文獻(xiàn)：[1] ：哈爾濱工業(yè)大學(xué)出版社，2000[2] ：上
試題試卷相關(guān)推薦
鄂ICP備17016276號(hào)-1