求助！！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行开始：
7 x7 Z2 t& N9 w8 d  %% find gens with violated Q constraints
" U) U( g- A9 a            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
( I  `4 X7 r- o( E% {            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
4 E/ H' U$ B) Y* J1 G2 fif ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
* b. Q3 ]0 D% z$ l( M; [5 h7 T                    if verbose
0 w+ P0 S/ c! G* a                        if ~isempty(mx)
- y7 l1 ~& e' i1 q# s  K% d                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
; @7 r: w4 x2 L7 r8 D6 H                        else
4 `$ Q! ]) A) B( H' |3 ]6 L                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
                        end
                    end
2 e6 E" n; h- y4 k, ?( t                    success = 0;
, Q  z& W3 {9 f7 \2 |# o! Z                    break;
- r8 @) \1 w7 \" Y) M. f                end

2 N1 \8 D2 j$ f; l4 |                %% one at a time?
) I' \4 w: E8 o8 f8 o                if qlim == 2    %% fix largest violation, ignore the rest
                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
                                     gen(mn, QMIN) - gen(mn, QG)]);
                    if k > length(mx)
                        mn = mn(k-length(mx));
                        mx = [];
8 {6 L5 u% H9 Q9 J                    else
                        mx = mx(k);
6 J1 k" l# X' R1 `( O) B                        mn = [];
                    end
                end

                if verbose && ~isempty(mx)
                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
, |2 {* Q8 W; e. r0 k' F& v( u                end
                if verbose && ~isempty(mn)
2 T! D( ^/ N4 N# O/ H9 a                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
  g5 `/ l7 d: }0 d0 D+ d                end
$ _! u9 E9 T2 G7 f               
2 H, i1 }# G4 c+ j                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
3 Z1 h) S6 t" P& P' {) i                fixedQg(mn) = gen(mn, QMIN);
% s7 s% j, Z4 P5 j                mx = [mx;mn];
               
6 i0 V9 A4 {* _/ r' Z                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
                for i = 1:length(mx)            %% (one at a time, since
                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
6 o; \- Y9 O% C3 s                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
3 d( a5 Q2 ~4 G2 H; X' Y( W                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
0 l- m( T" F/ y. ^0 W                end
                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
/ @3 y( B* g1 ?                    error('Sorry, MATPOWER cannot enforce Q limits for slack buses in systems with multiple slacks.');
- l/ a; H+ y: G1 B                end
                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
' N& I! g) X5 q0 E4 Q, Dhttp://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
8 R$ B0 q/ _7 D4 j; S$ H4 d楼主看看这个课件。
! H/ I" A& J6 H+ U希望对你又帮助~