%--------------------直角轉經緯度-------------------%
clear
clc

x=input（‘Please?input?the?parameter?（not?zero）:?x?=?‘）;
y=input（‘Please?input?the?parameter?（not?zero）:?y?=?‘）;
z=input（‘Please?input?the?parameter?（not?zero）:?z?=?‘）;

a=6378137;
e=1/298.257;
b=a-a*e;
e=sqrt（a*a-b*b）/a;
R=sqrt（x*x+y*y+z*z）;
B=atan（z/sqrt（x*x+y*y））;

%--------------------以下為迭代求fi-----------------%
fi1=acos（sqrt（x*x+y*y）/R）;
buchang=0.1/60/60/180*pi;
w=sqrt（1-e*e*sin（fi1）*sin（fi1））;
fi2=atan（tan（B）*（1+a*e*e/z*sin（fi1）/w））;
cha=abs（fi1-fi2）;
xuanze=（fi1switch?xuanze
????case?1
????????while（buchang>=0.0000001/60/60/180*pi）
????????????while（abs（fi1-fi2）<=cha）
????????????????cha=abs（fi1-fi2）;
????????????????fi1=fi1+buchang;
????????????????w=sqrt（1-e*e*sin（fi1）*sin（fi1））;
????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3301??2008-12-12?09:37??test2.m

?????文件???????1950??2008-11-12?20:10??xyztojwd.m

?????文件?????252928??2008-12-14?22:38??實驗報告.doc

?????文件?????302608??2006-11-07?21:15??GPS實驗1\GGA.gga

?????文件????????192??2008-12-11?22:22??GPS實驗1\GPS實驗1要求.txt

?????文件????????219??2006-11-29?10:01??GPS實驗2-DGPS\GPS實驗2要求.txt

?????文件???????2111??2008-10-27?22:16??assignment1.m

?????文件????????301??2008-10-21?19:33??assignment2.m

?????文件???????3241??2008-12-04?10:20??gps_position.m

?????文件???????2761??2008-12-11?23:20??test1.m

?????目錄??????????0??2009-04-19?13:18??GPS實驗2-DGPS\20061128-DGPS

?????目錄??????????0??2009-04-19?13:18??GPS實驗1

?????目錄??????????0??2009-04-19?13:18??GPS實驗2-DGPS

-----------?---------??----------?-----??----

???????????????569612????????????????????13