求助！！matpower算潮流PV节点转换问题

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
8 G1 R& f, a+ u9 u. G( a! v我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
, K( A9 h# j. w& z. J( p这是runpf.m文件中关于无功越界处理的语句，从234行开始：
  %% find gens with violated Q constraints
8 R3 m5 F6 p5 Z2 T0 D( K            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
# i! b/ J2 u/ `0 [1 |  z" B            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
% x% e+ t) n4 u* n  j4 Vif ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
( s& S' g5 F/ t! x- q& S9 F/ `                if isempty(pv)
                    if verbose
                        if ~isempty(mx)
, q6 h% U) y8 L                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
                        else
# ?$ C0 \- E. {1 q; W: D; b                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
7 y3 Z" K5 `) S! ?% A* m% t                        end
                    end
                    success = 0;
                    break;
                end

* ^8 u; l4 h2 Q  B8 r                %% one at a time?
! O5 j: T4 D' i8 X7 |& o. \6 c' S                if qlim == 2    %% fix largest violation, ignore the rest
                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
                                     gen(mn, QMIN) - gen(mn, QG)]);
                    if k > length(mx)
( [; m( X, W/ j3 Z/ }4 U- T5 L) A                        mn = mn(k-length(mx));
4 G$ H' x# D  ?7 w  j! N                        mx = [];
# S4 e- S9 ], S% l2 O                    else
  m6 f7 V1 M8 _( O/ g& Y                        mx = mx(k);
                        mn = [];
  L' s! G2 @) E6 B2 E1 p! Z                    end
                end

                if verbose && ~isempty(mx)
7 ~: L) E( A! y. a; T3 E. J                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
                end
5 ~2 z: \3 L) L9 F+ e" q! f                if verbose && ~isempty(mn)
                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
' q! e+ g( l# w" h3 i, o9 ?' x                end
" z2 [0 j8 h8 @4 {5 e( w; f               
8 s+ n7 F2 ^  `" F! j. T                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
                fixedQg(mn) = gen(mn, QMIN);
, f; c9 K9 [# H4 u+ F; ]7 V                mx = [mx;mn];
" z) Q/ \0 l3 m% v               
                %% convert to PQ bus
( i: e8 E7 o; U' ~4 q& Y0 z: F                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
! l5 `  L- `* n. d# E                for i = 1:length(mx)            %% (one at a time, since
                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
                end
+ g9 c4 G: K5 P$ @                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
                    error('Sorry, MATPOWER cannot enforce Q limits for slack buses in systems with multiple slacks.');
+ j7 Z) P3 k3 h0 z5 H                end
; U6 T' s3 x: E+ [: c' d                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
% Z+ T) G; K2 J" Phttp://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
楼主看看这个课件。
希望对你又帮助~