%**************************************************************************************************************************************************************
%   本程序為函數功能，目的是求取一列數據時間序列自回歸模型的特征值，以檢驗信號是否具有穩定性
%      
%   中國北方發動機研究所試驗測試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
function R=shijianxulie(data,n,err)   %f為采樣頻率,n為自回歸模型階數，err為時間序列模型的最終預測誤差
k=arcov(data,n);                      %該功能求解自回歸模型的前50階傳遞函數系數
R=poly2rc(k,err);                     %R為自回歸模型的特征值，模型穩定性的判據為所有特征值的絕對值小于1