已知道轨道初始状态矢量（即位置矢量、速度矢量）将上述轨道积分10个周期，用Matlab在直角坐标系画出该10圈轨道图。

好的，下面是一个示例代码，可以用来在 MATLAB 中绘制地球卫星的轨道图：

% 定义轨道参数
mu = 398600.4418; % 地球引力常数
t0 = 0; % 初始时间
tf = 2*pi*sqrt(a^3/mu)*10; % 终止时间
dt = (tf - t0)/1000; % 时间步长

% 定义初始状态矢量
r0 = [20000 0 0]'; % 初始位置矢量
v0 = [0 sqrt(mu/norm(r0))/2 0]'; % 初始速度矢量
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('地球卫星轨道图');

该代码可以通过数值积分计算出指定初始状态矢量下地球卫星的轨道数据，并绘制出轨道图。您可以根据需要修改初始状态矢量，比如位置矢量、速度矢量等等，来得到不同的轨道图。