%**************************************************************************************************************************************************************
%   本程序為主程序,目的是運用統計分析法計算WM1000發動機故障診斷所布各個測點加速度響應的相加平均結果
%   
%   中國北方發動機研究所試驗測試中心  蘇志霄   2006年10月
%**************************************************************************************************************************************************************
%
%輸入發動機轉速、采樣頻率
speed=1900;
fsamp=5000;
%讀入時間加速度數據,各測點在每個測試時間內時間數據列是一致的，因此統一標識
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);
subplot(2,2,1);R1=xiangjiapingjun(acc1,speed,fsamp);R1=R1';plot(t(1:size(R1,1)),R1);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 1st point');
subplot(2,2,2);R2=xiangjiapingjun(acc2,speed,fsamp);R2=R2';plot(t(1:size(R2,1)),R2);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 2nd point');
subplot(2,2,3);R3=xiangjiapingjun(acc3,speed,fsamp);R3=R3';plot(t(1:size(R3,1)),R3);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 3rd point');
subplot(2,2,4);R4=xiangjiapingjun(acc4,speed,fsamp);R4=R4';plot(t(1:size(R4,1)),R4);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 4th point');
figure(2);
subplot(2,2,1);R5=xiangjiapingjun(acc5,speed,fsamp);R5=R5';plot(t(1:size(R5,1)),R5);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 5th point');
subplot(2,2,2);R6=xiangjiapingjun(acc6,speed,fsamp);R6=R6';plot(t(1:size(R6,1)),R6);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 6th point');
subplot(2,2,3);R7=xiangjiapingjun(acc7,speed,fsamp);R7=R7';plot(t(1:size(R7,1)),R7);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 7th point');
subplot(2,2,4);R8=xiangjiapingjun(acc8,speed,fsamp);R8=R8';plot(t(1:size(R8,1)),R8);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 8th point');
figure(3);
subplot(2,2,1);R9=xiangjiapingjun(acc9,speed,fsamp);R9=R9';plot(t(1:size(R9,1)),R9);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 9th point');
subplot(2,2,2);R10=xiangjiapingjun(acc10,speed,fsamp);R10=R10';plot(t(1:size(R10,1)),R10);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 10th point');
subplot(2,2,3);R11=xiangjiapingjun(acc11,speed,fsamp);R11=R11';plot(t(1:size(R11,1)),R11);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 11th point');
subplot(2,2,4);R12=xiangjiapingjun(acc12,speed,fsamp);R12=R12';plot(t(1:size(R12,1)),R12);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 12th point');
figure(4);
subplot(2,2,1);R13=xiangjiapingjun(acc13,speed,fsamp);R13=R13';plot(t(1:size(R13,1)),R13);set(gca,'FontSize',7);grid;xlabel('time/Sec');ylabel('acceleration/g');title('The average response of 13th point');