求助！！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行开始：
# g* H2 U& \: d  %% find gens with violated Q constraints
2 _+ P' A8 w9 l8 [# ^  F+ A7 R            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
" {& p+ R& o' T$ F# H, [+ `/ x                if isempty(pv)
                    if verbose
                        if ~isempty(mx)
                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
                        else
2 M$ c' J' E7 N" ~1 n8 p                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
                        end
                    end
                    success = 0;
                    break;
) `: Y- g8 h. R: I- A5 c! l6 q                end
+ @  f" i- ]5 E# x
, \( @% B" H7 z                %% one at a time?
3 P- h2 ~8 L4 r2 a& Q                if qlim == 2    %% fix largest violation, ignore the rest
" k* R' L7 S  z% x                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
1 ]* I" m5 X/ m; j, d7 K; \% L* W                                     gen(mn, QMIN) - gen(mn, QG)]);
) B1 p( W7 J9 Z, U9 b                    if k > length(mx)
/ d! ~$ z; u* u7 m4 @                        mn = mn(k-length(mx));
% ]2 I) k& s6 V; T' E                        mx = [];
7 g; k9 t+ f+ E) M4 T9 [: B& a                    else
                        mx = mx(k);
2 ?2 F* M3 P6 f' r3 Y                        mn = [];
                    end
                end
7 [$ t, }! ^9 }, |
                if verbose && ~isempty(mx)
                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
                end
                if verbose && ~isempty(mn)
5 G. Y1 |- ^0 @, P                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
                end
6 n& z& e" y! e. A+ L               
                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
( o0 k+ K& ^5 z0 T1 g                fixedQg(mn) = gen(mn, QMIN);
  b" H2 O; j( @4 H4 x! H                mx = [mx;mn];
               
                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
0 e4 v! T) ^& d) ~                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,
  m, F7 p& U8 F" B                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
; {6 |4 z* Y' {; H  C9 ~                end
                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
  E' L$ N& _: m7 }5 U: Z                    error('Sorry, MATPOWER cannot enforce Q limits for slack buses in systems with multiple slacks.');
! R% j* L2 K) w# O! x                end
* P7 J' l* k5 y7 v" {                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
楼主看看这个课件。
# U3 _4 w, p! v6 `希望对你又帮助~