Ejemplo de calculo de contenido armonico en redes de AC

INSTITUTO TECNOLOGICO DE MORELIA

[pic 1]

Nombre: Erik Erudiel Santos González
Materia: Calidad de la energia
Maestro: Dr. Manuel Madrigal Martinez

Tarea 2

                        Lugar y fecha: Morelia, Michoacán 2 de abril del 2019

La tabla muestra los valores discretos de corriente y voltaje en terminales de una carga.

PARA LAS TRES CARGAS OBTENER:        

        1. CONTENIDO ARMÓNICO DE VOLTAJE Y CORRIENTE

        2. THD DE VOLTAJE Y CORRIENTE

        3. VALORES RMS DE VOLTAJE Y CORRIENTE

        4. POTENCIAS APARENTE, ACTIVA, REACTIVA Y FACTOR DE POTENCIA

1.- A partir de la transformada discreta de Fourier obtendremos la magnitud de contenido armónico en cada una de las cargas.

[pic 2]

Aplicando la ecuación anterior para cada una de los valores discretos de cada carga obtendremos la magnitud para cada armónico.  

2.-Distorsion armónica total

[pic 3]

[pic 4]

3.- Valores de RMS de voltaje y corriente

[pic 5]

[pic 6]

4.-Potencia Aparente, activa y reactiva.

[pic 7]

[pic 8]

[pic 9]

Resultados de contenido armónico para voltaje de la carga 1.

                                                      [pic 10]

 [pic 11]

Resultados de contenido armónico para corriente de la carga 1.[pic 12][pic 13]

Resultados de contenido armónico para voltaje de la carga 2.

[pic 14][pic 15]

Resultados de contenido armónico para corriente de la carga 2.[pic 16]

[pic 17]

Resultados de contenido armónico para voltaje de la carga [pic 18][pic 19]

Resultados de contenido armónico para corriente de la carga 3 [pic 20]

[pic 21]

Resultados finales para todas las cargas.

[pic 22]

Rutina en Matlab utilizada:

N=length(V);

x=zeros(N,1);

Fk=zeros(N,1);

for k=0:1:N-1

 for n=0:1:N-1

 x(n+1) = (V(n+1)*exp(-2*n*k*pi*j/N));

 end

 Fk(k+1) = sum(x)/N;

end

[theta, rho] = cart2pol(real(Fk), imag(Fk));

magnitud=zeros(N/2,1);

magnitud(1)=rho(1);

magnitud(2:N/2-1,1)=rho(2:N/2-1)*2

angulo=zeros(N,1);

theta=theta*180/pi;

angulo=theta(1:N)

figure(1)

bar(magnitud, 0.5, 'k')

grid on;

title('Contenido Armónico de voltaje ')

xlim([0 63])

ylabel('Magnitud');

xlabel('Armónicas');

N=length(I);

x=zeros(N,1);

Fk=zeros(N,1);

for k=0:1:N-1

 for n=0:1:N-1

 xI(n+1) = (I(n+1)*exp(-2*n*k*pi*j/N));

...