How to fix the discrepancy in the results obtained by ode45 versus Newton-Raphson method?

6 次查看（过去 30 天）
Minh Tran
Minh Tran 2021-7-12
Hello everyone,
I am trying to solve a 1D vibration problem: . Where .
With this problem, I have tried two ways:
  1. I discretized the ϕ in an N-loop and update it after each time step . Inside this N-loop, I use ode45 to solve the following equation: .
  2. I discretized the ϕ in an Newmark N-loop and update it after each time step . Inside this Newmark N-loop, I use Newton-Raphson to solve and update the above equation.
Theoretically, I expect there should be similar outputs on both of the above methods. But they are not. Is there any ideas why we have this discrepancy here? I am attaching the m-script here so that you can easily imagine the both algorithms and outputs.
Thanks in advance. Minh

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by