Hi ysm,
I understand that you are facing issues while utilizing geodetic2aer function. As per the documentation page, it is defined as follows
[az,elev,slantRange] = geodetic2aer(lat,lon,h,lat0,lon0,h0,spheroid)
In reference to your code, lat0, lon0 and h0 correspond to lat_receiver, lon_receiver and height_receiver while lat, lon and h correspond to lat_geosat, lon_geosat and height_geosat. For identification of the parameters, I have positioned the receiver at the origin of the local AER system and transformed the geodetic coordinates of the satellite to the local azimuth-elevation-range (AER) spherical coordinates specified by az, elev, and slantRange. Here is the code snippet for your reference
spheroid = referenceEllipsoid('GRS 80'); % referenceEllipsoid object
[az,elev,slantRange] = geodetic2aer(lat_geosat,lon_geosat,height_geosat, ...
lat_receiver,lon_receiver,height_receiver,spheroid)
Have a look at the below documentation page for better understanding
I hope this helps.
