qingxieshijianxulie.m

自己编写的比较全面的故障诊断matlab函数程序

%**************************************************************************************************************************************************************
%   本程序为主程序,目的是求取一列数据时间序列自回归模型的特征值，检验信号是否具有稳定性,以区别正常信号和故障信号的差异来进行故障诊断
%   
%   中国北方发动机研究所试验测试中心  苏志霄   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;