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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
7 G6 e; N6 ]4 R7 F! T. N这是runpf.m文件中关于无功越界处理的语句，从234行开始：
- Z6 W6 k. @0 d0 t  %% find gens with violated Q constraints
            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
" b* R. r! c4 gif ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
                    if verbose
                        if ~isempty(mx)
                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
                        else
                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
& L6 @8 B0 d' D+ F% R                        end
) x' {3 n  n3 e0 V                    end
                    success = 0;
2 n( V8 d( a) [# E, _  y% C                    break;
3 z  s& o/ i2 L; ]                end
$ [4 s: O) P8 q) v2 Q& t6 l
1 U8 p  F: k( L+ i) A9 W2 T                %% one at a time?
                if qlim == 2    %% fix largest violation, ignore the rest
                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
                                     gen(mn, QMIN) - gen(mn, QG)]);
8 _7 R1 I9 P$ b$ X. I# S- i                    if k > length(mx)
                        mn = mn(k-length(mx));
" h" G% g) `- Z. A- w                        mx = [];
                    else
                        mx = mx(k);
4 m3 I% A8 v* z- j                        mn = [];
7 M0 u6 a4 [* _$ s' a/ I                    end
                end
! d$ S% H; d: H4 b0 N. V: q
                if verbose && ~isempty(mx)
, o' d# u$ Y( q2 V3 u) b4 S! I' u                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
                end
                if verbose && ~isempty(mn)
                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
                end
               
$ s6 a; F+ {7 X9 U" X9 g& }6 t. H                %% save corresponding limit values
5 K" Q. c0 i. r7 G9 Q# l# B# q                fixedQg(mx) = gen(mx, QMAX);
2 L- B+ ~. W0 W, e% H' J/ d                fixedQg(mn) = gen(mn, QMIN);
                mx = [mx;mn];
8 ?9 g6 Q4 a; Q' x               
: v, Q8 ?- L# r& h  Z                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
1 d% N% z* Z% n# g) a: j1 n                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
; d8 \4 z0 o& b                for i = 1:length(mx)            %% (one at a time, since
6 u! G' f/ \$ M  q4 Y& @8 {+ @1 a; `                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
4 T1 `# s6 G. ?: K! j4 o                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
                end
( N& L8 H1 ]. T                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.');
                end
                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
5 }: i+ C/ t: r6 s! @8 ]; Ehttp://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
楼主看看这个课件。
  b2 d. W' y' i# `* I/ a1 s9 }希望对你又帮助~