如何使用脚本文件模拟微分方程中的时间延迟？ [MATLAB]

Question

因此，我正在编写一个Xdot = myfunc（t，X，zr）形式的函数来表示微分方程组。 编写此myfunc的最终目标是使用ode45或其他类似方法对其进行仿真。 myfunc主体中将使用一个变量，其时域形式应为zr_r（t）= zr（t-t_delay） 但我不知道如何在myfunc中表示时间延迟？

我已经知道一个解决方案：使用Xdot = myfunc（t，X，zr，zr_delayed）并设置zr_r = zr_delayed 这将需要我预先计算zr_delayed向量

上下文：我正在使用以下方程式模拟半车模型 在第三等式中，zr_r代表发送到后轮胎的道路轮廓（zr（t））。由于半车正在通过道路，因此前部将接收zr_f（t）= zr（t），而后部将接收zr_r（t）= zr（t-t_delay）时域