Re: [微積]轉換成離散的方程式

[chtm (臨危ψ)]

dx/dt=ax-bxy-cxz
dy/dt=-dy+exy-fyz
dz/dt=-gz+hyz

求：

Nx(t+1)=?
Ny(t+1)=?
Nz(t+1)=?

後項差分版

Nx(t+1)=Nx(t) + [aNx(t+1)-bNx(t+1)Ny(t+1)-cNx(t+1)Nz(t+1)]dt
Ny(t+1)=Ny(t) + [-dNy(t+1)+eNx(t+1)Ny(t+1)-fNy(t+1)]dt
Nz(t+1)=Nz(t) + [-gNz(t+1)+hNy(t+1)Nz(t+1)]dt

前項差分版

Nx(t+1)=Nx(t) + [aNx(t)-bNx(t)Ny(t)-cNx(t)Nz(t)]dt
Ny(t+1)=Ny(t) + [-dNy(t)+eNx(t)Ny(t)-fNy(t)]dt
Nz(t+1)=Nz(t) + [-gNz(t)+hNy(t)Nz(t)]dt


--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 140.112.43.3
→ chtm        :差異在速度項要選取t的值還是t+1的值