zhangxueying.zip

zhangxueying.zip

% 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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           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个文件信息