distributionnetwork-power2c.../theory/t.m

23 lines
453 B
Mathematica
Raw Permalink Normal View History

clc
clear
P=[12;16;15];
mP=P.*(1+normrnd(0,0.03,length(P),1));
Q=[7;15;10];
mQ=Q.*(1+normrnd(0,0.03,length(Q),1));
rVolt=[1;
1*exp(1j*-123/180*pi);
1*exp(1j*119/180*pi);
];
mVolt=[1*(1+normrnd(0,0.03));
1*(1+normrnd(0,0.03))*exp(1j*-120/180*pi);
1*(1+normrnd(0,0.03))*exp(1j*120/180*pi);
];
a=exp(1j*2*pi/3);
Tp2f=1/3*[1 1 1;
1 a a^2;
1 a^2 a];
rCurr=(P+1j*Q)./rVolt;
mCurr=(mP+1j*mQ)./mVolt;
rI=Tp2f*rCurr
mI=Tp2f*mCurr