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

wyz2727

(*求三序向量*)

(*A相量幅值*)fa = 10
9 m: d+ m8 Q; J2 K(*A相量角度*)ja = 0
4 B- u; v% G; ?/ N0 t(*B相量幅值*)fb = 10
- m0 ^8 A) {% {' u4 H& T(*B相量角度*)jb = 180
; a% Z# E/ G! @) ~1 @. k(*C相量幅值*)fc = 0
(*C相量角度*)jc = 0
% \- D0 c4 `* c9 g$ y
(*A相量正序实部*)rea1 = fa*Cos[ja*Pi/180] + fb*Cos[(jb + 120)*Pi/180] + jc*Cos[(jc + 240)*Pi/180]
) Z3 y! \1 w3 }# s(*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]
  V( L$ ]) E0 X/ ?(*A相量负序虚部*)ima2 = fa*Sin[ja*Pi/180] + fb*Sin[(jb + 240)*Pi/180] + jc*Sin[(jc + 120)*Pi/180]
4 }  ]( w* L) S0 N  y( S(*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]
2 m! J- W$ T  B0 N: [: b% Z# _3 h
(*A相量正序实部*)rea1 = N[rea1/3]
(*A相量正序虚部*)ima1 = N[ima1/3]
(*A相量负序实部*)rea2 = N[rea2/3]
/ h  k% n6 ^' ]" r(*A相量负序虚部*)ima2 = N[ima2/3]
/ G. e* I+ t2 U9 w4 M(*A相量零序实部*)rea0 = N[rea0/3]
% S6 ~, |7 m( z2 ^) }(*A相量零序虚部*)ima0 = N[ima0/3]

7 F' i. o1 |8 b3 q7 }; N(*作三序图*)sawline = Line[{{0,0},{rea1,ima1},
                            {rea1,0},{rea1,ima1},
% \/ Y6 e8 Q5 Q/ w5 m* ~; i7 X9 U                            {0,ima1},{rea1,ima1},
                            {0,0},{rea2,ima2},
7 n( K9 y! F! H4 ?0 a* m                            {rea2,0},{rea2,ima2},
5 D5 _5 P* L  N: I- \                            {0,ima2},{rea2,ima2},
0 a+ Q& P0 p3 x& e, a/ I( P                            {0,0},{rea0,ima0},
) Z0 s7 G1 A& O0 l. _                            {rea0,0},{rea0,ima0},
0 ?6 ?0 @. a1 A" }$ _% n                            {0,ima0},{rea0,ima0}}]
4 G$ u% f# H) E% A(*作三序图*)Show[Graphics[sawline],
' l: Z( O8 K% f            Graphics[{Text["    正序", {rea1, ima1}]}],
7 d8 |9 d+ F+ E- p! f            Graphics[{Text["    负序", {rea2, ima2}]}],
            Graphics[{Text["    零序", {rea0, ima0}]}],
7 T" z9 T- A9 g            Axes ->True]

1 T9 _3 V: M0 r. d" ^- W
" m; i/ n' z0 u* Y% r& B+ U以上是我用mathematica写的程序，不知对不对，还请大虾们指正。
最近我正在学习故障时的三序电流，苦于手里没有例子，不知道有谁可以贡献一些数据。
我想问一下：
( W" n4 g- G0 i' Q# H6 D. z假设有{A相：10安0度  B相：10安240度  C相：100安0度}，这组数据能大致表示大电流接地系统的C相单相接地吗？
假设有{A相：100安0度 B相：100安180度 C相：10安120度}，这组数据能大致表示系统的AB相间短路吗？
' y4 f1 B3 d) z/ P% D
5 n, O2 n! ~/ Z2 {: D- D) c[ 本帖最后由 wyz2727 于 2008-7-4 00:49 编辑 ]

sadever

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

damayi

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

& H1 J, s) v6 l9 c( B! T                               
登录/注册后可看大图

* r0 L6 Q7 G- M- {$ ~/ b, J(*求三序向量*)

(*A相量幅值*)fa = 10
(*A相量角度*)ja = 0
(*B相量幅值*)fb = 10
(*B相量角度*)jb = 180
(*C相量幅值*)fc = 0
( N" T4 s- D! g; I5 O* q1 D. J(*C相量角度*)jc = 0

3 Y/ ~$ d# y1 v5 a9 J(*A相量正序实部*)rea1 = fa*Cos[ja*Pi/180] + fb*Cos ...

7 E, a/ B- M: l- G4 d5 N写的太复杂了，算正负零序分量，用Mathematica简直是杀鸡用牛刀。

afandi

不错，支持一个。

wwwhjkhjk

非常好，感谢楼主

sdad

謝謝樓主提供!!!

kevinxjtu

感谢楼主分享

zouyanyun12

谢谢，赶紧分享下

paradisesdu

谢谢，赶紧分享下

mxj_123

感谢楼主分享啊