The graphical representation of rotor angle δ is given below,
The graphical representation of rotor speed ω is given below,
Codes:
1.
%main_function function
delta0=swing_equation(t,d)
H=2;Pm=0.8;Pmax=2;f=50;D=0.1;w=2*pi*f; delta0=zeros(size(d)); delta0(1)=d(2); if t<=0.2 delta0(2)=(0.5*w/H)*(Pm-Pmax*sin(d(1))+D); elseif t>=0.2 delta0(2)=(0.5*w/H)*(Pm-Pmax*sin(d(1))); end end
2.
close all;clear all;clc; d0=[0.4115 0];
[t,d]=ode45('swing_equation',[0 3],d0); figure(1) plot(t,d(:,1)*(180/pi)) title('Rotor angle vs time') grid on; hold on; ylabel('Angle,delta(degree)') xlabel('time (seconds)') figure(2) plot(t,d(:,2)+314.20) title('Rotor speed vs time') grid on xlabel('time(seconds)') ylabel('speed deviation,omega(rad/sec)')
إرسال تعليق