qingxiejunzhifangcha.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';
%**************************************************************************************************************************************************************
%故障诊断统计分析法：计算各个测点加速度信号数据的均值和方差
result=zeros(2,m);
for i=1:m;
    R=junzhifangcha(acc(i,:));result(1,i)=R(1);result(2,i)=R(2);
end;