Tarea potencia

clear all
clc
Z=[0.21248+1.3612*i 0.09528+0.5571*i 0.09528+0.493*i 0.09528+0.9135*i 0.09528+0.5517*i 0.09528+0.4892*i;
0.09528+0.5571*i 0.21248+1.3612*i 0.09528+0.5571*i 0.09528+0.5626*i 0.09528+0.9135*i 0.09528+0.5517*i;
0.09528+0.493*i 0.09528+0.5571*i 0.21248+1.3612*i 0.9529+0.4968*i 0.09528+0.5626*i 0.09528+0.9135*i;
0.09528+0.9135*i 0.09528+0.5626*i 0.9529+0.4968*i0.21248+1.3612*i 0.09528+0.5571*i 0.09528+0.493*i;
0.09528+0.5517*i 0.09528+0.9135*i 0.09528+0.5626*i 0.09528+0.5571*i 0.21248+1.3612*i 0.09528+0.5571*i;
0.09528+0.4892*i 0.09528+0.5517*i 0.09528+0.9135*i 0.09528+0.493*i 0.09528+0.5571*i 0.21248+1.3612*i]
Zs=Z(1:3,1:3)
Zm=Z(1:3,4:6)

Zk=(Zs-Zm)-(Zm.'-Zs);
Zk=(Zk*-1)
Zkin=inv(Zk)

ZN=(Zm-Zs)*(Zkin);
ZN=ZN*(Zm.'-Zs)
ZN=ZN*-1Znew=Zs-ZN


a=cos(120*pi/180)+i*sin(120*pi/180);
A=[1 1 1;1 a*a a;1 a a*a];
Ainv=(1/3)*(inv(A));
%Z012=Ainv*Zs*A


'SE HACE LA TRASPOSICION DE LINEAS'
Z1=Znew

'MATRIZ DE ROTACION'
Ro=[0 0 1;1 0 0;0 1 0];

'MATRIZ INVERSA DE LA DE ROTACION'
RINV=inv(Ro);

'CALCULAMOS Zcab:'
'Zcab=(RINV*(Ro*Z1))'
Z2=(RINV*(Z1*Ro))

'Zbca SE CALCULA DE LA MISMA MANERA QUE ELANTERIOR(SE REPITE EL PROCESO)'
Z3=(RINV*(Z2*Ro))

'FINALMENTE CALCULAMOS Zt:'

'Zt=(Zabc+Zcab+Zbca)'

Zt=(Z1+Z2+Z3)

Z012=Ainv*Zt*A

Z =

Columns 1 through 3

0.2125 + 1.3612i 0.0953 + 0.5571i 0.0953 + 0.4930i
0.0953 + 0.5571i 0.2125 + 1.3612i 0.0953 + 0.5571i
0.0953 + 0.4930i 0.0953 + 0.5571i 0.2125 + 1.3612i
0.0953 + 0.9135i 0.0953 + 0.5626i 0.9529 +0.4968i
0.0953 + 0.5517i 0.0953 + 0.9135i 0.0953 + 0.5626i
0.0953 + 0.4892i 0.0953 + 0.5517i 0.0953 + 0.9135i

Columns 4 through 6

0.0953 + 0.9135i 0.0953 + 0.5517i 0.0953 + 0.4892i
0.0953 + 0.5626i 0.0953 + 0.9135i 0.0953 + 0.5517i
0.9529 + 0.4968i 0.0953 + 0.5626i 0.0953 + 0.9135i
0.2125 + 1.3612i 0.0953 + 0.5571i 0.0953 + 0.4930i
0.0953 +0.5571i 0.2125 + 1.3612i 0.0953 + 0.5571i
0.0953 + 0.4930i 0.0953 + 0.5571i 0.2125 + 1.3612i

Zs =

0.2125 + 1.3612i 0.0953 + 0.5571i 0.0953 + 0.4930i
0.0953 + 0.5571i 0.2125 + 1.3612i 0.0953 + 0.5571i
0.0953 + 0.4930i 0.0953 + 0.5571i 0.2125 + 1.3612i

Zm =

0.0953 + 0.9135i 0.0953 + 0.5517i 0.0953 + 0.4892i
0.0953 + 0.5626i 0.0953 +0.9135i 0.0953 + 0.5517i
0.9529 + 0.4968i 0.0953 + 0.5626i 0.0953 + 0.9135i

Zk =

-0.2344 - 0.8954i 0 + 0.0001i 0.8576 + 0.0000i
0 + 0.0001i -0.2344 - 0.8954i 0 + 0.0001i
0.8576 + 0.0000i 0 + 0.0001i -0.2344 - 0.8954i

Zkin =

-0.0120 + 0.6007i 0.0000 + 0.0001i 0.5356 + 0.1517i
0.0000 + 0.0001i -0.2736 + 1.0452i 0.0000 + 0.0001i0.5356 + 0.1517i 0.0000 + 0.0001i -0.0120 + 0.6007i

ZN =

-0.0628 - 0.1135i -0.0010 - 0.0013i 0.1164 - 0.0296i
-0.0010 - 0.0013i -0.0586 - 0.2239i -0.0010 - 0.0013i
0.1164 - 0.0296i -0.0010 - 0.0013i -0.0628 - 0.1135i

ZN =

0.0628 + 0.1135i 0.0010 + 0.0013i -0.1164 + 0.0296i
0.0010 + 0.0013i 0.0586 + 0.2239i 0.0010 + 0.0013i
-0.1164 + 0.0296i 0.0010+ 0.0013i 0.0628 + 0.1135i

Znew =

0.1497 + 1.2477i 0.0942 + 0.5558i 0.2117 + 0.4634i
0.0942 + 0.5558i 0.1539 + 1.1373i 0.0942 + 0.5558i
0.2117 + 0.4634i 0.0942 + 0.5558i 0.1497 + 1.2477i

ans =

SE HACE LA TRASPOSICION DE LINEAS

Z1 =

0.1497 + 1.2477i 0.0942 + 0.5558i 0.2117 + 0.4634i
0.0942 + 0.5558i 0.1539 + 1.1373i 0.0942 + 0.5558i0.2117 + 0.4634i 0.0942 + 0.5558i 0.1497 + 1.2477i

ans =

MATRIZ DE ROTACION

ans =

MATRIZ INVERSA DE LA DE ROTACION

ans =

CALCULAMOS Zcab:

ans =

Zcab=(RINV*(Ro*Z1))

Z2 =

0.1539 + 1.1373i 0.0942 + 0.5558i 0.0942 + 0.5558i
0.0942 + 0.5558i 0.1497 + 1.2477i 0.2117 + 0.4634i
0.0942 + 0.5558i 0.2117 + 0.4634i 0.1497 + 1.2477i

ans =

Zbca SE...