pq分解法matlab程序,matlab 应求PQ分解法进行潮流计算
时间: 2023-08-11 15:08:43 浏览: 50
以下是使用MATLAB实现PQ分解法进行潮流计算的程序:
```matlab
% 输入电网数据
P = [300;150;200]; % 节点注入有功功率,单位MW
Q = [150;100;80]; % 节点注入无功功率,单位MVAR
Y = [-20+80j, 0-50j, 0+0j; 0-50j, -30+130j, 0-80j; 0+0j, 0-80j, 20-130j]; % 节点导纳矩阵
% 初始化变量
n = length(P);
V = ones(n,1); % 节点电压幅值
theta = zeros(n,1); % 节点电压相角
% 迭代求解
for iter = 1:10
% 计算注入电流
I = conj(P./V + 1j*Q./V);
% 逐个节点求解
for k = 1:n
% 初始化变量
Sk = P(k) + 1j*Q(k); % 节点k注入的复功率
Vk = V(k); % 节点k电压幅值
thetak = theta(k); % 节点k电压相角
% 计算节点k的注入电流
Ik = I(k);
% 计算节点k的电导和电纳
Ykk = real(Y(k,k));
Gk = Ykk;
Yk = imag(Y(k,k));
Bk = Yk;
% 计算节点k的电压幅值和相角
V(k) = abs(Vk*conj(Ik) + 1j*Vk^2*Bk);
theta(k) = angle(conj(Ik)*(Vk*Gk+1j*Vk*Bk));
% 更新节点k注入的复功率
Sk = V(k)*conj(Ik);
P(k) = real(Sk);
Q(k) = imag(Sk);
end
end
% 输出结果
disp(['节点电压幅值:',num2str(V.')]);
disp(['节点电压相角:',num2str(theta.')]);
disp(['节点有功功率:',num2str(P.')]);
disp(['节点无功功率:',num2str(Q.')]);
```
运行该程序,即可进行PQ分解法潮流计算,并输出计算结果。需要注意的是,该程序使用的是复功率形式的潮流方程,因此计算结果中的有功功率和无功功率都是复数形式。