将上述轨道积分10个周期，用Matlab在直角坐标系画出该10圈轨道图。

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

% 定义轨道参数
a = 20000; % 半长轴
e = 0.5; % 离心率
i = 30; % 轨道倾角
Omega = 60; % 升交点赤经
omega = 90; % 近地点幅角
M0 = 0; % 初始平近点角

% 定义轨道周期和时间步长
T = 2*pi*sqrt(a^3/398600.4418); % 轨道周期
dt = T/1000; % 时间步长

% 初始化轨道参数
E0 = kepler(M0,e); % 初始偏近点角
n = sqrt(398600.4418/a^3); % 平均角速度
t = 0:dt:T*10; % 时间范围

% 计算轨道参数随时间的变化
E = zeros(size(t));
for k = 1:length(t)
    M = M0 + n*t(k);
    E(k) = kepler(M,e);
end
theta = 2*atan(sqrt((1+e)/(1-e))*tan(E/2));
r = a*(1-e^2)./(1+e*cos(theta));

% 计算轨道在直角坐标系下的坐标
x = r.*cos(theta);
y = r.*sin(theta);

% 绘制轨道图
plot(x,y);
axis equal;
grid on;
xlabel('X轴');
ylabel('Y轴');
title('地球卫星轨道图');

% Kepler方程求解函数
function E = kepler(M,e)
    E = M;
    while abs(E-e*sin(E)-M) > 1e-8
        E = E - (E-e*sin(E)-M)/(1-e*cos(E));
    end
end

该代码可以计算出指定轨道参数下地球卫星的轨道数据，并绘制出轨道图。您可以根据需要修改轨道参数，比如半长轴、离心率、轨道倾角等等，来得到不同的轨道图。