写一个matlab代码，要求如下：已知高度，发射点经度，发射点纬度，地心方位角，速度大小和速度倾角，从而求出机械能，动量矩，半通径，半长轴和偏心率.请用中文详细注释代码，并且不要出现与上述参量无关的参数

以下是满足要求的MATLAB代码，已经加上了详细的中文注释：

% 已知参数
H = 40000; % 高度 （单位：米）
lon = 120.0; % 发射点经度 （单位：度）
lat = 30.0; % 发射点纬度 （单位：度）
az = 30.0; % 地心方位角 （单位：度）
v = 7800; % 速度大小 （单位：米/秒）
gamma = 20.0; % 速度倾角 （单位：度）

% 常数
G = 6.67430e-11; % 万有引力常数 （单位：m³/(kg·s²)）
M = 5.9722e24; % 地球质量 （单位：kg）
R = 6371000; % 地球半径 （单位：米）

% 将经纬度转换为弧度
lon = deg2rad(lon);
lat = deg2rad(lat);

% 计算速度向量在地心坐标系下的分量
vx = v*cosd(gamma)*cosd(az);
vy = v*cosd(gamma)*sind(az);
vz = v*sind(gamma);

% 计算机械能
E = 0.5*v^2 - G*M/(R+H);

% 计算动量矩
Lx = (R+H)*vy*M;
Ly = -(R+H)*vx*M;
Lz = 0;

% 计算半通径
a = -G*M/(2*E);

% 计算半长轴
h = norm([vx, vy, vz]);
r = R+H;
e = sqrt(1 - ((h^2)/(G*M))^2);
b = a*sqrt(1-e^2);
c = norm([r*cos(lat)*cos(lon), r*cos(lat)*sin(lon), r*sin(lat)]);
p = (a+b+c)/2;
A = sqrt(p*(p-a)*(p-b)*(p-c));
A1 = sqrt(p*(p-a)*(p-b)*(p-c-b));
A2 = sqrt(p*(p-a)*(p-b)*(p-c-a));
A3 = sqrt(p*(p-a)*(p-b)*(p-c-c));
S = A - A1 - A2 - A3;
c1 = S/b;

% 输出结果
fprintf('机械能为：%f J\n', E);
fprintf('动量矩为：(%f, %f, %f) kg*m/s\n', Lx, Ly, Lz);
fprintf('半通径为：%f m\n', a);
fprintf('半长轴为：%f m\n', b);
fprintf('偏心率为：%f\n', e);

这个代码的实现过程是，首先将经纬度转换为弧度。接下来，计算速度向量在地心坐标系下的分量，并利用万有引力定律计算机械能。再根据动量矩的定义公式计算动量矩。然后，根据机械能计算半通径。接着，根据速度向量的模长和地球半径计算半长轴，使用海伦公式计算椭圆面积并计算出椭圆的中心距和偏心率。最后输出结果。