已知卫星的目标轨道和当前轨道的速度和半径，都是圆轨道，当前轨道高度9800km，目标轨道高度10000km，霍曼变轨，MATLAB

实现霍曼变轨的 MATLAB 代码如下：

% 常数定义
G = 6.6742e-11; % 万有引力常数
M = 5.97e24; % 地球质量
R = 6378.137; % 地球半径
mu = G * M; % 地球标准引力参数

% 圆轨道速度计算
h = (R + 9800) * 1000; % 当前轨道的轨道高度
vc = sqrt(mu / (R + h)); % 当前轨道的圆轨道速度

% 转移轨道速度计算
h1 = (R + 10000) * 1000; % 目标轨道的轨道高度
va = sqrt(mu / (R + h1)); % 目标轨道的圆轨道速度

a = (h + h1) / 2 + R * 1000; % 椭圆轨道的长半轴
vp = sqrt(2 * mu / (R + h) - mu / a); % 第一个转移点的速度
vb = sqrt(2 * mu / (R + h1) - mu / a); % 第二个转移点的速度

Dv1 = vp - vc; % 第一个变轨点需要的速度增量
Dv2 = va - vb; % 第二个变轨点需要的速度增量

fprintf('第一个转移点速度：%f m/s\n', vp);
fprintf('第二个转移点速度：%f m/s\n', vb);
fprintf('第一个变轨点需要的速度增量：%f m/s\n', Dv1);
fprintf('第二个变轨点需要的速度增量：%f m/s\n', Dv2);

需要注意的是，此代码仅考虑了圆轨道情况，实际情况中卫星的轨道可能是椭圆形的，需要进行修正。同时，霍曼变轨的前提是卫星在同一轨道平面上运动，如果不在同一轨道平面上，则需要进行轨道倾角调整。