flow1.m

来自「前推潮流算法」· M 代码 · 共 34 行

% define values of variable:r,x,P,Q,U
r=[1.05 2.82 1.36 1.19 1.79 1.52 1.31 1.85 1.18 0.46 0.82] 
x=[0.42 1.17 0.57 0.83 1.23 1.04 0.89 1.26 0.49 0.21 0.34]
P=[0 10 8.96 56 12.01 14 35.28 10 35.28 8.96 14]
cosfi=0.7
Q=P*tan(acos(cosfi))
U0=10
U=[10 10 10 10 10 10 10 10 10 10 10]
dPkt=[0 0 0 0 0 0 0 0 0 0 0]
dQkt=[0 0 0 0 0 0 0 0 0 0 0]
Pkt=[0 0 0 0 0 0 0 0 0 0 0]
Qkt=[0 0 0 0 0 0 0 0 0 0 0]

% end to root,U to P and Q;root to end, P and Q to U
computing
% error examining
E=0.0001
g=11;
for c=1:g
    dU(c)=abs(U(c)-U1(c));
    if dU(c)>E
        computing;
    else 
        show;
    end
end

%result
Pkt
Qkt
U