clear
clc
NumberOfUnits=6;

%?Asynchronous?machine?data
Pnom=1.5e6/0.9*NumberOfUnits;?%?Nominal?power?（VA）
Vnom=575;???????%?Line-Line?voltage?（Vrms）
Fnom=60;????????%Hz
Rs=0.00706;?????%pu
Lls=0.171;??????%pu
Rr=0.005;???????%pu
Llr=0.156;??????%pu
Lm=2.9;?????????%pu
H=5.04;?????????%Inertia?constant?（s）
F=0.01;?????????%Friction?factor?（pu）
p=3;????????????%Number?of?pairs?of?poles


Ls=Lls+Lm;
Lr=Llr+Lm;
sig=1-（Lm^2/（Ls*Lr））;
%PWM?data
PWM_freq=27*Fnom;???%Not?used?with?average?model

%DC?link
Vdc_nom=1200;???????%Volts
C_DClink=10000e-6*NumberOfUnits;?%Farads

%Choke?data
R_RL=0.30/100;??%pu
L_RL=0.30;??????%pu

%Turbine?data
Pmec=1.5e6*NumberOfUnits;??%?Nominal?power?（W）
power_C=0.73;???%pu
speed_C=1.2;????%pu

%Tracking?characteristic?speeds:?[speed_A?...?speed_D]
speed_ABCD=[0.7?0.71?1.2?1.21];
speed_A=speed_ABCD（1）;?%pu?of?synchronous?speed
speed_B=speed_ABCD（2）;?%pu?of?synchronous?speed
speed_C=speed_ABCD（3）;?%pu?of?synchronous?speed
speed_D=speed_ABCD（4）;

power_A=0;?%pu?of?mechanical?power
power_B=power_C*（speed_B/speed_C）^3;?%pu?of?mechanical?power
power_D=1;



Kp_grid_side_cur_reg=2.5;
Ki_grid_side_cur_reg=500;
Kp_dc_reg=0.002;
Ki_dc_reg=0.05;

Kp_Q=0.05;
Ki_Q=5;
Kp_rotor_side_cur_reg=0.3;
Ki_rotor_side_cur_reg=8;

%Imax_grid_conv=Pmax_pu/Vnom_pu=Pmax_pu/1;
Pmax=500;
Imax_grid_conv=Pmax;




%Ts_Control=1.000e-04;
%Ts_Power=5.0e-06;
%B=0;

%%%%%??????????????//Pnom=Pmec;
%%
	%Wind?turbine?typical?characteristic?data
????c2=116;c3=0.4;c4=5;c5=21;
????%Computation?of?c1?and?c6?for?lambda_nom?et?cpmax
????Pelec_base=Pnom;
????wind_base=11;
????P_wind_base=power_C;
????speed_nom=speed_C;
????pitch_angle=0;
????cp_nom=0.48;
????lambda_nom=8.1;
????%cp_nom=0.4382;
????%lambda_nom=6.325;
????cp_max=cp_nom;

????lambda_i_nom=1/（1/lambda_nom-0.035）;
????k_nom=-（c2*c5/lambda_i_nom-c4*c5-c2）*exp（-c5/lambda_i_nom）/lambda_nom^2;
????c1=cp_max/（（c2/lambda_i_nom-c4）*exp（-c5/lambda_i_nom）+k_nom*lambda_nom）;
????c6=k_nom*c1;
????

	wr_elec_pu=1e-3:0.01:1.5;
????wr_mec_pu=wr_elec_pu/speed_nom;
????wind_max=ceil（（1.2/P_wind_base）^（1/3）*wind_base）;
????wind_min=round（（（0.1）/P_wind_base）^（1/3）*wind_base）;
????wind1=wind_min:（wind_base-wind_min）/5:wind_base;
????wind2=wind_base+（wind_base-wind_min）/5:（wind_base-wind_min）/5:wind_max;
????wind_pu=[wind1?wind2]/wind_base;
????
???	fig_num=figure;
	clf（fig_num）;?
????for?k=1:length（wind_pu）
????????for?i=1:length（wr_mec_pu）;
????????????lambda_pu=wr_mec_pu（i）/wind_pu（k）;
????????????lambda=lambda_pu*lambda_nom;
????????????lambda_i=1/（1/（lambda+0.08*pitch_angle）-0.035/（pitch_angle^3+1））;
????????????cp=c1*（c2/lambda_i-c3*pitch_angle-c4）*exp（-c5/lambda_i）+c6*lambda;
????????????cp_pu=cp/cp_nom;
????????????Pwind_pu（ki）=wind_pu（k）^3*cp_pu*P_wind_base;
????????end
????end
????for?k=1:length（wind_pu）
????????plot（wr_elec_puPwind_pu（k:））
????????hold?on
????end
????
????plot（[speed_Dspeed_D0

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????5093??2010-09-27?10:02??Initialisation.m

?????文件?????836552??2010-11-12?17:29??IMC_cb.mdl

-----------?---------??----------?-----??----

???????????????841645????????????????????2