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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
! c3 @2 w0 T2 K8 a; K这是runpf.m文件中关于无功越界处理的语句，从234行开始：
  %% find gens with violated Q constraints
            mx = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) > gen(:, QMAX) );
' X0 d' r' j3 Z! B            mn = find( gen(:, GEN_STATUS) > 0 & gen(:, QG) < gen(:, QMIN) );
% y5 a0 W" Y, zif ~isempty(mx) || ~isempty(mn)  %% we have some Q limit violations
                if isempty(pv)
                    if verbose
                        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
: l' t2 m& S" M# q6 L                    end
                    success = 0;
                    break;
; O! f- @9 B0 a, W                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)]);
: T5 w* Z7 e* J                    if k > length(mx)
                        mn = mn(k-length(mx));
                        mx = [];
+ r, o( V! C& f! \) S, c" h7 _                    else
5 w4 m6 j+ d5 s9 c  `$ x% z2 a% \                        mx = mx(k);
                        mn = [];
+ X; Y+ D& ?$ a! G0 B' [6 b$ v# T                    end
3 ~) x+ K# e/ h+ X                end
" W+ ~* U# E+ |+ S- z2 {
  B; q# x7 Y4 d/ b/ s' y                if verbose && ~isempty(mx)
  O) @. M1 |7 n* x, W5 E                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
                end
9 g& S2 L3 M+ C. a7 l4 A" o                if verbose && ~isempty(mn)
! @  c$ F' e! C4 {. V- ~) A/ S                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
                end
               
                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
6 U: B1 ^6 S" N# A                fixedQg(mn) = gen(mn, QMIN);
                mx = [mx;mn];
1 G. O& e& G* M: U% i               
5 Z+ J4 n4 s" m) w  ?$ h8 ~                %% convert to PQ bus
                gen(mx, QG) = fixedQg(mx);      %% set Qg to binding limit
! f" Y& v; J; k; D2 S! l9 Y2 p/ I                gen(mx, GEN_STATUS) = 0;        %% temporarily turn off gen,
( G3 d7 v' v$ j/ `0 u                for i = 1:length(mx)            %% (one at a time, since
                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
8 U5 X5 z) ^3 q2 f9 o: Q0 c                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
8 y, z- }5 l0 e' b                end
                if length(ref) > 1 && any(bus(gen(mx, GEN_BUS), BUS_TYPE) == REF)
                    error('Sorry, MATPOWER cannot enforce Q limits for slack buses in systems with multiple slacks.');
; ^0 F8 c3 x7 j! ~: m: }                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
楼主看看这个课件。
* [& e0 D" D' _) p  I0 B希望对你又帮助~