EDA365电子论坛网

标题: 发电机暂态摇摆曲线仿真程序 [打印本页]

作者: 小白的白    时间: 2020-10-20 17:29
标题: 发电机暂态摇摆曲线仿真程序
对这个曲线程序仿真总是不对,不知道是哪里出问题了1 `) a) I1 z' D/ l% d7 Q2 @
程序如下
) T+ m/ \3 ~; kfunction dy=power(t,y)
* r9 l" }! }* ]  `" eglobal w0 Pt E U X1 Tj
+ l* `8 O( m2 W* W2 g( {7 c0 Udy=[(y(2)-1)*w0;(Pt-sin(y(1))*E*U/X1)/Tj]' x, a& j+ u3 u0 y1 Z

2 U7 W: u0 k& Z  U>> clear
; F$ V0 |- o" z! D; a' b' V- w>> global w0 E U Pt Tj X1
9 u- a5 S, `$ D, M, L  {" G- W>> w0=2*pi*50;
; \5 ?3 s0 z0 c>> E=1.47;
. C, c2 {; W; d( |1 X>> U=1;
% N9 ?/ C2 D- P. t6 |' H>> Pt=1;
) R: A+ g  _; m1 O. m$ K' G# r# A>> Tj=11.28;
4 ^* C: `4 |. n3 ^>> X1=2.82;9 X8 Y; p% b9 K* Z. Y0 j
>> tspan=[0,0.3];& l6 K: @+ k0 W2 I
>> y1=[31.54*pi/180,1];+ k5 Z3 h  i0 w2 m) ]' F
>> [t,dy]=ode45('power',tspan,y1);
- D3 z( y1 F% @0 Y- T5 _>> x=dy(:,1);
! k+ l7 h6 M1 k' ^( A' n( x>> plot(t,x*180/pi)
; U# P% V8 B$ E4 d
作者: kekek    时间: 2020-10-20 18:12
帮你顶一下
作者: llbnmo    时间: 2020-10-20 19:52
来学习一下
作者: IBB-EUT    时间: 2020-10-21 13:39
楼主,程序弄好了没,我也是这种情况
作者: zaiyiaaaa    时间: 2020-10-21 13:44

作者: 梁山乐    时间: 2020-10-22 19:12
学习一下。




欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2