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

edwardsy

现在用matpower工具算潮流，发现一个问题。执行runpf('case××')算得某一个PV节点的无功出力是24.3，然后我修改generator data数据中的Qmax和Qmin，把Qmax改成了24，然后重新用runpf函数来算，可是得到最后的结果，这个PV节点的无功出力还是24.3。
4 j5 F0 h/ P5 F2 Z5 n! ?我查看runpf.m文件，发现里面有关于无功限制的语句，就是当无功越界时，将PV节点转为PQ节点。那么为什么执行不了呢？求教各位朋友~~谢谢~~
这是runpf.m文件中关于无功越界处理的语句，从234行开始：
  %% find gens with violated Q constraints
            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
                if isempty(pv)
                    if verbose
3 a( B! @- D. w9 a$ W9 v                        if ~isempty(mx)
                            fprintf('Gen %d (only one left) exceeds upper Q limit : INFEASIBLE PROBLEM\n', mx);
+ T0 u: w2 ^3 [1 O                        else
                            fprintf('Gen %d (only one left) exceeds lower Q limit : INFEASIBLE PROBLEM\n', mn);
7 G) s2 m- |. K                        end
                    end
                    success = 0;
( H% @. k/ I& z/ K" C0 c! C# y" S                    break;
                end

2 K+ b# o3 e) M! i- s" o4 u                %% one at a time?
6 ?' l  O# |. Q: G2 p% f                if qlim == 2    %% fix largest violation, ignore the rest
4 ^% F5 j3 u" F6 {, C                    [junk, k] = max([gen(mx, QG) - gen(mx, QMAX);
  ]) E+ Y: z$ y                                     gen(mn, QMIN) - gen(mn, QG)]);
$ m0 L) S2 n5 j( b2 a5 ~6 n" P0 L                    if k > length(mx)
( ?2 |4 G. T: A8 W) X( D                        mn = mn(k-length(mx));
                        mx = [];
                    else
( T+ p, `( x# I/ }' e. q& L                        mx = mx(k);
                        mn = [];
                    end
9 p, Q4 i! ^' L9 p5 t                end
+ h  U& d; S! }# [$ Q
                if verbose && ~isempty(mx)
8 o: n8 {7 `; k- F                    fprintf('Gen %d at upper Q limit, converting to PQ bus\n', mx);
6 v1 g6 b) I6 \$ ]6 V1 @% ~& ^* ?. n                end
                if verbose && ~isempty(mn)
  Y2 B/ F5 p0 Y) f: g                    fprintf('Gen %d at lower Q limit, converting to PQ bus\n', mn);
( V/ y) @) M) m: n                end
: e9 G+ A7 m6 r2 w4 E$ s9 f               
7 n* e7 V# O" v4 e8 o  E/ s                %% save corresponding limit values
                fixedQg(mx) = gen(mx, QMAX);
                fixedQg(mn) = gen(mn, QMIN);
$ F1 p  R2 L+ I9 z3 `                mx = [mx;mn];
/ n, ~' y* {( Q3 X3 j. e6 V) J               
                %% convert to PQ bus
" ~, ^- w# z+ P' g                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
4 V* }4 p. \& @4 c; L                    bi = gen(mx(i), GEN_BUS);   %%  they may be at same bus)
                    bus(bi, [PD,QD]) = ...      %% adjust load accordingly,
' f( C6 ^2 P3 d5 K                        bus(bi, [PD,QD]) - gen(mx(i), [PG,QG]);
! K- e1 M4 R3 ?+ y5 r1 X                end
2 N. c: }6 a- f! J                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.');
  e! k: _, e9 a1 O% W                end
                bus(gen(mx, GEN_BUS), BUS_TYPE) = PQ;   %% & set bus type to PQ

快乐哈哈

要想无功限制，必须设置mpoption中的qlim值，当qlim为1时才能无功限制，默认值为0，不进行无功限制。楼主可以看看mpoption.m文件。
1 y. @7 S6 [) ?+ D% p6 P# ghttp://wenku.baidu.com/view/1df27b718e9951e79b89278d.html
6 M+ B3 h3 ]) C楼主看看这个课件。
希望对你又帮助~