马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!
您需要 登录 才可以下载或查看,没有账号?立即加入
×
附件前半部分和本文相同,后半部分为部分代码,在原网站上能找得到,只是综合了一下。
+ ?- R5 B# C5 a+ `
, I$ j. ]6 E2 v9 M$ Q: W1 q
1 V% X2 [% O" Z& | 在PSSE中使用Python脚本:简要介绍 translated by Jackain
1 s8 q0 U0 H* a3 Y- E 9 d' n9 C% v7 x U6 C; Q3 [
PSSE30 提供了 Python 的 API 接口,使在 PSSE中使用Python脚本语言 成为了可能,这不仅使 PSSE 批量运行,同时可以根据你的 Python 代码,使处理功能更强大,更灵活。 1 q8 g1 m5 i& J5 I
我没有找到很好的关于在 PSSE 中使用 Pyhton 的资料,我把我能找到的东西写在这里供初学者入门。 % `# [- v# b$ B/ W
前面的一些资料能在使用说明书 USERS.pdf 的第 12 章找到,其他资料是通过 API.pdf 和 IPLAN.pdf 综合而来的,可以从中找到。 . S( l! k0 `9 N, s& L1 H
我发现使用 Adobe 的搜索功能非常有用,你可以根据一个关键词同时搜索在一个文件夹下面的所有的 pdf 文件,这样寻找某个函数的 IPLAN 和 API 的描述遍得很简单。 IPLAN 的描述一般在查找详细描述“某些传递数据的参数”的解释时非常有用。 IPLAN 函数的名称同 API 文档中 Python 的那些函数相同,只是后者是大写字母而已。 2 z. q" ]6 J3 U' q9 U9 i$ O
我使用的是 30 版本的 PSSE ,和 2.3 版的 Pyhton ,这篇文章的前提是你了解一些基本的 Python , Mack Lutz 和 David Ascher 写的《 Learning Pyhon 》非常不错。
0 o" H% m4 K& w. y# K Python + K& z0 `* |# \3 K, w1 h
Python.org homepage
- n# g2 V$ d+ R5 Z& B- C PSSE
9 C g, q, G+ _ ?( X$ i Power System Simulator for Engineering.
! M8 |6 A3 `# X4 m, D
( f* i$ @) Z9 q# j- V! B) b% ^ 在 PSSE 中使用 Python 简介 $ W8 u4 F3 ~) o- s" _& V
PSSE 说明文档说:根据你在 PSSE 中的操作可以生成 Python 代码。这花了我一段时间才理解,因为我是 PSSE 和 Python 的初学者。根据以下步骤,可以生成你的第一个 python 代码,用来进行潮流计算分析。前提是你有一个包含有网络信息的 .raw 文件,格式符合 PSSE30 要求。
' X8 w4 I6 H5 K. u% [$ s& m; Y6 f 1、
: S. V+ K- S% w' h# ]+ k 打开 PSSE 。 0 P5 p. |1 p& K6 e* \9 R& d
2、
1 H6 g5 d3 L+ l G* Q: q2 j( m8 e- \ 选择“ I/O Control ”中的“ Start Recording ”菜单,输入一个 PSSE 写入的文件名。 PSSE 会自动地添加 .py 后缀。 4 k' I6 i4 z" s) E+ F; f
3、
0 K2 ]# S& g" v+ z$ s% g 使用“ File ”菜单中的“ Open ”按钮,读入网络信息。将文件类型设定为 .raw 告诉 PSSE 你在给他的是 raw 文件(不知道是原文傻还是我翻译的傻 :’()
# V% d; u* ^: t2 Y) U( ~# s9 E 4、/ ^- U* e4 {" }) L% y2 J
选择“ Power Flow ”中“ Solution ”,选择“ Solve ”,运行潮流计算。选择所需要的设置模式,点击“ Solve ”按钮。 1 ~$ \, J3 C0 I' M
5、
d2 ?' {+ [1 g: k 选择“ File ”菜单中的“ Save or Show ”保存运算结果。选择“ Power Flow Raw Data ”选项卡,点击“ Data File ”单选按钮,并在旁边的文本框中输入文件名。点击 OK 按钮。
5 Y# n+ u* T" q3 S3 \2 K 使用你最喜欢的 python 浏览器查看生成的 python 代码。
0 j* y# d" W! j7 `" X" u& ^1 m0 r # File:"C:\Documents and Settings\jbo\MySensibleFolderName\pssetest.py", generated on TUE, FEB 08 2005 14:11 psspy.read(0,r"""C:\Documents and Settings\jbo\MySensibleFolderName\testCase_30EV30test.raw""") 7 [3 u7 s& V, K( a
psspy.fdns([0,0,0,1,1,0,99,0])
5 ?$ a% W6 W6 p9 _& `' N: Q* P psspy.rawd(0,1,[1,1,1,0,0,0],0,r"""C:\Documents and Settings\jbo\MySensibleFolderName\SOLVEDtestCase_30EV30test.raw""")
注意,所有的方法都是通过一个叫做 psspy 的对象调用的。这个对象包括所有可用的方法和对象。他们看起来只不过是通过一个列表来组织的。我还没有写任何接口对象来组织和调用 psspy 方法时,就是这么想的。(瞎翻译的) " P% k7 R) ^/ S: k( U
注意声称路径字符串的方法,这不是 python 提倡的格式,下面的写法更好: " Y- F0 Q& e1 [, N. j
"C:\\Documents and Settings\\jbo\\MySensibleFolderName\\testCase_30EV30test.raw".
5 Z" i# X8 q5 \6 b. ]- _% S 怎样来运行 python 代码 ; r6 \. P) T% _9 J# W/ i
我试过三种方式
7 ]+ t. r0 c7 \& X0 V' U 1、) G0 c/ W4 Y- ?! |3 G7 G1 m1 K
在 PSSE 中打开命令输入窗口( command line input, CLI )任务。
4 W, r( g( V) |# v; I' d 2、0 ?, x+ O% \) ` ]" C1 L
在 PSSE 中选择“ Run Automation File ”按钮运行一个 python 文件。
' d3 f. E' W( ^8 {* x 3、
[ u. x& R# B9 [ 在 PSSE30 中 python 文件作为一个命令行变量。你可以指定图形界面只有在完成之后才出现。这可以为机器节省内存。在这种情况下,你可以指定仿真的母线数量。可以参照 PSSE Users Manual 中的 PSSE Overview 。 - C& L t9 h5 Q% P- x
在 PSSE 中打开 CLI . e }" E' G- h2 r: j& Y
1、
. u D# ?* u3 q6 _ 在“ I/O Control ”菜单中选择“ Begin command line input session ”。 6 o* I8 v# x! D) a. [/ ]4 M
2、+ n, c0 E" D' f- b
使用“ Command Language ”选择框将语言更改为 python 。 7 D9 G. B( C9 s6 n
3、
8 Z4 M! U* {) y: C 选择“ RunAuto ”按钮并选择一个 python 文件, PSSE 会运行那个文件。 ) ^1 b0 K; z3 s- B4 d2 H1 H+ o
注意 python 命令可以直接在控制台输入。控制台动作同编译 python 的交互式平台不同,它不自动显示调用方法所返回的结果。要想显示的话必须使用 python 的 print 命令,例如: print(psspy.aMethod(myParameters)) 。
: y+ R) c* w4 j 一个重要的点是所有列出的方法都是 psspy 方法。所以使用要使用 totbus 方法返回母线数量时,输入 CLI 时不能使用 ival=totbus() ,而应该使用 ival=psspy.totbus() 。在 CLI 中使用 print(ival) 来显示 ival 。 Python 语法也可以使其简化,如: from psspy import * 。采用何种编程风格由用户决定。
! p: v: ~( H6 w ' n5 E7 I0 c. |' F2 B! {- y. h
在 IDLE 编译平台测试 python 代码
n% u9 |4 `6 I" s$ ~$ N$ t6 I 我花了让人崩溃的长时间来寻找下面这个加入到 python 路径的方法。
) Z1 l& S* R- @: X) }' f# ^& o 当你想将自己写的一个模块导入 python ,或者要把一个处在其他文件夹下的模块加入到 python 文件中去,那么下面的方法就大有用啦: 4 N3 h) a' e( z+ ]3 S4 z7 S" e
1、
; i; J, m/ E1 Z- i 创建一个扩展名为 .pth 的文件。我们起个很漂亮比尔盖茨又很喜欢的名字 myconfig.pth 。当你在使用记事本之类的编辑器时就要小心了,不要保存成 myconfig.pth.txt 哦。 1 x* |' j- m! t9 q( L" p- F, u4 v
2、3 L. d+ i6 u7 x$ R9 V" u
把这个文件放到 python 根目录下面去,比如我的就是 C:\Python23 。
3 Q: k: w# m$ c 3、
6 M6 z- @9 j, P8 ]! |9 _" N 把每个你想加入 PYTHONPATH 的文件夹全路径,单独一行的加入到这个文件中去。例如,第一行可能就是: C:\aFolder\work\Python 。
4 y# R8 z$ o0 O" Y 这样可以解决下面的错误,当你导入 Spam 这个模块的时候就不会出现: ImportError: No module named Spam.
) D1 {1 `7 }& a 0 r+ }* {' Q2 r
获取和设置网络数据
: F! x" k2 r, ^% d: Q3 C, r: _# d7 B 直接获取数据是很费力费时的事情。在网络中的每个部分都有几个不同的参量,比如负荷。每个参量均必须分别通过函数调用来访问。所以如果要找出每个母线的电压,那么必须对每个母线通过 psspy 对象调用方法 busdat 。它的参数通过一个字符串来解释,就如一开始所说的,这些参数可以在 IPLAN 说明书里找到,同 Python 里的命令名字相同。
( Z e* E$ [) c9 C, ~ 生成
N; o C$ ?5 t* S h6 t' w 字符 / 文件
$ z9 n% f+ {% U. U) _ 格式的报告是可以的,我曾经试图寻找转存这些变量的方法,这样就可以不通过文件系统直接可以通过 python 代码来访问了,但是没找到。如果有谁找到了这种方法请告诉我。如果西门子能改进一下 API 的话这将变得很简单,但是我们现在只能根据已有的东西来做。
% Y7 h6 H; ?; ^- M 反过来,也令人郁闷的是,要设置数据,每个部分的所有参量都要提供出来。幸运的是,有默认的对象 _i 和 _f ,当你不想改变某个值时,插入这些对象即可。当运行和导入我自己的主模块时,我把这些对象存储在 psspy 对象中,这样每当 psspy 对象使用的时候,他们就能使用。对于整形和浮点数参数分别有个默认值(对象),我想他们与各自能取到的最大值相等。
% |+ G4 z. B* M* R# ?9 @! h. u 3 j/ O7 w- x3 D1 A) D. b4 t! o: f/ s
在 Python 中组织 PSSE 数据 , R. ~8 M' {' B5 H4 O
当使用 Python 扩展 PSSE 时,经常没有必要抽取所有的 PSSE 使用的数据。我创建了两个主要的类,一个与一个母线有关,叫做 Bus ,另一个与两条母线有关,叫做 Branch 。我所使用的所有层次如下所示:
& R# }3 n* V U/ ~9 V o
4 W# x3 P- ]$ x1 v' I3 O Machine
· ]/ t! }/ e4 T5 ~
Branch
o
& \: K$ V# n. F Transformer
我把这些对象做一个列表来标识整个的网络。 Python 在快速和简化处理对象列表方面非常擅长。选择或者增加方法可以用于大批 bus 和 branch 类,这些类包含有用于对象中所有参量的处理方法。这样可以提供一种便利:一个部分的任何参量,通过大量的方法均可以对其他部分可见。可是呢,我正在有目的的在我的工作中避免出现这种情况。
1 r, j: e; [/ L o' a 注意当在 Python 中创建一个 PSSE 模型的拷贝时,这些拷贝还处于同步状态。因为潮流计算的原因,一些参数会改变而其他参数不会,例如,改变线路参数并重新计算,貌似改变了复杂(复数)的母线电压,但是不改变负荷。这样必须把从 PSSE 模型里观察到的数据同存储在 python 模型中但是在 PSSE 模型中使用的数据区分开来。
0 u0 Z7 L' u8 C- G* G0 m# l7 \
* Z+ l( ], `4 t% {+ N8 l0 x& ^1 `* t 编译更复杂的 Python 程序
* K, q. z3 m. V4 y 我保留了一个叫做 PSSPYStub.py 的 python 文件,它包含一个仿冒的 psspy 类,其中对应于每个所使用的 psspy 方法均创建一个短小的仿冒方法。这样当我进行简单测试编译 python 代码时把这个 PSSPYStub 类代替 psspy 类,这样就避免了开启 PSSE ,可以允许 python 调试工具的简单使用。
- | ~- h% l+ _0 q4 t# w! A 这个短小的仿冒方法对于 psspy 函数来说是微不足道的,但是在使用一些函数时需要注意,例如 nxtbus 函数,这个函数用来返回 8 A/ v" p2 l1 L( [
调用 inibus 函数之后还没有输出的下一个载入 PSSE 的母线节点号。为在这个仿冒方法中实现这个函数,需要在 inibus 函数中创建一个 count 变量,并且在 nxtbus 函数中对这个变量进行递减,并返回一个母线节点号直到其达到 0 。 . Z6 @( U4 B+ }
注意!更难的 python 内容 ) W8 g2 x$ \# A1 ?9 W* L* F
调试过后的代码如果能在 PSSE 中运行同样可以使用 IDLE 界面。打开 PSSE 和 CLI 界面并且输入并运行 Runldle.py ,不出意外地话会打开 IDLE 编程环境。通过打开这个调试器运行你的 python 代码。 & M* M9 T9 v- v0 P( `! M$ R( z
使用这个方法有点不同,当你使用 import 命令时,主要的 psspy 对象和伴随的 _i 和 _f 将在你运行中的模块中不能访问。在 PSSE 中的“ Run Auto ”按钮不同于标准 python 中的 import 命令。“”按钮允许运行一个 python 模块时把 psspy 作为一个全局变量。但是在你得到的 IDLE 解释器中的 import 命令不会把 psspy 作为导入模块的全局变量。在 python 中一个全局变量只能在本模块中可见。
; S K! x0 ?( N9 R 为了解决这个问题,我在本应在“ Run Auto ”按钮运行的代码之外创建了一个类, __init__ 的定义会将psspy对象作为一个参数分配到一个类变量中。这样他在整个类中都会可见。 ) r5 s4 s$ p: O2 g: t8 G$ ^1 A* ]
包含这个类的模块被导入后,一个这种类的对象会被创建,这样这个类的所有函数都能被调用。例如,这个类可能包含如下的函数: : z+ V5 M# e' U8 s U+ Z1 U/ }
·
0 u4 v7 c8 j4 o loadRawFile 装载一个 raw 文件到 PSSE 中 .
·; n" u! q, `7 P6 T8 q" Q2 g
parseDataFile
/ ^* g% d. a9 G: B. Z m! S8 C. n 装载额外的仿真参数 .
·
9 G! b* P v. e runSimulation - O* P: d- j( k& z
运行一个潮流计算 .
·
; r! J' T6 I/ d saveResults
' f, D& n9 y O- ` 存储重要的数据到文件中 .
你可以从Python解释器中手动调用这些函数。 e X5 a/ ~: `( T2 U! Q* n
同样的,你可以使用__init__ - Q% n$ t5 ]4 ?( m; c2 K7 B
函数仅按顺序包含所有你想执行的脚本。这些命令就像是初始化 PSSE 和你想使用的任何 Python 对象。最终这些命令会使 PSSE 进行一些例如潮流的计算。
3 G3 s) }9 _; g6 A9 B - P9 v* X6 ^0 ^3 |" e
处理数组
7 N$ \ k/ ^! i) n" d$ F 我逐渐能熟练处理数组和数字矩阵了。不在重复强调使用列表类型,而使用数字类型,这样会比使用列表快的多,而且你的程序读起来也更直观。它同在 MATLAB 里面使用一样方便,数字类型不仅仅是一个 Python 模块,也与最优化并且与编译的矩阵运算代码相关。
- T! S$ P7 o; Y0 F 与数字类型相关的信息可以在 这里 找到
% y4 s+ u+ ?5 x% y/ S, j9 H 在导入数字类型(Numeric)之前必须进行安装。使用在Python中导入模块的标准方法 进行导入,你可以通过dir(Numeric)命令观察可用的方法及类结构。
& I5 D a- q* `: i. P8 M$ a 我使用的最多的方法是: ) M+ s; l: K* W" t/ I
Array 9 D W0 n" Y* B
使用这种方法从一个列表里创建数组可以包含一个数字或者一列数字。加些参数就可以使数组包含浮点数等等。例如:Numeric.array([1,2,3],'f')可以返回一个包含有浮点数1,2和3的数组。 - v8 \5 K0 B; A& z& T" d X
Arrange
T8 \% r$ I. \: \* l- T/ ] 这个方法根据你的描述来产生一个范围的等差数字,例如:Numeric.arange(0,0.5,0.1) 返回 & A3 F f& G9 k. c
array([ 0. , 0.1, 0.2, 0.3, 0.4]) . 注意数组中不包括最终值。 6 q5 J: U: [4 g
Transpose
3 c/ G2 h% ~2 Z- D& H" h 给定 a=array([array([0,1,2]),array([3,4,5])]) , 那么 transpose(a) 可以返回: 4 B5 e6 n0 L$ I9 l$ h& `3 {0 p
, Q% l8 O2 G3 K* t/ Y0 o# O6 O ' E3 F' M3 d7 v1 r
5 C+ u. `+ S. T v9 n. F# _3 t$ A array([[0, 3],
' ^. g, B. {- S" b9 f; \ 9 U( q8 q# n1 C0 y, d z4 ?' B' @- L
[1, 4],( k8 @: k9 ]" H$ N$ f7 L% r
[2, 5]])数字数组中最有用的功能是数字定义了标准的操作如何执行例如加法,乘法和逻辑运算。 保存数据以便在MATLAB中使用 导入数据到MATLAB中有多种方法,广义来说主要有以下三种: l% M/ h9 O5 G% k. O
MATLAB 格式.mat文件 l; {+ I% F1 F7 d6 `4 q2 O1 u
ASCII 编码数据文件 l: ^5 @: d% O4 i+ D$ V
文本文件 MATLAB 格式.mat文件 这个格式需要各种文件头。我不知道有没有pythond的模块可以执行保存为.mat文件,但是在这个文件中细节描述了这种格式。 ASCII 编码数据文件 就我所注意到的,这些文件只能存储数字信息。这些值按照行列来组织,每一行的列数必须相等。列之间可以使用各种字符隔开,比较常用的选择就是使用tab建或者空格键键入的空格。行使用转行分开,这样他们在文本编辑器中显示新的一行。 下面是一个例子: from Numeric import * m = array([1,10,2,20]) m = reshape(m,(2,2)) f = open('test.dat','w') for r in m:
) q5 M! E+ F* x- ~" |, M% J& m for c in r:
& \- }6 ] T$ P% K f.write("%i ")
: M7 V8 \' o6 K) I f.write("\n")f.close()
在MATLAB里面,下面的语句就可以装载test.dat文件并把它作为一个变量test。 文本文件 最强大的,看起来可以被容易执行的方法,就是生成一个文本字符串,当执行的时候可以生成需要的对象或者在MATLAB中的对象。但是我想这可能不是最有效的办法。 例如,设想一个bus类含有fields, voltage, load 和 generation over time,如下所示: class Bus: 3 }1 u" L7 O* P7 Z/ O
number=None
, G0 k8 y& P# y Pgen=[]
( J, @7 t/ a) ~+ B6 y Pload=[]
+ }3 e1 m; f. p; l V=[]
2 Z7 O b9 {/ N def __init__(self,number):3 m1 e. o: f+ o& P Z
4 X: s9 C I. M* p: y/ u1 P; ]
self.number = number $ D. _( X3 P3 S0 Z
def storeStepData(self,g,l,v):5 V! n3 |7 I8 W* [, o! R+ L) J
/ ?$ j/ N9 r0 Y
self.Pgen.append(g)
3 }% J/ d& |0 N* v self.Pload.append(l). `$ P! C" k1 ~$ D4 }& s1 U
self.V.append(v)
在类的定义中加入一种方法,这个方法可以产生适用于MATLAB的自我解释的字符串描述。 ! W1 @! Z& B! Z5 L+ |. \
def matlabString(self):
) C) W9 q* A: q2 A s1="bus%i.Pgen=%s"%(self.number,str(self.Pgen))3 e0 N& {; | [
s2="bus%i.Pload=%s"%(self.number,str(self.Pload)), K' N4 `- f& j5 L& e/ z7 F0 [
s3="bus%i.V=%s"%(self.number,str(self.V)); x8 Q: T2 Z0 w3 Y
return("%s;%s;%s;")%(s1,s2,s3)
在 Python 中,在主程序中插入退出字符来建立包含变量的字符串。 紧接着一个包含需要插入主程序的值的 tuple( 元组 ) 。最后一样连接三个其他的字符串,并在每个字符串后面放置“ ; ”字符,这同时把三个命令分开,并能组织 MATLAB 输出这些结构内容。 4 L, M% s9 P: f) i5 d+ ~
下面的测试代码在 Python 中运行会产生一个示例保存文件。 i! c! {2 s8 ?& F
b=Bus(10) b.storeStepData(1,2,3) b.storeStepData(1.5,2.3,3.1) f=open('test.dat','w') f.write(b.matlabString()) f.close()
% K/ X0 |8 b, h
下面一段 MATLAB 中的代码可以产生一个结构,包含有从 bus 工程中保存下来的数据
: K c, T+ M! m; T( H fhandle=fopen('test.dat') fstring=fscanf(fhandle,'%s') eval(fstring)
结构“ bus10 ”现在可以在键入的 matlab 代码中起到作用。其中“ 10 ”部分是根据 bus 和母线节点号在 matlabSring() 函数中生成的。 4 `/ x$ Y0 H2 ]# `8 R
8 A: b, I6 R0 p+ Q& j" c" t 异常处理 ! @3 X4 y! R+ t7 J. J0 h8 k& P
异常处理,是现代编程语言中使用的一种错误处理和诊断的方法。例如,没有必要在其他返回值再增加一个错误代码,而是如果有错误发生的话直接引发一个异常。这样做的好处是,错误不会被默认忽略。如果不检查错误代码的话,这个方法的失败就可能不被发现,可能会造成跟踪程序错误的困难。
- d& p) F6 M. ~( {! T; g# o: d; s 在 python.org有一个异常处理的手册 ,它同样包含有用户自定义异常。
# n2 O) X- [( P/ B! n 9 l. y- Q S, x& Y
但是 PSSE 提供的 API 接口不适用异常处理,这样你在每次调用 API 时,都要严格检查错误代码变得很重要。我建立了一个与 PSSE 交互的特殊类来解决这个问题,所有 API 的调用都要适用这些类。每个 API 调用的时候都跟随错误代码检查。如果错误代码只是的是失败,则就引发一个异常,根据具体的错误代码使用合适的,定制的异常和信息。
k( v+ x4 y9 M, K
5 g G2 V5 c( w' ~, f+ w! s! M5 h5 ~ 使用 Python 模块 random
6 k5 A- f6 }4 E6 \8 T/ M Random模块 可以根据不同的分布状态产生随机数,并允许控制产生的状态。使用在 Python 中引入 模块的标准方法 ,导入 random 。你可以通过 dir(random) 函数查看可用的方法和类的构造。
) L/ x5 O; W+ \+ m 我用的最多的方法是:
4 E9 v7 i: E/ ]- S/ f. ?) h# Q9 j Getstate ' D4 P" a5 l6 F/ K3 f+ L
这个方法返回产生随机数的当前状态。这个状态可以保存并在稍后用来设置模块的状态回到以前状态,具体方法可以查看下面的 setstate ,例如: savedState = random.getstate() .
. H; R; S& ~, p/ a) Y* a Jumpahead : `3 N$ q8 H: r$ K
这个方法用来改变当前状态,只要模块不会产生一系列数字覆盖在执行 jumping ahead 命令之前的状态。例如 random.jumpahead(42) 9 ~) `0 e, D& ]/ y
normalvariate 和 uniform
4 D& K* d( q* C 这些方法是少数返回随即数的方法的两个,他们分别返回符合正态分布和平均分布的随机数。方法 mormalvariate 需要平均值和标准偏差两个参数,方法 uniform 需要最大值最小值得范围 rang[min, max) 。 3 \ f% i# D9 X8 Y) Y# S! \" e, G
例如: random.normalvariate(0,1) 和random.normalvariate(0,10) . 1 |8 p, c+ z8 P- x* [2 a5 A, M
Seed
3 ~0 r1 d# T$ n+ ~$ n% o 这种方法可以使用任何不规则格式设置随机数产生器的状态,例如字符串。这样,不必要首先在文件系统中存储一个状态类就可以设置状态了。例如: random.seed("My! I do like coding, it's such fun!") . * A4 w; p# l* s3 T; O" ?3 ^0 D
Setstate
& R5 ]( P, S3 X9 u8 {# C* L 这种方法使随机数产生器的状态设置为之前使用 getstate 保存的状态,例如: random.setstate(saveState) 。
- x! [: b W% x 我发现在使用 random 模块时创建自己的类是十分有用的,在其中,产生随机数的方法保证了 random 模块的状态不发生变化。这可能被忽略了,但是一次只有一个 random 模块导入。那样我就可以创建多个自定义的 random 模块来产生随机数字,没有必要干涉每个产生器的状态。
8 k# J6 _) Q9 o P3 x: ?! u
, S" s& n8 }9 u$ V5 U 其他注意事项
' d, q7 J) K; w3 w" z4 G l
! h" t; X2 }3 Z- k6 f' e" v 我也写了一些实现以下功能的 python 代码: c+ j. H# m- c) Y; y/ e
²7 o1 Z6 C( b5 D s
潮流运算和获取网络所有的数据。 6 P, n, b; L* Y% n W
²
: p5 _$ V9 `9 e5 `% U' a3 [ 基本的低 / 高 / 带通滤波器。 $ g# _6 l9 S4 P6 w
²0 e2 c z4 L u/ g
梯形三角形模糊逻辑 / V" @! J* L: ]) O* ~7 r
², |8 V( C2 p8 a
如下所说的我不能使用 scipy ,基本的插补类 0 [6 z8 K1 c% X) P
l! x+ M6 A: s P6 [ t5 i
我在 PSSE CLI 或者在 CLI 运行代码时,使用 import scipy 时出现了问题,当我从 Fortran machine 里导入时,出现了浮点数溢出的现象,所以我尝试了很多次重建 wheel 因为我不能使用。
4 f J) |/ L8 G d9 i" w l o+ z3 ^% O; r, J$ {
我又碰到浮点数溢出的现象了,这次是当我把自己的模糊逻辑模块应应用为 PSSE 扩展脚本的一部分时出现的,我想这是甚至因为我再 scipy 只引入了 interp 部分造成的。如果只在 python 解释器下测试的话就不会出现这个问题。
2 n) Z% U* S2 ^% e7 a. `) @; _ l7 ^: a. b# u Y8 p- k7 v+ ^
我碰到了一个 Fortran 的错误,可能是由 Python 引起的。我忘记了它是不是被零除过,还是因为我试图当作一个 list 来读取数字,或者这两个方面都有。这里的课程可能会在 Python 没有抛出异常的地方出现问题,结果造成了不良的运行时错误,但是原因可能很简单而且错误能确定。
* X$ ~ ]+ E8 W+ L% j l
8 H/ ]' N1 i( r2 M4 _ 在潮流运行完毕之后需要使用 solved ()来检验潮流是不是收敛了。如果潮流没有收敛,一些网络可能从步参数比默认参数小 100 或者更多( -_- )受益。当先减小补偿时,代数限制应该也要放宽。
- _- e( m* @/ I6 W6 e6 z' K, { |# z+ N l7 R# K$ E4 S& k/ T4 @# |! ~
读一下说明书。他有不同求解方法的正面和反面介绍。 9 F3 e! u, ~, c( D) l5 T+ F5 l
l& U3 Y7 J, c- H4 y* e Z9 E0 N
在我的版本中 GSLV 会让电脑崩溃,更郁闷的是他在我的工程里运行比在 MSLV 中运行的好,而我只能呢个使用 MSLV 进行扩展仿真。
& a, g2 W) E3 Y! W | l
8 Y8 r/ K" j6 U7 @, ^/ D 警告。逐步备份你的 .raw 文件。 PSSE 有个坏习惯,覆盖的时候没有警告,即使那个模式已经关掉。 PSSE 还可能把不能自动调节的变压器,做成能自动调节的变压器模型。我不确定什么时候会发生这种情况,但是我已经被这个问题浪费了一些时间了。
0 o* T2 l, Z- {+ _1 h/ ~7 m' t l: C7 K2 q" A/ u
一个不错的网站 ,关于 Python 的。 " y3 E( x- C) O7 D7 _) h& U
4 Q- l: W- p+ }* M! ^4 k4 K1 x, c 建议和评论 ; }# j M5 z% F+ L$ ^& J6 e, R$ f
如果有任何建议或者改进欢迎发 Email 给我,我可以分别把具体的 PSSE 功能性的使用和详细地执行信息联系起来。 Email jbo at ee ed ac uk 。
, _ |/ g8 V# _3 a& O7 Y 发邮件给我的任何问题或者在 AQ网页 里的问题我都会回答的。 , ^! G' u, o! D6 n2 o g
" J* S4 `* t$ \2 R
+ a. K Y8 A7 {4 w: C
[ 本帖最后由 dsyljh 于 2007-9-24 16:27 编辑 ]
使用Python扩展PSSE.rar
21.65 KB, 下载次数: 106, 下载积分: 威望 -2 点, 学分 -5 点
介绍在PSSE中使用Python脚本语言的初级入门方法和一些代码
评分
查看全部评分
楼主热帖