qingxiexiangjiapingjun.m

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

%**************************************************************************************************************************************************************
%   本程序为主程序,目的是运用统计分析法计算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');