Karen Iter10
Karen Iter10
Karen Iter10
http://umss-ie-sistemasdepotencia2.jimdo.com/
PRACTICA 1
FLUJO DE CARGA RESUELTO POR EL METODO GAUSS-SEIDEL
Determinar las tensiones de barra, considera un error maxima (Epsilon) igual a 0.0001
Datos de Barra
Datos de Linea
Datos de Demanda
Factor de Potencia
Grupos
B1 B3
P (MW) P (MW)
grupo 7 grupo 8
40 45
30 25
20 20
90 90
FP grupo 7 FP grupo 8
0.95 0.95
0.97 0.97
0.88 0.88
Grupo 4
Determinar las tensiones de barra, considera un error maxima (Epsilon) igual a 0.0001
Datos de Barra
Datos de Linea
CARGA
Linea Z pu B100 Ysh/2 pu B100 1 15+4,93026157768295i 0,15+0,0493026157768295i
LIN 1-2 0,008+0,08i 0,0045i 2 55+13,7842993394408i 0,55+0,137842993394408i
LIN 2-3 0,01+0,1i 0,006i 3 20+10,7948564427617i 0,2+0,107948564427617i
P (MW)
Barra grupo 2 Q[MVAr] Barra FP grupo 2
1 15 4.930 1 0.95 alfa
2 55 13.784 2 0.97 1.20E+00
3 20 10.795 3 0.88
Total 90
1. Lectrura de datos
2. Conversión a por unidad
3. Formacion del YBUS
Ylin12= 1,23762376237624-12,3762376237624i
YLlin23= 0,99009900990099-9,9009900990099i
1,23762376237624-12,3717376237624i -1,23762376237624+12,3762376237624i
0 -0,99009900990099+9,9009900990099i 0,99009900990099-9,89499
1 SL
2 PV
3 PQ
5. Calculo de las potencias especificadas
1 SL Sesp1= Pgen1+Qgen1-0,15-0,0493026157768295i
2 PV Sesp2= 0,2+Qgen2 -0,55-0,137842993394408i Sinyesp2= -0,35+Qiny2'
3 PQ Sesp3= 0 -0,2-0,107948564427617i
V1= 1,2
VE= V2= 1,05
V3= 1
7. PROCESO ITERATIVO 1
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,33910891089109-23,3800641089109i
Y23*V3= -0,99009900990099+9,9009900990099i
Y21*V1+Y22*V2+Y23*V3= -0,13613861386139+1,3724111386139i
V2*= 1,05
V2*{ }= -0,14294554455446+1,44103169554459i
Qiny2= 0,14294554455446-1,44103169554459i
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,20705445544554-1,44103169554459i
Sesp2*= -0,20705445544554+1,44103169554459i
V2*= 1,05
Sesp2*/V2*= -0,197194719471943+1,3724111386139i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,33910891089109-23,3800641089109i
Y23*V3= -0,99009900990099+9,9009900990099i
[-Y21V1-Y22V2-Y23V3]= -2,47524752475248+24,7524752475248i
{-[ ]}= 2,47524752475248-24,7524752475248i
[ + ]= 2,27805280528054-23,3800641089109i
V2= 1,04972838608404-0,00271485895266421i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
1-alfa)*V2^k= -0.21
alfa*V2^k+1= 1,25967406330085-0,00325783074319705i
V2acel= 1,04967406330085-0,00325783074319705i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.049679119
|V2esp|/|V2acel|= 1,00030569448043
V2correg= 1,04999494286825-0,00325882664407342i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1
Sesp3*/V3*= -0,0018348623853211+0,000990353802088229i
Y31xV1= 0
Y32xV2= -1,00733334299754+10,3992160943827i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 1,00733334299754-10,3992160943827i
[ + ]= 1,00549848061222-10,3982257405806i
V3= 1,05050765360405-0,00349753832679661i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)*V3^k -0.2
alfa*V3^k+1= 1,26060918432486-0,00419704599215593i
V3acel= 1,06060918432486-0,00419704599215593i
|V3acel|= 1.060617489
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -5,05713175003031E-06-0,00325882664407342i
(|V3^(k+1) - V3^k| ) ( 0,06060918432486-0,00419704599215593i
7. PROCESO ITERATIVO 2
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,29950495049505+12,9950495049505i
Y22*V2= 2,26653423942423-23,3872112654611i
Y23*V3= -1,00855319247852+10,5052365240007i
Y21*V1+Y22*V2+Y23*V3= -0,0415239035493402+0,113074763490101i
V2*= 1,04999494286825+0,00325882664407342i
V2*{ }= -0,04396837978699+0,118592610627377i
Qiny2= 0,04396837978699-0,118592610627377i
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,30603162021301-0,118592610627377i
Sesp2*= -0,30603162021301+0,118592610627377i
V2*= 1,04999494286825+0,00325882664407342i
Sesp2*/V2*= -0,291106739974726+0,113849383599246i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,26653423942423-23,3872112654611i
Y23*V3= -1,00855319247852+10,5052365240007i
[-Y21V1-Y22V2-Y23V3]= -2,49370170733001+25,3567216725156i
{-[ ]}= 2,49370170733001-25,3567216725156i
[ + ]= 2,20259496735528-25,2428722889164i
V2= 1,13222244651204-0,0143569708132774i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,20999898857365+0,000651765328814684i
alfa*V2^k+1= 1,35866693581445-0,0172283649759329i
V2acel= 1,1486679472408-0,0165765996471182i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.148787551
|V2esp|/|V2acel|= 0,914007119373709
V2correg= 1,04989068157448-0,0151511300924737i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,06060918432486+0,00419704599215593i
Sesp3*/V3*= -0,188165163378337+0,102524383042804i
Y31xV1= 0
Y32xV2= -0,889484535296775+10,4099583622151i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,889484535296775-10,4099583622151i
[ + ]= 0,701319371918438-10,3074339791723i
V3= 1,03837763674727-0,0330245199700664i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V3^k-0,212121836864972+0,000839409198431186i
alfa*V3^k+1= 1,24605316409672-0,0396294239640797i
V3acel= 1,03393132723175-0,0387900147656485i
|V3acel|= 1.034658714
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,000104261293770058-0,0118923034484003i
(|V3^(k+1) - V3^k| ) ( -0,0266778570931099-0,0345929687734926i
7. PROCESO ITERATIVO 3
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,00149929128372-23,4113824628275i
Y23*V3= -0,639634831262639+10,2753497891912i
Y21*V1+Y22*V2+Y23*V3= -0,123284054830409+1,7154524748786i
V2*= 1,04989068157448+0,0151511300924737i
V2*{ }= -0,155425823967505+1,79916967530586i
Qiny2= 0,155425823967505-1,79916967530586i
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,194574176032495-1,79916967530586i
Sesp2*= -0,194574176032495+1,79916967530586i
V2*= 1,04989068157448+0,0151511300924737i
Sesp2*/V2*= -0,160564317898009+1,71599047195347i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,00149929128372-23,4113824628275i
Y23*V3= -0,639634831262639+10,2753497891912i
[Y21V1+Y22V2+Y23V3]= -2,12478334611413+25,1268349377061i
{-[ ]}= 2,12478334611413-25,1268349377061i
[ + ]= 1,96421902821612-23,4108444657526i
V2= 1,04970091473136-0,0168064031801335i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209978136314896+0,00303022601849474i
alfa*V2^k+1= 1,25964109767763-0,0201676838161602i
V2acel= 1,04966296136273-0,0171374577976655i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.04980285
|V2esp|/|V2acel|= 1,00018779669018
V2correg= 1,04986008459268-0,017140676155518i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03393132723175+0,0387900147656485i
Sesp3*/V3*= -0,189253067336274+0,111506142301261i
Y31xV1= 0
Y32xV2= -0,869755765383663+10,4116252693884i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,869755765383663-10,4116252693884i
[ + ]= 0,680502698047389-10,3001191270871i
V3= 1,0374373000179-0,0350341882170638i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V3^k-0,20678626544635+0,0077580029531297i
alfa*V3^k+1= 1,24492476002148-0,0420410258604766i
V3acel= 1,03813849457513-0,0342830229073469i
|V3acel|= 1.038704414
7. PROCESO ITERATIVO 4
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 1,95713044921288-23,4151333252358i
Y23*V3= -0,688424025249169+10,3125425432264i
Y21*V1+Y22*V2+Y23*V3= -0,216442090887779+1,7488943665055i
V2*= 1,04986008459268+0,017140676155518i
V2*{ }= -0,257211143815341+1,8323844237768i
Qiny2= 0,257211143815341-1,8323844237768i
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,092788856184659-1,8323844237768i
Sesp2*= -0,092788856184659+1,8323844237768i
V2*= 1,04986008459268+0,017140676155518i
Sesp2*/V2*= -0,0598703024062671+1,74633807699528i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 1,95713044921288-23,4151333252358i
Y23*V3= -0,688424025249169+10,3125425432264i
[Y21V1+Y22V2+Y23V3]= -2,17357254010066+25,1640276917413i
{-[ ]}= 2,17357254010066-25,1640276917413i
[ + ]= 2,11370223769439-23,417689614746i
V2= 1,05067027452481-0,0101900854873814i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209972016918536+0,0034281352311036i
alfa*V2^k+1= 1,26080432942977-0,0122281025848577i
V2acel= 1,05083231251123-0,0087999673537541i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.050869159
|V2esp|/|V2acel|= 0,999172914565602
V2correg= 1,04996318441156-0,00879268902893263i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03813849457513+0,0342830229073469i
Sesp3*/V3*= -0,189012512060244+0,110224681298601i
Y31xV1= 0
Y32xV2= -0,952511182299241+10,4043807258857i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,952511182299241-10,4043807258857i
[ + ]= 0,763498670238997-10,2941560445871i
V3= 1,03767159392011-0,0266699556910215i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V3^k= -0,207627698915026+0,00685660458146938i
alfa*V3^k+1= 1,24520591270413-0,0320039468292258i
V3acel= 1,0375782137891-0,0251473422477564i
|V3acel|= 1.037882912
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( 0,000103099818880015+0,00834798712658537i
(|V3^(k+1) - V3^k| ) ( -0,000560280786029876+0,0091356806595905i
7. PROCESO ITERATIVO 5
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,1432424834081-23,3988320198065i
Y23*V3= -0,778321575555976+10,2979499803354i
Y21*V1+Y22*V2+Y23*V3= -0,120227606999366+1,7506031090438i
V2*= 1,04996318441156+0,00879268902893263i
V2*{ }= -0,141627069850141+1,83701169105137i
Qiny2= 0,141627069850141-1,83701169105137i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209992636882312+0,00175853780578653i
alfa*V2^k+1= 1,25958808678078-0,0139400991913073i
V2acel= 1,04959544989847-0,0121815613855208i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.049666137
|V2esp|/|V2acel|= 1,00031806602238
V2correg= 1,04992929054833-0,0121854359262971i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,0375782137891+0,0251473422477564i
Sesp3*/V3*= -0,190123324194575+0,108646904139154i
Y31xV1= 0
Y32xV2= -0,918886070579563+10,4074042984253i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,918886070579563-10,4074042984253i
[ + ]= 0,728762746384988-10,2987573942861i
V3= 1,03778422592714-0,0301916813667382i
(1-alfa)= -2.00E-01
1-alfa)*V3^k= -0,20751564275782+0,00502946844955128i
alfa*V3^k+1= 1,24534107111257-0,0362300176400858i
V3acel= 1,03782542835475-0,0312005491905345i
|V3acel|= 1.03829432
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,0000338938632300501-0,00339274689736447i
(|V3^(k+1) - V3^k| ) ( 0,0002472145656498-0,0060532069427781i
7. PROCESO ITERATIVO 6
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06762160588125-23,4056354139063i
Y23*V3= -0,718633600444955+10,306390923503i
Y21*V1+Y22*V2+Y23*V3= -0,136160509415195+1,7522406581116i
V2*= 1,04992929054833+0,0121854359262971i
V2*{ }= -0,164310723317866+1,83806961587788i
Qiny2= 0,164310723317866-1,83806961587788i
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,185689276682134-1,83806961587788i
Sesp2*= -0,185689276682134+1,83806961587788i
V2*= 1,04992929054833+0,0121854359262971i
Sesp2*/V2*= -0,156519665303353+1,75247694563289i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06762160588125-23,4056354139063i
Y23*V3= -0,718633600444955+10,306390923503i
[Y21V1+Y22V2+Y23V3]= -2,20378211529644+25,1578760720179i
{-[ ]}= 2,20378211529644-25,1578760720179i
[ + ]= 2,04726244999309-23,405399126385i
V2= 1,0498282143777-0,0130896543840686i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V2^k= -0,209985858109666+0,00243708718525942i
alfa*V2^k+1= 1,25979385725324-0,0157075852608823i
V2acel= 1,04980799914357-0,0132704980756229i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.049891871
|V2esp|/|V2acel|= 1,00010299042765
V2correg= 1,04991611931835-0,0132718648098948i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03782542835475+0,0312005491905345i
Sesp3*/V3*= -0,189412426425283+0,109708562774469i
Y31xV1= 0
Y32xV2= -0,908116308138022+10,4083495623697i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,908116308138022-10,4083495623697i
[ + ]= 0,718703881712739-10,2986409995952i
V3= 1,03767187007868-0,0311970003370854i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
1-alfa)*V3^k= -0,20756508567095+0,0062401098381069i
alfa*V3^k+1= 1,24520624409442-0,0374364004045025i
V3acel= 1,03764115842347-0,0311962905663956i
|V3acel|= 1.038110005
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -0,0000131712299800135-0,0010864288835977i
(|V3^(k+1) - V3^k| ) ( -0,000184269931279957+4,25862413890032E-06i
7. PROCESO ITERATIVO 7
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,04340104789105-23,407762396079i
Y23*V3= -0,718493319563875+10,3045622522783i
Y21*V1+Y22*V2+Y23*V3=-0,160240786524315+1,7482850047142i
V2*= 1,04991611931835+0,0132718648098948i
V2*{ }= -0,191442386975862+1,83342591355621i
Qiny2= 0,191442386975862-1,83342591355621i
V2*= 1,04991611931835+0,0132718648098948i
Sesp2*/V2*= -0,128924456136964+1,74788913870586i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,04340104789105-23,407762396079i
Y23*V3= -0,718493319563875+10,3045622522783i
[Y21V1+Y22V2+Y23V3]= -2,20364183441536+25,1560474007932i
{-[ ]}= 2,20364183441536-25,1560474007932i
[ + ]= 2,0747173782784-23,4081582620873i
V2= 1,0500730351814-0,0118811456287278i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
0 (1-alfa)*V2^k= -0,20998322386367+0,00265437296197896i
alfa*V2^k+1= 1,26008764221768-0,0142573747544734i
99009900990099+9,9009900990099i V2acel= 1,05010441835401-0,0116030017924944i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.05016852
|V2esp|/|V2acel|= 0,999839531106601
V2correg= 1,04993590926004-0,0116011398716367i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
V3*= 1,03764115842347+0,0311962905663956i
Sesp3*/V3*= -0,189445920926263+0,109728275039168i
Y31xV1= 0
Y32xV2= -0,924677733211557+10,4068913192792i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,924677733211557-10,4068913192792i
[ + ]= 0,735231812285294-10,29716304424i
V3= 1,03768946451991-0,0295284276383251i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,207528231684694+0,00623925811327912i
alfa*V3^k+1= 1,24522735742389-0,0354341131659901i
V3acel= 1,0376991257392-0,029194855052711i
|V3acel|= 1.03810973
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( 0,0000197899416900515+0,00
(|V3^(k+1) - V3^k| ) ( 0,0000579673157301563+0,00
7. PROCESO ITERATIVO 8
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,08064671169466-23,4044811413311i
Y23*V3= -0,738366906150584+10,3031545667769i
Y21*V1+Y22*V2+Y23*V3=-0,142868709307414+1,7501585739607i
V2*= 1,04993590926004+0,0116011398716367i
V2*{ }= -0,17030682262555+1,83589689382073i
Qiny2= 0,17030682262555-1,83589689382073i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V2^k= -0,209987181852008+0,00232022797432734i
alfa*V2^k+1= 1,25987010267799-0,0143977990657559i
V2acel= 1,04988292082598-0,0120775710914286i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.04995239
|V2esp|/|V2acel|= 1,00004534768056
V2correg= 1,0499305305813-0,0120781187812644i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
) V3*= 1,0376991257392+0,029194855052711i
) Sesp3*/V3*= -0,189657258906037+0,10936270234422i
) Y31xV1= 0
Y32xV2= -0,919949844325402+10,4073103213805i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,919949844325402-10,4073103213805i
[ + ]= 0,730292585419365-10,2979476190363i
V3= 1,0377185170381-0,0300304990588845i
(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,20753982514784+0,0058389710105422i
alfa*V3^k+1= 1,24526222044572-0,0360365988706614i
V3acel= 1,03772239529788-0,0301976278601192i
|V3acel|= 1.03816168
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -5,37867873995346E-06-0,0004
(|V3^(k+1) - V3^k| ) ( 0,0000232695586799281-0,001
7. PROCESO ITERATIVO 9
Qiny2= -Im{V2*[Y21*VI+Y22*V2+Y23*V3]}
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,07001396997936-23,405423952535i
Y23*V3= -0,728461501679889+10,3043778028108i
Y21*V1+Y22*V2+Y23*V3=-0,143596046552019+1,7504389987907i
V2*= 1,0499305305813+0,0120781187812644i
V2*{ }= -0,17190788349249+1,83610497664374i
Qiny2= 0,17190788349249-1,83610497664374i
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,17809211650751-1,83610497664374i
Sesp2*= -0,17809211650751+1,83610497664374i
V2*= 1,0499305305813+0,0120781187812644i
Sesp2*/V2*= -0,149485402607074+1,7505067483607i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,07001396997936-23,405423952535i
Y23*V3= -0,728461501679889+10,3043778028108i
[Y21V1+Y22V2+Y23V3]= -2,21361001653138+25,1558629513257i
{-[ ]}= 2,21361001653138-25,1558629513257i
[ + ]= 2,06412461392431-23,405356202965i
V2= 1,04990131873476-0,0123396875322646i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V2^k= -0,20998610611626+0,00241562375625288i
alfa*V2^k+1= 1,25988158248171-0,0148076250387175i
V2acel= 1,04989547636545-0,0123920012824646i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.04996861
|V2esp|/|V2acel|= 1,00002990021246
V2correg= 1,04992686846325-0,0123923718059357i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
) V3*= 1,03772239529788+0,0301976278601192i
) Sesp3*/V3*= -0,189542162322847+0,109540170496768i
) Y31xV1= 0
Y32xV2= -0,916834802380092+10,4075852043945i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,916834802380092-10,4075852043945i
[ + ]= 0,727292640057245-10,2980450338977i
V3= 1,03769822882881-0,0303316472153118i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,207544479059576+0,00603952557202384i
alfa*V3^k+1= 1,24523787459457-0,0363979766583742i
V3acel= 1,03769339553499-0,0303584510863504i
|V3acel|= 1.03813738
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( -3,66211805014061E-06-0,0003
(|V3^(k+1) - V3^k| ) ( -0,0000289997628899474-0,00
7. PROCESO ITERATIVO 10
V2=1/Y22[((Sesp2*=-0,35+Qiny2)/V2*)-Y21*V1-Y22*V2-Y23*V3]
1/Y22= 0,004448595488407+0,0444649871706679i
Sesp2= -0,17189845626565-1,83520302669839i
Sesp2*= -0,17189845626565+1,83520302669839i
V2*= 1,04992686846325+0,0123923718059357i
Sesp2*/V2*= -0,143073278579984+1,74962285387627i
Y21*V1= -1,48514851485149+14,8514851485149i
Y22*V2= 2,06300842525917-23,4060424777688i
Y23*V3= -0,726840479872758+10,3042499073626i
[Y21V1+Y22V2+Y23V3]= -2,21198899472425+25,1557350558775i
{-[ ]}= 2,21198899472425-25,1557350558775i
[ + ]= 2,06891571614427-23,4061122020012i
V2= 1,04995624789793-0,0121300143674224i
V2acel=(1-alfa)*V2^(K)+ALFA*V2^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V2^k= -0,20998537369265+0,00247847436118714i
alfa*V2^k+1= 1,25994749747752-0,0145560172409069i
V2acel= 1,04996212378487-0,0120775428797198i
V2correg=(|V2esp|/|V2acel|)xV2acel
|V2esp|= 1.05
|V2acel|= 1.05003158
|V2esp|/|V2acel|= 0,999969920438937
V2correg= 1,04993054138505-0,0120771795925313i
|V2correg=| 1.05000000
Estimar el valor de V3
V3=1/Y33[((Sesp3*)/V3*)-Y31xV1-Y32xV2-Y33xV3]
1/Y33= 0,01001201077256+0,100059434940318i
Sesp3*= -0,2+0,107948564427617i
) V3*= 1,03769339553499+0,0303584510863504i
) Sesp3*/V3*= -0,189529544482794+0,109572238122028i
) Y31xV1= 0
Y32xV2= -0,919959153920532+10,4073094984584i
Y33xV3= 0
[-Y21V1-Y22V2-Y23V3]= 0,919959153920532-10,4073094984584i
[ + ]= 0,730429609437738-10,2977372603364i
V3= 1,03769884055141-0,0300146823994633i
V3acel=(1-alfa)*V3^(K)+ALFA*V3^(k+1)
(1-alfa)= -2.00E-01
(1-alfa)*V3^k= -0,207538679106998+0,00607169021727008i
alfa*V3^k+1= 1,24523860866169-0,036017618879356i
V3acel= 1,03769992955469-0,0299459286620859i
|V3acel|= 1.03813193
(|V1^(k+1) - V1^k| ) ( 0
MAX ( |V2^(k+1) - V2^k| ) < E ( 3,67292180003531E-06+0,0003
(|V3^(k+1) - V3^k| ) ( 6,53401969996636E-06+0,0004
)
)
)
)
)
)
)
)
)
. PROCESO ITERATIVO 7
2850047142i
0474007932i
8913192792i
)
,0000197899416900515+0,0016707249382581i )
,0000579673157301563+0,0020014355136846i )
. PROCESO ITERATIVO 8
1585739607i
6397152918i
3103213805i
)
5,37867873995346E-06-0,000476978909627699i )
,0000232695586799281-0,0010027728074082i )
. PROCESO ITERATIVO 9
4389987907i
8629513257i
5852043945i
)
3,66211805014061E-06-0,0003142530246713i )
0,0000289997628899474-0,000160823226231201i )
. PROCESO ITERATIVO 10
6925781087i
7350558775i
3094984584i
)
,67292180003531E-06+0,000315192213404399i ) < 0,00001
,53401969996636E-06+0,000412522424264501i )