Code covered by the BSD License
 Detection_pitch(t,x,u,fla...end csfunc
 FDI_DriveTrain_c(t,x,u,fl...
 FD_fault2_c(t,x,u,flag,A,...end csfunc
 FD_pitch1_a(t,x,u,flag,A,...end csfunc
 FD_pitch3_a(t,x,u,flag,A,...end csfunc
 FDIBenchMarkData.mclose all
 initialization.mFDI design initialization
 plot_fault_index.m
 FDIBenchMark_BLW_test_b_1025

View all files
Code of An Estimationbased Approach of Fault Detection and Isolation of the Wind TurbineBenchmark
by
Qi
17 Jan 2012
The matlab code of the solution of “FDI Design Competition for Wind Turbines" (winning 2nd place)

initialization.m 
%% FDI design initialization
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Np_vm=0;
% Np_omega_r=0;
% Np_omega_g=0;
% Np_tau_g=0;
% Np_P_g=0;
% Np_beta=0;
%pitch system enhanced
Ppb = [2; 2.5];
A_f = .5*Bpb*Cpb;
Apb_E = Apb+A_f;
%Apb_E = Apb;
Lpb = place(Apb_E',Cpb',Ppb);
Lpb = Lpb';
Apb_E_bar = Apb_ELpb*Cpb;
%pitch system original
Lpb_FD = place(Apb',Cpb',Ppb);
Lpb_FD = Lpb_FD';
Apb_bar = ApbLpb_FD*Cpb;
%drive train system
A_DT = Addt(2:3,2:3);
B_DT = [Bddt(2:3,2) Addt(2:3,1)];
C_DT = [1 0 ];
P_DT = [2; 2.5];
L_DT = place(A_DT',C_DT',P_DT);
L_DT = L_DT';
A_DT_bar = A_DTL_DT*C_DT;


Contact us