%**************************************************************************************************************************************************************
%   本程序為主程序,目的是運用統計分析法計算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');