%?
clear?all;?clc;?close?all;

fle=‘snn27.wav‘;????????????????????????????%?指定文件名
[xxfs]=wavread（fle）;???????????????????????%?讀入一幀語音信號
u=filter（[1?-.99]1xx）;????????????????????%?預加重
wlen=length（u）;?????????????????????????????%?幀長
p=12;???????????????????????????????????????%?LPC階數
a=lpc（up）;?????????????????????????????????%?求出LPC系數
U=lpcar2pf（a255）;??????????????????????????%?由LPC系數求出功率譜曲線
freq=（0:256）*fs/512;????????????????????????%?頻率刻度
df=fs/512;??????????????????????????????????%?頻率分辨率
U_log=10*log10（U）;??????????????????????????%?功率譜分貝值
subplot?211;?plot（u‘k‘）;???????????????????%?作圖
axis（[0?wlen?-0.5?0.5]）;
title（‘預加重波形‘）;
xlabel（‘樣點數‘）;?ylabel（‘幅值‘）
subplot?212;?plot（freqU_log‘k‘）;
title（‘聲道傳遞函數功率譜曲線‘）;
xlabel（‘頻率/Hz‘）;?ylabel（‘幅值/dB‘）;

n_frmnt=4;??????????????????????????????????%?取四個共振峰
const=fs/（2*pi）;????????????????????????????%?常數??
rts=roots（a）;???????????????????????????????%?求根
k=1;????????????????????????????????????????%?初始化
yf?=?[];
bandw=[];
for?i=1:length（a）-1?????????????????????
????re=real（rts（i））;????????????????????????%?取根之實部
????im=imag（rts（i））;????????????????????????%?取根之虛部
????formn=const*atan2（imre）;???????????????%?計算共振峰頻率
????bw=-2*const*log（abs（rts（i）））;???????????%?計算帶寬
????
????if?formn>150?&?bw?<700?&?formn????????yf（k）=formn;
????????bandw（k）=bw;
????????k=k+1;
????end
end

[y?ind]=sort（yf）;??????????????????????????%?排序
bw=bandw（ind）;
F?=?[NaN?NaN?NaN?NaN];??????????????????????%?初始化
Bw?=?[NaN?NaN?NaN?NaN];
F（1:min（n_frmntlength（y）））?=?y（1:min（n_frmntlength（y）））;???%?輸出最多四個
Bw（1:min（n_frmntlength（y）））?=?bw（1:min（n_frmntlength（y）））;?%?輸出最多四個
F0?=?F（:）;??????????????????????????????????%?按列輸出
Bw?=?Bw（:）;
p1=length（F0）;??????????????????????????????%?在共振峰處畫線
for?k=1?:?p1
????m=floor（F0（k）/df）;
????P（k）=U_log（m+1）;
????line（[F0（k）?F0（k）][-10?P（k）]‘color‘‘k‘‘linestyle‘‘-.‘）;
end
fprintf（‘F0=%5.2f???%5.2f???%5.2f???%5.2f\n‘F0）;
fprintf（‘Bw=%5.2f???%5.2f???%5.2f???%5.2f\n‘Bw）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-10-10?15:39??源碼\
?????目錄???????????0??2016-04-15?20:40??源碼\9550\
?????目錄???????????0??2016-04-15?19:53??源碼\9550\pro\
?????目錄???????????0??2016-04-15?20:17??源碼\9550\pro\basic_tbx\
?????文件????????1050??2013-10-23?04:28??源碼\9550\pro\basic_tbx\add_noisedata.m
?????文件????????1010??2013-10-23?04:31??源碼\9550\pro\basic_tbx\add_noisefile.m
?????文件????????2895??2009-11-02?10:08??源碼\9550\pro\basic_tbx\enframe.m
?????文件????????4242??2010-10-28?23:41??源碼\9550\pro\basic_tbx\findpeaks.m
?????文件?????????753??2013-06-12?09:17??源碼\9550\pro\basic_tbx\findSegment.m
?????文件?????????110??2012-06-12?21:40??源碼\9550\pro\basic_tbx\frame2time.m
?????文件?????????234??2004-05-27?03:14??源碼\9550\pro\basic_tbx\freqz_m.m
?????文件????????2535??2010-08-01?21:41??源碼\9550\pro\basic_tbx\frq2mel.m
?????文件?????????645??2013-05-02?05:20??源碼\9550\pro\basic_tbx\Gnoisegen.m
?????文件?????????136??2004-05-27?03:05??源碼\9550\pro\basic_tbx\ideal_lp.m
?????文件????????2931??1999-10-30?05:01??源碼\9550\pro\basic_tbx\instfreq.m
?????文件?????????564??2013-10-24?03:54??源碼\9550\pro\basic_tbx\linsmoothm.m
?????文件????????1315??2007-05-04?21:01??源碼\9550\pro\basic_tbx\lpcar2ff.m
?????文件????????1441??2007-05-04?21:01??源碼\9550\pro\basic_tbx\lpcar2pf.m
?????文件????????1307??2007-05-04?21:01??源碼\9550\pro\basic_tbx\lpcar2zz.m
?????文件????????2525??2010-08-01?21:41??源碼\9550\pro\basic_tbx\mel2frq.m
?????文件???????11383??2010-01-03?09:58??源碼\9550\pro\basic_tbx\melbankm.m
?????文件????????1585??2013-10-23?23:52??源碼\9550\pro\basic_tbx\Mtmpsd_ssb.m
?????文件????????1663??2008-01-14?10:30??源碼\9550\pro\basic_tbx\OverlapAdd2.m
?????文件????????1352??2013-10-24?03:58??源碼\9550\pro\basic_tbx\pitch_vad1.m
?????文件????????2358??2013-08-02?10:28??源碼\9550\pro\basic_tbx\pitch_vads.m
?????文件?????????398??2013-01-03?06:33??源碼\9550\pro\basic_tbx\pitfilterm1.m
?????文件?????????363??2013-06-08?23:35??源碼\9550\pro\basic_tbx\polydetrend.m
?????文件????????1888??2009-06-04?00:57??源碼\9550\pro\basic_tbx\rfft.m
?????文件????????2206??2016-04-15?19:43??源碼\9550\pro\LPC.m
?????文件?????????684??2012-11-23?00:30??源碼\9550\pro\snn27.wav
?????文件??????268288??2018-10-10?18:36??源碼\9550\基于LPC法的共振峰檢測20160415.doc
............此處省略6個文件信息