%?a_lsf_conversion.m
function?lsf?=?a_lsf_conversion（a）
%如果a不是實數，輸出錯誤信息：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
%判斷線性預測多項式的根是否都在單位圓內，如果不在，則輸出錯誤信息
if?（max（abs（roots（a）））?>=?1.0）
????error（‘The?polynomial?must?have?all?roots?inside?of?the?unit?circle.‘）;
end
%?求對稱和反對稱多項式的系數
?p?=?length（a）-1;??%?求對稱和反對稱多項式的階次
a1?=?[a;0];????????%給行矩陣a再增加一個元素為0的行??????
a2?=?a1（end:-1:1）;?%a2的第一行為a1的最后一行，最后一行為a1的第一行
P1?=?a1+a2;????????%?求對稱多項式的系數
Q1?=?a1-a2;????????%?求反對稱多項式的系數?
%如果階次p為偶數次，從P1取掉實數根z?=-1，從Q1取掉實數根z?=1
%如果階次為奇數次，從Q1取掉實數根z?=?1及z?=-1
if?rem（p2）??%?求解p除以2的余數，如果p為奇數次，余數為1，否則為0
????Q?=?deconv（Q1[1?0?-1]）;%奇數階次，從Q1取掉實數根z?=?1及z?=-1
????P?=?P1;
else??????????%?p為偶數階次執行下面操作
????Q?=?deconv（Q1[1?-1]）;%從Q1取掉實數根z?=?1
????P?=?deconv（P1[1??1]）;%從P1取掉實數根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

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

?????文件???????1537??2010-01-13?14:58??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\a_lsf_conversion.m

?????文件?????223941??2007-11-06?17:25??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\beijing.txt

?????文件??????79564??2010-01-13?11:58??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\beijing.wav

?????文件???????3511??2010-01-14?17:02??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter11_1.m

?????文件????????679??2010-01-13?22:34??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter2_1.m

?????文件????????853??2010-08-17?16:42??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter2_2.asv

?????文件????????853??2010-08-17?16:43??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter2_2.m

?????文件???????1551??2010-01-12?15:57??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_1.m

?????文件???????2909??2010-01-13?09:59??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_10.m

?????文件????????500??2010-01-13?23:16??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_11.m

?????文件???????1357??2010-01-13?10:33??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_12.m

?????文件???????1416??2010-01-13?22:41??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_2.m

?????文件????????988??2010-01-13?22:44??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_3.m

?????文件???????1386??2010-01-13?08:45??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_4.m

?????文件????????778??2010-01-13?09:02??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_5.m

?????文件????????445??2010-01-13?22:56??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_6.m

?????文件???????1123??2010-05-02?19:03??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_7.m

?????文件???????1235??2010-01-13?09:26??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_8.m

?????文件???????1107??2010-01-13?09:39??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter3_9.m

?????文件???????4774??2010-01-13?10:46??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter4_1.m

?????文件???????2261??2010-05-02?18:40??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter4_2.m

?????文件???????2099??2010-05-02?19:02??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter4_3.m

?????文件???????1062??2010-01-13?23:44??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter4_4.m

?????文件????????854??2010-01-13?11:57??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter5_1.m

?????文件???????1435??2010-01-13?12:03??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter5_2.m

?????文件???????1778??2010-01-14?00:50??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter5_3.m

?????文件???????1977??2010-01-13?15:38??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter6_1.m

?????文件????????670??2010-01-13?15:51??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter6_2.m

?????文件???????1279??2010-01-15?17:00??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter6_3.m

?????文件????????929??2010-01-13?16:19??數字語音處理及MATLAB仿真%5B張雪英%5D%5B程序源代碼%5D\chapter6_4.m

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