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

小白的白

对这个曲线程序仿真总是不对，不知道是哪里出问题了
. y9 e8 }. Q- K: e0 d1 ^% J3 J4 \程序如下
function dy=power(t,y)
global w0 Pt E U X1 Tj
dy=[(y(2)-1)*w0;(Pt-sin(y(1))*E*U/X1)/Tj]
. V  W# g# f+ q, H, r
9 g( o4 K; j) u' L0 E>> clear
>> global w0 E U Pt Tj X1
>> w0=2*pi*50;
>> E=1.47;
9 s; [0 F4 b. D# G5 R>> U=1;
>> Pt=1;
>> Tj=11.28;
>> X1=2.82;
>> tspan=[0,0.3];
>> y1=[31.54*pi/180,1];
' S% a, A( y3 ?0 C% m>> [t,dy]=ode45('power',tspan,y1);
>> x=dy(:,1);
>> plot(t,x*180/pi)
" F" O8 W) X5 l9 I

kekek

帮你顶一下

llbnmo

来学习一下

IBB-EUT

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

zaiyiaaaa

梁山乐

学习一下。