diff --git a/smooth_resampled_traj.py b/smooth_resampled_traj.py
index d49f5a4..e87a38b 100644
--- a/smooth_resampled_traj.py
+++ b/smooth_resampled_traj.py
@@ -336,6 +336,10 @@ def resample_trajectory(times: np.ndarray, joint_data: dict, target_dt: float, j
     """三次样条重采样，夹爪做限位。返回 (new_times, resampled_joints)。"""
     t_start, t_end = float(times[0]), float(times[-1])
     new_times = np.arange(t_start, t_end, target_dt)
+    if new_times.size == 0 or not np.isclose(new_times[-1], t_end):
+        new_times = np.append(new_times, t_end)
+    else:
+        new_times[-1] = t_end
     resampled = {}
     for name in joint_names:
         cs = CubicSpline(times, joint_data[name], bc_type='natural')
@@ -687,7 +691,7 @@ def visualize_end_effector(
 ):
     """
     绘制轨迹优化前后的机器人末端位置对比图（3D）。
-    需要 7 轴臂关节与 roboticstoolbox；若不满足则退化为关节空间对比图。
+    需要 7 轴臂关节与 roboticstoolbox；若不满足则退化为关节空间图。
     show_interactive=True 时保存后不关闭窗口，并调用 plt.show() 以便旋转/拖动。
     """
     q_before = _get_arm_q_matrix(clean_joints, joint_names)
@@ -925,4 +929,4 @@ def main():
 
 
 if __name__ == "__main__":
-    main()
+    main()
\ No newline at end of file
diff --git a/test_resample_endpoint.py b/test_resample_endpoint.py
new file mode 100644
index 0000000..0867345
--- /dev/null
+++ b/test_resample_endpoint.py
@@ -0,0 +1,33 @@
+import unittest
+
+import numpy as np
+
+from smooth_resampled_traj import resample_trajectory
+
+
+class ResampleTrajectoryEndpointTests(unittest.TestCase):
+    def setUp(self):
+        self.times = np.array([0.0, 0.5, 1.0])
+        self.joint_names = ["joint1"]
+        self.joint_data = {"joint1": np.array([0.0, 0.5, 1.0])}
+
+    def test_keeps_final_timestep_when_dt_divides_span(self):
+        new_times, resampled = resample_trajectory(
+            self.times, self.joint_data, 0.1, self.joint_names
+        )
+
+        self.assertAlmostEqual(new_times[-1], 1.0)
+        self.assertEqual(len(new_times), 11)
+        self.assertAlmostEqual(resampled["joint1"][-1], 1.0)
+
+    def test_keeps_exact_endpoint_with_shorter_final_interval(self):
+        new_times, _ = resample_trajectory(
+            self.times, self.joint_data, 0.3, self.joint_names
+        )
+
+        np.testing.assert_allclose(new_times, [0.0, 0.3, 0.6, 0.9, 1.0])
+        self.assertAlmostEqual(new_times[-1], self.times[-1])
+
+
+if __name__ == "__main__":
+    unittest.main()