% 定义轨道参数 a=314*exp(2); mu = 398600.4418; % 地球引力常数 t0 = 0; % 初始时间 tf = 2*pi*sqrt(a^3/mu); % 终止时间 dt = (tf - t0)/1000; % 时间步长 r0 = [14.4601 2.3155 12.9992]'; % 初始位置矢量 v0 = [-108.9716 20.9281 117.4898]'; % 初始速度矢量 y0 = [r0; v0]; % 定义ODE函数 odefun = @(t,y) [y(4:6); -mu/norm(y(1:3))^3.*y(1:3)]; % 解ODE方程 [t,y] = ode45(odefun, [t0, tf], y0); % 绘制轨道图 plot3(y(:,1), y(:,2), y(:,3)); axis equal; grid on; xlabel('X轴'); ylabel('Y轴'); zlabel('Z轴'); title('地球卫星轨道图');优化代码使得只在10个周期内积分只展示这10个轨道

可以通过修改终止时间 tf 的值来实现只展示10个周期内的轨道。假设要展示的是前10个周期，则tf的值应该为：tf = 10 * 2*pi*sqrt(a^3/mu)。同时，由于时间步长 dt 已经在代码中定义，不需要再修改。修改后的代码如下：

% 定义轨道参数 a=314*exp(2); mu = 398600.4418; % 地球引力常数 t0 = 0; % 初始时间 tf = 10*2*pi*sqrt(a^3/mu); % 终止时间，只展示前10个周期 dt = (tf - t0)/1000; % 时间步长    r0 = [14.4601 2.3155 12.9992]'; % 初始位置矢量 v0 = [-108.9716 20.9281 117.4898]'; % 初始速度矢量 y0 = [r0; v0];  % 定义ODE函数 odefun = @(t,y) [y(4:6); -mu/norm(y(1:3))^3.*y(1:3)];  % 解ODE方程 [t,y] = ode45(odefun, [t0, tf], y0);  % 绘制轨道图 plot3(y(:,1), y(:,2), y(:,3)); axis equal; grid on; xlabel('X轴'); ylabel('Y轴'); zlabel('Z轴'); title('地球卫星轨道图，展示前10个周期');```

运行修改后的代码，即可展示只有10个周期的轨道。