求助！！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行开始：
( p* ?. {$ X( P  %% 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) );
if ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
1 `' s& y9 R, r: [. F* j1 T5 `5 ~8 R                    if verbose
# d0 g0 t8 L! C. E                        if ~isempty(mx)
' s( D% T3 e7 s, R+ q                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
                        else
8 V9 T* X  r- E- n7 d5 y! J                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
                        end
# _7 |* r& H2 v, a2 a3 n2 C                    end
                    success = 0;
                    break;
3 r2 y3 o4 h+ G                end

$ h8 a8 l% Y* N/ q% b4 Z: B+ V                %% 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)]);
9 r) n' p  t2 P( Y$ v0 l: v6 H                    if k > length(mx)
                        mn = mn(k-length(mx));
, M8 m: a8 C8 w6 U+ Q                        mx = [];
                    else
                        mx = mx(k);
                        mn = [];
                    end
                end

                if verbose && ~isempty(mx)
# Y) w; e  i5 B/ d4 `8 V- c                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
) G5 P3 u9 b4 S6 ^3 W1 m# P: y$ c7 }                end
$ g& r* u' j0 d% L- E                if verbose && ~isempty(mn)
                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
0 P8 u0 Y0 i( O& l8 d4 l- `  I                end
3 e+ m: }6 l+ J) x, s7 A% Y               
                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
; k4 k) t7 b8 D/ E1 h                fixedQg(mn) = gen(mn, QMIN);
. T. ]! m) T: K( s# i                mx = [mx;mn];
% O) J- V7 A1 `( B0 _               
                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
8 f1 |# N7 s$ U) _! M- C$ _                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
                for i = 1:length(mx)            %% (one at a time, since
6 J& L4 A. t, Y                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
) R2 `) A+ @* |. r" V" b, l, h* c9 P                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
. g/ w& C$ S- ]: X( n% T                end
                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文件。
( P2 j  o) a( E1 t; b2 {http://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
6 p0 }5 j) D+ T! C# d楼主看看这个课件。
+ M- L' P* Q# P% e希望对你又帮助~