Matlab - Ayuda con iteraciones y while loop

   
Vista:

Ayuda con iteraciones y while loop

Publicado por Luis enrique luizzrosado@hotmail.com (1 intervención) el 29/05/2016 22:53:04
* Hola buenas tardes , lo siento , pero tengo el siguiente código que he estado haciendo , pero no puedo hacer iteraciones . Todo es correcto , las ecuaciones me dan los resultados que quiero , pero deben cumplir con el mayor número de iteraciones como sea posible hasta que el resultado de " increq = 0,001 " . En cada iteración los valores de v2, s2, s3, q20 ; Deben cambiar de la siguiente manera : v2 = v21 , s2 = s21 , s31 = s3 , q20 q21 = . si me pueden ayudar por favor. Estoy desesperado.*


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
clc
  disp ('===========================================')
  disp ('Modelado de sistemas electricos de potencia')
  disp ('Unidad 2: Flujo de potencias; 3 nodos')
  disp ('===========================================')
  z12=input('Coloque impedancia de linea 1-2 ====> ');
  z13=input('Coloque impedancia de linea 1-3 ====> ');
  z23=input('Coloque impedancia de linea 2-3 ====> ');
  y12=-(z12)^-1;
  y13=-(z13)^-1;
  y23=-(z23)^-1;
  disp ('===========================================')
 
  B11=y12+y13;,B12=-y12;,B13=-y13;
  B21=-y12;,B22=y12+y23;,B23=-y23;
  B31=-y13;,B32=-y23;,B33=y13+y23;
 
  matriz=[B11 B12 B13;B21 B22 B23;B31 B32 B33]
  disp ('¿EN QUE BUS SE ENCUENTRA EL CONDENSADOR? 1, 2 o 3')
  bus=input(' ')
  if bus == 1
      ybus=-[B22 B23;B32 B33]
      elseif bus == 2
          ybus=-[B11 B13;B31 B33]
      else bus == 3
          ybus=-[B11 B12;B21 B22]
      end
  ybusinversa=inv(ybus)
  %ITERACIONES%
  disp('===================================')
  disp ('Iteraciones')
  disp('===================================')
  p20=input ('ingrese la potencia "P" de la carga (recuerde que las cargas se representan con signo negativo)===> ');
  q20=input ('ingrese la potencia "Q" de la carga (recuerde que las cargas de representan con signo negativo)===> ');
  p30=input ('ingrese la potencia "P" del generador====> ');
      v1=1;, s1=0;
  v2=1;, s2=0;
  v3=1;, s3=0;
            %CALCULANDO POTENCIA ACTIVA P%
  p21=v2*((B21*v1*sin(s2-s1))+(B23*v3*sin(s2-s3)));
  p31=v3*((B31*v1*sin(s3-s1))+(B32*v2*sin(s3-s2)));
            %CALCULAR INCREMENTO DE P%
  incp2v2=p20-p21;
  incp3v3=p30-p31;
            %CALCULAR INCREMENTO DE  Q%
        matrizincp=[incp2v2;incp3v3];
  incdeq=(ybusinversa)*(matrizincp);
                  %s21 y s31%
  des2=[incdeq];
  s21=des2(1:2:end);
  s31=des2(2:2:end);
  %CALCULO DE Q%
  singrad1=sin(acos(s21-s1));
  singrad2=sin(acos(s21-s31));
  q21=((-B22*(v2^2))-((v2*(B21*v1*(singrad1)))+(B23*v3*(((singrad2))))));
  %INCREMENTO DE Q%
  increq=(q20-q21)/v2;
  %INCREMENTO DE V2%
  increv2=-(1/B22)*(increq);
  %CALCULO V2%
  v21=v2+increv2;
%FIN%
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder