qingxiezixiefangcha.m

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

%**************************************************************************************************************************************************************
%   本程序为主程序,目的是运用统计分析法计算WM1000发动机故障诊断所布各个测点加速度响应的自协方差函数
%   
%   中国北方发动机研究所试验测试中心  苏志霄   2006年10月
%**************************************************************************************************************************************************************
%
m=13;                                                             %m取值为测点数目
%读入时间加速度数据,各测点在每个测试时间内时间数据列是一致的，因此统一标识
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';
%**************************************************************************************************************************************************************
%故障诊断统计分析法：计算各个测点加速度信号数据的自协方差函数
for i=1:m;
    figure(i);
    R=zixiefangcha(acc(i,:));
    length=round(2/3*max(size(R)));
    subplot(2,1,1);
    plot(t(1:length),R(1:length));set(gca,'FontSize',7);
    grid;
    xlabel('time/Sec');
    ylabel('Self-Convariance Coefficient');
    title(['The Self-Convariance Curve and its Frequency Spectrum of ',num2str(i),' point']);
    subplot(2,1,2);
    f=5000*(0:511)/1024;
    A=fft(R,1024);
    pu=sqrt(A.*conj(A))/1024;
    plot(f,pu(1:512));xlabel('Frequcncy/Hz');set(gca,'FontSize',7);
    hold on;
    [a n]=max(pu);
    plot(f(n),pu(n),'r.','MarkerSize',15,'EraseMode','none');
    N=num2str(n/512*2500,'%.1f');
    text(f(n)+30,pu(n),[N,'Hz'],'FontSize',8);
    hold off;
end;