?? dirprob.m
字號:
echo on%--------------------------------------------------% DirProb% Example direct geodetic problem. Computes% lat,lon,ht of 2nd point from distance,azimuth,% vertical angle from 1st point.% 7 Jul 96%% M-files: dms2rad, rad2dms, inverse, direct,% refell%--------------------------------------------------clear%---------- Define lat,lon,ht of 1st (from) point(s)lat1=dms2rad([50 15 31.67214]); % Enter latitudelon1=dms2rad([-95 51 58.18951]); % Enter longitudeh1=248.39; % Enter ell. heightn=max(size(lat1));%---------- Define az,va,dist to 2nd point(s)azdms=[100 30 10]; % Enter azimuthvadms=[10 45 05]; % Enter vertical angled=1500; % Enter distanceaz=dms2rad(azdms);va=dms2rad(vadms);%---------- Compute lat,lon,ht of 2nd point(s)[a,b,e2,finv]=refell('NAD83');[lat2,lon2,h2]=direct(lat1,lon1,h1,az,va,d,a,e2);latdms=rad2dms(lat2);londms=rad2dms(lon2);%---------- List resultsfor i=1:n fprintf('\nPoint %4.0f\n',i); fprintf('Latitude %4.0f %2.0f %9.6f\n',latdms(i,1),latdms(i,2),latdms(i,3)); fprintf('Longitude %4.0f %2.0f %9.6f\n',londms(i,1),londms(i,2),londms(i,3)); fprintf('Height %9.4f\n',h2(i));end
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -