对称分量法求出正序、负序、零序的mathematica程序

wyz2727

(*求三序向量*)

(*A相量幅值*)fa = 10
. ?4 T9 r) i/ j. n1 E6 T) O(*A相量角度*)ja = 0
: W! d9 \: w2 H% C3 Z2 w(*B相量幅值*)fb = 10
: `+ K1 H4 \! o: ^$ @(*B相量角度*)jb = 180
: t, I/ b; {0 u) b. [- R( e(*C相量幅值*)fc = 0
(*C相量角度*)jc = 0

( L& l3 e, D; L8 D3 l2 d(*A相量正序实部*)rea1 = fa*Cos[ja*Pi/180] + fb*Cos[(jb + 120)*Pi/180] + jc*Cos[(jc + 240)*Pi/180]
(*A相量正序虚部*)ima1 = fa*Sin[ja*Pi/180] + fb*Sin[(jb + 120)*Pi/180] + jc*Sin[(jc + 240)*Pi/180]
(*A相量负序实部*)rea2 = fa*Cos[ja*Pi/180] + fb*Cos[(jb + 240)*Pi/180] + jc*Cos[(jc + 120)*Pi/180]
(*A相量负序虚部*)ima2 = fa*Sin[ja*Pi/180] + fb*Sin[(jb + 240)*Pi/180] + jc*Sin[(jc + 120)*Pi/180]
: I4 \5 ^; `6 R% W: q* n(*A相量零序实部*)rea0 = fa*Cos[ja*Pi/180] + fb*Cos[jb*Pi/180] + jc*Cos[jc*Pi/180]
(*A相量零序虚部*)ima0 = fa*Sin[ja*Pi/180] + fb*Sin[jb*Pi/180] + jc*Sin[jc*Pi/180]
( A* O% Z/ w6 ~# ?1 W7 K
( ^4 |  A; D' s, Y) n0 V(*A相量正序实部*)rea1 = N[rea1/3]
(*A相量正序虚部*)ima1 = N[ima1/3]
; p. R  y$ s* L+ I(*A相量负序实部*)rea2 = N[rea2/3]
(*A相量负序虚部*)ima2 = N[ima2/3]
(*A相量零序实部*)rea0 = N[rea0/3]
(*A相量零序虚部*)ima0 = N[ima0/3]

) V0 h6 T* n$ B& J(*作三序图*)sawline = Line[{{0,0},{rea1,ima1},
                            {rea1,0},{rea1,ima1},
                            {0,ima1},{rea1,ima1},
                            {0,0},{rea2,ima2},
5 L5 J1 T3 L9 ~6 L                            {rea2,0},{rea2,ima2},
" b' b! J5 h. N" X8 H, L/ y                            {0,ima2},{rea2,ima2},
& {/ @5 d$ G3 l+ Q$ |                            {0,0},{rea0,ima0},
, v2 \0 S$ A7 b$ p5 W# m                            {rea0,0},{rea0,ima0},
                            {0,ima0},{rea0,ima0}}]
(*作三序图*)Show[Graphics[sawline],
            Graphics[{Text["    正序", {rea1, ima1}]}],
3 x& J: w0 d! C2 j6 V( J' F  X' q            Graphics[{Text["    负序", {rea2, ima2}]}],
            Graphics[{Text["    零序", {rea0, ima0}]}],
0 f* i4 n6 |0 |6 g; _4 P1 K2 z; E            Axes ->True]

$ y8 A( H" _% R
; S' Y) _! ^2 E% _8 I3 s5 v以上是我用mathematica写的程序，不知对不对，还请大虾们指正。
) f! d, L( k& Q( T, J最近我正在学习故障时的三序电流，苦于手里没有例子，不知道有谁可以贡献一些数据。
& ]- v, S1 P+ S( t3 u, c+ m- Y$ _我想问一下：
. G5 `! x3 O3 B0 S6 ~1 E3 m假设有{A相：10安0度  B相：10安240度  C相：100安0度}，这组数据能大致表示大电流接地系统的C相单相接地吗？
假设有{A相：100安0度 B相：100安180度 C相：10安120度}，这组数据能大致表示系统的AB相间短路吗？

' `/ O( Q$ Q% n; x# K; w7 k[ 本帖最后由 wyz2727 于 2008-7-4 00:49 编辑 ]

sadever

相信这里的人用Matlab比较多，要数据，你可以用matlab仿

damayi

原帖由 wyz2727 于 2008-7-4 00:27 发表

                               
登录/注册后可看大图

(*求三序向量*)

$ O9 }; j) q3 k# v9 v; J- b(*A相量幅值*)fa = 10
+ E# o/ }4 r: [% G5 v(*A相量角度*)ja = 0
( e. O) t: O; e/ P) Y/ p(*B相量幅值*)fb = 10
4 |; v+ a, N0 K& C# a; t8 m* Y(*B相量角度*)jb = 180
& }+ d# [' L. R* O1 j(*C相量幅值*)fc = 0
(*C相量角度*)jc = 0

; P" n% Y( c" h(*A相量正序实部*)rea1 = fa*Cos[ja*Pi/180] + fb*Cos ...

. K( W1 B  {' p  q1 y$ D写的太复杂了，算正负零序分量，用Mathematica简直是杀鸡用牛刀。

afandi

不错，支持一个。
( i+ w' u6 q* ]& T7 `

wwwhjkhjk

非常好，感谢楼主

sdad

謝謝樓主提供!!!

kevinxjtu

感谢楼主分享

zouyanyun12

谢谢，赶紧分享下

paradisesdu

谢谢，赶紧分享下

mxj_123

感谢楼主分享啊