风力发电机组故障诊断系统源程序.doc

(完整版)风力发电机组故障诊断系统源程序(完整版)风力发电机组故障诊断系统源程序(完整版)风力发电机组故障诊断系统源程序x=x（1：25600);N=length（x)；%读取时域信号fs=2560；%采样频率df=fs/N;%频率分辨率t=(0:N-1)/fs;y=（abs（fft（x）））。^2。/N；％功率谱y=y（1：1+N/2）;%取一半的有用信号f=(0:N/2)*2560/N；［d,e]=max（y（ceil(12/df）:ceil(33/df)））；％找到频域中12Hz到33Hz中最大值v=（ceil（12/df)+e-2)＊df；％最大值对应频率，即电机转频f0=27*v/(102＊5.25*5.25）；f1=27*v/(102＊5.25)；f2=f0；f2z=27＊24＊24＊v/(126＊39*126）；f4=27＊v/102;f5=f1;f5z=24*27*v/（126＊39);f7=f4；f8=v;ff3=27*v；ff2=27＊v/5.25；ff1=27＊v/（5.25*5.25）；ff11=v;ff12=f7；ff21=f4；ff22=f5z；ff31=f1；ff32=f2z；ffr=［ff32，ff31,ff22,ff21，ff12，ff11］;％根据电机转频计算齿轮箱各个齿轮转频和三级啮合频率a1=ceil（(ff1—5)/df)；a2=ceil((ff1+5）/df）;b1=ceil（（ff2—20)/df）;b2=ceil（(ff2+20）/df）;c1=ceil（(ff3—40）/df);c2=ceil（（ff3+40）/df）;％计算三级啮合频率附近区间的频谱对应点数范围gxxx1=f(a1：a2);gyyy1=y(a1:a2);gxxx2=f（b1:b2)；gyyy2=y（b1:b2）；gxxx3=f（c1：c2)；gyyy3=y（c1：c2）；％频谱和时域对应赋值,Labview中画出频谱中各啮合频率出放大图％变频程序x=x（1：25600）；N=length(x);%读取时域信号fs=2560;%采样频率df=fs/N;％频率分辨率t=(0:N-1)/fs;%时间序列y=(abs（fft（x)））。^2./N;%功率谱y=y（1:1+N/2）;%取一半的有用信号f=(0：N/2）*2560/N；［d,e］=max(y（ceil(12/df):ceil（33/df))）；%找到频域中12Hz到33Hz中最大值v=（ceil（12/df）+e-2）＊df；％最大值对应频率，即电机转频f0=27＊v/（102＊5。25*5.25）；f1=27*v/(102＊5。25）;f2=f0;f2z=27＊24＊24*v/(126＊39*126)；f4=27*v/102;f5=f1;f5z=24＊27*v/(126*39）；f7=f4;f8=v；ff3=27*v；ff2=27＊v/5。25;ff1=27＊v/（5.25*5.25);ff11=v；ff12=f7；ff21=f4；ff22=f5z;ff31=f1；ff32=f2z；ffr=[ff32,ff31，ff22,ff21,ff12,ff11］；％根据电机转频计算齿轮箱各个齿轮转频和三级啮合频率a11=ceil((ff1—5）/df)；a12=ceil（（ff1+5)/df）；a21=ceil（(ff2—5）/df）;a22=ceil(（ff2+5）/df);a31=ceil（(ff3-5）/df);a32=ceil(（ff3+5）/df);％找到三级啮合频率对应谱线y1=zeros(1，9）;fr1=zeros(1,9）；y2=zeros(1,9);fr2=zeros（1,9)；y3=zeros（1，9）;fr3=zeros(1,9)；[y1(1)，fr1（1)]=max(y（a11：a12））；fx1=(fr1(1)+a11-2)＊df;%找到频谱中第一级啮合频率值b111=（ceil(（fx1-f2z-0.2)/df))：(ceil（（fx1—f2z+0.2)/df))；b112=(ceil(（fx1+f2z-0.2)/df））：（ceil（（fx1+f2z+0.2）/df)）;b121=(ceil(（fx1-2*f2z-0。4)/df))：(ceil（(fx1-2＊f2z+0.4）/df）);b122=（ceil((fx1+2＊f2z—0.4）/df）)：（ceil(（fx1+2＊f2z+0.4）/df）);c111=（ceil((fx1—f1—0。5)/df))：(ceil(（fx1-f1+0。5）/df))；c112=（ceil（（fx1+f1—0。5)/df)）:(ceil（(fx1+f1+0。5）/df)）；c121=（ceil