GPS卫星坐标计算

1、卫星定位技术与方法-根据广播星历参数计算卫星坐标 作业报告指导教师：能 八、永良班 级：一测%经二班学生姓名：一东学生学号： 2 0 0 8 0 7 8 3作业日期：2010年12月08日错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。错误!未定义书签。已知数据 计算步骤 平均角速度(mean angular speed) 规化时刻 (normal time) 平近点角 (mean anomaly) 偏近点角(eccentric anomaly) 真近点角(true anomaly

2、) 升交距角(argument of ascending node)轨道向径( Orbital radius ) 扰动改正( Perturbed correction ) 10. 卫星在升交点轨道直角坐标系中的坐标错误!未定义书签。11. 升交点经度( Longitude of ascending node ) 错误 !未定义书 签。源程序错误!未定义书签。程序运行结果错误!未定义书签。七 . 作业体会错误!未定义书签。根据广播星历参数计算卫星坐标一.已知数据：根据以下的广播星历参数计算 UTC200侔1月30日8点0分00秒一20分00秒，每隔一分钟的PRN7勺卫星坐标。Computethe

3、 coordinate of PRN7 with interval of 1 minute.Navigation data:之NAVITIQII MTARIREX VERSION / TEPEEFHTOt海23-L 版本30-JAN-04 碌;41 PCM / RUN BY / DATEHMD OF HEADEB4 i i sc， 3 o a a- or 5/ ：as 1332157D-04-0. saz 17 55 it? 23c-11 a. oa maoa oaaoo choo D. 13SDOOD0OOQ(nm)3 QB6先町50叩0加02 0.14755430 13TD-O0 (J.

4、100130&339()1 D, 3&5W号g硒BD05 Q 时的26聊州6alM)2 Q. 5Q53Q9995W-05 口口1535TB空门八附！ 口，4b口Booooaoomm)曰 o.sieeosFsaegc-aT o.305054077as9tHJL a.slsesss65662d-ot Q. 96172S74996H)O 0. 2794fi6750000(I+03-0. lOS70ia71CHM-0. 3 L8S 19831)-08 D, 6O83921ia0D-U 0. lODOOOOOOOOOMJl 0. L255MOOOQD+00-0. 0569676000-00 G 1390

5、0QQQQMQM4S D. QEXJDOQEXIQOODDHJQ7 4 1 a0 9 0 Oi 0 Q62932S3&7313M3-0 54K96鸵1064D-11 QQOQQWQIXKKKJtHOO D,210WKX)0(KDM)3 Q 5肝312印用0财02 Q. 4dKJ05S部77mgF. WL49924两/？眄1 0. 2754B5ZL7571-0S O. 125S56591351D-01 O. 90375542640715-05 (J. 5 L5SF3Q 38319(04 口, 4网瞅M妙计用。.额67限611删&HK 0.193198315097MI-0.726431606200

6、0-07 D. 93?曲0945网口 0.19ia7500000a(I+03-0. L85137234909DKIL-0. 7SlflS109M2IH)8 -D. 13U34M4204D-09 0. lOOOOOOOOOOaCHOl 0.125EOOOOOOOOOHM 0. 0000900000001X)0 D, ODODOOOOOWOQHJO 0. OOOOOOOOOODOJWQ-0. LaG2M51W3M8 0.210000000X01X)3 D. OQOOCOO3QOOCiIWQ卫星导航文件格式:PRN弱比a结H M skit (iicMuij)工：；父也比0登:财：3K)IX出比小）r

7、 radian 曲M/jsifirs)弧侬艰援加回：:运娘d sec cfGPStek篇孰你孙唠）四 raiersjiOfiirfinsj|羸L工日蚂cn=?i Ci&nii)OMGA DC 7 (dEi5 式；HXJT(ridiam senCodfi ctt L2 chsnnriGPSAVk-LiP Jitafla?SV anzrzcy .底:gr/htlrETGDfirafc33C cfDi：iTntfm：tcondisparespar*spine计算步骤：The steps for satellite coordinates.平均角速度 (mean angular speed)nn0n由

8、广播星历获得，GM=+14.规化时刻(normal time)tk t tot0已知(由广播星历获得)，t为GPSW秒.平近点角(mean anomaly):Mk Mo n tkM0已知(由广播星历获得).偏近,电角(eccentric anomaly):Ek Mk e sin Ek迭代求解：初始值取E=M以弧度为单位.真近点角(true anomaly):1 e2 sin EkVkarctancosEk e.升交品巨角(argument of ascending node)近地点角品巨(argument of perigee).轨道向径(Orbital radius ):rka (1 eco

9、s Ek). 扰动改正(Perturbed correction ):?升交角距(Argument of ascending node?轨道向径u= CuC cos2CuS sin 2 0(Orbital radius )r C rC CoS 2CrS sin 2 0?轨道顷角(Orbital inclination )i C iC cos 2CiS sin 2.改正后升交角距、轨道向径、轨道倾角改正后升交角品巨( Corrected argument of ascending nodu k 0 u改正后的轨道向径(Corrected orbital radius)rk a (1 e cosE

10、k)r改正后的轨道倾角(Corrected orbital inclination )ik i0 i (IDOT)tk10.卫星在开交点轨道直角坐标系中的坐标：如下图所示xk rkcosukyk rksinuk.升交点经度(Longitude of ascending node ):如下图所示t 0 (e)t ket 0e_ -5e 7.2921151467 10 rad/s.在地固坐标系中的卫星位置(Expressed in spheric coordinate system)k xk cos k yk cosik sin k Yk xk sin k yk cosi k cos k Zk y

11、ksinik三.源程序：using System;using ;using ;using ;using ;using ;using Testpublic partial class Form1 : Formpublic Form1()InitializeComponent();序号, 40);星历内容,130);导航数据,130);= true; oString();Mi = (arrsi);double t0e=;double t;double XK = new double21;double YK = new double21;double ZK = new double21;for (i

12、nt l = 0; l 21;l+ )ukukukoString();四.程序运行结果：星历内容的读取、显示:4显示卫星星历文件内容:130日口 口4 5 6 123455T89O1-23 1 Qw 1- 1- 1- 1- JJ 1- 1- JJ _2 _2_序号星力内容导航数据FYMdHM Ll 稿 a a 工 CJdMLJ wc sic Q c i- 545696321。联eT 1 O.OOOOOOOOOOOOeWO O.ElOOOOOOOOOOeWS 0.5375125 OOOCOH)2 0 第53口5923T7X-0B0. 2149475672e+01 0.275485217571e-

13、05 O,1Z5U5S6913514-01 0.903755426407-05 0 51537303981 5+4 0 q&口如WOQUQ口*打6 D.E06753611565-06 0.193193315tNTl -C.T264316032 OOe-07 0 8376s0094门 S刀口卫星坐标计算结果:讦SOM里粽(in rii iiirfin! mriin llll，F，i-时间支坐标建标第咐钥1-13113645.4&3S0051609ZCE7.2&E309ZieSZ60Z4.39269T&第1分钥1-13W56BT,T53557110556,南1303165(ng 31Kx4第渤钥1

14、T3557口16,270341316329,1953Z2216385603.3302314第3分钥1-13力7在乙5T 6556916H6S96.633TOOT16 烟 543.7533164第吗钥1-13396496.3567713IB1&5239.29E3028161302&3,475666&第5分钥1-1014335.424OCK16163EL7,46401116012771.1636663第凶朝1-14132225,721007516202471.5293033l5S8607B.75&04第俏智1-1 1243651.321409416221241.635Z33T15756138.32

15、09013第吩钓1-14X4134.430894716240116.579861915 陷 116.4 10132第涉智1-1447363T,333356G16259091.9QEB50Z154 兆 B69.35 比 364第 1 口分钟 -1459216. 8&7口子曰16278151,6&923S515357457.1625691第 1L 分钟 -14704702.3T5622216297236 47124T3152348S8.9691003第1工分钟 -14S葭印.彩934491631&49J.6500163 117% &114355第 13分钟 -14926799.201Q522163

16、35751.137534914966323.8101145第 14分钟 -15036341. EE2E51216353056.S561E9T14B303银41总咨第 15 分钟 -15144869,931133163796.63057614日93252.羽 6792第 16 分钟 -15E523T6.1rB35 国 9163937S5.6&9765514555051.472E第17分钟 T5358朋1,592。1516413146.2&85T021441s7弭 02 笔 543第归分钟 -151&4310,4211616432535.359609914ZT53T0,1200T4B第1日分钟-1

17、5%仃级电721，16451-9L6.61T5113KL33910.02=36661第匆分钟 -15&T20S5,4436巩816411231.4580113991334. J7S2333运行界面:品聚II早星丽立丽西1磅里是逑TME CWECtjs s:ar t (b.1 心 CicCi iO星l：W脾 旧 Yer Mtn.d印1 HE c 口弓 2 点Al诙racrCL。Q. 战把SMW 鸵 3U -(1 写 4cin4e-ii C . OOOCOCCCDMD*a O.LODOmCOMD r.53mi 25rtm 4853JS厂TE 一0. 4m3rG /在帕 i n.sm如 T*ia Q

18、JHH期欣Id o n3refi4HTflr6 C153rC098nBKkl C.460aXCC0MDeKfl 0 206?53615*5#6 0. 193198315W7cKil F25431 的道口 5T I WdLLHl咖司观标港而画里巾肿7 Ml初5 U39jO516092527,26232l&6BDB4.69aS1B第方钟-335307 7535271&1LQS5E 251,383lfi50T433. 3303944费好肿-13SSUC1S_27DB4131&12B&W. 3216JSB&J3.3302314笔也神-13TTO22 5TTT91614WKC 6937007102ft2

19、543.75331M第一钟-J第I&W目强STTU151后号酬.工艇3目虎质】脂狗再卬钟5第美即钟7中虱435铝在陶16163017.44011IGOt ETTl.锢聊6革自上曾-14132225 72IK75162024?!. 5M3SBIsaasnre-. fiirs&n4豺I丑肿7m:：段1.切叫16S2124L 83ES39?iST&Biae.eamDia串竹科-141M134 43甘通仃1fii?4rn w 6Tqmig139110.99tQl3e第9分制1-J417SE.1T-3B535S51K59OS1. gCEBSJZ15493S9. 5563&1强0加-J459S1G2.66709l&7B15L06Se3eG153674S?. 165B91第由钟-14T04702 8T982221K9TOT 4TP4T.3l5254W,mcD0B芾、冽r-HJ1625Q 75J331&3LB4Q1.&3V048SlElOim.GlL