EDA365电子论坛网
标题:
发电机暂态摇摆曲线仿真程序
[打印本页]
作者:
小白的白
时间:
2020-10-20 17:29
标题:
发电机暂态摇摆曲线仿真程序
对这个曲线程序仿真总是不对,不知道是哪里出问题了
1 `) a) I1 z' D/ l% d7 Q2 @
程序如下
) T+ m/ \3 ~; k
function dy=power(t,y)
* r9 l" }! }* ] `" e
global w0 Pt E U X1 Tj
+ l* `8 O( m2 W* W2 g( {7 c0 U
dy=[(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