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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
这是runpf.m文件中关于无功越界处理的语句，从234行开始：
  %% find gens with violated Q constraints
            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
5 k3 i4 A# P6 |, F( L# K( `% b1 ^* l            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
if ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
) Q- B, G1 d! K" P7 A4 j0 j                    if verbose
# I' K/ d9 b1 x                        if ~isempty(mx)
                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
6 F0 \( {3 ]% J                        else
                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
0 [1 b, n7 Y% n                        end
* C, }0 r0 c4 }/ ?* f6 d. p0 R                    end
8 Q) d9 [9 w! @$ j                    success = 0;
                    break;
6 r& D% d' y. u; t7 _                end
1 h2 D  }8 d. E; ]. b6 V
) `# o$ q0 S) F( p. a; E; m3 {                %% one at a time?
9 T/ D$ M- W) @, X                if qlim == 2    %% fix largest violation, ignore the rest
" f# s2 F* S" }& @- \                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
                                     gen(mn, QMIN) - gen(mn, QG)]);
                    if k > length(mx)
                        mn = mn(k-length(mx));
+ g7 g; v% t* Y( S# z                        mx = [];
                    else
                        mx = mx(k);
                        mn = [];
                    end
                end
4 ^. T5 u" }) ]6 U
4 R: H: {% L! B8 B% {                if verbose && ~isempty(mx)
                    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);
: @7 w2 [2 Q: x! x4 |( O" C                end
6 ~4 }& `$ q! K$ D               
                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
9 i, l+ P8 h2 |4 R: V- \, L1 y6 ]                fixedQg(mn) = gen(mn, QMIN);
  G0 J: p& Q: v; B7 g4 T                mx = [mx;mn];
) D8 g' P, K) {& k4 V3 v, r               
$ K5 z* k, h& n2 n                %% convert to PQ bus
9 _& w; \" h& W1 U                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
/ R: r8 t% @9 Q8 v4 N+ y4 `                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
                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,
8 }' l5 B& }+ m& n                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
5 G. K* z* x; E) Y, |                end
& c6 C" u* v1 i7 h: R. n. a                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
( i5 W! ]* G2 N% G' o% }5 A                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
. A  J% N7 I8 i* _8 b+ Jhttp://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
* K! b6 t9 [) g楼主看看这个课件。
希望对你又帮助~