qingxiejiegoucanshu.m

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

%**************************************************************************************************************************************************************
%   本程序为主程序,目的是运用Taylor变换对WM1000发动机的前13阶质量、刚度及阻尼等“伪”结构参数进行识别,以区别正常信号和故障信号的差异来进行故障诊断
%   
%   中国北方发动机研究所试验测试中心  苏志霄   2006年10月
%**************************************************************************************************************************************************************
%
speed=1500;                                                       %发动机转速
m=13;                                                             %m取值为测点数目
H=1;                                                              %Taylor变换中的常数，该数值对计算结果影响很小，故赋值为1
%**************************************************************************************************************************************************************
%读入时间加速度数据,各测点在每个测试时间内时间数据列是一致的，因此统一标识
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';
%将加速度数据由单精度格式转换为双精度格式
acc=double(acc);
%输入“伪”外力，假设仅在第11号传感器部位作用了一个转速频率的正弦外力
omg=2*pi*speed/60;
%F11=sin(omg*t);
%**************************************************************************************************************************************************************
%取时间范围从1/5开始为实际有效计算时间
t=t(round(len/5):len);length=max(size(t));acc=acc(:,round(len/5):len);%F11=acc(F11(len/5):len);
%Taylor变换中的两个系数变量kc(k)和yt(k)，由于使用时间区域均匀分割的识别方法，因此这里只需计算kc(0)和yt(0)
kc=(0+1)*(0+2)/H^2;yt=(0+1)/H;
%将计算时间分为3倍测点数目个小时间区域
dengfen=round(length/(3*m)-3*m);                        %为防止取整后会造成丢点问题而造成点数不够，故每段区域的点数适当减少(-3*m)
t=t(1:dengfen+1)';t=t-t(1);length=max(size(t));         %由于时间区域均匀分割，因此每子段内的时间点数是相同的，故统一代表
%计算Taylor变换中的系数矩阵
xishujuzhen=zeros(3*m);
for i=1:m;
    i,
    for j=1:3*m;
        for k=1:length;
            denfenxinhao(k)=acc(i,(j-1)*dengfen+k);
        end;
        xishujuzhen(i,j)=kc*taylor(t,denfenxinhao,2,H);
    end;
end;
for i=m+1:2*m;
    i,
    for j=1:3*m;
        for k=1:length;
            denfenxinhao(k)=acc(i-m,(j-1)*dengfen+k);
        end;
        xishujuzhen(i,j)=yt*taylor(t,denfenxinhao,1,H);
    end;
end;
for i=2*m+1:3*m;
    i,
    for j=1:3*m;
        for k=1:length;
            denfenxinhao(k)=acc(i-2*m,(j-1)*dengfen+k);
        end;
        xishujuzhen(i,j)=taylor(t,denfenxinhao,0,H);
    end;
end;
%输入“伪”外力的Taylor变换项
%FT=zeros()