?? powerflow.m
字號:
function V=powerflow(A0)
A0T = A0'; % 形成 A0 的轉置
S = [0;-NodeData(:,2) - i*NodeData(:,3)]; % 形成 S
ZL = [0;BranchData(:,3) + i*BranchData(:,4)]; % 形成 ZL
V = ones(NN,1);
V(1) = 1.05; % 各個節點電壓賦初值
IL(NN,1) = -conj(S(NN) / V(NN)); % 最末支路電流賦初值
Delta = 1; % 收斂判據賦初值
TempV = V; % 賦初值,用于記憶上次迭代結果
while Delta > 1e-8
IN = conj(S ./ V); % 節點注入電流
for n = 1:NN-1
IL(NN-n) = A0(NN-n,NN-n+1:end) * IL(NN-n+1:end) - IN(NN-n);
end % 電流回代過程
for n = 2:NN
V(n) = A0T(n,1:n-1) * V(1:n-1) - ZL(n) * IL(n);
end % 電壓前推過程
Delta = max(abs(V-TempV)); % 更新收斂判據
TempV = V; % 記憶迭代結果
end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -