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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
, C% N5 G, `$ ~2 {. [我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
( V/ D! ]/ ~9 e2 a) T这是runpf.m文件中关于无功越界处理的语句，从234行开始：
  %% find gens with violated Q constraints
% G1 b0 D8 G, _3 C: y            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
* J: M+ Z+ s9 U8 n            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
, N; X1 r8 U- Y6 a& xif ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
                    if verbose
, b, |2 V+ b$ o( @                        if ~isempty(mx)
                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
                        else
6 K9 M3 Y- A* J: R  c                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
8 N/ ^5 K' l+ w: e                        end
2 Z' s& J! b# `( G                    end
                    success = 0;
2 k3 [% S1 u& `                    break;
                end

: W$ W; X7 H- v0 G& N6 R                %% one at a time?
                if qlim == 2    %% fix largest violation, ignore the rest
                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
4 ]: G* k+ J2 m/ V                                     gen(mn, QMIN) - gen(mn, QG)]);
' @0 h% j& m; j5 Q9 Z' c: o1 u                    if k > length(mx)
/ s: O8 ?, J5 V' U' r+ s                        mn = mn(k-length(mx));
                        mx = [];
                    else
                        mx = mx(k);
- r4 x# X! y; D3 V1 o8 E2 i                        mn = [];
: R$ W) f: o* `& }! ^+ x                    end
3 u3 y" F3 t% v7 F. t                end

5 {8 X0 H' p/ v6 m9 h/ G9 y+ Q                if verbose && ~isempty(mx)
1 O% W5 E1 G2 K3 u4 R. p                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
                end
4 f2 ~  {  {2 L- d1 E/ @6 x                if verbose && ~isempty(mn)
                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
                end
               
* `  j2 j6 _1 }; O( V                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
                fixedQg(mn) = gen(mn, QMIN);
* E- R3 H/ g# K# X7 \                mx = [mx;mn];
               
) O" X8 w* d7 Q2 [2 L( [# }: w' d                %% convert to PQ bus
2 U! M6 ~% S8 _+ ?3 d                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
3 F  K1 T; E  z, d- h4 H                for i = 1:length(mx)            %% (one at a time, since
+ D7 a* ]. Z/ N; v3 C                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
& {" |  w+ ^; G9 k, |( D( F: w) F                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
0 D  C2 {  C% N- X9 T1 J                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
! A( ]' g- _- o6 S: d: U                end
% |7 }' _/ T5 M, I6 j% T2 Q                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
5 O; D6 r; ]0 ^$ T) \* L4 l  c                    error('Sorry, MATPOWER cannot enforce Q limits for slack buses in systems with multiple slacks.');
                end
* V( H1 T5 K: C/ }( j& c. l: V8 }                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
8 Y. ^2 g% R9 p; X; dhttp://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
楼主看看这个课件。
5 f& ^; Y3 ~$ ]! J2 b& Z/ y希望对你又帮助~