Matlab Circuitos

En este laboratorio vimos el funcionamiento interno de transformadores, se nos dio un circuito en el que debíamos hallar el voltaje resultante, que dependía de la impedancia Z colocada al final delcircuito.
Esta impedancia variaba en cada inciso. Estos son los códigos que se utilizaron:
t=[0:0.0001:0.04]
V1=220
f=50
w=pi*2*f
R=1.2
Xl=j*3.8
Rc=890
Xc=-j*5000
R=50
L=0.255
C=0.0000354Z=R
%Verificar el funcionamiento de transformadores para carga R
%Circuito Equivalente
Z1=1.2+j*3.8
Z2=890*j*5/(890+j*5)
Z1max=abs(Z1)
desfase1=angle(Z1)
Z2max=abs(Z2)
desfase2=angle(Z2)Z3=Z*Z2/(Z+Z2)
V2=V1*Z3/(Z3+Z1)
V2max=abs(V2)
DV=angle(V2)
I1=V1/(Z1+Z3)
I2=I1*Z3/Z
I2max=abs(I2)
DI=angle(I2)
GV2=V2max*cos(w*t+DV)
GI2=I2max*cos(w*t+DI)
plot(t,GV2,t,GI2*30)

La gráficaresultante es la siguiente:

En el segundo ejercicio la impedancia conectada era de jwL, se utilizó el siguiente código:
t=[0:0.0001:0.04]
V1=220
f=50
w=pi*2*f
R=1.2
Xl=j*3.8
Rc=890Xc=-j*5000
R=50
L=0.255
C=0.0000354
Z=R-j*w*L
%Verificar el funcionamiento de transformadores para carga R-L
%Circuito Equivalente
Z1=1.2+j*3.8
Z2=890*j*5/(890+j*5)
Z1max=abs(Z1)desfase1=angle(Z1)
Z2max=abs(Z2)
desfase2=angle(Z2)
Z3=Z*Z2/(Z+Z2)
V2=V1*Z3/(Z3+Z1)
V2max=abs(V2)
DV=angle(V2)
I1=V1/(Z1+Z3)
I2=I1*Z3/Z
I2max=abs(I2)
DI=angle(I2)
GV2=V2max*cos(w*t+DV)GI2=I2max*cos(w*t+DI)
plot(t,GV2,t,GI2*30)

Y finalmente en el tercero utilizamos una impedancia capacitiva, se utilizó el siguiente código:
t=[0:0.0001:0.04]
V1=220
f=50
w=pi*2*f
R=1.2
Xl=j*3.8
Rc=890Xc=-j*5000
R=50
L=0.255
C=0.0000354
Z=R-j*(1/(w*C))
%Verificar el funcionamiento de transformadores para carga R,R-L,R-C
%Circuito Equivalente
Z1=1.2+j*3.8
Z2=890*j*5/(890+j*5)
Z1max=abs(Z1)desfase1=angle(Z1)
Z2max=abs(Z2)
desfase2=angle(Z2)
Z3=Z*Z2/(Z+Z2)
V2=V1*Z3/(Z3+Z1)
V2max=abs(V2)
DV=angle(V2)
I1=V1/(Z1+Z3)
I2=I1*Z3/Z
I2max=abs(I2)
DI=angle(I2)
GV2=V2max*cos(w*t+DV)...