classical - wd=x(1); delta=x(2); d % Now write differential...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function xdot=classical(t,x) global ws H KD Pm Pemax; % Global parameters used for convenience % You can pass them to this routine or % cut-and-paste % called as [t,x]=ode45('classical', [ 0 tend], [wd0 delta0]'); % Name state variables for readability
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: wd=x(1); delta=x(2); d % Now write differential equations xdot(1)=ws/(2*H)*(Pm-Pemax*sin(delta)-KD*wd/ws); % Acceleration xdot(2)=wd; % Speed deviation xdot=xdot'; x return; r...
View Full Document

Ask a homework question - tutors are online