发电机暂态摇摆曲线仿真程序

小白的白

对这个曲线程序仿真总是不对，不知道是哪里出问题了
程序如下
+ m1 ?" o% [9 G) S" s# x, P7 Tfunction dy=power(t,y)
global w0 Pt E U X1 Tj
dy=[(y(2)-1)*w0;(Pt-sin(y(1))*E*U/X1)/Tj]

! h  ~$ O5 ~: n& |# |>> clear
>> global w0 E U Pt Tj X1
>> w0=2*pi*50;
# S% P, K% k6 ^1 I1 v- ^0 D>> E=1.47;
>> U=1;
>> Pt=1;
>> Tj=11.28;
& r3 W7 U& X! y/ ]) t& V>> X1=2.82;
>> tspan=[0,0.3];
>> y1=[31.54*pi/180,1];
1 W/ m4 g% p- [  K2 j8 Y( s>> [t,dy]=ode45('power',tspan,y1);
>> x=dy(:,1);
' I$ ?0 g# s! C, s>> plot(t,x*180/pi)
( D# c/ Y  T) W" X% j: ]

kekek

帮你顶一下

llbnmo

来学习一下

IBB-EUT

楼主，程序弄好了没，我也是这种情况

zaiyiaaaa

梁山乐

学习一下。