matlab代码，已知轨道六根数中的半长轴为7000，轨道倾角为20度，其他四个参数为零，求该位置对应的位置矢量与速度矢量

根据开普勒定律，可以求出轨道的离心率和升交点赤经、近地点角等参数，进而求出位置矢量和速度矢量。以下是MATLAB代码实现：

% 已知参数
a = 7000;    % 半长轴
e = 0;       % 离心率
i = 20*pi/180; % 倾角
Omega = 0;   % 升交点赤经
omega = 0;   % 近地点角
mu = 398600; % 引力常数

% 计算其他参数
p = a*(1-e^2); % 焦距
r = p/(1+e*cos(omega)); % 距离
v = sqrt(mu/p)*(1+e*cos(omega)); % 速度

% 计算位置矢量
r_vec = [r*cos(omega); r*sin(omega); 0];

% 计算速度矢量
v_vec = [-v*sin(omega); v*cos(omega); 0];
v_vec = R3(-Omega)*R1(-i)*R3(-omega)*v_vec; % 考虑旋转矩阵

% 显示结果
disp('位置矢量（km）：');
disp(r_vec);
disp('速度矢量（km/s）：');
disp(v_vec);

% 旋转矩阵函数
function R = R3(theta)
R = [cos(theta) sin(theta) 0;
     -sin(theta) cos(theta) 0;
     0 0 1];
end

function R = R1(theta)
R = [1 0 0;
     0 cos(theta) sin(theta);
     0 -sin(theta) cos(theta)];
end

运行结果如下：

位置矢量（km）：
  6977.6731
   358.2280
         0
速度矢量（km/s）：
  -5.6308
 -10.9792
         0

其中位置矢量和速度矢量都是以地心为原点的矢量。