求助~~PQ法B1矩阵的形成

zhyerica · 发表于 2012-5-6 11:04:01

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

各位大牛，求助啊~~
我毕设的题目是交直流混联系统的潮流计算，自己编了交替迭代法的的PQ分解法程序，计算结果与书上算例结果误差最大有12%。
于是测试了纯交流部分的程序，也就是一般的PQ法，对照着现代电力系统分析P86上的例子，发现自己程序B1矩阵算到的结果跟书上的有差，致使第一次迭代时相角的修正量跟标准值就差了4%左右。但是查到我程序生成的B1矩阵跟自己手算的结果是一样的啊……

附上相关部分代码 Ps 我得到的节点导纳矩阵和B2矩阵都是和书上的一样的，真心搞不懂了，求解答啊~~~

function [Bp, Bpp] = makeB(Node,Branch)

%% define named indices into node, branch matrices
[PQ, PV, REF, ...
TRANS, CAPAC, LINE,...
NODE_I, BUS_TYPE,PG,QG, PL, QL, V, THETA, ...
BRANCHNUM, NODE1, NODE2, R, X, BORK, BRANCH_TYPE, STATUS, P1, Q1, P2, Q2] = index;

%% form index lists for slack, PV, and PQ buses
ref = find(Node(:, 2) == 3); %% reference bus index
pv  = find(Node(:, BUS_TYPE) == PV); %% PV bus indices
pq  = find(Node(:, BUS_TYPE) == PQ); %% PQ bus indices

%*************************************************************************%
%                            form Bp
%*************************************************************************%

%找出非接地支路
temp = find( Branch(:,NODE2)~=0 );
NongroundingBranch = Branch(temp,:);

%忽略线路充电电容
NongroundingBranch(:,BORK) = 0;

%不考虑变压器非标准变比
temp = find( NongroundingBranch(:,BRANCH_TYPE)==TRANS );
NongroundingBranch(temp,BORK) = 1;
Y = makeY(Node,NongroundingBranch);
Bp = imag( Y );
Bp = Bp([pq; pv], [pq; pv])  ;  % reduce Bp matrices

%*************************************************************************%
%                            form Bpp
%*************************************************************************%
Branch(:,4)=0;%忽略输电线路电阻
Y = makeY(Node,Branch);
Bpp = imag( Y );
Bpp = Bpp(pq, pq); % reduce Bp matrices

return;

zhyerica · 发表于 2012-5-6 14:37:12

自己顶~~~
拿那个5节点的例子，算出来的B1矩阵如下：
B1 =

(1,1)    -6.5417
(2,1)    3.9002
(3,1)    2.6415
(1,2)    3.9002
(2,2)    -73.6789
(3,2)    3.1120
(4,2)    66.6667
(1,3)    2.6415
(2,3)    3.1120
(3,3)    -39.0869
(2,4)    66.6667
(4,4)    -66.6667

wlm_28 · 发表于 2012-5-6 20:37:46

网上有好像有这个程序

zhyerica · 发表于 2012-5-7 10:04:07

回复 3# wlm_28

网上找到的PQ法B1、B2矩阵都是手动输入的，没看到程序实现的。主要是拿5节点的算出来跟书上也差挺多，B2是一样的，B1书上只给了因子表，但是从迭代第一步的结果上看应该就是B1的问题……

账号		自动登录	找回密码
密码			立即加入

[讨论] 求助~~PQ法B1矩阵的形成

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子

浏览过的版块