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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
( `+ |! j0 P- z这是runpf.m文件中关于无功越界处理的语句，从234行开始：
; _7 \* [2 V( Z0 G+ @  %% find gens with violated Q constraints
            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
1 x# j) ]( Y' _$ A/ G            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
if ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
5 R% F, Z; B6 J. g0 ~5 z                if isempty(pv)
8 V9 g5 W% R* W! f" O% U+ B5 H                    if verbose
. g9 r$ m5 T8 P7 z( ^                        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);
$ ]  R1 b4 z8 I- \' ~                        end
& ^- {6 X' u9 l! z8 f                    end
                    success = 0;
                    break;
                end
7 D0 Y! t, G8 r) S. k# w
7 B6 |2 @! n, L; b1 D. \& j                %% one at a time?
. r( C2 n* O1 N1 ~# d$ d8 c) {. ]                if qlim == 2    %% fix largest violation, ignore the rest
6 g9 i3 D6 X, n: d; j) `                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
) C/ S7 Z5 v% B5 u                                     gen(mn, QMIN) - gen(mn, QG)]);
                    if k > length(mx)
                        mn = mn(k-length(mx));
                        mx = [];
9 E: W/ }7 k5 f, x/ W; I( T                    else
                        mx = mx(k);
* g* G8 L+ b( A  l1 F, Z                        mn = [];
                    end
8 q% L" B5 L4 V                end

: Q: V$ X1 k0 J) K                if verbose && ~isempty(mx)
0 d% d. @; c" T/ g, 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);
2 g0 i; @) Z& u0 i9 S3 H3 T                end
- C- x' a( }' [+ g! \% f               
                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
                fixedQg(mn) = gen(mn, QMIN);
                mx = [mx;mn];
$ f0 y( [1 t) _; {7 `6 I2 Q               
) R! `# c/ \# m0 o) }! B3 t, S3 K                %% convert to PQ bus
( M9 C& i% R( I0 L' n; ~' o                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
, ~  h8 w  ]3 a5 A, b                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
2 v" v7 `6 }' H                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]);
/ c* {9 A1 X  o                end
0 w# B0 G0 _( g5 Q/ p' g                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
! M+ t1 \$ g! I# ~& [. |( R                    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文件。
http://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
楼主看看这个课件。
希望对你又帮助~