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

小白的白

对这个曲线程序仿真总是不对，不知道是哪里出问题了
程序如下
% S+ h1 S, ]" \4 F8 q0 c. _$ Efunction dy=power(t,y)
global w0 Pt E U X1 Tj
dy=[(y(2)-1)*w0;(Pt-sin(y(1))*E*U/X1)/Tj]
9 v' g4 p; U; I! T+ D
" _# r( s: _8 F! C>> clear
>> global w0 E U Pt Tj X1
$ v. e( a# ]# x) `/ V( z; b/ c>> w0=2*pi*50;
8 e* p( S3 ?- w>> E=1.47;
0 G6 S) W! G9 A/ S>> U=1;
>> Pt=1;
* m. p; u  U. O  W2 T5 _2 b5 d>> Tj=11.28;
% V* q; t6 S& g7 l+ ^2 G) t) m>> X1=2.82;
>> tspan=[0,0.3];
/ l5 E5 K9 F0 g4 J>> y1=[31.54*pi/180,1];
>> [t,dy]=ode45('power',tspan,y1);
! N% I% D5 d# B. d- c& s' u>> x=dy(:,1);
>> plot(t,x*180/pi)

kekek

帮你顶一下

llbnmo

来学习一下

IBB-EUT

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

zaiyiaaaa

梁山乐

学习一下。