matpower潮流计算发电机功率及节点负荷修改问题

陌上花开为君故

matpower30节点算例中，我需要修改发电机的出力，使其在一定范围内波动，比如说50%~150%Pg之间, 0.05步长，自动循环运行程序，怎么修改runpf.m或者case30.m中的程序？在调用runpf('case30')的主程序中修改有效吗？
* h8 Z) c& ?% O! ?30节点算例中支路的线路长度数据有吗，哪里可以找到~紧急求助啊，谢谢大家

凉风

yicqu 发表于 2012-3-7 20:22
8 O1 J' u. N  j可以在runpf.m中修改，如果仅修改发电机出力而不修改负荷大小，那么不平衡的功率完全由平衡节点承担？如果 ...

具体怎么加循环

sunshine0703

matlab痴 发表于 2012-5-23 22:39
回复 8# heavenying

: g2 y/ c. q$ T请问33节点建模时，为什么支路阻抗有六组数据，但是输入时只有一组呢？麻烦你了，谢谢！！！

JIUSHIGUANG1

同问 MATPOWER除了算潮流还有什么具体的用法吗

yicqu

可以在runpf.m中修改，如果仅修改发电机出力而不修改负荷大小，那么不平衡的功率完全由平衡节点承担？如果同时修改发电机和负荷功率的话，修改比较简单
/ N( K# O6 R0 u: [- g# srunpf函数默认的是采用牛顿-拉夫逊法进行潮流计算，可以对下面这句话增加循环（可以ctrl+f找到这句话）：
5 b( V+ U" S7 K$ G[V, success, iterations] = newtonpf(Ybus, Sbus, V0, ref, pv, pq, mpopt);
0 K4 e/ k+ K" {3 f每次调用这句话之前需要根据循环控制变量修改Sbus
4 m* h5 d, Z  d& h  w如果确实需要仅修改发电机出力而不修改负荷大小，可以在[baseMVA, bus, gen, branch] = loadcase(casename);即读入数据以后加入循环。
) {2 U9 m/ k& S# J4 j6 M, S8 W% X30节点系统更多的信息可以参考Alsac, O. & Stott, B., "Optimal Load Flow with Steady State Security",IEEE Transactions on Power Apparatus and Systems, Vol. PAS 93, No. 3,1974, pp. 745-751.

542764435

回复 2# yicqu
7 `  s% L2 J; j. {- a
1 _$ n/ g3 \3 l) ]* ILz 我最近看http://www.ee.washington.edu/research/pstca/ 里面的IEEE common data format 和Data in Common Data Format (file ieeeXcdf.txt)  但是这两个东西好像对不上，貌似ieeexcdf.txt 是为结果的？我用matpower运行IEEE的算例时候 好像结果和它的部分是一样，希望能向你请教下！

yicqu

回复 3# 542764435
3 M  C+ A5 O( x
) Y  k+ v# S* ?# _8 d' k8 P是一致的，比如节点数据的data cards 部分，以前面几列为例说明如下
Columns  1- 4   Bus number (I) *节点编号
* _$ v! `7 {/ ~! z" i! H9 cColumns  7-17   Name (A) (left justify) *节点名称
Columns 19-20   Load flow area number (I) Don't use zero! *潮流区域编号
0 i8 B: |" y2 r8 N6 EColumns 21-23   Loss zone number (I)损耗分区编号
4 h& `4 U/ [/ ]1 RColumns 25-26   Type (I) *节点类型
                 0 - Unregulated (load, PQ)
; y) m* \, u+ \1 q7 t% ~2 o                 1 - Hold MVAR generation within voltage limits, (PQ)
                 2 - Hold voltage within VAR limits (gen, PV)
: n' `4 a. ?7 p# ]9 a                 3 - Hold voltage and angle (swing, V-Theta) (must always
. T* q6 O/ J# ^- s' j                      have one)
: X# j$ T! u0 t8 r# X/ [Columns 28-33   Final voltage, p.u. (F) *最终电压幅值
Columns 34-40   Final angle, degrees (F) *最终电压相角
) U! g* T1 j# Y1 k, \% r3 f* n' i5 d1 y. . .
0 U6 V" v# |; s: [' e7 N- F一般程序都可以不理会Columns多少到多少，根据空格区分字段就可以了。

yicqu

回复 2# yicqu
" g( C' y% v+ O/ M

    另外一种情况，通过增加负荷，描述电网电压变化的情况，类似于电压稳定中研究的问题，但一般来说牛顿法无法达到电压崩溃点，无法画出完整地PV曲线。

kamui17@163.com

回复 5# yicqu
+ d; |4 \+ w" G+ i6 Z% E

    你好，我想问一下怎样用编程语句修改case14或30原有的参数，比如问我要将全系统的负荷改变成弹性负荷，该怎样编程？

matlab痴

回复 4# yicqu

$ _. u+ Y) v' |1 n* a( T
    请教一下，能否用matpower来计算ieee33配电网系统呢？可以的话，编写case33.m文件里的发电机参数怎么设置呢？

heavenying

回复 7# matlab痴
! g& Z  p- h# p2 k- x, C
& A# t7 y* g' `- W- `
   我最近也想用matpower编配电网潮流。昨天试了一下，跑了个33节点的，收敛到是收敛了，但结果和paper里的不太相同

matlab痴

回复 8# heavenying

* o; `$ o5 [9 x3 k7 x7 V: w/ Q
% l/ B2 N3 [. x8 R* x3 I, U9 r    我这几天又重新修改了配网33节点输入的数据，结果收敛，感觉计算的还可以，但是之前没有计算出来的正确的结果，现在的结果也不知道到底对不对！

大号

请教：
) U9 Z5 x( t8 ~7 B! ]5 ?我用下面的语句修改基值：
3 ?' K% P6 X9 Ympc=loadcase('case30');
! h% O" C, S( Q1 R6 i+ Ympc.baseMVA=157;
3 ?4 U. A$ A' f但是各台发电机所在母线基值仍然是100MVA。要怎么处理啊。
; x2 W2 {  i8 C4 H* L( W: ]另外，我主要是想模拟机组增减有功无功的发电机机端的波动情况。该如何做？好像直接幅值有功无功有点问题。