MATLAB仿真[张雪英][程序源代码].rar

MATLAB仿真[张雪英][程序源代码].rar

% 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个文件信息