%**************************************************************************************************************************************************************
%   本程序為主程序,目的是求取一列數據時間序列自回歸模型的特征值，檢驗信號是否具有穩定性,以區別正常信號和故障信號的差異來進行故障診斷
%   
%   中國北方發動機研究所試驗測試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
m=13;                                                             % m取值為測點數目
n=40;                                                             % n為自回歸模型階數err為時間序列模型的最終預測誤差
err=0.1;                                                          % err為時間序列模型的最終預測誤差     
%**************************************************************************************************************************************************************
%讀入時間加速度數據,各測點在每個測試時間內時間數據列是一致的，因此統一標識
t=acceleration_X;
len=max(size(t));
acc=zeros(m,len);
acc(1,:)=acceleration_Y';acc(2,:)=acceleration2_Y';acc(3,:)=acceleration3_Y';acc(4,:)=acceleration4_Y';acc(5,:)=acceleration5_Y';acc(6,:)=acceleration6_Y';
acc(7,:)=NONE_Y';acc(8,:)=NONE2_Y';acc(9,:)=NONE3_Y';acc(10,:)=NONE4_Y';acc(11,:)=NONE5_Y';acc(12,:)=NONE6_Y';acc(13,:)=NONE7_Y';
%將加速度數據由單精度格式轉換為雙精度格式
acc=double(acc);
%將序列值換算為[0 2*pi]之內的序列
detangle=2*pi/n;
angle=0:detangle:(n-1)*detangle;
len=max(size(angle));
%*************************************************************************************************************************************************************
figure(1);
for i=1:9;
    subplot(3,3,i);
    set(gca,'FontSize',7);
    R=abs(shijianxulie(acc(i,:),n,err));
    R=R';
    polar(angle,R,'-..r');
    title(['The ',num2str(i),' point']);
    set(get(gca,'Title'),'Color','k','Fontsize',9);
end;
figure(2);
for i=10:13;
    subplot(3,3,i-9);
    set(gca,'FontSize',7);
    R=abs(shijianxulie(acc(i,:),n,err));
    R=R';
    polar(angle,R,'-..r');
    title(['The ' num2str(i) ' point']);
    set(get(gca,'Title'),'Color','k','Fontsize',9);
end;