qingxiegailvmidu.m

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

%**************************************************************************************************************************************************************
%   本程序为主程序,目的是运用统计分析法计算WM1000发动机故障诊断所布各个测点加速度响应幅值的概率密度函数
%   
%   中国北方发动机研究所试验测试中心  苏志霄   2006年10月
%**************************************************************************************************************************************************************
%
%读入时间加速度数据,各测点在每个测试时间内时间数据列是一致的，因此统一标识
t=acceleration_X;
acc1=acceleration_Y;acc1=double(acc1);    %将加速度数据由单精度格式转换为双精度格式
acc2=acceleration2_Y;acc2=double(acc2);
acc3=acceleration3_Y;acc3=double(acc3);
acc4=acceleration4_Y;acc4=double(acc4);
acc5=acceleration5_Y;acc5=double(acc5);
acc6=acceleration6_Y;acc6=double(acc6);
acc7=NONE_Y;acc7=double(acc7);
acc8=NONE2_Y;acc8=double(acc8);
acc9=NONE3_Y;acc9=double(acc9);
acc10=NONE4_Y;acc10=double(acc10);
acc11=NONE5_Y;acc11=double(acc11);
acc12=NONE6_Y;acc12=double(acc12);
acc13=NONE7_Y;acc13=double(acc13);
%**************************************************************************************************************************************************************
%故障诊断统计分析法：计算各个测点加速度信号幅值的概率密度函数
figure(1);
R=gailvmidu(acc1);subplot(2,2,1);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 1st point');
R=gailvmidu(acc2);subplot(2,2,2);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 2nd point');
R=gailvmidu(acc3);subplot(2,2,3);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 3rd point');
R=gailvmidu(acc4);subplot(2,2,4);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 4th point');
figure(2);
R=gailvmidu(acc5);subplot(2,2,1);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 5th point');
R=gailvmidu(acc6);subplot(2,2,2);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 6th point');
R=gailvmidu(acc7);subplot(2,2,3);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 7th point');
R=gailvmidu(acc8);subplot(2,2,4);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 8th point');
figure(3);
R=gailvmidu(acc9);subplot(2,2,1);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 9th point');
R=gailvmidu(acc10);subplot(2,2,2);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 10th point');
R=gailvmidu(acc11);subplot(2,2,3);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 11th point');
R=gailvmidu(acc12);subplot(2,2,4);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 12th point');
figure(4);
R=gailvmidu(acc13);subplot(2,2,1);hist(R,220);set(gca,'FontSize',7);grid;xlabel('Acc. Amp./g');ylabel('Prob. of 13th point');