> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 1

Taller 1 de Electrónica de Potencia

Fonseca Díaz Cristian David, Rincón Miranda Adriana

Universidad Distrital Francisco José de Caldas

3) Empiezan a conducir los diodos D3 y D4:

I. EJERCICIO 1.

Debido a la periodicidad de 𝐼𝑥 de periodo π, los diodos

empiezan a conducir nuevamente en un tiempo π + α.

𝑑𝑉𝑐 (𝜃)
𝑉𝑐 = −𝑉𝑝 sin 𝜃; 𝐼𝑐 (𝜃) = 𝐶𝜔 ; 𝐼𝑐 (𝜃) = −𝐶𝜔 cos 𝜃 ;
𝑑𝜃

𝐼𝑥 (𝜃) = −𝐶𝜔 cos 𝜃 + 𝐼𝑜 (3)

4) Utilizando las conclusiones del análisis anterior, se

obtienen los valores de las variables solicitadas:

- Teniendo en cuenta que en β 𝐼𝑥 = 0, entonces se

despeja la ecuación (1)

𝐼𝑥 (𝛽) = 𝐶𝜔 cos 𝛽 + 𝐼𝑜 = 0
Fig. 1. Circuito del primer ejercicio.
−𝐼0
Para el ejercicio, se asignan los siguientes valores: 𝑉𝑝 = 120 ∗ 𝛽 = 𝐶𝑜𝑠 −1 ( )
𝜔𝐶𝑉𝑝
√2 ∗ 𝑠𝑒𝑛(𝜃), 𝜃 = 𝜔𝑡 y 𝜔 = 2 ∗ 𝜋 ∗ 60. 𝐶 = 100𝜇𝐹 𝑒 𝐼𝑜 =
3𝐴. −3
𝛽 = 𝐶𝑜𝑠 −1 ( )
1) Conducen D1 y D2: 2𝜋60 ∗ 100 ∗ 10−6 ∗ 120 ∗ √2

𝛽 = 2.05886 𝑟𝑎𝑑.
θ = α 𝑉𝑐 = 𝑉𝑔 = 𝑉𝑝 sin 𝜃; 𝐼𝑥 (𝜃) = 𝐼𝑐 (𝜃) + 𝐼𝑥 (𝜃);
Se utilizan métodos numéricos para encontrar el valor de α por
𝑑𝑉𝑐 (𝜃)
𝐼𝑐 (𝜃) = 𝐶𝜔 ; medio de la ecuación (2). Esto teniendo en cuenta que
𝑑𝜃
𝑉𝑐 (𝜋 + 𝛼) = −𝑉𝑝 𝑆𝑒𝑛(𝜋 + 𝛼) = 𝑉𝑝 𝑆𝑒𝑛(𝛼):
𝐼𝑐 (𝜃) = 𝐶𝜔 cos 𝜃 ;
−𝐼𝑜
𝑉𝑐 (𝜃) = (𝜃 − 𝛽) + 𝑉𝑐 (𝛽) (4)
𝜔𝐶
𝐼𝑥 (𝜃) = 𝐶𝜔 cos 𝜃 + 𝐼𝑜 (1);
−𝐼𝑜
𝐼𝑥 es cero en β >
𝜋 𝑉𝑝 𝑆𝑒𝑛(𝛼) = (𝜋 + 𝛼 − 𝛽) + 𝑉𝑐 (𝛽)
2 𝜔𝐶

2) No hay conducción de ningún diodo (𝜃 = 𝛽):

120 ∗ √2 ∗ 𝑆𝑒𝑛(𝛼)
−3
𝑑𝑉𝑐 (𝜃) 𝐼𝑜 𝑑𝑉𝑐 (𝜃) = (𝜋 + 𝛼
𝑉𝑐 ≠ 𝑉𝑝 ; 𝐼𝑐 = −𝐼𝑜 ; −𝐼𝑜 = 𝐶𝜔 ;− = ; 2 ∗ 60 ∗ 𝜋 ∗ 100 ∗ 10−6
𝑑𝜃 𝐶𝜔 𝑑𝜃
− 2.0589) + 120 ∗ √2 ∗ 𝑆𝑒𝑛(2.0589)
𝜃 𝐼
𝑜 𝜃
∫𝛽 − 𝐶𝜔 𝑑𝜃 = ∫𝛽 𝑑𝑉𝑐 (𝜃) 𝑑𝜃;
𝛼 = 0.257586 rad.

Haciendo uso de estos datos, se puede encontrar 𝐼𝑔𝑅𝑀𝑆 de la

Integro y despejo de la ecuación anterior para encontrar 𝑉𝑐 (𝜃) siguiente forma:
−𝐼𝑜
𝑉𝑐 (𝜃) = (𝜃 − 𝛽) + 𝑉𝑐 (𝛽) (2)
𝜔𝐶
> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 2

2𝜋+∝ 1) Graficas de corriente de generador ig y corriente ix

1 2
𝐼𝑔𝑅𝑀𝑆 = √ ∗∫ (𝐶𝜔𝑉𝑝 cos 𝜃 + 𝐼𝑜 ) 𝑑𝜃
2𝜋 ∝

𝛽
1 2
𝐼𝑔𝑅𝑀𝑆 = √ ∗ ∫ (𝐶𝜔𝑉𝑝 cos 𝜃 + 𝐼𝑜 ) 𝑑𝜃
𝜋 ∝

𝐼𝑔𝑅𝑀𝑆
2.089
1 2
= √ ∗∫ (100 ∗ 10−6 ∗ 2 ∗ 𝜋 ∗ 60 ∗ 120 ∗ √2 cos 𝜃 + 3) 𝑑𝜃
𝜋 0.25759

𝐼𝑔𝑅𝑀𝑆 = 4.5023 𝐴

- Para encontrar P, se necesita de la armónica

fundamental de 𝐼𝑔 , por lo que es necesario recurrir al
uso de la serie de Fourier sobre la función:

𝐼𝑥(𝜃) , ∝< 𝜃 < 𝛽 Fig. 2. Gráfica Ig(roja) e Ix(azul).

𝐼𝑔 (𝜃) = {
−𝐼𝑥 (𝜃 − 𝜋), 𝜋+∝< 𝜃 < 𝜋 + 𝛽
2) Código Matlab
𝐶𝜔 cos( 𝜃) + 𝐼𝑜 , ∝< 𝜃 < 𝛽
𝐼𝑔 (𝜃) = { clc; clear;
−𝐶𝜔 cos(𝜃 − 𝜋) − 𝐼𝑜 , 𝜋+∝< 𝜃 < 𝜋 + 𝛽
Vp=120*sqrt(2);
w=120*pi; Io=3;
Al realizar la serie de Fourier para obtener el término c=100*10^(-6);
fundamental, se obtiene lo siguiente: b=acos(-Io/(w*c*Vp));
syms a;
𝐼𝑔1 = 3.5234 Cos(𝜃) + 4.989𝑆𝑒𝑛(𝜃) sola=vpasolve(Vp*sin(a)+(Io/(w*c))*(pi+a-
b)-Vp*sin(b)==0,a);
theta0=[0 sola];
𝐼𝑔1 = 5.4814𝑆𝑒𝑛(𝜃 + 40°) Ip0=[0 0];
theta1=sola:pi/200:b;
Igp1=(w*c*Vp*cos(theta1)+Io);
Para obtener 𝐼𝑔1𝑅𝑀𝑆 se divide sobre raíz de dos la magnitud de theta2=[b (pi+sola)];
Igp2=[0 0];
𝐼𝑔1 , al igual que para obtener 𝑉𝑔𝑅𝑀𝑆 con 𝑉𝑔. Con esto,
theta3=(sola+pi):pi/200:(b+pi);
encuentro el valor de las potencias activas y aparentes.
Igp3=(w*c*Vp*cos(theta3)-Io);
Ig=[Ip0,Igp1,Igp2,Igp3];
𝑃 = 𝑉𝑔𝑅𝑀𝑆 ∗ 𝐼𝑔1𝑅𝑀𝑆 ∗ cos(𝜃 − 𝛼)
x=[theta0,theta1,theta2,theta3];
5.4814
𝑃 = 120 ∗ ∗ cos(−40°) = 356.121 𝑊 theta0x=[0 sola];
√2 Ix0=[0 0];
theta1x=sola:pi/200:b;
𝑆 = 𝑉𝑔𝑅𝑀𝑆 ∗ 𝐼𝑔𝑅𝑀𝑆 = 120 ∗ 4.5023 = 540.276 𝑉𝐴 Ixp1=(w*c*Vp*cos(theta1)+Io);
theta2x=[b (pi+b)];
𝑃 356.121 Ixp2=[0 0];
𝐹𝑃 = =
𝑆 540.276 Ix=[Ix0,Ixp1,Ixp2];
x0=[theta0x,theta1x,theta2x];
𝐹𝑃 = 0.6591
subplot(2,1,1);
plot(x,Ig,'r','LineWidth',2);
title('Ig')

subplot(2,1,2);
plot(x0,Ix,'b','LineWidth',2);
title('Ix')
> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 3

II. EJERCICIO 2
Fig. 3. Circuito del segundo ejercicio
𝑣𝑐2 (𝜃) 𝜙 𝑃0
{ =− (𝜑 − 𝛽)
2 𝛽 𝑤∗𝑐

2 ∗ 𝑃0
𝑣𝑐2 (𝜑) = (𝜑 − 𝛽) + 𝑣𝑐2 (𝛽)
𝑤∗𝑐
2 ∗ 𝑃0
𝑣𝑐2 (𝜑) = − (𝜑 − 𝛽) + 𝑣𝑐2 (𝛽)
𝑤∗𝑐
𝛽 <𝜙 <𝜋+𝛼

𝑃𝑎𝑟𝑎 𝑒𝑛𝑐𝑜𝑛𝑡𝑟𝑎𝑟 𝑎 𝑎𝑙𝑝ℎ𝑎 𝑠𝑒 𝑑𝑒𝑏𝑒 𝑒𝑣𝑎𝑙𝑢𝑎𝑟 𝑒𝑛 𝑒𝑠𝑒 𝑝𝑢𝑛𝑡𝑜

𝑠𝑢𝑝𝑜𝑛𝑖𝑒𝑛𝑑𝑜 𝑞𝑢𝑒 𝑝𝑎𝑟𝑎 𝑒𝑠𝑒 𝑣𝑎𝑙𝑜𝑟 𝑉𝑐 = 𝑉𝑝, 𝑞𝑢𝑒 𝑒𝑠 𝑑𝑜𝑛𝑑𝑒
𝑙𝑜𝑠 𝑑𝑖𝑜𝑑𝑜𝑠 𝑣𝑢𝑙𝑣𝑒𝑛 𝑎 𝑐𝑜𝑛𝑑𝑢𝑐𝑖𝑟, 𝑝𝑜𝑟 𝑞𝑢𝑒 𝑒𝑠𝑡𝑎𝑛 𝑝𝑜𝑙𝑎𝑟𝑖𝑧𝑎𝑑𝑜𝑠
𝑑𝑖𝑟𝑒𝑐𝑡𝑎𝑚𝑒𝑛𝑡𝑒.
2 ∗ 𝑃0
El análisis se empieza con la siguiente configuración para los 𝑣𝑐2 (𝜋 + 𝛼) = − (𝜋 + 𝛼 − 𝛽) + 𝑣𝑐2 (𝛽)
diodos: 𝑤∗𝑐
𝐷1 𝑦 𝐷2 𝑂𝑁
{ 2𝑃0
𝐷3 𝑦 𝐷4 𝑂𝐹𝐹 𝑉𝑝2 ∗ 𝑠𝑒𝑛2 (𝜋 + 𝛼) + (𝜋 + 𝛼 − 𝛽) − 𝑉𝑝2 𝑠𝑒𝑛(𝛽)2 = 0
Se considera que se ha conectado la fuente justo en cero, con 𝑤∗𝑐
vc=0, entonces: 𝑆𝑒 𝑑𝑒𝑏𝑒 𝑟𝑒𝑠𝑜𝑙𝑣𝑒𝑟 𝑚𝑎𝑡𝑒𝑚á𝑡𝑖𝑐𝑎𝑚𝑒𝑛𝑡𝑒 𝑝𝑎𝑟𝑎 𝑒𝑛𝑐𝑜𝑛𝑡𝑟𝑎𝑟
−𝑣𝑔 + 𝑣𝑐 = 0 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑎𝑙𝑝ℎ𝑎
𝑉𝐶 = 𝑉𝑔
𝐼𝑔 𝑠𝑒 𝑒𝑣𝑎𝑙𝑢𝑎𝑟𝑎 𝑐𝑜𝑛 𝑟𝑒𝑠𝑝𝑒𝑐𝑡𝑜 𝑎 𝐼𝑋 , 𝑎𝑠í, 𝑡𝑒𝑛𝑑𝑟𝑎
𝑣𝑐 = 𝑉𝑝 ∗ 𝑆𝑒𝑛(𝜃)
𝑑𝑣𝐶 𝑙𝑎 𝑠𝑖𝑔𝑢𝑖𝑒𝑛𝑡𝑒 𝑓𝑜𝑟𝑚𝑎
𝑖𝑐 = 𝐶 ∗ 𝐼𝑥 ∝< 𝜃 < 𝛽
𝑑𝑡 {
𝑑𝑣𝑐 −𝐼𝑋 (𝜃 − 𝜋) 𝜋+∝< 𝜃 < 𝜋 + 𝛽
𝑖𝑐 = 𝑤 ∗ 𝐶 ∗
𝑑𝜃 𝑃𝑎𝑟𝑎 𝑒𝑛𝑐𝑜𝑛𝑡𝑟𝑎𝑟 𝑎 𝐵𝑒𝑡𝑎, 𝑒𝑛 𝑢𝑛 𝑝𝑢𝑛𝑡𝑜 𝑙𝑎 𝑐𝑜𝑟𝑟𝑖𝑒𝑛𝑡𝑒 𝑖 = 0
𝑑𝑉𝑝 ∗ 𝑆𝑒𝑛(𝜃)
𝑖𝑐 = 𝑤 ∗ 𝐶 ∗ 𝑝𝑎𝑟𝑎 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝐵𝑒𝑡𝑎, 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠,
𝑑𝜃
𝑖𝑐 = 𝑤 ∗ 𝐶 ∗ 𝑉𝑝 ∗ cos(𝜃)
𝑠𝑖 𝜃 = 𝛽 𝐼𝑥 = 0
𝑖𝑥 = 𝑖𝑐 + 𝑖𝑜
𝑃𝑜
𝑖𝑜 = 𝑐𝑜𝑠𝑒𝑐𝑎𝑛𝑡𝑒(𝜃) 𝑃𝑜
𝑉𝑝 0 = 𝑤 ∗ 𝑐 ∗ 𝑉𝑝 ∗ 𝑐𝑜𝑠(𝐵) +
𝑐𝑜𝑠𝑒𝑐(𝛽)
𝑃𝑜 𝑉𝑝
𝑖𝑥 = 𝑖𝑐 + 𝑐𝑜𝑠𝑒𝑐𝑎𝑛𝑡𝑒(𝜃) 𝑆𝑒 𝑟𝑒𝑠𝑢𝑒𝑙𝑣𝑒 𝑛𝑢𝑚𝑒𝑟𝑖𝑐𝑎𝑚𝑒𝑛𝑡𝑒 𝑝𝑎𝑟𝑎 𝑒𝑛𝑐𝑜𝑛𝑡𝑟𝑎𝑟 𝑎 𝐵
𝑉𝑝

𝑃𝑜 III. EJERCICIO 3
𝑖𝑥 = 𝑤 ∗ 𝐶 ∗ 𝑉𝑝 ∗ cos(𝜃) + 𝑐𝑜𝑠𝑒𝑐𝑎𝑛𝑡𝑒(𝜃)
𝑉𝑝

Cuando 𝑖𝑥 = 0, se abrirán los diodos y quedara solo la fuente

de potencia con el condensador.

Periodo en el cual Vc se descarga

𝐶𝑑𝑣𝑐
𝐼𝑐 =
𝑑𝑡
𝑑𝑣𝑐
𝐼𝑐 = 𝑤 ∗ 𝐶 ∗
𝑑𝜃
−𝐼0 = 𝐼𝑐
𝑃𝑜 𝑑𝑣𝑐
− =𝑤∗𝑐∗ Fig. 4. Circuito del tercer ejercicio.
𝑣𝑐 𝑑𝜃
𝑃𝑜
𝑣𝑐 𝑑𝑣𝐶 = − 𝑑𝜃 Para el ejercicio, se asignan los siguientes valores: 𝑉𝑝 = 120 ∗
𝑤∗𝑐
√2 ∗ 𝑠𝑒𝑛(𝜃), 𝜃 = 𝜔𝑡 y 𝜔 = 2 ∗ 𝜋 ∗ 60. 𝐿 = 500𝑚𝐻, 𝐶 =
∅ 𝜙
𝑃𝑜 100𝜇𝐹 𝑒 𝐼𝑜 = 1 𝐴. Para este caso particular, los diodos se
∫ 𝑣𝑐 𝑑𝑣𝑐 = − ∫ 𝑑𝜃 encuentran en modo de conducción continuo. Al llevar a cabo
𝛽 𝑤∗𝑐 𝛽
el desarrollo del ejercicio, se obtiene que el comportamiento del
> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 4

circuito se encuentra descrito por una ecuación de segundo Para los valores de rizado se debe derivar la función a la cual se
orden. le va a encontrar el rizado, igualar la derivada a cero y encontrar
las raíces. Por último estos valores se evalúan y se clasifican
𝐼𝑜 𝑑𝐼𝑙2 (𝜃) 𝐼𝐿 (𝜃) según su valor.
𝑉𝑝 cos(𝜃) + = +
𝜔𝐶 𝑑𝜃 𝜔𝐶
Para Il, la derivada es:
La respuesta de esta ecuación diferencial, está compuesta de 𝐼𝑙 ′ = −0.393015𝑠𝑒𝑛(0.3751𝜃) − 0.5879254 cos(0.3751𝜃)
una parte natural y una forzada. + 1.0478𝑠𝑒𝑛(𝜃) = 0 0<𝜃<𝜋

𝐼𝐿 (𝜃) = 𝐴1 cos(𝑘𝜃) + 𝐴2 𝑠𝑒𝑛(𝑘𝜃) + 𝐴𝑘 cos(𝑘𝜃) + 𝐼𝑜 Usando Texas Nspire los puntos son los siguientes:

0 < 𝑉𝑐 < 𝜋 𝜃 = 0.692727 𝐼𝑙 = 0.804 𝑀𝑖𝑛𝑖𝑚𝑜

𝜃 = 2.4488332 𝐼𝑙 = 1.2 𝑀á𝑥𝑖𝑚𝑜
𝜔𝑜 1 𝜔𝐶𝑉𝑝 𝑘 2 ∆𝐼𝑙 = 1.2 − 0.804 = 0.396
𝑘= ; 𝜔𝑜 = ; 𝐴𝑘 =
𝜔 √𝐿𝐶 𝑘 2 −1 ∆𝐼𝑙 = 0.396

Para obtener los valores de los coeficientes de la ecuación, se Para vc, la derivada es:
resuelve el siguiente sistema de ecuaciones: 𝑣𝑐 ′ = 28.0895 cos(0.37513𝜃) − 41.57615𝑆𝑒𝑛(0.3751𝜃)
− 27.792733 cos(𝜃) = 0
2 ∗ 𝐴𝐾 = 𝐴1 [cos(𝑘𝜋) − 1] + 𝐴2 𝑠𝑒𝑛(𝑘𝜋)
0 = −𝐴1 𝑠𝑒𝑛(𝑘𝜋) + 𝐴2 [cos(𝑘𝜋) − 1] 𝜃 = 0.01931467 𝑣𝑐 = 110.82784 𝑀á𝑥𝑖𝑚𝑜
𝜃 = 1.5430 𝑣𝑐 = 105.52417 𝑀𝑖𝑛𝑖𝑚𝑜
El valor de la tensión del condensador se obtiene del análisis ∆𝑣𝑐 = 110.82784 − 105.52417 = 5.3036648
circuital. ∆𝑣𝑐 = 5.3036648
1
𝑉𝑐 (𝜃) = [𝐴1 𝑠𝑒𝑛(𝑘𝜃) − 𝐴2 𝑐𝑜𝑠(𝑘𝜃) + 𝐾𝐴𝑘 𝑠𝑒𝑛(𝜃)]
𝑘𝜔𝐶
Para encontrar Vc e Il, se debe usar el teorema de valor medio:
0 < 𝑉𝑐 < 𝜋
1 𝜋
𝑉𝑐 = ∫ (74.087915𝑠𝑒𝑛(𝑘𝜃) + 110.83072𝑐𝑜𝑠(𝑘𝜃)
Al resolver el sistema, se obtiene lo siguiente: 𝜋 0
− 27.792733𝑠𝑒𝑛(𝜃)) 𝑑𝜃
𝐴1 = 1.04776
𝑉𝑐 = 108.03795 𝑣
𝐴2 = −1.56738

𝐼𝑙 0 < 𝜃 < 𝜋
𝑖𝐷1 = {
Al reemplazar, se obtienen las expresiones de 𝑣𝑐 e 𝑖𝐿 , para 0 < 0 𝜋 < 𝜃 < 2𝜋
𝜃<𝜋
1 𝜋
𝐼𝐷1 = ∫ 𝐼𝑙 𝑑𝜃
𝑣𝑐 (𝜃) = 74.087915𝑠𝑒𝑛(𝑘𝜃) + 110.83072𝑐𝑜𝑠(𝑘𝜃) 𝜋 0
− 27.792733𝑠𝑒𝑛(𝜃) 1 𝜋
𝐼𝐷1 = ∫ (1.0477613 cos(𝑘Ɵ) − 1.56738𝑠𝑒𝑛(𝑘𝜃)
2𝜋 0
𝐼𝐿 (𝜃) = 1.0477613 cos(𝑘Ɵ) − 1.56738𝑠𝑒𝑛(𝑘𝜃) − 1.0478 cos(𝜃) + 1) 𝑑𝜃
− 1.0478 cos(𝜃) + 1
Donde; 𝐼𝐷1 = 0.499999 𝐴
𝑘 = 0.375132
𝐴𝑘 = −1.04776 Para el factor de potencia se necesita 𝐼𝑔𝑟𝑚𝑠 e 𝐼𝑔1 𝑟𝑚𝑠

Como 𝐼𝐿 (𝜃) tiene un período de π, la corriente del diodo D1,

será la siguiente: 𝐼𝑙(𝜃), 0<𝜃<𝜋
𝐼𝑔 (𝜃) = {
−𝐼𝑙(𝜃 − 𝜋), 𝜋 < 𝜃 < 2𝜋
𝐼𝐷1 (𝜃) = 𝑆𝑒 𝑡𝑜𝑚𝑎 𝑒𝑙 𝑎𝑟𝑒𝑎 𝑑𝑒 𝑢𝑛𝑎 𝑦 𝑠𝑒 𝑚𝑢𝑙𝑡𝑖𝑝𝑙𝑖𝑐𝑎 𝑝𝑜𝑟 2
1.0477613 cos(𝑘Ɵ) − 1.56738𝑠𝑒𝑛(𝑘𝜃) − 1.0478 cos(𝜃) + 1;
{ 0<𝜃<𝜋 𝜋
0 ; 𝜋 < 𝜃 < 2𝜋 1
𝐼𝑔𝑟𝑚𝑠 = √ (2 ∗ ∫ (𝐼𝑙)2 𝑑𝜃
2∗𝜋 0
𝐼𝑔𝑟𝑚𝑠 = 1.0098566
> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 5

𝑃𝑎𝑟𝑎 𝐼𝑔1 𝑟𝑚𝑠 𝑠𝑒 ℎ𝑎𝑐𝑒 𝑑𝑒 𝑢𝑠𝑜 𝑑𝑒 𝑙𝑜𝑠 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒𝑠 𝑑𝑒

𝑓𝑜𝑢𝑟𝑖𝑒𝑟 𝑝𝑎𝑟𝑎 𝑙𝑎 𝑓𝑢𝑛𝑑𝑎𝑚𝑒𝑛𝑡𝑎𝑙
𝜋 2𝜋
1
𝑎1 = (∫ 𝐼𝑙(𝜃) cos(Ɵ) 𝑑𝜃 + ∫ −𝐼𝑙(𝜃 − 𝜋) cos(Ɵ) 𝑑𝜃 )
𝜋 0 𝜋 2) Código Matlab
𝑎1 = −0.17652
𝜋 2𝜋
clear; clc
1 theta=0:pi/100:pi;
𝑏1 = (∫ 𝐼𝑙(𝜃)𝑆𝑒𝑛(𝜃) 𝑑𝜃 + ∫ −𝐼𝑙(𝜃 − 𝜋)𝑆𝑒𝑛(𝜃) 𝑑𝜃)
𝜋 0 𝜋 pha=0;
𝑏1 = 1.2732395 theta1=theta;
theta2=theta+pi;
𝐼𝑔1 = −0.17652 cos(𝜃) + 1.2732385𝑠𝑒𝑛(𝜃) Vp=120*sqrt(2);
w=120*pi;
|𝐼𝑔1 | = 1.2854181 l=500*10^(-3);
c=100*10^(-6);
𝐷𝑒𝑠𝑝𝑙𝑎𝑧𝑎𝑚𝑖𝑒𝑛𝑡𝑜 Io=1;
0.17652 k=(1/sqrt(l*c))/(w);
arctan (− ) = −0.13776 𝑟𝑎𝑑 Ak=(w*c*Vp*(k^2))/(k^2-1);
1.2732395
syms x y;
𝐼𝑔1 = 1.2854181 𝑆𝑒𝑛(𝜃 − 0.13776) [A1,A2]=vpasolve([(cos(k*pi)-
1)*x+y*sin(k*pi)==2*Ak,(cos(k*pi)-1)*y-
𝑃 = 𝑉𝑔𝑟𝑚𝑠 𝐼𝑔1 𝑟𝑚𝑠 cos(𝜃1 − 𝜙1 ) x*sin(k*pi)==0],[x,y]);
Vc=(A1/(w*c*k))*sin(k*theta)-
1.2854181 ((A2)/(w*c*k))*cos(k*theta)+(Ak/(w*c))*si
𝑃 = 120 ∗ cos(0.13776) n(theta);
√2 Il=A1*cos(k*theta)+A2*sin(k*theta)+Ak*cos
(theta)+Io;
𝑃 = 108.03796 𝑤 Iln=-
1*(A1*cos(k*theta)+A2*sin(k*theta)+Ak*cos
𝑆 = 𝑉𝑔𝑟𝑚𝑠 ∗ 𝐼𝑔𝑟𝑚𝑠 = 120 ∗ 1.0098566 = 121.18279 (theta)+Io);
fVg=120*sqrt(2)*sin(theta);
𝑃 108.03796
𝑃𝐹 = =
𝑆 121.18279 subplot(2,1,1);
plot(theta1,Il,'r','LineWidth',2);
𝑃𝐹 = 0.8915289 hold on
plot(theta2,Il,'r','LineWidth',2);
1) Graficas Corriente en el inductor y tensión en el title('Il')
condensador.
subplot(2,1,2);
plot(theta1,Vc,'b','LineWidth',2)
hold on
plot(theta2,Vc,'b','LineWidth',2)
title('Vc')

IV. EJERCICIO 4

Fig. 5. Grafica Il(rojo) y grafica Vc(azul)

> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 6

Fig. 6. Circuito del tercer ejercicio. Primero se deben encontrar los valores de k y 𝐴𝑘 para resolver
el sistema de ecuaciones. Esto se hace de la siguiente manera:
Para el ejercicio, se asignan los siguientes valores: 𝑉𝑝 = 120 ∗
√2 ∗ 𝑠𝑒𝑛(𝜃), 𝜃 = 𝜔𝑡 y 𝜔 = 2 ∗ 𝜋 ∗ 60. 𝐿 = 40𝑚𝐻, 𝐶 = 𝜔𝑜 1
𝑘= = = 1.3263
100𝜇𝐹 𝑒 𝐼𝑜 = 1 𝐴. 𝜔 √𝐿 ∗ 𝑐 ∗ 120 ∗ 𝜋
Para el desarrollo del ejercicio, se hace uso del circuito de la
figura 3, pero se tiene en cuenta el hecho de que el modo de ∗ 120 ∗ √2 ∗ (1.3263)2
funcionamiento es discontinuo, por lo que se añaden un par de 𝐴𝑘 = = 393.28247
(1.3263)2 − 1
ecuaciones con el fin de encontrar los puntos de α y β de
conducción y no conducción del puente de diodos. El resto de Para resolver el sistema de ecuaciones se necesario usar un
las ecuaciones, se encuentran desarrolladas en el punto 3. programa, en este caso se usó un código en Matlab, con
condiciones para los ángulos 𝛼 𝑦 𝐵, se muestran los valores
𝑑𝐼𝑙 (𝜃) encontrados.
𝑉𝑝 𝑠𝑒𝑛(𝜃) = 𝜔𝐿 ∗ + 𝑉𝑐 (1)
𝑑𝜃

𝑑𝑉𝑐 (𝜃) 𝐴1 = 123.298496

𝐼𝑙 = 𝜔𝐶 ∗ + 𝐼𝑜 (2)
𝑑𝜃
𝐴2 = −243.1454727
Al derivar la ecuación (2) y hacer el reemplazo en (1), se
obtiene: 𝛽 = 2.8817061 𝑟𝑎𝑑.

𝑉𝑐 = 𝐴1 cos(𝑘𝜃) + 𝐴2 𝑠𝑒𝑛(𝑘𝜃) + 𝐴𝑘 𝑠𝑒𝑛(𝜃) (3) ∝ = 0.85901219 𝑟𝑎𝑑.

Lo anterior, aplica para α <θ< β. Al utilizar la ecuación (4) del Haciendo uso de estos valores, se procede a encontrar 𝐼𝐿 𝑦 𝑣𝑐
punto 1, se deduce que esta es aplicable durante β <θ< α + π. por medio de las ecuaciones (2) y (4).
𝐼𝑜
𝑉𝑐 (𝜃) = (𝛼 + 𝜋 − 𝜃) + 𝑉𝑐 (𝛼) 𝑉𝑐
𝜔𝐶
𝐴1 ∗ cos(𝑘𝜃) + 𝐴1 ∗ 𝑠𝑒𝑛(𝑘𝜃) + 𝐴𝑘 ∗ 𝑠𝑒𝑛(𝜃); 𝛼 < 𝜃 < 𝛽
𝐼𝑜 𝐼𝑜
𝑉𝑐 (𝜃) = (𝛼 + 𝜋 − 𝜃) + 𝑉𝑝 𝑆𝑒𝑛(𝛼) (4) = {
𝜔𝐶 ∗ (𝜋 + 𝛼 − 𝜃) + 𝑉𝑝 ∗ 𝑠𝑒𝑛(𝛼); 𝛽 < 𝜃 < 𝜋 + 𝛼
𝜔𝐶
Al igualar las ecuaciones (3) y (4) para los puntos que cumplen
las condiciones y derivarla, se obtienen las dos ecuaciones 𝐼𝐿
faltantes para conocer los valores de 𝐴1 , 𝐴2 , ∝ 𝑦 𝛽. Las 𝜔𝐶[−𝐴1 ∗ 𝑘 ∗ 𝑠𝑒𝑛(𝑘𝜃) + 𝐴2 ∗ 𝑘 ∗ 𝑐𝑜𝑠(𝑘𝜃) + 𝐴𝑘 cos(𝜃)] + 𝐼𝑜 ; 𝛼 < 𝜃
= {
condiciones son las siguientes: 0; 𝛽 <𝜃 <𝜋+𝛼

𝑉𝑐 (∝) = 𝑉𝑐 (∝ +𝜋) Para encontrar Vc se debe usar el teorema de valor medio:

𝛽
𝑉𝑐 (𝛽) = 𝑉𝑐 (𝛽) 1
𝑉𝐶 = (∫ 𝐴1 𝑐𝑜𝑠(𝑘𝜃) + 𝐴2 𝑠𝑒𝑛(𝑘𝜃) + 𝐴𝑘 𝑠𝑒𝑛(𝜃) 𝑑𝜃
𝜋 𝛼
𝑑𝑉𝑐 (𝛼) 𝑑𝑉𝑐 (𝛼 + 𝜋) 𝜋+𝛼
𝐼𝑜
= +∫ (𝜋 + 𝛼 − 𝜃) + 𝑣𝑝
𝑑 𝑑 𝛽 𝑤 ∗𝑐
𝑑𝑉𝑐 (𝛽) 𝑑𝑉𝑐 (𝛽) ∗ 𝑠𝑒𝑛(𝛼)𝑑𝜃)
=
𝑑𝜃 𝑑𝜃
𝑉𝑐 = 138.54108 𝑣
De esta forma, se obtiene el siguiente sistema de ecuaciones:
Para encontrar 𝐼𝐿𝑅𝑀𝑆 , se hace lo siguiente:
𝐴1 cos(𝑘𝛼) + 𝐴2 𝑠𝑒𝑛(𝑘𝛼) + 𝐴𝑘 𝑠𝑒𝑛(𝛼) = 𝑉𝑝 𝑠𝑒𝑛(𝛼) (5)
𝐼𝐿𝑅𝑀𝑆
𝛽
𝐴1 cos(𝑘𝛽) + 𝐴2 𝑠𝑒𝑛(𝑘𝛽) + 𝐴𝑘 𝑠𝑒𝑛(𝛽) = 𝑉𝑝 𝑆𝑒𝑛(𝛼) + 1
= √ ∗ ∫ (𝜔𝐶[−𝐴1 ∗ 𝑘 ∗ 𝑠𝑒𝑛(𝑘𝜃) + 𝐴2 ∗ 𝑘 ∗ 𝑐𝑜𝑠(𝑘𝜃) + 𝐴𝑘 cos(𝜃)] + 𝐼𝑜 )2 𝑑𝜃
𝐼𝑜 2𝜋 ∝
(𝜋 + 𝛼 − 𝛽) (6)
𝜔𝐶

−𝐼𝑜 De lo anterior, se obtiene que

−𝐴1 ∗ 𝑘 sen(𝑘𝛼) + 𝐴2 ∗ 𝑘 𝑐𝑜𝑠(𝑘𝛼) + 𝐴𝑘 𝑐𝑜𝑠(𝛼) = (7)
𝜔𝐶
𝐼𝐿𝑅𝑀𝑆 = 1.4608404 𝐴
−𝐼𝑜
−𝐴1 ∗ 𝑘 sen(𝑘𝛽) + 𝐴2 ∗ 𝑘 𝑐𝑜𝑠(𝑘𝛽) + 𝐴𝑘 𝑐𝑜𝑠(𝛽) = (8)
𝜔𝐶
> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 7

Para encontrar el 𝑓𝑝, se necesita de 𝑉𝑔1𝑅𝑀𝑆 𝑒 𝐼𝑔1𝑅𝑀𝑆 . 𝐼𝑔1𝑅𝑀𝑆 1) Graficas de vc e il

se obtiene del valor RMS de la serie de Fourier de:

𝐼𝐿(𝜃) , ∝< 𝜃 < 𝛽

𝐼𝑔 (𝜃) = {
−𝐼𝐿 (𝜃 − 𝜋), 𝜋+∝< 𝜃 < 𝜋 + 𝛽

𝑃𝑎𝑟𝑎 𝐼𝑔1 𝑟𝑚𝑠 𝑠𝑒 ℎ𝑎𝑐𝑒 𝑑𝑒 𝑢𝑠𝑜 𝑑𝑒 𝑙𝑜𝑠 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒𝑠 𝑑𝑒

𝑓𝑜𝑢𝑟𝑖𝑒𝑟 𝑝𝑎𝑟𝑎 𝑙𝑎 𝑓𝑢𝑛𝑑𝑎𝑚𝑒𝑛𝑡𝑎𝑙
𝜋 2𝜋
1
𝑎1 = (∫ 𝐼𝑙(𝜃) cos(Ɵ) 𝑑𝜃 + ∫ −𝐼𝑙(𝜃 − 𝜋) cos(Ɵ) 𝑑𝜃 )
𝜋 0 𝜋
𝑎1 = −0.856591
𝜋 2𝜋
1
𝑏1 = (∫ 𝐼𝑙(𝜃)𝑆𝑒𝑛(𝜃) 𝑑𝜃 + ∫ −𝐼𝑙(𝜃 − 𝜋)𝑆𝑒𝑛(𝜃) 𝑑𝜃)
𝜋 0 𝜋
𝑏1 = 1.6327223

𝐼𝑔1 = −0.856591 cos(𝜃) + 1.6327223𝑠𝑒𝑛(𝜃) Fig. 7. Grafica vc (rojo) y grafica Il (azul)

|𝐼𝑔1 | = 1.8437816

𝐷𝑒𝑠𝑝𝑙𝑎𝑧𝑎𝑚𝑖𝑒𝑛𝑡𝑜
0.856591
arctan (− ) = −0.483164𝑟𝑎𝑑
1.6327223

𝐼𝑔1 = 1.8437816𝑆𝑒𝑛(𝜃 − 0.483164) 2) Código Matlab

clear; clc
𝑃 = 𝑉𝑔𝑟𝑚𝑠 𝐼𝑔1 𝑟𝑚𝑠 cos(𝜃1 − 𝜙1 )
Vp=120*sqrt(2);
1.8437816 w=120*pi;
𝑃 = 120 ∗ cos(−0.483164)
√2 l=40*10^(-3);
c=100*10^(-6);
𝑃 = 138.54108𝑤 Io=1;
k=(1/sqrt(l*c))/(w);
𝑆 = 𝑉𝑔𝑟𝑚𝑠 ∗ 𝐼𝑔𝑟𝑚𝑠 = 120 ∗ 1.0098566 = 175.3008 ak=(Vp*(k^2))/(k^2-1);
syms a1 a2 a b;
𝑃 138.54108𝑤 [SA1, SA2, Sa, Sb]=vpasolve([
𝑃𝐹 = = a1*cos(k*a)+a2*sin(k*a)+ak*sin(a)-
𝑆 175.3008
Vp*sin(a)==0
𝑃𝐹 = 0.7903049
a1*cos(k*b)+a2*sin(k*b)+ak*sin(b)-
Vp*sin(a)-(Io/(w*c))*(pi+a-b)==0
-
a1*k*sin(k*a)+a2*k*cos(k*a)+ak*cos(a)+Io/
(w*c)==0
-
a1*k*sin(k*b)+a2*k*cos(k*b)+ak*cos(b)+Io/
(w*c)==0],[a1 a2 a b],[-inf inf;-inf
inf;0 pi;0 pi])
theta=Sa:pi/2000:Sb;
Vc=SA1*cos(k*theta)+SA2*sin(k*theta)+ak*s
in(theta);
theta1=Sb:pi/2000:(pi+Sa);
Vc2=Vp*sin(Sa)+(Io/(w*c))*(pi+Sa-theta1);

theta0=[0 Sa];
Il0=[0 0];
Il1=(-
SA1*w*c*k*sin(k*theta)+SA2*w*c*k*cos(k*th
eta)+ak*w*c*cos(theta))+Io;
> Fonseca Díaz Cristian David 20161007013, Rincón Miranda Adriana 20161007040< 8

theta2=[Sb (pi+Sa)];
Il2=[0 0];
Ig=[Il0,Il1,Il2];
x=[theta0,theta,theta2];

subplot(2,1,1);
plot(theta,Vc,'r','LineWidth',2);
hold on
plot(theta1,Vc2,'r','LineWidth',2);
title('Vc')

subplot(2,1,2);
plot(x,Ig,'b','LineWidth',2);
title('Il')