用matlab实现将ECEF坐标转化为WGS-84坐标，请给出完整且可以直接运行的代码（未安装工具箱

在MATLAB中，如果你不使用特定的坐标转换工具箱，可以手动编写代码来实现从ECEF（Earth-Centered Earth-Fixed）坐标到WGS84坐标系统的转换。这里有一个简单的示例，但请注意，它假设你知道经纬度（φ, λ）和海拔高度（h）。由于没有特定的工具箱，我们不能直接调用现成的函数，但你可以按照这个公式进行计算：

% 假设你有ECEF坐标(x, y, z)
[x, y, z] = ...; % 这里应替换为你实际的ECEF坐标

% WGS84椭球参数
a = 6378137; % 赤道半径
f = 1/298.257223563; % 反离心率
b = a*(1-f); % 极半径
e_squared = f*(2-f); % 第二偏心率的平方
N = a/sqrt(1-e_squared*sin(phi).^2); % 地理坐标系下的大地水准面距离

% 经纬度和海拔已知，代入转换公式
lat_rad = atan2(sqrt(1-e_squared)*z, sqrt(1-e_squared*y + e_squared*z*cos(lat_rad*pi/180))); % φ in radians
lon_rad = atan2(y, x); % λ in radians

% 计算WGS84坐标
lat_degrees = lat_rad*180/pi; % 转换为度
lon_degrees = lon_rad*180/pi;

% 如果需要考虑海拔高度，高度h可以添加到相应的分量上
% WGS84坐标 (lat_degrees, lon_degrees, h)

注意：这个代码没有考虑到精度问题，如果需要高精度的结果，可能需要引入一些额外的库或算法。同时，确保输入的φ和λ是以度为单位的。