牛頓拉夫遜法&PQ分解法計算潮流MATLAB程序計算模型。
本程序模型基于MATLAB進行潮流計算,建議先安裝matpower插件(MATLAB中非常重要的潮流計算的插件)。
本程序可進行牛拉法和PQ分解法潮流計算的切換,對比潮流計算的結果。很適合新手學習使用。
部分程序實例:
.rtcContent { padding: 30px; } .lineNode {font-size: 10pt; font-family: Menlo, Monaco, Consolas, "Courier New", monospace; font-style: normal; font-weight: normal; }
tic %開始計時
clear %清除工作空間變量
clc %清屏
%%?********************數據輸入(全局標幺值p.u.)*********************
%% 基本參數設置
SB=100; %基準容量,單位MVA
mode=2; % 1_極坐標下, 2_PQ分解法
Tmax=100; %最大迭代次數
limit=1.0e-5; %要求精度
%% 母線參數:母線編號 節點類型 初始電壓 初始相角(度) 負荷有功 負荷無功 基準電壓(有名值)
% 注意:1.母線填寫順序遵循:1_PQ節點-->2_PV節點-->3_平衡節點
% ? ? ? 2.注入功率為+,消耗功率為-
Bus=[1 1 1 0 -0.3 -0.18 115;
? ? 2 1 1 0 -0.55 -0.13 10.5;
? ? 3 2 1.1 0 0.5 0 115;
? ? 4 3 1.05 0 0 0 115;];
%% 線路參數:I側母線 J側母線 阻抗 1/2接地導納
Line=[1 2 0.1+0.4i 0.01528i;
? ? ? 1 4 0.12+0.5i 0.0192i;
? ? ? 2 4 0.08+0.4i 0.01413i;];
%% 變壓器參數:I側母線 J側母線 阻抗(歸算I側復數) 變比(I->J為k:1)
Trans=[3 1 0.3i 1.1];
%% 加接地無功補償裝置: 節點 導納
% 靜止電容器,設置單組容量和總組數,以投入組數為自變量進行無功優化求解
Cap=[2 0];
%%?**************************參數運算*******************************
%% 變壓器π型等效導納參數
%Trans_pi:1_I側母線 2_J側母線
Trans_pi(:,1:2)=Trans(:,1:2);
Trans_pi(:,3)=1/(Trans(:,3)./Trans(:,4));
Trans_pi(:,4)=1/(Trans(:,3)./(1-Trans(:,4)));
Trans_pi(:,5)=1/(Trans(:,3)./(Trans(:,4).^2-Trans(:,4)));
%% 計算PQ節點數量m(去除PV節點和平衡節點),其中包含浮游節點,其初始P&Q=0
n=size(Bus,1); %總節點數
m=n-1;% 去除1個平衡節點數
for i=1:n
? ? if Bus(i,2)==2 %除去PV節點數
? ? ? ? m=m-1;
? ? end
end
%% 形成P,Q,U,θ向量,以節點順序排列
for i=1:n
? ? P(Bus(i,1))=Bus(i,5);
? ? Q(Bus(i,1))=Bus(i,6);
? ? U(Bus(i,1))=Bus(i,3);
? ? cita(Bus(i,1))=Bus(i,4);
潮流計算結果: