www.pudn.com > matlabsimulink.rar > sy_01.m, change:2010-10-23,size:545b


% 求解发电机转子摇摆曲线的微分方程 
% 在主程序中定义全局变量传递参数 
global y0 Tj Pt E U X1 
y0=2*pi*50;Tj=11.28;Pt=1;E=1.47;U=1; 
% 系统转移电抗 
% X1=2.82;   %故障时 
X1=1.062;       %故障切除后 
% 指定解算微分方程的时间区间 
tspan=[0.0 0.6]; 
% tspan=linspace(0,10,3000); 
% 给定初值向量 
% y1=[31.54*pi/180;1]     %故障时 
y1=[37.307*3.14/180;1.0063] %故障切除后 
% 求解微分方程 
[t,YY]=ode45('power_tra',tspan,y1); 
% 输出求解结果 
x=YY(:,1); 
y=YY(:,2); 
% 曲线绘制 
plot(t+0.1,x*180/pi); 
xlabel('t/s'); 
ylabel('\delta/deg'); 
grid on