写一个matlab代码，要求如下：已知高度，发射点经度，发射点纬度，地心方位角，速度大小和速度倾角，从而求出机械能，动量矩，半通径，半长轴和偏心率

以下是一个可能的 Matlab 代码，需要输入一些参数，包括高度 h、发射点经度 lon、发射点纬度 lat、地心方位角 azimuth、速度大小 v、速度倾角 elevation。该代码使用 WGS 84 椭球模型计算。

% Constants
mu = 3.986004418e14; % Earth gravitational parameter
R = 6378137.0; % Earth radius
f = 1/298.257223563; % Earth flattening
omega = 7.292115e-5; % Earth rotation rate
g0 = 9.80665; % Standard gravity

% Convert latitude and longitude to ECEF coordinates
lat = deg2rad(lat);
lon = deg2rad(lon);
N = R / sqrt(1 - f*(2-f)*sin(lat)^2);
x = (N + h) * cos(lat) * cos(lon);
y = (N + h) * cos(lat) * sin(lon);
z = (N*(1-f)^2 + h) * sin(lat);
r_ecef = [x; y; z];

% Calculate velocity in ECEF coordinates
azimuth = deg2rad(azimuth);
elevation = deg2rad(elevation);
v_horizontal = v * cos(elevation);
v_vertical = v * sin(elevation);
v_ecef = [v_horizontal*cos(azimuth); v_horizontal*sin(azimuth); v_vertical];

% Calculate specific mechanical energy and specific angular momentum
v_ecef_norm = norm(v_ecef);
r_ecef_norm = norm(r_ecef);
v_r = dot(v_ecef, r_ecef) / r_ecef_norm;
v_theta = sqrt(v_ecef_norm^2 - v_r^2);
specific_mechanical_energy = v_ecef_norm^2/2 - mu/r_ecef_norm;
specific_angular_momentum = cross(r_ecef, v_ecef);

% Calculate semi-latus rectum, specific kinetic energy, specific potential energy
semi_latus_rectum = specific_angular_momentum^2 / mu;
specific_kinetic_energy = v_ecef_norm^2/2;
specific_potential_energy = -mu / r_ecef_norm;

% Calculate semi-major axis, eccentricity, perigee radius, apogee radius
semi_major_axis = -mu / (2*specific_mechanical_energy);
eccentricity = sqrt(1 + 2*specific_mechanical_energy*(specific_angular_momentum^2)/(mu^2));
perigee_radius = semi_latus_rectum / (1 + eccentricity);
apogee_radius = semi_latus_rectum / (1 - eccentricity);

% Calculate mechanical energy, momentum and eccentricity vector
mechanical_energy = specific_mechanical_energy * mu;
momentum = r_ecef_norm * v_theta;
eccentricity_vector = (1/mu) .* ((v_ecef_norm^2 - mu/r_ecef_norm) .* r_ecef - v_r .* v_ecef);

此代码计算了许多值，其中包括特定机械能、特定动量矩、半通径、半长轴和偏心率等。请根据需要选择所需的值。