مساعدة


(system) #1

[CENTER]السلام عليكم
رجاءا لمن لديه ادنى فكرة في هذا التمرين ان يتقدم لكي يشرح لي كيفية تطبيقه على الماتلاب
لديا قاعدة بيانات هذا التطبيق ولكن لا اعرف تطبيقها في الماتلاب الرجاء المساعدة وفقكم الله

global ng ng=4;%-Nombre de générateurs du systéme npv=3;%- Nombre de noeuds PV (code : 2) npq=9;%- Nombre de noeuds PQ (code : 0) ns=1;%- Numéro du noeud sluck bus % (noeud balancie ou de régulation)(code : 1) nbus=13;%- Nombre total de noeuds du systéme nbr=14;%- Nombre total de branches du systéme Hz=60;%-fréquence du systéme accel=1.4;%- Facteur d'accélération % (pour le calcul de lodad flow par la méthode de gauss saidel) accuracy=0.0001;%- Tolérance (accurucy) relative au load flow (méthode de guss saidel) basemva=900;%- Puissance de base du systéme maxiter=100;%- Nombre d'itérations maximal Nsys=1;%- Identificateur du systéme % DONNEES DES NOEUDS: % bus data format % bus: % col1 number % col2 voltage magnitude (pu) % col3 voltage angle (degree) % col4 p_gen (pu) % col5 p_gen (pu) % col6 p_load (pu) % col7 q_load (pu) % col8 G shunt (pu) % col9 B shunt (pu) % col10 bus_type  %    bus_type -1,swing bus %         -2,generateur bus (PV bus) %         -3,load bus (PQ bus) % col11 q_gen_max (pu) % col12 q_gen_min (pu) % col13 v_rated (kv) % col14 v_max pu % col15 v_min pu  % Pour charge nominale  bus=     [     1  1.03   0.0 7.00 1.85 0.00  0.00 0.00 0.00 1 99.0 -99.0 22.0  1.1  .9;     2  1.01   0.0 7.00 2.35 0.00  0.00 0.00 0.00 2 99.0 -99.0 22.0  1.1  .9;     3  1.03   0.0 7.19 1.76 0.00  0.00 0.00 0.00 2 99.0 -99.0 22.0  1.1  .9;     4  1.01   0.0 7.00 2.02 0.00  0.00 0.00 0.00 2 99.0 -99.0 22.0  1.1  .9;     5  1.0103 0.0 0.00 0.00 0.00  0.00 0.00 0.00 3 0.0    0.0 230.0 1.5  .5;     6  0.9781 0.0 0.00 2.00 0.00  0.00 0.00 0.00 3 0.0    0.0 500.0 1.5  .5;     7  0.95   0.0 0.00 0.00 9.76  1.00 0.00 0.00 3 0.0    0.0 115.0 1.05 .95;     8  0.9899 0.0 0.00 3.50 0.00  0.00 0.00 0.00 3 0.0    0.0 500.0 1.5  .5;     9  0.95   0.0 0.00 0.00 17.67 1.00 0.00 0.00 3 0.0    0.0 115.0 1.05 .95;     10 0.9879 0.0 0.00 0.00 0.00  0.00 0.00 0.00 3 99.0   0.0 230.0 1.5  .5;     11 1.0    0.0 0.00 0.00 0.00  0.00 0.00 1.60 3 0.0  -99.0 500.0 1.5  .5;     12 1.0125 0.0 0.00 0.00 0.00  0.00 0.00 0.00 3 0.0    0.0 230.0 1.5  .5;     13 0.9938 0.0 0.00 0.00 0.0   0.00 0.00 0.00 3 0.0    0.0 230.0 1.5  .5]; % line data format : from bus, to bus, resistance (pu), reactance (pu), % line charging (pu), tap ratio, tap phase, tapmax, tapmin, tapsize % operating condition 1-Tow lines between busses 6- 11 and 8- 11  linedata=[      1  5  0.0    0.0167 0.00   1.0 0. 0.  0.  0.;       2  10 0.0    0.0167 0.00   1.0 0. 0.  0.  0.;      5  10 0.0025 0.25 0.0437   1.0 0. 0.  0.  0.;      10 6  0.001  0.0100 0.0175 1.0 0. 0.  0.  0.;      6  7  0.0    0.005  0.00   1.0 0. 1.2 0.8 0.5;      6  11 0.011  0.110  0.1925 1.0 0. 0.  0.  0.;      6  11 0.011  0.110  0.1925 1.0 0. 0.  0.  0.;      11 8  0.011  0.11   0.1925 1.0 0. 0.  0.  0.;      11 8  0.011  0.11   0.1925 1.0 0. 0.  0.  0.;      8  9  0.0    0.005  0.00   1.0 0. 1.2 0.8 0.5;      8  13 0.001  0.01   0.0175 1.0 0. 0.  0.  0.;      13 12 0.0025 0.025  0.0437 1.0 0. 0.  0.  0.;      3  12 0.0    0.0167 0.0    1.0 0. 0.  0.  0.;      4  13 0.0    0.0167 0.0    1.0 0. 0.  0.  0.]; % Machine data format % 1.machine number, % 2.bus number, % 3.base mva, % 4.leakage reactance x_l(pu), % 5.resystance r_a(pu), % 6.d-axis synchronous reactance x_d (pu), % 7.d-axis transient reactance x'_d (pu), % 8.d-axis subtransient reactance x"_d (pu), % 9.d-axis open-circuit time constant T'_do (sec), % 10.d-axis open-circuit subtransient time constant T"_do (sec), % 11.q-axis synchronous reactance x_q (pu), % 12.q-axis transient reactance x'_q (pu), % 13.q-axis subtransient reactance x"_q (pu), % 14.q-axis open-circuit time constant T'_qo (sec), % 15.q-axis open-circuit subtransient time constant T"_qo (pu), % 16.inertia constant H (sec), % 17.M(s²/rad), % 18.D #,  mac_con=[     1 1  900 0.200 0.0025 1.8 0.30 0.25 8.00 0.03 1.7 0.55 0.25 0.4 0.05 6.5 0.1254 2;     2 2  900 0.200 0.0025 1.8 0.30 0.25 8.00 0.03 1.7 0.55 0.25 0.4 0.05 6.5 0.0340 2;     3 11 900 0.200 0.0025 1.8 0.30 0.25 8.00 0.03 1.7 0.55 0.25 0.4 0.05 6.5 0.0160 2;     4 12 900 0.200 0.0025 1.8 0.30 0.25 8.00 0.03 1.7 0.55 0.25 0.4 0.05 6.5 0.0053 2];  % Exiter data format  % column data % 1  exciter type  % 2  machine number % 3  input filter time constant  % 4  voltage regulator gain K_A  % 5  voltage regulator time constant T_A (sec) % 6  voltage regulator time constant T_B (sec) % 7  voltage regulator time constant T_C (sec) % 8  maximum voltage regulator output VR_max % 9  minimum voltage regulator output VR_min % 10 exciter constant K_e % 11 exciter time constant T_E % 12 E_1  % 13 saturation function S_E (E_1)  % 14 E_2 % 15 saturation function S_E (E_2) % 16 stabilizer gain K_f % 17 stabilizer time constant (T_f)  exc_con=[     0 1 0.01 200.0 0.5 0 0 10.0 -10.0 0 0 0 0 0 0;      0 2 0.01 200.0 0.5 0 0 10.0 -10.0 0 0 0 0 0 0;      0 3 0.01 200.0 0.5 0 0 10.0 -10.0 0 0 0 0 0 0;      0 4 0.01 200.0 0.5 0 0 10.0 -10.0 0 0 0 0 0 0];

[/center]