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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
) E# J8 E6 N, ^# t我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
这是runpf.m文件中关于无功越界处理的语句，从234行开始：
0 Y+ c* @2 h; D/ q, `+ O$ b  %% find gens with violated Q constraints
9 ~. `2 S- E/ ~5 S            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
8 F4 t9 x( z% G+ V            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
9 i* y- ]% g1 t: bif ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
                    if verbose
9 t+ ^( c: g9 g( u, H7 m                        if ~isempty(mx)
                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
                        else
: a4 e9 ]- J3 L, d: d. L9 i' m2 O% S! m                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
                        end
# B! Q( A* u8 e, e1 T. R$ ]- k                    end
                    success = 0;
                    break;
/ J2 U7 u7 ]1 j/ A% d2 \3 P                end

                %% 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)]);
) [: C( c$ o2 L7 w" S2 K8 z5 j% L7 b                    if k > length(mx)
0 H1 M6 q! G& v  O; ?. I                        mn = mn(k-length(mx));
                        mx = [];
) O/ j# p4 R9 ?) E                    else
. g  `5 ]) D" T                        mx = mx(k);
& s' H  l/ ~( R" l; T                        mn = [];
                    end
                end
; N. v$ \. z4 ~3 ?. E0 g
                if verbose && ~isempty(mx)
% d6 u* K+ \- A8 T% H0 H                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
0 |2 Y3 o5 z, S9 q                end
                if verbose && ~isempty(mn)
                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
                end
( u# D- l: Z% V               
                %% save corresponding limit values
) Z4 T* Y3 @* t& i6 Y- V                fixedQg(mx) = gen(mx, QMAX);
! H6 t0 i; i( k& ]                fixedQg(mn) = gen(mn, QMIN);
                mx = [mx;mn];
) ^; C( y1 q7 B) |( X3 R/ o               
                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
5 n- v% e5 X) D& \' P                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
                for i = 1:length(mx)            %% (one at a time, since
2 u; m6 _% V. R' X& [                    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]);
                end
/ _; f3 C- j4 k" {/ n& a8 K                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
) O$ |$ Y: S# L                    error('Sorry, MATPOWER cannot enforce Q limits for slack buses in systems with multiple slacks.');
" x% x5 y8 L4 X* K3 l                end
                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
7 x3 ^8 }1 a+ e6 x  ~http://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
楼主看看这个课件。
希望对你又帮助~