Saltar al contenido principal

# 5.5: Ejemplo

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$

( \newcommand{\kernel}{\mathrm{null}\,}\) $$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\id}{\mathrm{id}}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\kernel}{\mathrm{null}\,}$$

$$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$

$$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$

$$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\AA}{\unicode[.8,0]{x212B}}$$

$$\newcommand{\vectorA}[1]{\vec{#1}} % arrow$$

$$\newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow$$

$$\newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vectorC}[1]{\textbf{#1}}$$

$$\newcommand{\vectorD}[1]{\overrightarrow{#1}}$$

$$\newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}}$$

$$\newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}}$$

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\avec}{\mathbf a}$$ $$\newcommand{\bvec}{\mathbf b}$$ $$\newcommand{\cvec}{\mathbf c}$$ $$\newcommand{\dvec}{\mathbf d}$$ $$\newcommand{\dtil}{\widetilde{\mathbf d}}$$ $$\newcommand{\evec}{\mathbf e}$$ $$\newcommand{\fvec}{\mathbf f}$$ $$\newcommand{\nvec}{\mathbf n}$$ $$\newcommand{\pvec}{\mathbf p}$$ $$\newcommand{\qvec}{\mathbf q}$$ $$\newcommand{\svec}{\mathbf s}$$ $$\newcommand{\tvec}{\mathbf t}$$ $$\newcommand{\uvec}{\mathbf u}$$ $$\newcommand{\vvec}{\mathbf v}$$ $$\newcommand{\wvec}{\mathbf w}$$ $$\newcommand{\xvec}{\mathbf x}$$ $$\newcommand{\yvec}{\mathbf y}$$ $$\newcommand{\zvec}{\mathbf z}$$ $$\newcommand{\rvec}{\mathbf r}$$ $$\newcommand{\mvec}{\mathbf m}$$ $$\newcommand{\zerovec}{\mathbf 0}$$ $$\newcommand{\onevec}{\mathbf 1}$$ $$\newcommand{\real}{\mathbb R}$$ $$\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}$$ $$\newcommand{\laspan}[1]{\text{Span}\{#1\}}$$ $$\newcommand{\bcal}{\cal B}$$ $$\newcommand{\ccal}{\cal C}$$ $$\newcommand{\scal}{\cal S}$$ $$\newcommand{\wcal}{\cal W}$$ $$\newcommand{\ecal}{\cal E}$$ $$\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}$$ $$\newcommand{\gray}[1]{\color{gray}{#1}}$$ $$\newcommand{\lgray}[1]{\color{lightgray}{#1}}$$ $$\newcommand{\rank}{\operatorname{rank}}$$ $$\newcommand{\row}{\text{Row}}$$ $$\newcommand{\col}{\text{Col}}$$ $$\renewcommand{\row}{\text{Row}}$$ $$\newcommand{\nul}{\text{Nul}}$$ $$\newcommand{\var}{\text{Var}}$$ $$\newcommand{\corr}{\text{corr}}$$ $$\newcommand{\len}[1]{\left|#1\right|}$$ $$\newcommand{\bbar}{\overline{\bvec}}$$ $$\newcommand{\bhat}{\widehat{\bvec}}$$ $$\newcommand{\bperp}{\bvec^\perp}$$ $$\newcommand{\xhat}{\widehat{\xvec}}$$ $$\newcommand{\vhat}{\widehat{\vvec}}$$ $$\newcommand{\uhat}{\widehat{\uvec}}$$ $$\newcommand{\what}{\widehat{\wvec}}$$ $$\newcommand{\Sighat}{\widehat{\Sigma}}$$ $$\newcommand{\lt}{<}$$ $$\newcommand{\gt}{>}$$ $$\newcommand{\amp}{&}$$ $$\definecolor{fillinmathshade}{gray}{0.9}$$

Considera el sistema que se muestra en la Figura 1. Este sencillo sistema cuenta con cinco autobuses (numerados del 1 al 5) y cuatro líneas. Dos de los autobuses están conectados a generadores, dos a cargas y el bus 5 es el “bus oscilante”, representado como un “bus infinito”, o suministro de voltaje.

A los efectos de este ejercicio, supongamos que las impedancias de línea son:

\ [\\ begin {array} {l}
\ mathbf {Z} _ {0} =.05+j .1\
\ mathbf {Z} _ {1} =.05+j .05\
\ mathbf {Z} _ {2} =.15+j .2\
\ mathbf {Z} _ {3} =.04+j .12
\ end {array}\ label {12}\]

También especificamos la potencia real y la magnitud de voltaje para los generadores y la potencia real y reactiva para las cargas:

• Bus 1: La potencia real es 1, el voltaje es de 1.05 por unidad
• Bus 2: La potencia real es 1, el voltaje es 1.00 por unidad
• Bus 3: La potencia real es de -.9 por unidad, la potencia reactiva es 0
• Bus 4: La potencia real es -1, la potencia reactiva es -.2 por unidad.

Tenga en cuenta que la potencia de carga se toma como negativa, ya que este programa simple asume que toda la potencia se mide en la red.

Ejemplo de flujo de carga de% de mente simple

% Primero, impedancias

z1=.05+j*.1;

z2=.05+j*.05;

z3=.15+j*.2;

z4=.04+j*.12;

% Esta es la Matriz de incidencia de nódulos

NI= [1 0 0 0; 0 0 0 1; -1 1 1 0; 0 0 -1 -1; 0 -1 0 0 0];

% Este es el vector de magnitudes de voltaje “conocidas”

VNM = [1.05 1 0 0 1] ';

% Y el vector de ángulos de voltaje conocidos

VNA = [0 0 0 0 0] ';

% y esta es la “clave” a la que realmente se conocen

KNM = [1 1 0 0 1] ';

KNA = [0 0 0 0 1] ';

% y que van a ser manipulados por el sistema

KUM = 1 - KNM;

KUA = 1 - KNA;

% Aquí están las cargas conocidas (positivo es la red INTO

% Usa ceros para incógnitas

P= [1 1 -.9 -1 0] ';

Q= [0 0 0 -.2 0] ';

% y aquí están los vectores correspondientes para indicar

% qué elementos deben verificarse en la comprobación de errores

PC = [1 1 1 1 0] ';

QC = [0 0 1 1 0] ';

Comprobación = KNM + KNA + PC + QC;

% de vectores P y Q desconocidos

PU = 1 - PC;

QU = 1 - QC;

fprintf ('Aquí está la matriz de admisión de líneas:\ n');

Y= [1/Z1 0 0 0; 0 1/Z2 0 0; 0 0 1/Z3 0; 0 0 0 1/Z4]

fprintf ('Y aquí está la matriz de admisión de bus\ n')

YN=NI*Y*NI'

% Ahora: aquí hay algunas magnitudes y ángulos de voltaje de arranque

VM = [1.05 1 .993 .949 1] ';

VA = [.0965 .146 .00713 .0261 0] ';

% Aquí comienza un bucle

Error = 1;

TOL=1e-10;

N = longitud (VNM);

% Construir una tensión candidata a partir de lo que tenemos hasta ahora

VMAG = VNM .* KNM + VM .* KUM;

VANG = VNA .* KNA + VA .* KUA;

V = VMAG .* exp (j .* VANG);

% y calcular la potencia para comenzar

I = (YN*V);

PI = real (V .* conj (I));

QI = imag (V .* conj (I));

%pausa

mientras que (Error>Tol);

para i=1:N,% Corre por todos los autobuses

% ¡Lo que hagamos depende de qué autobús!

si (KUM (i) == 1) & (KUA (i) == 1),% no conoce la magnitud o el ángulo del voltaje

pvc= (P (i) -j*q (i)) /conj (V (i));

para n=1:N,

si n ~=i, pvc = pvc - (YN (i, n) * V (n)); final

final

V (i) = PVC/Yn (i, i);

elseif (KUM (i) == 0) & (KUA (i) == 1),% conoce magnitud pero no ángulo

% primero debe generar una estimación para Q

Qn = imag (V (i) * conj (YN (i,:) *V));

pvc= (P (i) -j*qn) /conj (V (i));

para n=1:N,

si n ~=i, pvc = pvc - (YN (i, n) * V (n)); final

final

PV=PVC/Yn (i, i);

V (i) = VM (i) * exp (j*ángulo (pv));

% final probablemente debería tener más casos

% final a través de la lista de voltaje: error de verificación

% Ahora calcula las corrientes indicadas por esta expresión de voltaje

I = (YN*V);

% Para fines de comprobación de errores, calcule la potencia indicada

PI = real (V .* conj (I));

QI = imag (V .* conj (I));

% Ahora nos enteramos de lo cerca que estamos de las condiciones deseadas

PERR = (P-PI) .* PC;

QERR = (Q-QI) .* QC;

Error = suma (abs (PERR) .^2 + abs (QERR) .^2);

final

fprintf ('Aquí están los voltajes\ n')

V

fprintf ('Poder real\ n')

P

fprintf ('Potencia reactiva\ n')

Q

This page titled 5.5: Ejemplo is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by James Kirtley (MIT OpenCourseWare) via source content that was edited to the style and standards of the LibreTexts platform.