%F1F2F3F4

clear;

[yfsnbits]=wavread（‘F:\語音庫\實驗語音\dress6.wav‘）;

y=filter（[1?-0.975][1]y）;%預加重，系數等于0.975
all_length=length（y）;
figure（100）;set（gcf‘Color‘‘White‘）;
specgram（y（:1）1024fs4016）;
%本實驗中采樣率為8kHz幀長取256持續時間t=256/8000=31.75ms幀移為1/4幀長
N=256;%幀長
fm=N*1/4;?%幀移
franum=floor（（length（y）-N+fm）/fm）;%計算總幀數
D_lpcc=16;?p=13;?LSFArray=[zeros（1p）‘];
xl?=?0:fs/（2*N-1）:fs;
for?num=1:franum
????for?q=1:N
????????h（q）=0.54-0.46*cos（2*pi*q/N）;%漢明窗
????????x（q）=y（（num-1）*fm+q）*h（q）;
????end%加窗后的短時信號
????%subplot（211）;
????%plot（x）;axis（[1?256?0?0.045]）;
????x_fft?=?abs（fft（xN））.^2;
????A?=?lpc（xp）;
????A_fft?=?abs（fft（AN））.^2;
????lpcc=lpc_to_lpcc（Ap）;%lpc?to?lpcc?
????%figure（2）;plot（A_fft（1:128）‘.‘）;
????[A_fft_decind_a_fft]?=?sort（A_fft（1:128））;%找出極值點
????lsf?=?poly2lsf（A）;
????freq?=?lsf/2/pi*fs;
????freq_array?=?[0];
????LSFArray?=?[LSFArray?lsf];
????for?k?=?1?:?p-1
????????delta_lsf（k）?=?lsf（k+1）?-?lsf（k）;?%差分線譜對
????end
????[d_lsf_decind_lsf]?=?sort（delta_lsf）;%升序排序后的差分線譜對
????%figure（3）;plot（delta_lsf‘.‘）;
????%?找出相鄰頻帶
????[ind_lsf_dind_ind_lsf]?=?sort（ind_lsf（1:8））;%對差分線譜對的標號進行排序找出標號所對應的相鄰頻帶
????k?=?1;????????sum?=?1;??????index?=?1;
????for?i?=?k+1?:?8
????????if?abs（?ind_lsf_d（k）?-?ind_lsf_d（i）?）==sum%相鄰頻段
???????????sum?=?sum?+?1;
????????else
????????????if?sum?>?1?&&?i?==?k?+?sum%在相鄰的（一段元素）區域中選擇最小差分線譜對
???????????????if?sum?>?3%對長相鄰區域的處理??選取出兩個共振峰
???????????????????[minlsfmin_indx]?=?min（delta_lsf（ind_lsf_d（k:i-1）））;?
???????????????????%[maxlsfmax_indx]?=?max（delta_lsf（ind_lsf_d（k:i-1）））;
???????????????????[lsf1ind_lfs1]?=?sort（?delta_lsf（ind_lsf_d（k:i-1））?）;
???????????????????for?j?=?2?:?sum
???????????????????????if?abs（ind_lfs1（j）?-?ind_lfs1（1））?>?1
???????????????????????????max_indx?=?ind_lfs1（j）;
???????????????????????????break;
???????????????????????end
???????????????????end
???????????????????freq_array?=?[freq_array?ind_lsf_d（min_indx+k-1）];
???????????????????freq_array?=?[freq_array?ind_lsf_d（max_indx+k-1）];
???????????????????index?=?index?+?2;
???????????????????sum?=?1;
???????????????????k?=?i;
???????????????else???????%短相鄰區域的處理?選取出一個共振峰
???????????????????[minlsfmin_indx]?=?min（delta_lsf（ind_lsf_d（k:i-1）））;?
???????????????????freq_array?=?[freq_array?ind_lsf_d（min_indx+k-1）];
???????????????????index?=?index?+?1;
???????????????????sum?=?1;
???????????????????k?=?i;
???????????????end
????????????else
????????????????freq_array?=?[freq_array?ind_lsf_d（k）];
????????????????index?=?index?+?1;
????????????????k?=?k?+?1;
????????????end
????????end
????end
????if?sum?<=?3
????????if?k?==?8%對數組末尾元素的處理
????????????freq_array?=?[freq_array?ind_lsf_d（k）];
????????????index?=?index?+?1;
????????else
????????????[minlsfmin_indx]?=?min（delta_lsf（ind_lsf_d（k:i）））;?
????????????freq_array?=?[freq_array?ind_lsf_d（min_indx+k-1）];
????????????index?=?index?+?1;
????????end
????else?%對長相鄰區域的處理
????????[minlsfmin_indx

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4947??2008-01-04?15:51??共振峰\21840260lsp\lsp.m

?????文件???????9194??2007-12-24?14:56??共振峰\39709610formant~\formant~.c

?????文件???????2797??2009-07-07?15:27??共振峰\39709610formant~.rar

?????文件???????1141??2009-02-18?12:19??共振峰\5956457lpc\新建文件夾\LPCfrequency.m

?????文件??????44389??1994-04-16?22:41??共振峰\650375klatt.3.04\klatt.3.04\example1.par

?????文件?????172430??1994-04-16?22:41??共振峰\650375klatt.3.04\klatt.3.04\example2.par

?????文件??????18994??1994-04-17?02:00??共振峰\650375klatt.3.04\klatt.3.04\getopt.c

?????文件???????3950??1994-04-17?00:03??共振峰\650375klatt.3.04\klatt.3.04\getopt.h

?????文件???????1908??1994-05-08?02:47??共振峰\650375klatt.3.04\klatt.3.04\INSTALL

?????文件???????8358??1994-05-09?15:51??共振峰\650375klatt.3.04\klatt.3.04\klatt.c

?????文件??????11436??1994-05-08?02:47??共振峰\650375klatt.3.04\klatt.3.04\klatt.man

?????文件????????500??1994-05-09?15:59??共振峰\650375klatt.3.04\klatt.3.04\Makefile

?????文件??????24263??1994-05-09?15:56??共振峰\650375klatt.3.04\klatt.3.04\parwave.c

?????文件???????6565??1994-04-16?22:41??共振峰\650375klatt.3.04\klatt.3.04\parwave.h

?????文件????????227??1994-04-16?22:41??共振峰\650375klatt.3.04\klatt.3.04\proto.h

?????文件??????10806??1994-04-21?17:06??共振峰\650375klatt.3.04\klatt.3.04\README

?????文件?????312940??2003-05-14?21:28??共振峰\9927426speech_analysis\Feature

?????文件???????5905??2003-05-14?16:52??共振峰\9927426speech_analysis\Feature.dsp

?????文件????????539??2003-05-14?16:52??共振峰\9927426speech_analysis\Feature.dsw

?????文件?????435200??2009-07-06?20:20??共振峰\9927426speech_analysis\Feature.ncb

?????文件??????59904??2009-07-06?20:20??共振峰\9927426speech_analysis\Feature.opt

?????文件???????4106??2009-07-06?20:19??共振峰\9927426speech_analysis\Feature.plg

?????文件???????4984??2003-04-04?13:21??共振峰\9927426speech_analysis\FeatureLib\FeatureLib.dsp

?????文件???????1446??2003-04-04?12:38??共振峰\9927426speech_analysis\FeatureLib\FeatureLib.plg

?????文件???????9372??2003-05-15?08:22??共振峰\9927426speech_analysis\FE_basic.cpp

?????文件???????4935??2003-05-15?08:22??共振峰\9927426speech_analysis\FE_common.h

?????文件???????8985??2003-05-14?18:23??共振峰\9927426speech_analysis\FE_complex.h

?????文件???????2337??2003-05-14?16:51??共振峰\9927426speech_analysis\FE_delta.cpp

?????文件??????21618??2003-05-14?16:51??共振峰\9927426speech_analysis\FE_endpoint.cpp

?????文件???????4127??2003-05-14?16:51??共振峰\9927426speech_analysis\FE_endpoint.h

............此處省略74個文件信息