ep3_p2.m

来自「此代码展示了Euler计算方法与差分方法」· M 代码 · 共 25 行

%  Ep3_p2: << 实验三 >>  Euler法  (原式)
%  Designed by FGH

h= 0.05;
H= 120; Vw= 450; Ve= 90;

clear x y;
tk= 0; 
k= 1;
% (3.23)
x(1)= 0; y(1)= 0;
while y(k) < H
   % (3.21)
   x(k+1)= x(k) + Vw*h*(Ve*tk-x(k))/sqrt((Ve*tk-x(k))^2+(H-y(k))^2);
   % (3.22)
   y(k+1)= y(k) + Vw*h/sqrt(1+((Ve*tk-x(k))/(H-y(k)))^2);
   
   k= k+1; tk= tk+h;
end
x
y
sprintf(' k = %d , tk = %7.4f \n',k-1,tk)
sprintf(' L = %8.4f , T = %8.4f \n',x(k),x(k)/Ve)