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

小白的白

对这个曲线程序仿真总是不对，不知道是哪里出问题了
# `6 o, m) w* C9 Q: t+ V程序如下
function dy=power(t,y)
global w0 Pt E U X1 Tj
: T$ @" \1 ]3 ]* @8 Z1 g  ~dy=[(y(2)-1)*w0;(Pt-sin(y(1))*E*U/X1)/Tj]

>> clear
>> global w0 E U Pt Tj X1
>> w0=2*pi*50;
>> E=1.47;
. D/ C2 g& B/ a8 X4 W8 R>> U=1;
>> Pt=1;
>> Tj=11.28;
2 A1 ]; [7 |" ?>> X1=2.82;
>> tspan=[0,0.3];
7 H, I* r/ _% x& p6 }$ |* v>> y1=[31.54*pi/180,1];
, z$ Y. j, E3 c: X$ [4 S>> [t,dy]=ode45('power',tspan,y1);
>> x=dy(:,1);
>> plot(t,x*180/pi)

kekek

帮你顶一下

llbnmo

来学习一下

IBB-EUT

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

zaiyiaaaa

梁山乐

学习一下。