科里奥利用有限差分法Matlab

Question

我试图创建一个算法来绘制旋转参考帧上移动宗地的轨迹。在包裹上只影响科里奥利的虚拟力量。该议案由

描述

其中

以有限差分形式制定上述方程式
所以包裹的坐标是

计算和绘制坐标的代码部分是

x=zeros(1,T/Dt)
y=zeros(1,T/Dt)

x(1)=x0
y(1)=y0
i = 0
for t = 1:Dt:T
    i = i + 1
    u1 = u0 + Dt*f*u0;
    v1 = v0 - Dt*f*v0;

    x(i+1) = x(i) + Dt*u1;
    y(i+1) = y(i) + Dt*v1;

    u0 = u1; v0 = v1;        

end

plot(x,y,'-*')

结果图是

我对最终结果不满意，并且我想知道我的代码实现是否正确完成。我希望更像是螺旋形或圆形运动 示例结果应该是

任何帮助，将不胜感激。谢谢。