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

小白的白

对这个曲线程序仿真总是不对，不知道是哪里出问题了
) ^! {1 b- T- ~$ y& A' \9 E程序如下
function dy=power(t,y)
( M0 `8 [/ M' q' e; [global w0 Pt E U X1 Tj
dy=[(y(2)-1)*w0;(Pt-sin(y(1))*E*U/X1)/Tj]

  u& Y9 X  C! P>> clear
>> global w0 E U Pt Tj X1
>> w0=2*pi*50;
( @) u8 Q( A, M>> E=1.47;
>> U=1;
8 n# B2 p: M  m5 R( O1 ]>> Pt=1;
>> Tj=11.28;
/ i! W& G: e8 j# J3 C: Z>> X1=2.82;
>> tspan=[0,0.3];
' d& p" H* x* B' L>> y1=[31.54*pi/180,1];
>> [t,dy]=ode45('power',tspan,y1);
" J" E7 P4 W! ?8 S>> x=dy(:,1);
3 A. d: ^8 Q% k8 y4 @$ C  c( S>> plot(t,x*180/pi)
- g& v" i& R  f+ b

kekek

帮你顶一下

llbnmo

来学习一下

IBB-EUT

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

zaiyiaaaa

梁山乐

学习一下。