www.pudn.com > ship.zip > f.asv, change:2016-08-31,size:254b


function dy = f(t,x) 
global L Lw B d Cwa Cp Cb m Iz rou V n Dp Cpp delta H wp wr P AR Xuu Xvr Xvv Xrr Yv Yvv Yrr Yvrr Yvvr Nv  Nr Nvv Nvvr m mx my Jz Iz 
u = x(1); 
v = x(2); 
r = x(3); 
dy(1) = ((m+my)*r*v+Xuu*u^2+Xvv*v^2+Xvr*v*r+Xrr*r^2+Xp+Xr)/(m+mx)