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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
+ I( ]1 A$ Z+ u! z) w  s; [这是runpf.m文件中关于无功越界处理的语句，从234行开始：
2 g) T, q1 ~, \* i' {: j. c  %% find gens with violated Q constraints
7 e6 i+ T% s1 g6 w/ [            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
. q3 c4 L9 Q1 ~* h( G: L            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
, Y/ t! H# T  s5 \if ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
* s/ m/ V4 w: Y/ [2 J                if isempty(pv)
                    if verbose
3 n, K' p0 ^0 u! [. r                        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);
                        end
                    end
                    success = 0;
                    break;
                end

7 K% f! L+ k: U, R2 A' v                %% one at a time?
$ B% H7 V: N: ?% N; w                if qlim == 2    %% fix largest violation, ignore the rest
$ y  ~# V% C- M% }+ ~                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
2 m4 a3 A4 C; _1 l' X                                     gen(mn, QMIN) - gen(mn, QG)]);
                    if k > length(mx)
7 Z9 M" }9 H( G) v4 J. h, G                        mn = mn(k-length(mx));
, \) y7 D' i* s6 @$ \+ O$ h                        mx = [];
9 |" }& i. f1 o5 c( \+ g& U                    else
                        mx = mx(k);
                        mn = [];
. l3 l9 z, q0 b2 z- y) Z                    end
& l- E$ ?8 N1 r7 R                end

! p6 O& L" M. ~                if verbose && ~isempty(mx)
                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
                end
                if verbose && ~isempty(mn)
- z( n5 d" P" x* A8 K1 @; s0 r0 K0 B                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
                end
( f/ R( w3 e; ?" u               
                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
7 q3 V" Y4 n! O8 @- e                fixedQg(mn) = gen(mn, QMIN);
: l; i8 u  d3 [  Z* K7 M                mx = [mx;mn];
               
6 O% l, P7 p% W* k4 M6 V                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
( @+ y' J* _8 S! ]! c0 w                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
                for i = 1:length(mx)            %% (one at a time, since
% `2 |4 X- g9 P                    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]);
& \$ f) k! h* x                end
. T# ~) m! _  B0 B# n6 N                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
7 l, F4 Y8 R9 ?' T1 h- w4 n- `' s4 H                    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 U" e; r& V# t+ X! Y' @" P+ ]http://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
楼主看看这个课件。
希望对你又帮助~