%?a_lsf_conversion.m
function?lsf?=?a_lsf_conversion（a）
%如果a不是實數(shù)，輸出錯誤信息：LSF不適用于復多項式的求解
if?~isreal（a）
????error（‘Line?spectral?frequencies?are?not?defined?for?complex?polynomials.‘）;
end
%?如果a（1）不等于1，將其歸一化為1
if?a（1）?~=?1.0
????a?=?a./a（1）;%將矩陣a的每個元素除以a（1）再賦給矩陣a
end
%判斷線性預測多項式的根是否都在單位圓內(nèi)，如果不在，則輸出錯誤信息
if?（max（abs（roots（a）））?>=?1.0）
????error（‘The?polynomial?must?have?all?roots?inside?of?the?unit?circle.‘）;
end
%?求對稱和反對稱多項式的系數(shù)
?p?=?length（a）-1;??%?求對稱和反對稱多項式的階次
a1?=?[a;0];????????%給行矩陣a再增加一個元素為0的行??????
a2?=?a1（end:-1:1）;?%a2的第一行為a1的最后一行，最后一行為a1的第一行
P1?=?a1+a2;????????%?求對稱多項式的系數(shù)
Q1?=?a1-a2;????????%?求反對稱多項式的系數(shù)?
%如果階次p為偶數(shù)次，從P1取掉實數(shù)根z?=-1，從Q1取掉實數(shù)根z?=1
%如果階次為奇數(shù)次，從Q1取掉實數(shù)根z?=?1及z?=-1
if?rem（p2）??%?求解p除以2的余數(shù)，如果p為奇數(shù)次，余數(shù)為1，否則為0
????Q?=?deconv（Q1[1?0?-1]）;%奇數(shù)階次，從Q1取掉實數(shù)根z?=?1及z?=-1
????P?=?P1;
else??????????%?p為偶數(shù)階次執(zhí)行下面操作
????Q?=?deconv（Q1[1?-1]）;%從Q1取掉實數(shù)根z?=?1
????P?=?deconv（P1[1??1]）;%從P1取掉實數(shù)根z?=-1
end
rP??=?roots（P）;%求去掉實根后的多項式P的根
rQ??=?roots（Q）;%求去掉實根后的多項式Q的根
aP??=?angle（rP（1:2:end））;%將多項式P的根轉換為角度（為歸一化角頻率）賦給ap
aQ??=?angle（rQ（1:2:end））;%將多項式Q的根轉換為角度（為歸一化角頻率）賦給aQ
lsf?=?sort（[aP;aQ]）;%%將P、Q的根（歸一化角頻率）按從小到大順序排序后即為lsf
%?EOF?a_lsf_conversion.m

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-03-16?23:12??源代碼\
?????文件????????1537??2010-01-13?14:58??源代碼\a_lsf_conversion.m
?????文件??????223941??2007-11-06?17:25??源代碼\beijing.txt
?????文件???????79564??2010-01-13?11:58??源代碼\beijing.wav
?????文件????????3511??2010-01-14?17:02??源代碼\chapter11_1.m
?????文件?????????679??2010-01-13?22:34??源代碼\chapter2_1.m
?????文件?????????853??2010-08-17?16:42??源代碼\chapter2_2.asv
?????文件?????????853??2010-08-17?16:43??源代碼\chapter2_2.m
?????文件????????1551??2010-01-12?15:57??源代碼\chapter3_1.m
?????文件????????2909??2010-01-13?09:59??源代碼\chapter3_10.m
?????文件?????????500??2010-01-13?23:16??源代碼\chapter3_11.m
?????文件????????1357??2010-01-13?10:33??源代碼\chapter3_12.m
?????文件????????1416??2010-01-13?22:41??源代碼\chapter3_2.m
?????文件?????????988??2010-01-13?22:44??源代碼\chapter3_3.m
?????文件????????1386??2010-01-13?08:45??源代碼\chapter3_4.m
?????文件?????????778??2010-01-13?09:02??源代碼\chapter3_5.m
?????文件?????????445??2010-01-13?22:56??源代碼\chapter3_6.m
?????文件????????1123??2010-05-02?19:03??源代碼\chapter3_7.m
?????文件????????1235??2010-01-13?09:26??源代碼\chapter3_8.m
?????文件????????1107??2010-01-13?09:39??源代碼\chapter3_9.m
?????文件????????4774??2010-01-13?10:46??源代碼\chapter4_1.m
?????文件????????2261??2010-05-02?18:40??源代碼\chapter4_2.m
?????文件????????2099??2010-05-02?19:02??源代碼\chapter4_3.m
?????文件????????1062??2010-01-13?23:44??源代碼\chapter4_4.m
?????文件?????????854??2010-01-13?11:57??源代碼\chapter5_1.m
?????文件????????1435??2010-01-13?12:03??源代碼\chapter5_2.m
?????文件????????1778??2010-01-14?00:50??源代碼\chapter5_3.m
?????文件????????1979??2014-03-10?11:55??源代碼\chapter6_1.m
?????文件?????????670??2010-01-13?15:51??源代碼\chapter6_2.m
?????文件????????1279??2010-01-15?17:00??源代碼\chapter6_3.m
?????文件?????????929??2010-01-13?16:19??源代碼\chapter6_4.m
............此處省略39個文件信息