谐振子系统的仿真速度

Question

我正在编写一个模拟以获得谐波振荡器系统的真实速度

其中 P=[p1 p2;p2 p3] 可以使用 Rung-Kutta 积分方法找到 P(0)=[1 0; 0 1]

现在，我想编写 matlab 代码来获取系统的真实位置z 并估计其中 x=[x1 x2]' 的 x。这是我查找 x 和 z 的代码。然而，结果并不正确。你能帮我修改一下吗。这是我的matlab代码

function  z=getPos(x)
  %% Function to get velocity
  v=0+1*randn; %r=1;
  z=[0 1]*x+v; 
return
    

这是我的结果

让我看看帮助我是解决方案的正确实施。谢谢

Answer 1

您的代码可以返回一个负值，因为 randn 是一个标准差为 1 的正态分布随机数。理论上这可以返回任何数字，并且它有大约 15% 的机会小于 -1。

你也是正确的，你应该在每一步更新 x。我会将 x 作为 getPos 的输入，将 xdot 作为输出。然后添加另一个步骤来更新位置。