Scribd
Upload
20 views18 pages

GREYHOUND

Uploaded by

Esteban Mauricio Viellard Corrales
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
20 views18 pages

GREYHOUND

Uploaded by

Esteban Mauricio Viellard Corrales
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 18

f//Caso conexion Tipo T

function ZRT=zrt(Za,Zb,Zc)
Z11=Za+Zb
Z12=Zb
Z21=Zb
Z22=Zb+Zc
ZRT=[Z11,Z12;Z21,Z22];
endfunction
//Caso conexion Tipo Pi
function YRT=yrt(Ya,Yb,Yc)
Y11=Ya+Yb
Y12=-Yb
Y21=-Yb
Y22=Yb+Yc
YRT=[Y11,Y12;Y21,Y22]
endfunction
//Caso transformador lineal
function TransformerT=TranT(L1,L2,M,V1,V2)
s=poly(0,"s");
Z=[s*L1,s*M;s*M,s*L2];
V=[V1,V2];
I=inv(Z)*V;
TransformerT=I
endfunction
//Configuración T de un transformador lineal
function Tparam=Tres(L1, L2, M)
Tparam=[L1-M,M,L2-M];
endfunction
//Caso Doble-Fuente
function GAMMA=doublesource(Imp,Vf1,Impf1,Vf2,Impf2)
// GAMMA=[I1;I2;V1;V2]
L=[0;0;Vf1;Vf2]
Z=[Imp,-eye(Imp);Impf1,0,1,0;0,Impf2,0,1]
GAMMA=inv(Z)*L
endfunction
//Transformaciones de parametros dado Z
function Z=Z2Z(Z)
endfunction
function Y=Z2Y(Z)
Y=inv(Z)
endfunction
function H=Z2H(Z)
h11=(det(Z))/(Z(2,2))
h12=Z(1,2)/Z(2,2)
h21=-Z(2,1)/Z(2,2)
h22=1/Z(2,2)
H=[h11,h12,h12,h21,h22]
endfunction
function G=Z2G(Z)
g11=1/(Z(1,1))
g12=-Z(1,2)/Z(1,1)
g21=Z(2,1)/Z(1,1)
g22=det(Z)/Z(1,1)
G=[g11,g12,g12,g21,g22]
endfunction
function T=Z2T(Z)
A=Z(1,1)/(Z(2,1))
B=det(Z)/Z(2,1)
C=1/Z(2,1)
D=Z(2,2)/Z(2,1)
T=[A,B,C,D]
endfunction
function T=Z2T(Z)
A=Z(1,1)/(Z(2,1))
B=det(Z)/Z(2,1)
C=1/Z(2,1)
D=Z(2,2)/Z(2,1)
T=[A,B,C,D]
endfunction
function t=Z2t(Z)
a=Z(2,2)/(Z(1,2))
b=det(Z)/Z(1,2)
c=1/Z(1,2)
d=Z(1,1)/Z(1,2)
t=[a,b,c,d]
endfunction
///////////////////
///////////////////
//Transformaciones de parametros dado Y
function Z0=Y2Z(Y)
z11=Y(2,2)/det(Y)
z12=-Y(1,2)/det(Y)
z21=-Y(2,1)/det(Y)
z22=Y(1,1)/det(Y)
Z0=[z11,z12,z21,z22]
endfunction
function Y0=Y2Y(Y)
endfunction
function H0=Y2H(Y)
h11=1/(Y(1,1))
h12=-Y(1,2)/Y(1,1)
h21=Y(2,1)/Y(1,1)
h22=det(Y)/Y(1,1)
H0=[h11,h12,h12,h21,h22]
endfunction
function G0=Y2G(Y)
g11=det(Y)/(Y(2,2))
g12=Y(1,2)/Y(2,2)
g21=-Y(2,1)/Y(2,2)
g22=1/Y(2,2)
G=[g11,g12,g12,g21,g22]
endfunction
function T0=Y2T(Y)
A=-Y(2,2)/(Y(2,1))
B=-1/Y(2,1)
C=-det(Y)/Y(2,1)
D=-Y(1,1)/Y(2,1)
T0=[A,B,C,D]
endfunction
function t0=Y2t(Y)
a=-Y(1,1)/(Y(1,2))
b=-1/Y(1,2)
c=-det(Y)/Y(1,2)
d=-Y(2,2)/Y(1,2)
t0=[a,b,c,d]
endfunction
///////////////////
///////////////////
//Transformaciones de parametros dado H
function Z1=H2Z(H)
z11=det(H)/H(2,2)
z12=H(1,2)/H(2,2)
z21=-H(2,1)/H(2,2)
z22=1/H(2,2)
Z1=[z11,z12,z21,z22]
endfunction
function Y1=H2Y(H)
y11=1/H(1,1)
y12=-H(1,2)/H(1,1)
y21=H(2,1)/H(1,1)
y22=det(H)/H(1,1)
Y1=[y11,y12,y21,y22]
endfunction
function H1=H2H(H)
endfunction
function G1=H2G(H)
g11=H(2,2)/det(H)
g12=-H(1,2)/det(H)
g21=-H(2,1)/det(H)
g22=H(1,1)/det(H)
G1=[g11,g12,g12,g21,g22]
endfunction
function T1=H2T(H)
A=-det(H)/(H(2,1))
B=-1/Y(2,1)
C=-H(2,2)/H(2,1)
D=-1/H(2,1)
T1=[A,B,C,D]
endfunction
function t1=H2t(H)
a=1/(H(1,2))
b=H(1,1)/H(1,2)
c=H(2,2)/H(1,2)
d=det(H)/H(1,2)
t1=[a,b,c,d]
endfunction
///////////////////
///////////////////
//Transformaciones de parametros dado G
function Z2=G2Z(G)
z11=1/G(1,1)
z12=-G(1,2)/G(1,1)
z21=G(2,1)/G(1,1)
z22=det(G)/G(1,1)
Z2=[z11,z12,z21,z22]
endfunction
function Y2=G2Y(G)
y11=det(G)/G(2,2)
y12=G(1,2)/G(2,2)
y21=-G(2,1)/G(2,2)
y22=1/G(2,2)
Y2=[y11,y12,y21,y22]
endfunction
function H2=G2H(G)
h11=G(2,2)/det(G)
h12=-G(1,2)/det(G)
h21=-G(2,1)/det(G)
h22=G(1,1)/det(G)
H2=[h11,h12,h21,h22]
endfunction
function G2=G2G(G)
endfunction
function T2=G2T(G)
A=1/(G(2,1))
B=G(2,2)/G(2,1)
C=G(1,1)/G(2,1)
D=det(G)/G(2,1)
T2=[A,B,C,D]
endfunction
function t2=G2t(G)
a=-det(G)/(G(1,2))
b=-G(2,2)/G(1,2)
c=-G(1,1)/G(1,2)
d=-1/G(1,2)
t2=[a,b,c,d]
endfunction
///////////////////
///////////////////
//Transformaciones de parámetros dado T
function Z3=T2Z(T)
z11=T(1,1)/T(2,1)
z12=det(T)/T(2,1)
z21=1/T(2,1)
z22=T(2,2)/T(2,1)
Z3=[z11,z12;z21,z22]
endfunction
function Y3=T2Y(T)
y11=T(2,2)/T(1,2)
y12=-det(T)/T(1,2)
y21=-1/T(1,2)
y22=T(1,1)/T(1,2)
Y3=[y11,y12;y21,y22]
endfunction
function H3=T2H(T)
h11=T(1,2)/T(2,2)
h12=det(T)/T(2,2)
h21=-1/T(2,2)
h22=T(2,1)/T(2,2)
H3=[h11,h12;h21,h22]
endfunction
function G3=T2G(T)
g11=T(2,1)/(T(1,1))
g12=-det(T)/T(1,1)
g21=1/T(1,1)
g22=T(1,2)/T(1,1)
G3=[g11,g12;g21,g22]
endfunction
function T3=T2T(T)
endfunction
function t3=T2t(T)
a=T(2,2)/det(T)
b=T(1,2)/det(T)
c=T(2,1)/det(T)
d=T(1,1)/det(T)
t3=[a,b;c,d]
endfunction
///////////////////
///////////////////
//Transformaciones de parametros dado t
function Z4=t2Z(t)
z11=t(2,2)/t(2,1)
z12=1/t(2,1)
z21=det(t)/T(2,1)
z22=t(1,1)/T(2,1)
Z4=[z11,z12,z21,z22]
endfunction
function Y4=t2Y(t)
y11=t(1,1)/T(1,2)
y12=-1/T(1,2)
y21=-det(t)/T(1,2)
y22=t(2,2)/T(1,2)
Y4=[y11,y12,y21,y22]
endfunction
function H4=t2H(t)
h11=t(1,2)/t(1,1)
h12=1/t(1,1)
h21=det(t)/t(1,1)
h22=T(2,1)/t(1,1)
H4=[h11,h12,h21,h22]
endfunction
function G4=t2G(t)
g11=t(2,1)/(t(2,2))
g12=-1/t(2,2)
g21=det(t)/t(2,2)
g22=-t(1,2)/det(t)
G4=[g11,g12,g21,g22]
endfunction
function T4=t2T(t)
A=T(2,2)/det(t)
B=T(1,2)/det(t)
C=T(2,1)/det(t)
D=T(1,1)/det(t)
T4=[A,B,C,D]
endfunction
function t4=t2t(t)
endfunction

SIMPLIFICADOR DEL MEJOR PROFE


function [xsim]=simpFrac(num, den)
a=coeff(num)
b=coeff(den)
divisor=a(size(a,'*'))
a=a/divisor
b=b/divisor
nsim=poly(a,"s","coeff")
dsim=poly(b,"s","coeff")
xsim=nsim/dsim
endfunction
function [Xsim]=miSimp(X)
nf=size(X,1)
nc=size(X,2)
Xsim=zeros(nf,nc)
for f=[1:nf]
for c=[1:nc]
num=X('num')(f,c)
den=X('den')(f,c)
Xsim(f,c)=simpFrac(num,den)
end
end
endfunction

// VERTICAL
function V = vertical(z)
V = [1,0;1/z,1]
endfunction

// HORIZONTAL
function H = horizontal(z)
H = [1,z;0,1]
endfunction

// TRANSFORMADOR
function TRANS = transformer(a)
TRANS = [1/a,0;0,a]
endfunction

################
function YRT=yrt(Ya,Yb,Yc)
Y11=Ya+Yb
Y12=-Yb
Y21=-Yb
Y22=Yb+Yc
YRT=[Y11,Y12;Y21,Y22]
endfunction

function ZRT=zrt(Za,Zb,Zc)
Z11=Za+Zb
Z12=Zb
Z21=Zb
Z22=Zb+Zc
ZRT=[Z11,Z12;Z21,Z22];
endfunction

function T=Z2T(Z)
A=Z(1,1)/(Z(2,1))
B=det(Z)/Z(2,1)
C=1/Z(2,1)
D=Z(2,2)/Z(2,1)
T=[A,B;C,D]
endfunction

function T0=Y2T(Y)
A=-Y(2,2)/(Y(2,1))
B=-1/Y(2,1)
C=-det(Y)/Y(2,1)
D=-Y(1,1)/Y(2,1)
T0=[A,B;C,D]
endfunction

function GAMMA=doublesource(Imp,Vf1,Impf1,Vf2,Impf2)
// GAMMA=[I1;I2;V1;V2]
L=[0;0;Vf1;Vf2]
Z=[Imp,-eye(Imp);Impf1,0,1,0;0,Impf2,0,1]
GAMMA=inv(Z)*L
endfunction

function Z3=T2Z(T)
z11=T(1,1)/T(2,1)
z12=det(T)/T(2,1)
z21=1/T(2,1)
z22=T(2,2)/T(2,1)
Z3=[z11,z12;z21,z22]
endfunction

s = poly(0,"s")
y_a = yrt(1/(10*s),s/10,1/(10*s))
z_b = zrt(10/s,10*s,10/s)
z_c = zrt(0,10*s,0)

t_a = Y2T(y_a)
t_b = Z2T(z_b)
t_c = Z2T(z_c)

t_r = t_a*t_b*t_c

z_r = T2Z(t_r)

v0 = doublesource(z_r,1,0,0,1)

EJERCICIO 2 CLASE 21 NOV MÉTODO 1


VF=50;
z1=1j*5;
z2=3;
z3=-1j*4;
z4=1j*10;
z5=5;
m=1j*6;

eqn1=-VF+z1*I1-m*I2+z2*(I1-I2)+z3*(I1-I2)==0;
eqn2= z5*I2+z3*(I2-I1)+z2*(I2-I1)+z4*I2-m*I1==0;

a = solve(eqn1,eqn2,I1,I2);

L1 = 1j*1;
L2 = 4*1j;
M = 0.9*sqrt(L1*L2);
R1 = 0.1;
R2 = 1;

imp1 = L1-M;
imp2 = L2-M;

op1 = R2+imp2;

paralelo = ((op1*M)/(op1+M));

final = paralelo+R1+imp1;

disp(final);

PASAR A ADMITANCIAS

function Y=ad(a)
Y=1/a
endfunction

function YRT=yrt(Ya,Yb,Yc)
Y11=Ya+Yb
Y12=-Yb
Y21=-Yb
Y22=Yb+Yc
YRT=[Y11,Y12;Y21,Y22]
endfunction
function H = horizontal(z)
H = [1,z;0,1]
endfunction
function T0=Y2T(Y)
A=-Y(2,2)/(Y(2,1))
B=-1/Y(2,1)
C=-det(Y)/Y(2,1)
D=-Y(1,1)/Y(2,1)
T0=[A,B;C,D]
endfunction
function [xsim]=simpFrac(num, den)
a=coeff(num)
b=coeff(den)
divisor=a(size(a,'*'))
a=a/divisor
b=b/divisor
nsim=poly(a,"s","coeff")
dsim=poly(b,"s","coeff")
xsim=nsim/dsim
endfunction

function [Xsim]=miSimp(X)
nf=size(X,1)
nc=size(X,2)
Xsim=zeros(nf,nc)
for f=[1:nf]
for c=[1:nc]
num=X('num')(f,c)
den=X('den')(f,c)
Xsim(f,c)=simpFrac(num,den)
end
end
endfunction
function T=Z2T(Z)
A=Z(1,1)/(Z(2,1))
B=det(Z)/Z(2,1)
C=1/Z(2,1)
D=Z(2,2)/Z(2,1)
T=[A,B;C,D]
endfunction
function Y=ad(a)
Y=1/a
endfunction
function res=paralelo(x,y)
res=(x*y)/(x+y)
endfunction
function Z3=T2Z(T)
z11=T(1,1)/T(2,1)
z12=det(T)/T(2,1)
z21=1/T(2,1)
z22=T(2,2)/T(2,1)
Z3=[z11,z12;z21,z22]
endfunction
j=%i;
ZA=-2*j;
ZB=paralelo(1,2*j);
ZC=-2*j;
PI=yrt(ad(ZA),ad(ZB),ad(ZC));
A=horizontal(1);
PIT=Y2T(PI);
l=A*A;
res=l*PIT;
resf=T2Z(res);
disp(res);

EJERCICIO DE CLASE

function Tparam=Tres(L1, L2, M)


Tparam=[L1-M,M,L2-M];
endfunction

function ZRT=zrt(Za, Zb, Zc)


Z11=Za+Zb
Z12=Zb
Z21=Zb
Z22=Zb+Zc
ZRT=[Z11,Z12;Z21,Z22];
endfunction

function ZTH=zth(R1, R2, ZL1, ZL2, ZL, M)


ZTH=((R1*ZL2)+(R1*R2)+(R1*ZL)+(ZL1*ZL2)+(R2*ZL1)+(ZL1*ZL)-(M*M))/(ZL2+R2+ZL)
endfunction
function GAMMA=doublesource(Imp, Vf1, Impf1, Vf2, Impf2)
// GAMMA=[I1;I2;V1;V2]
L=[0;0;Vf1;Vf2]
Z=[Imp,-eye(Imp);Impf1,0,1,0;0,Impf2,0,1]
GAMMA=inv(Z)*L
endfunction

///DATOS/////
j=%i;
ZL1=1*j;
ZL2=4*j;
k=0.9;
M=k*sqrt(ZL1*ZL2);
port=Tres(ZL1,ZL2,M);
Z=zrt(port(1),port(2),port(3));
vf=10*exp(j*0/180);
res=doublesource(Z,vf,0.1,0,1);
zth=vf/res(1);
disp(zth);
zth2=zth(0.1,0,ZL1,ZL2,1,M);
disp(zth2);

s = poly(0,"s")
function [xsim]=simpFrac(num, den)
a=coeff(num)
b=coeff(den)
divisor=a(size(a,'*'))
a=a/divisor
b=b/divisor
nsim=poly(a,"s","coeff")
dsim=poly(b,"s","coeff")
xsim=nsim/dsim
endfunction

function Z3=T2Z(T)
z11=T(1,1)/T(2,1)
z12=det(T)/T(2,1)
z21=1/T(2,1)
z22=T(2,2)/T(2,1)
Z3=[z11,z12;z21,z22]
endfunction

function GAMMA=doublesource(Imp, Vf1, Impf1, Vf2, Impf2)


// GAMMA=[I1;I2;V1;V2]
L=[0;0;Vf1;Vf2]
Z=[Imp,-eye(Imp);Impf1,0,1,0;0,Impf2,0,1]
GAMMA=inv(Z)*L
endfunction

function T=Z2T(Z)
A=Z(1,1)/(Z(2,1))
B=det(Z)/Z(2,1)
C=1/Z(2,1)
D=Z(2,2)/Z(2,1)
T=[A,B;C,D]
endfunction

function [Xsim]=miSimp(X)
nf=size(X,1)
nc=size(X,2)
Xsim=zeros(nf,nc)
for f=[1:nf]
for c=[1:nc]
num=X('num')(f,c)
den=X('den')(f,c)
Xsim(f,c)=simpFrac(num,den)
end
end
endfunction

function ZRT=zrt(Za,Zb,Zc)
Z11=Za+Zb
Z12=Zb
Z21=Zb
Z22=Zb+Zc
ZRT=[Z11,Z12;Z21,Z22];
endfunction

a=zrt(0,10*s,10/s)
b=zrt(0,10*s,0)

T_a=Z2T(a)
T_b=Z2T(b)

mult=((T_a)^3)*(T_b)
ParZ=T2Z(mult)

a=doublesource(ParZ, 1, 0, 0, 1)

#####################################################

function YRT=yrt(Ya, Yb, Yc)


Y11=Ya+Yb
Y12=-Yb
Y21=-Yb
Y22=Yb+Yc
YRT=[Y11,Y12;Y21,Y22]
endfunction

function ZRT=zrt(Za, Zb, Zc)


Z11=Za+Zb
Z12=Zb
Z21=Zb
Z22=Zb+Zc
ZRT=[Z11,Z12;Z21,Z22];
endfunction

function T=Z2T(Z)
A=Z(1,1)/(Z(2,1))
B=det(Z)/Z(2,1)
C=1/Z(2,1)
D=Z(2,2)/Z(2,1)
T=[A,B;C,D]
endfunction

function T0=Y2T(Y)
A=-Y(2,2)/(Y(2,1))
B=-1/Y(2,1)
C=-det(Y)/Y(2,1)
D=-Y(1,1)/Y(2,1)
T0=[A,B;C,D]
endfunction

function GAMMA=doublesource(Imp, Vf1, Impf1, Vf2, Impf2)


// GAMMA=[I1;I2;V1;V2]
L=[0;0;Vf1;Vf2]
Z=[Imp,-eye(Imp);Impf1,0,1,0;0,Impf2,0,1]
GAMMA=inv(Z)*L
endfunction

function Z3=T2Z(T)
z11=T(1,1)/T(2,1)
z12=det(T)/T(2,1)
z21=1/T(2,1)
z22=T(2,2)/T(2,1)
Z3=[z11,z12;z21,z22]
endfunction

function [xsim]=simpFrac(num, den)


a=coeff(num)
b=coeff(den)
divisor=a(size(a,'*'))
a=a/divisor
b=b/divisor
nsim=poly(a,"s","coeff")
dsim=poly(b,"s","coeff")
xsim=nsim/dsim
endfunction

function [Xsim]=miSimp(X)
nf=size(X,1)
nc=size(X,2)
Xsim=zeros(nf,nc)
for f=[1:nf]
for c=[1:nc]
num=X('num')(f,c)
den=X('den')(f,c)
Xsim(f,c)=simpFrac(num,den)
end
end
endfunction

s = poly(0,"s");
y_a = yrt(1/(10*s),s/10,1/(10*s));
z_b = zrt(10/s,10*s,10/s);
z_c = zrt(0,10*s,0);
t_a = Y2T(y_a);
t_b = Z2T(z_b);
t_c = Z2T(z_c);
t_r = t_a*t_b*t_c;

// APLICAR MISIMP EN EL T EQUIVALENTE


z_r = miSimp(T2Z(t_r));
v0=doublesource(z_r,1,0,0,1);
H=syslin("c",v0(4,1));
bode(H);

APLICAR “miSimp” a la matriz equivalente.

s=poly(0,"s");

//Simplificación dada por el profesor Duarte:

function [xsim]=simpFrac(num, den)


a=coeff(num)
b=coeff(den)
divisor=a(size(a,'*'))
a=a/divisor
b=b/divisor
nsim=poly(a,"s","coeff")
dsim=poly(b,"s","coeff")
xsim=nsim/dsim
endfunction

function [Xsim]=miSimp(X)
nf=size(X,1)
nc=size(X,2)
Xsim=zeros(nf,nc)
for f=[1:nf]
for c=[1:nc]
num=X('num')(f,c)
den=X('den')(f,c)
Xsim(f,c)=simpFrac(num,den)
end
end
endfunction

// Caso conexión Tipo T


function Impedancia_T = zct(Za, Zb, Zc)
// Calcula la matriz de impedancia equivalente para una conexión en T.
Z11 = Za + Zb;
Z12 = Zb;
Z21 = Zb;
Z22 = Zb + Zc;
Impedancia_T = [Z11, Z12; Z21, Z22];
endfunction

a=zct(4*s+1.6*s, -1.6*s, 1*s+1.6*s)

function Matriz_Parametros_T = ZaT(Z)


// Calcula la matriz de parámetros T a partir de la matriz de impedancia Z.
A = Z(1, 1) / Z(2, 1);
B = det(Z) / Z(2, 1);
C = 1 / Z(2, 1);
D = Z(2, 2) / Z(2, 1);
Matriz_Parametros_T = [A, B; C, D];
endfunction

b=ZaT(a)
function V = vertical(z)
// Matriz de transformación vertical.
V = [1, 0; 1/z, 1];
endfunction

function H = horizontal(z)
// Matriz de transformación horizontal.
H = [1, z; 0, 1];
endfunction

c=horizontal(1)
d=vertical(4)

//Multiplicamos los T
x=c*b*d
//Pasamos a Zeq

// Transformaciones de parámetros dado T


function Impedancia_T = TaZ(T)
// Calcula la matriz de impedancia a partir de la matriz de parámetros T.
Z11 = T(1, 1) / T(2, 1);
Z12 = det(T) / T(2, 1);
Z21 = 1 / T(2, 1);
Z22 = T(2, 2) / T(2, 1);
Impedancia_T = [Z11, Z12; Z21, Z22];
endfunction

Zeq=TaZ(x)
Simp=miSimp(Zeq)

//Utilizamos doble fuente:


// Caso Doble-Fuente
function Corrientes_DobleFuente = doblefuente(Zeq, Vf1, Zf1, Vf2, Zf2)
// Calcula corrientes en un circuito con dos fuentes.
// Corrientes_DobleFuente = [I1; I2; V1; V2]
L = [0; 0; Vf1; Vf2];
Z = [Zeq, -eye(Zeq); Zf1, 0, 1, 0; 0, Zf2, 0, 1];
Corrientes_DobleFuente = inv(Z)*L;
endfunction

DOBLE=doblefuente(Simp, 0, 1, 0, 4)

########## DESDE AQUÍ EMPIEZA ########


// PUNTO 2
L1 = 1;
L2 = 4;
ZA = 1;
ZB = 2;
M = 1.6;

gama = L1+M;
betha = L2+M;
epsilon = -M;

ZTH = (ZA+gama)+((epsilon*(betha+epsilon))/(epsilon+(betha+epsilon)));
disp(ZTH)

l1=4
l2=1
m=1.6

// Caso conexión Tipo T


function Impedancia_T = zct(Za, Zb, Zc)
// Calcula la matriz de impedancia equivalente para una conexión en T.
Z11 = Za + Zb;
Z12 = Zb;
Z21 = Zb;
Z22 = Zb + Zc;
Impedancia_T = [Z11, Z12; Z21, Z22];
endfunction
s=poly(0,"s")
zequ=zct(l1*s+(m*s), -m*s,l2*s+(m*s))

function [xsim]=simpFrac(num, den)


a=coeff(num)
b=coeff(den)
divisor=a(size(a,'*'))
a=a/divisor
b=b/divisor
nsim=poly(a,"s","coeff")
dsim=poly(b,"s","coeff")
xsim=nsim/dsim
endfunction
function [Xsim]=miSimp(X)
nf=size(X,1)
nc=size(X,2)
Xsim=zeros(nf,nc)
for f=[1:nf]
for c=[1:nc]
num=X('num')(f,c)
den=X('den')(f,c)
Xsim(f,c)=simpFrac(num,den)
end
end
endfunction

function DobleFuente=doublesource(Imp,Vf1,Impf1,Vf2,Impf2)
// DobleFuente=[I1;I2;V1;V2]
L=[0;0;Vf1;Vf2]
Z=[Imp,-eye(Imp);Impf1,0,1,0;0,Impf2,0,1]
DobleFuente=inv(Z)*L
endfunction

d=doublesource(zequ,1,1,0,4)

//Hallar Zth: V=IR


ZTH=1/d(1,1)

You might also like