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

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:a1