数值积分与鲁棒性分析实战

最新推荐文章于 2025-11-15 12:51:27 发布

原创

最新推荐文章于 2025-11-15 12:51:27 发布 · 758 阅读

24 ·

CC 4.0 BY-SA版权

文章标签：

#数值积分 # 鲁棒性分析 # MATLAB

1、（MATLAB）针对以下三种容差情况进行数值积分并绘制图形：（i）相对容差（RelTol） = 1e - 3，绝对容差（AbsTol） = 1e - 6；（ii）相对容差 = 1e - 6，绝对容差 = 1e - 9；（iii）相对容差 = 1e - 9，绝对容差 = 1e - 12。同时，在图形中添加图例和坐标轴标签。提示：使用odeset函数设置容差，并将选项传递给ode45。

可按以下步骤实现：

使用 odeset 函数分别设置三种容差情况的选项；
调用 ode45 函数进行数值积分；
绘制图形；
添加图例和坐标轴标签。

示例代码如下：

% 情况（i）
ode_options1 = odeset('RelTol', 1e-3, 'AbsTol', 1e-6);
[tout1, qout1] = ode45(@(time, state) dqdt_attitude_kinematics(time, state), time_interval, q0, ode_options1);

% 情况（ii）
ode_options2 = odeset('RelTol', 1e-6, 'AbsTol', 1e-9);
[tout2, qout2] = ode45(@(time, state) dqdt_attitude_kinematics(time, state), time_interval, q0, ode_options2);

% 情况（iii）
ode_options3 = odeset('RelTol', 1e-9, 'AbsTol', 1e-12);
[tout3, qout3] = ode45(@(time, state) dqdt_attitude_kinematics(time, state), time_interval, q0, ode_options3);

% 绘制图形
figure;
plot(tout1, qout1, tou