CN115032670B - 一种纠正gps观测得到的分潮误差的方法及设备 - Google Patents

一种纠正gps观测得到的分潮误差的方法及设备 Download PDF

Info

Publication number
CN115032670B
CN115032670B CN202210953013.5A CN202210953013A CN115032670B CN 115032670 B CN115032670 B CN 115032670B CN 202210953013 A CN202210953013 A CN 202210953013A CN 115032670 B CN115032670 B CN 115032670B
Authority
CN
China
Prior art keywords
tide
point
gps
amplitude
correcting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210953013.5A
Other languages
English (en)
Other versions
CN115032670A (zh
Inventor
徐晓庆
魏泽勋
王永刚
李淑江
潘海东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
First Institute of Oceanography MNR
Original Assignee
First Institute of Oceanography MNR
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by First Institute of Oceanography MNR filed Critical First Institute of Oceanography MNR
Priority to CN202211260841.7A priority Critical patent/CN115616630A/zh
Priority to CN202210953013.5A priority patent/CN115032670B/zh
Publication of CN115032670A publication Critical patent/CN115032670A/zh
Application granted granted Critical
Publication of CN115032670B publication Critical patent/CN115032670B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/40—Correcting position, velocity or attitude
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明属于海洋监测领域，提供了一种纠正GPS观测得到的分潮误差的方法及设备，方法包括以下步骤：确定分潮对为确定N2对M2、S2对M2、K2对M2、K2对S2、P1对K1、O1对K1、Q1对K1和Q1对O1的分潮对；建立各分潮对的相对关系，获取各个网格点的振幅比及迟角差；计算平衡潮调和常数，得到各分潮经地潮校正后的平衡潮振幅；计算GPS数据调和常数，获取由GPS数据调和分析得到的各分潮的振幅；根据步骤2、步骤3和步骤4的结果并利用相对导纳纠正K1和K2分潮估值误差。本发明能够有效改善GPS观测得到的K1和K2分潮的误差并提高K1和K2分潮的准确度，为提高研究海潮负荷效应的准确度提供一种新的方法。

Description

一种纠正GPS观测得到的分潮误差的方法及设备
技术领域
本发明属于海洋监测技术领域，尤其涉及一种纠正GPS观测得到的分潮误差的方法及设备。
背景技术
随着GNSS技术与数据处理方法的快速发展，可以通过GNSS准确地获取观测站位由海潮负荷引起的位移形变，并进而得到海潮负荷位移调和常数，为研究海潮负荷效应提供一种新的方法。一般在利用GNSS研究海潮负荷位移时，主要使用GPS观测数据，对GPS观测数据进行调和分析时需要提前估算主要分潮的周期，而GPS卫星轨道的周期（约11.967小时）与K2分潮的周期（约11.9672小时）十分相近，GPS星座的重复周期（约23.9319小时）与K1分潮的周期（约23.9345小时）非常相近，容易产生耦合现象，在利用GNSS技术提取海潮负荷位移参数时，容易造成周期相近的误差会被对应相同周期的K1/K2分潮吸收，导致调和分析估算的K1/K2分潮误差较大。因此如何解决K1/K2分潮估值误差过大一直是GNSS研究海潮负荷位移的难点和热点问题。
发明内容
本发明利用S2对M2，N2对M2，K2对S2，O1对K1，P1对K1和Q1对O1分潮的相对导纳纠正GPS观测得到K1和K2分潮估值误差过大问题。
本发明第一方面提供了一种纠正GPS观测得到的K1和K2分潮误差的方法，包括以下步骤：
步骤1，确定分潮对；所述确定分潮对为确定N2对M2、S2对M2、K2对M2、K2对S2、P1对K1、O1对K1、Q1对K1和Q1对O1的分潮对；
步骤2，建立各分潮对的相对关系，获取各个网格点的振幅比
Figure 654690DEST_PATH_IMAGE001
及迟角差
Figure 368568DEST_PATH_IMAGE002
；
步骤3，计算平衡潮调和常数，得到各分潮经地潮校正后的平衡潮振幅
Figure 210622DEST_PATH_IMAGE003
；
步骤4，计算GPS数据调和常数，获取由GPS数据调和分析得到的各分潮的振幅
Figure 904034DEST_PATH_IMAGE004
；
步骤5，根据步骤2、步骤3和步骤4的结果并利用相对导纳纠正K1和K2分潮估值误差；具体为：
令
Figure 515144DEST_PATH_IMAGE005
式中，
Figure 399923DEST_PATH_IMAGE004
为由GPS数据调和分析得到的各分潮的振幅，
Figure 463694DEST_PATH_IMAGE003
为各分潮平衡潮振幅，
Figure 459332DEST_PATH_IMAGE006
为模数，将各半日分潮N2、M2、S2和K2的
Figure 423484DEST_PATH_IMAGE007
值进行曲线拟合，然后将K2分潮的异常
Figure 479164DEST_PATH_IMAGE007
值纠正到正常的拟合曲线上，从而完成K2分潮误差的纠正；同理，将各全日分潮K1、P1、O1和Q1的
Figure 30231DEST_PATH_IMAGE007
值进行曲线拟合，然后将K1分潮的异常
Figure 298402DEST_PATH_IMAGE007
值纠正到正常的拟合曲线上，从而完成K1分潮误差的纠正。
在一种可能的设计中，所述步骤1具体过程为：
对某一分潮i，导纳Mi为
Figure 884104DEST_PATH_IMAGE008
（1）
其中
Figure 612151DEST_PATH_IMAGE009
，
Figure 650514DEST_PATH_IMAGE010
和
Figure 722375DEST_PATH_IMAGE011
为分潮的调和常数，振幅和迟角，
Figure 162584DEST_PATH_IMAGE012
为引潮力系数，在同一潮族中，主要分潮n相对于次要分潮m的关系可表示为
Figure 560067DEST_PATH_IMAGE013
（2）
其中
Figure 590121DEST_PATH_IMAGE014
（3）
Figure 465673DEST_PATH_IMAGE015
（4）
式（3）中
Figure 760388DEST_PATH_IMAGE016
即是传统潮汐学中的振幅比，
Figure 328772DEST_PATH_IMAGE017
为迟角差，
Figure 76149DEST_PATH_IMAGE018
为主要分潮n相对于次要分潮m的振幅关系，
Figure 256856DEST_PATH_IMAGE019
为主要分潮n相对于次要分潮m的迟角关系；如果实际分潮的振幅比与引潮力之比相同，则
Figure 406078DEST_PATH_IMAGE020
等于1；实际海洋中
Figure 145364DEST_PATH_IMAGE020
不会正好等于1，但对大部分海域它接近于1；如果两个分潮的迟角相等，则
Figure 380036DEST_PATH_IMAGE021
等于零；实际海洋中
Figure 597391DEST_PATH_IMAGE021
一般不等于零，对大部分海域它大于零或小于零；
确定N2对M2、S2对M2、K2对M2、K2对S2、P1对K1、O1对K1、Q1对K1和Q1对O1分潮对；
Figure 630813DEST_PATH_IMAGE022
（5）
Figure 541000DEST_PATH_IMAGE023
（6）
令振幅比分别为
Figure 262968DEST_PATH_IMAGE024
迟角差分别为
Figure 284014DEST_PATH_IMAGE025
（7）
对某两个分潮，设振幅比为
Figure 407828DEST_PATH_IMAGE026
，迟角差为
Figure 223337DEST_PATH_IMAGE027
,记
Figure 934066DEST_PATH_IMAGE028
，
Figure 758803DEST_PATH_IMAGE029
；
所述步骤2的具体过程为：
已知有N个观测点，分别用
Figure 2702DEST_PATH_IMAGE030
代表，各点的经纬度为
Figure 254692DEST_PATH_IMAGE031
，其中
Figure 709508DEST_PATH_IMAGE032
对南纬用负值表示；
Figure 603514DEST_PATH_IMAGE033
，
Figure 436341DEST_PATH_IMAGE032
以度为单位；
Figure 593653DEST_PATH_IMAGE010
和
Figure 43089DEST_PATH_IMAGE011
为各验潮站观测点的振幅比和迟角差；
对于观测站
Figure 711093DEST_PATH_IMAGE034
点，计算它对其它各站点j的距离为
Figure 788639DEST_PATH_IMAGE035
（8）
求
Figure 208863DEST_PATH_IMAGE034
点的加权平均值
Figure 614436DEST_PATH_IMAGE036
（9）
Figure 584666DEST_PATH_IMAGE037
（10）
其中，
Figure 392085DEST_PATH_IMAGE038
代表对所有观测点求和，但不包括
Figure 658244DEST_PATH_IMAGE034
点；
Figure 551113DEST_PATH_IMAGE039
是修正后的开尔文波波长，其根据
Figure 59455DEST_PATH_IMAGE034
点附近的水深确定的，方法如下：
设已知海区内及周围各网格点的水深值，即对点
Figure 986960DEST_PATH_IMAGE040
其经纬度为
Figure 656976DEST_PATH_IMAGE041
，
Figure 541536DEST_PATH_IMAGE042
,水深为Dm；对
Figure 587989DEST_PATH_IMAGE034
点，计算网格点
Figure 635580DEST_PATH_IMAGE040
与验潮站点
Figure 476497DEST_PATH_IMAGE034
的距离为
Figure 343959DEST_PATH_IMAGE043
（11）
选取所有
Figure 961147DEST_PATH_IMAGE044
的点，对其水深值进行平均，记为
Figure 332086DEST_PATH_IMAGE045
；
计算相应的开尔文波波速
Figure 875062DEST_PATH_IMAGE046
Figure 964241DEST_PATH_IMAGE047
（12）
开尔文波波长
Figure 618076DEST_PATH_IMAGE048
（13）
其中，ui为开尔文波波速，T为潮波周期，即
Figure 873215DEST_PATH_IMAGE049
（14）
Figure 55935DEST_PATH_IMAGE050
是潮波角速率；
最后选取
Figure 897989DEST_PATH_IMAGE051
Figure 355515DEST_PATH_IMAGE052
（15）
对所有
Figure 701045DEST_PATH_IMAGE053
各点算出
Figure 87290DEST_PATH_IMAGE054
和
Figure 885481DEST_PATH_IMAGE055
，并计算
Figure 146698DEST_PATH_IMAGE056
（16）
及均方根偏差
Figure 346736DEST_PATH_IMAGE057
（17）
若
Figure 667996DEST_PATH_IMAGE058
（18）
则认为第
Figure 953483DEST_PATH_IMAGE034
点数据可能异常，考虑是否舍去；
为验证去除的异常值后的精确度，由
Figure 268925DEST_PATH_IMAGE054
和
Figure 589048DEST_PATH_IMAGE055
反算各观测点的振幅比
Figure 815630DEST_PATH_IMAGE059
及迟角差
Figure 853993DEST_PATH_IMAGE060
Figure 958478DEST_PATH_IMAGE061
（19）
Figure 398687DEST_PATH_IMAGE062
Figure 796170DEST_PATH_IMAGE063
（20）
再由
Figure 321829DEST_PATH_IMAGE064
（21）
Figure 197381DEST_PATH_IMAGE065
（22）
计算振幅比与迟角差的均方差
Figure 256211DEST_PATH_IMAGE066
（23）
Figure 559016DEST_PATH_IMAGE067
（24）
按上述方法先舍弃不合理的观测值，按预先给定的分辨率设定各网格点，记各网格点的序号为
Figure 571971DEST_PATH_IMAGE068
，其对应的经纬度为
Figure 251215DEST_PATH_IMAGE069
，计算点
Figure 400436DEST_PATH_IMAGE070
到各观测点
Figure 641187DEST_PATH_IMAGE034
的距离为
Figure 875859DEST_PATH_IMAGE071
（25）
求
Figure 358793DEST_PATH_IMAGE070
点的加权平均值
Figure 628100DEST_PATH_IMAGE072
（26）
Figure 272708DEST_PATH_IMAGE073
（27）
其中
Figure 499071DEST_PATH_IMAGE074
计算方法与前述方法相同，
由
Figure 520117DEST_PATH_IMAGE075
和
Figure 909510DEST_PATH_IMAGE076
反算各个网格点的振幅比
Figure 990598DEST_PATH_IMAGE001
及迟角差
Figure 465442DEST_PATH_IMAGE002
；
Figure 791643DEST_PATH_IMAGE077
（28）
Figure 769964DEST_PATH_IMAGE078
Figure 21953DEST_PATH_IMAGE079
（29）
在一种可能的设计中，所述步骤3的具体过程为：
对任意观测点
Figure 718514DEST_PATH_IMAGE080
，其对应的平衡潮潮高
Figure 845477DEST_PATH_IMAGE081
可表示为
Figure 678303DEST_PATH_IMAGE082
（30）
式中，
Figure 101195DEST_PATH_IMAGE083
和
Figure 285051DEST_PATH_IMAGE084
为分潮的交点因子和交点订正角，
Figure 186011DEST_PATH_IMAGE085
为平衡潮第
Figure 905968DEST_PATH_IMAGE034
个分潮在
Figure 499760DEST_PATH_IMAGE086
时刻的位相；
Figure 905334DEST_PATH_IMAGE087
是分潮角速率；
Figure 875564DEST_PATH_IMAGE088
为族号，对全日潮
Figure 948562DEST_PATH_IMAGE089
，对半日潮
Figure 205931DEST_PATH_IMAGE090
；
Figure 364380DEST_PATH_IMAGE033
为经度；
Figure 607143DEST_PATH_IMAGE091
为北京标准时区号；
Figure 534647DEST_PATH_IMAGE092
为各分潮经地潮校正后的平衡潮振幅，
Figure 735822DEST_PATH_IMAGE093
Figure 351873DEST_PATH_IMAGE094
Figure 663906DEST_PATH_IMAGE095
Figure 711496DEST_PATH_IMAGE096
Figure 817992DEST_PATH_IMAGE097
Figure 685454DEST_PATH_IMAGE098
Figure 299713DEST_PATH_IMAGE099
Figure 936231DEST_PATH_IMAGE100
（31）
式中，
Figure 948049DEST_PATH_IMAGE032
为纬度。
在一种可能的设计中，所述步骤4的具体过程为：
根据潮汐理论，假设
Figure 568386DEST_PATH_IMAGE101
为海潮负荷对观测站中某一个方向的瞬时位移，则主要分潮的调和常数可表示为，
Figure 222221DEST_PATH_IMAGE102
（32）
式中，
Figure 214710DEST_PATH_IMAGE103
为位移的平均位置；
Figure 928589DEST_PATH_IMAGE026
和
Figure 505063DEST_PATH_IMAGE027
为各分潮的调和常数；
Figure 697010DEST_PATH_IMAGE083
和
Figure 308120DEST_PATH_IMAGE084
为分潮的交点因子和交点订正角，
Figure 697294DEST_PATH_IMAGE085
为平衡潮第
Figure 761065DEST_PATH_IMAGE034
个分潮在
Figure 491124DEST_PATH_IMAGE086
时刻的位相；
Figure 956740DEST_PATH_IMAGE087
是分潮角速率；这里
Figure 278000DEST_PATH_IMAGE104
、
Figure 330532DEST_PATH_IMAGE105
及
Figure 129861DEST_PATH_IMAGE106
与时间有关，与地点无关；而
Figure 449984DEST_PATH_IMAGE026
和
Figure 410986DEST_PATH_IMAGE027
与地点有关，与时间无关；故上式可写作
Figure 449349DEST_PATH_IMAGE107
(33)
其中，
Figure 816483DEST_PATH_IMAGE108
(34)
Figure 991113DEST_PATH_IMAGE109
(35)
按最小二乘法原理，解算出
Figure 388596DEST_PATH_IMAGE110
和
Figure 648676DEST_PATH_IMAGE111
，再按照下式计算出调和常数
Figure 291272DEST_PATH_IMAGE112
和
Figure 851567DEST_PATH_IMAGE113
，
Figure 888793DEST_PATH_IMAGE114
(36)
本发明第二方面还提供了一种纠正GPS观测得到的分潮误差的设备，所述设备包括至少一个处理器和至少一个存储器，所述处理器和存储器相耦合；所述存储器中存储有计算机程序或指令；所述处理器执行所述计算机程序或指令时，可以实现如第一方面所述的纠正GPS观测得到的K1和K2分潮误差的方法。
本发明第三方面还提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序或指令，所述计算机程序或指令被处理器执行时，可以实现如第一方面所述的纠正GPS观测得到的K1和K2分潮误差的方法。
与现有技术相比，本发明提供了一种纠正GPS观测得到的2分潮误差的方法及设备，并产生如下的有益效果：
1.针对GPS轨道周期与K2分潮周期一致、GPS星座的重复周期与K1分潮的周期一致从而产生耦合作用，导致K2分潮和K1分潮估值误差过大的难题。本发明提出了一种基于相对导纳纠正GPS观测得到的K1和K2分潮误差的方法及设备，能够有效改善GPS观测得到的K1和K2分潮的误差并提高K1和K2分潮的准确度。
2.通过GNSS准确地获取观测站位由海潮负荷引起的位移形变，并进而得到海潮负荷位移调和常数，本发明为提高研究海潮负荷效应的准确度提供一种新的方法。
3.为提高研究海潮负荷效应的准确度，本发明采用建立8个分潮对相对关系的方法，能有效得到各潮族的相对导纳关系。结合本发明采用的GPS数据分析调和常数与相应平衡潮调和常数的比值，能有效的纠正K1和K2分潮的误差并提高K1和K2分潮的准确度。
附图说明
图1为本发明纠正GPS观测得到的K1和K2分潮误差方法的流程框图。
图2为本发明纠正K2分潮误差的示意图。
图3为本发明纠正GPS观测得到的K1和K2分潮误差的设备结构简易示意图。
具体实施方式
下面结合具体实施例对发明进行进一步说明。
实施例1：
为实现基于相对导纳纠正GPS观测得到的K1和K2分潮误差的推算技术，其包括相对导纳计算和纠正GPS观测得到的K1和K2分潮误差的计算步骤，大体流程如图1所示。
（一）确定分潮对
对某一分潮i，导纳Mi为
Figure 901748DEST_PATH_IMAGE008
（1）
其中
Figure 846570DEST_PATH_IMAGE009
；
Figure 730213DEST_PATH_IMAGE010
和
Figure 985612DEST_PATH_IMAGE011
为分潮的调和常数，振幅和迟角；
Figure 220284DEST_PATH_IMAGE012
为引潮力系数。在同一潮族中，主要分潮n相对于次要分潮m的关系可表示为
Figure 703218DEST_PATH_IMAGE013
（2）
其中
Figure 706946DEST_PATH_IMAGE014
（3）
Figure 617133DEST_PATH_IMAGE015
（4）
式（3）中
Figure 840567DEST_PATH_IMAGE016
即是传统潮汐学中的振幅比，
Figure 127192DEST_PATH_IMAGE017
为迟角差。如果实际分潮的振幅比与引潮力之比相同，则
Figure 251005DEST_PATH_IMAGE020
等于1。实际海洋中
Figure 332094DEST_PATH_IMAGE020
不会正好等于1，但对大部分海域它接近于1。如果两个分潮的迟角相等，则
Figure 275779DEST_PATH_IMAGE021
等于零。实际海洋中
Figure 864630DEST_PATH_IMAGE021
一般不等于零，对大部分海域它大于零或小于零。
确定N2对M2，S2对M2，K2对M2，K2对S2，P1对K1，O1对K1，Q1对K1和Q1对O1分潮对。
Figure 842950DEST_PATH_IMAGE022
（5）
Figure 94940DEST_PATH_IMAGE023
（6）
令振幅比分别为
Figure 57080DEST_PATH_IMAGE024
迟角差分别为
Figure 154349DEST_PATH_IMAGE025
（7）
对某两个分潮，设振幅比为
Figure 488640DEST_PATH_IMAGE026
，迟角差为
Figure 177111DEST_PATH_IMAGE027
，记
Figure 360967DEST_PATH_IMAGE028
，
Figure 527507DEST_PATH_IMAGE029
。
（二）建立各分潮对的相对关系
已知有N个观测点，分别用
Figure 984814DEST_PATH_IMAGE030
代表，各点的经纬度为
Figure 578606DEST_PATH_IMAGE031
，其中
Figure 984180DEST_PATH_IMAGE032
对南纬用负值表示。
Figure 688830DEST_PATH_IMAGE033
，
Figure 761829DEST_PATH_IMAGE032
以度为单位。
Figure 260943DEST_PATH_IMAGE010
和
Figure 655278DEST_PATH_IMAGE011
为各验潮站观测点的振幅比和迟角差。
对于观测站
Figure 429199DEST_PATH_IMAGE034
点，计算它对其它各站点j的距离为
Figure 356703DEST_PATH_IMAGE035
（8）
求
Figure 557877DEST_PATH_IMAGE034
点的加权平均值
Figure 170999DEST_PATH_IMAGE036
（9）
Figure 483032DEST_PATH_IMAGE037
（10）
其中，
Figure 265043DEST_PATH_IMAGE038
代表对所有观测点求和，但不包括
Figure 371539DEST_PATH_IMAGE034
点。
Figure 504580DEST_PATH_IMAGE039
是根据
Figure 590610DEST_PATH_IMAGE034
点附近的水深确定的，方法如下：
设已知海区内及周围各网格点的水深值，即对点
Figure 227128DEST_PATH_IMAGE040
其经纬度为
Figure 770105DEST_PATH_IMAGE041
，
Figure 390442DEST_PATH_IMAGE042
，水深为Dm。对
Figure 44277DEST_PATH_IMAGE034
点，计算网格点
Figure 559136DEST_PATH_IMAGE040
与验潮站点
Figure 741855DEST_PATH_IMAGE034
的距离为
Figure 318330DEST_PATH_IMAGE115
（11）
选取所有
Figure 775856DEST_PATH_IMAGE044
的点，对这些点的水深值进行平均，记为
Figure 154010DEST_PATH_IMAGE045
。
计算相应的开尔文波波速
Figure 773210DEST_PATH_IMAGE046
Figure 836981DEST_PATH_IMAGE047
（12）
开尔文波波长
Figure 832619DEST_PATH_IMAGE048
（13）
其中，T为潮波周期，即
Figure 298236DEST_PATH_IMAGE049
（14）
这里
Figure 88337DEST_PATH_IMAGE050
是潮波角速率。
最后选取
Figure 137939DEST_PATH_IMAGE116
Figure 202847DEST_PATH_IMAGE052
（15）
对所有
Figure 788549DEST_PATH_IMAGE053
各点算出
Figure 313334DEST_PATH_IMAGE054
和
Figure 351697DEST_PATH_IMAGE055
，并计算
Figure 423558DEST_PATH_IMAGE056
（16）
及均方根偏差
Figure 598188DEST_PATH_IMAGE057
（17）
若
Figure 261250DEST_PATH_IMAGE058
（18）
则认为第
Figure 25725DEST_PATH_IMAGE034
点数据可能异常，考虑是否舍去，这里
Figure 166856DEST_PATH_IMAGE040
可取2。
为验证去除的异常值后的精确度，由
Figure 930413DEST_PATH_IMAGE054
和
Figure 764377DEST_PATH_IMAGE055
反算各观测点的振幅比
Figure 511753DEST_PATH_IMAGE059
及迟角差
Figure 190996DEST_PATH_IMAGE060
Figure 841682DEST_PATH_IMAGE061
（19）
Figure 315389DEST_PATH_IMAGE062
Figure 815640DEST_PATH_IMAGE063
（20）
再由
Figure 298574DEST_PATH_IMAGE064
（21）
Figure 302302DEST_PATH_IMAGE065
（22）
计算振幅比与迟角差的均方差
Figure 711025DEST_PATH_IMAGE066
（23）
Figure 698572DEST_PATH_IMAGE067
（24）
按上述方法先舍弃不合理的观测值，按预先给定的分辨率设定各网格点，记各网格点的序号为
Figure 719618DEST_PATH_IMAGE068
，其对应的经纬度为
Figure 577852DEST_PATH_IMAGE069
，计算点
Figure 658941DEST_PATH_IMAGE070
到各观测点
Figure 369670DEST_PATH_IMAGE034
的距离为
Figure 194407DEST_PATH_IMAGE117
（25）
求
Figure 172727DEST_PATH_IMAGE070
点的加权平均值
Figure 424717DEST_PATH_IMAGE072
（26）
Figure 590119DEST_PATH_IMAGE073
（27）
Figure 239DEST_PATH_IMAGE074
计算方法与前述方法类似。
由
Figure 98645DEST_PATH_IMAGE075
和
Figure 255957DEST_PATH_IMAGE076
反算各个网格点的振幅比
Figure 174234DEST_PATH_IMAGE001
及迟角差
Figure 606353DEST_PATH_IMAGE002
。
Figure 60730DEST_PATH_IMAGE077
（28）
Figure 654522DEST_PATH_IMAGE078
Figure 60096DEST_PATH_IMAGE079
（29）
（三）计算平衡潮调和常数
对任意观测点
Figure 764747DEST_PATH_IMAGE080
，其对应的平衡潮潮高
Figure 837745DEST_PATH_IMAGE081
可表示为，
Figure 100974DEST_PATH_IMAGE082
（30）
式中，
Figure 993843DEST_PATH_IMAGE083
和
Figure 236606DEST_PATH_IMAGE084
为分潮的交点因子和交点订正角，
Figure 164111DEST_PATH_IMAGE085
为平衡潮第
Figure 99706DEST_PATH_IMAGE034
个分潮在
Figure 479871DEST_PATH_IMAGE086
时刻的位相；
Figure 293369DEST_PATH_IMAGE087
是分潮角速率；
Figure 809801DEST_PATH_IMAGE088
为族号，对全日潮
Figure 916297DEST_PATH_IMAGE089
，对半日潮
Figure 783759DEST_PATH_IMAGE090
；
Figure 633903DEST_PATH_IMAGE033
为经度；
Figure 774815DEST_PATH_IMAGE091
为北京标准时区号；
Figure 52213DEST_PATH_IMAGE092
为各分潮经地潮校正后的平衡潮振幅，
Figure 672550DEST_PATH_IMAGE093
Figure 60806DEST_PATH_IMAGE094
Figure 551830DEST_PATH_IMAGE095
Figure 767173DEST_PATH_IMAGE096
Figure 609227DEST_PATH_IMAGE097
Figure 801174DEST_PATH_IMAGE098
Figure 943443DEST_PATH_IMAGE099
Figure 61178DEST_PATH_IMAGE100
（31）
式中，
Figure 124949DEST_PATH_IMAGE032
为纬度。
（四）计算GPS数据调和常数
利用传统调和分析方法分析GPS观测数据，根据潮汐理论，假设
Figure 120587DEST_PATH_IMAGE101
为海潮负荷对观测站东向、北向和径向中某一个方向的瞬时位移，则主要分潮的调和常数可表示为，
Figure 586203DEST_PATH_IMAGE102
（32）
式中，
Figure 641884DEST_PATH_IMAGE103
为位移的平均位置；
Figure 694416DEST_PATH_IMAGE026
和
Figure 228165DEST_PATH_IMAGE027
为各分潮的调和常数；
Figure 548288DEST_PATH_IMAGE083
和
Figure 774870DEST_PATH_IMAGE084
为分潮的交点因子和交点订正角，
Figure 78813DEST_PATH_IMAGE085
为平衡潮第
Figure 908929DEST_PATH_IMAGE034
个分潮在
Figure 83558DEST_PATH_IMAGE086
时刻的位相；
Figure 215462DEST_PATH_IMAGE087
是分潮角速率；这里
Figure 741122DEST_PATH_IMAGE104
、
Figure 616674DEST_PATH_IMAGE105
及
Figure 911389DEST_PATH_IMAGE106
与时间有关，与地点无关；而
Figure 981238DEST_PATH_IMAGE026
和
Figure 728614DEST_PATH_IMAGE027
与地点有关，与时间无关；故上式可写作
Figure 673437DEST_PATH_IMAGE107
(33)
其中，
Figure 557079DEST_PATH_IMAGE108
(34)
Figure 296365DEST_PATH_IMAGE109
(35)
按最小二乘法原理，解算出
Figure 29572DEST_PATH_IMAGE110
和
Figure 246927DEST_PATH_IMAGE111
，再按照下式计算出调和常数
Figure 516234DEST_PATH_IMAGE112
和
Figure 426422DEST_PATH_IMAGE113
，
Figure 148390DEST_PATH_IMAGE114
(36)
在实际计算过程中，由于GPS卫星轨道的周期（约为11.967小时）、GPS星座的重复周期（约为23.9319小时）与部分分潮的周期(K2周期约为11.9672小时，K1周期约为23.9345小时)非常相近，容易产生耦合现象。GPS卫星轨道的误差极易被这些分潮吸收，导致K1/K2分潮负荷位移的估算结果误差较大。
（五）利用相对导纳纠正K1和K2分潮估值误差
对于同一个潮族中，都有主要分潮和次要分潮，主要分潮占比大，次要分潮占比小，但同一个潮族中分潮的潮汐特征基本一致。
令
Figure 936480DEST_PATH_IMAGE005
(37)
式中，
Figure 794714DEST_PATH_IMAGE004
为由GPS数据调和分析得到的各分潮的振幅；
Figure 141382DEST_PATH_IMAGE003
为各分潮平衡潮振幅，
Figure 350646DEST_PATH_IMAGE006
为模数。
对同一个潮族，各分潮的相对导纳特征具有类似特征，对半日潮族，S2对M2，N2对M2，K2对S2相对导纳特征相似，所以M2、N2、S2和K2分潮对应的R值也具有相似性，将各半日分潮的
Figure 175383DEST_PATH_IMAGE007
值进行曲线拟合，将K2分潮的异常
Figure 658098DEST_PATH_IMAGE007
值纠正到正常的拟合曲线上，这样就解决了K2分潮负荷位移的估算结果误差较大的问题，如图2所示。同理，对全日潮族，O1对K1，P1对K1和Q1对O1分潮的相对导纳特征相似，所以K1、P1、O1和Q1分潮对应的R值也具有相似性，将各半日分潮的
Figure 644508DEST_PATH_IMAGE007
值进行曲线拟合，将K1分潮的异常
Figure 341069DEST_PATH_IMAGE007
值纠正到正常的拟合曲线上，这样就解决了K1分潮负荷位移的估算结果误差较大的问题。
实施例2：
如图3所示，本发明同时提供了一种纠正GPS观测得到的K1和K2分潮误差的设备，设备包括至少一个处理器和至少一个存储器，同时还包括通信接口和内部总线；所述存储器中存储有计算机程序或指令；所述处理器执行所述计算机程序或指令时，可以实现如实施例1所述的纠正GPS观测得到的K1和K2分潮误差的方法。其中内部总线可以是工业标准体系结构(Industry Standard Architecture，ISA)总线、外部设备互连(PeripheralComponent，PCI)总线或扩展工业标准体系结构(.XtendedIndustry StandardArchitecture，EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本申请附图中的总线并不限定仅有一根总线或一种类型的总线。其中存储器可能包含高速RAM存储器，也可能还包括非易失性存储NVM，例如至少一个磁盘存储器，还可以为U盘、移动硬盘、只读存储器、磁盘或光盘等。
设备可以被提供为终端、服务器或其它形态的设备。
图3是为示例性示出的一种设备的框图。设备可以包括以下一个或多个组件：处理组件，存储器，电源组件，多媒体组件，音频组件，输入/输出(I/O)的接口，传感器组件，以及通信组件。处理组件通常控制电子设备的整体操作，诸如与显示，电话呼叫，数据通信，相机操作和记录操作相关联的操作。处理组件可以包括一个或多个处理器来执行指令，以完成上述的方法的全部或部分步骤。此外，处理组件可以包括一个或多个模块，便于处理组件和其他组件之间的交互。例如，处理组件可以包括多媒体模块，以方便多媒体组件和处理组件之间的交互。
存储器被配置为存储各种类型的数据以支持在电子设备的操作。这些数据的示例包括用于在电子设备上操作的任何应用程序或方法的指令，联系人数据，电话簿数据，消息，图片，视频等。存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。
电源组件为电子设备的各种组件提供电力。电源组件可以包括电源管理系统，一个或多个电源，及其他与为电子设备生成、管理和分配电力相关联的组件。多媒体组件包括在所述电子设备和用户之间的提供一个输出接口的屏幕。在一些实施例中，屏幕可以包括液晶显示器(LCD)和触摸面板(TP)。如果屏幕包括触摸面板，屏幕可以被实现为触摸屏，以接收来自用户的输入信号。触摸面板包括一个或多个触摸传感器以感测触摸、滑动和触摸面板上的手势。所述触摸传感器可以不仅感测触摸或滑动动作的边界，而且还检测与所述触摸或滑动操作相关的持续时间和压力。在一些实施例中，多媒体组件包括一个前置摄像头和/或后置摄像头。当电子设备处于操作模式，如拍摄模式或视频模式时，前置摄像头和/或后置摄像头可以接收外部的多媒体数据。每个前置摄像头和后置摄像头可以是一个固定的光学透镜系统或具有焦距和光学变焦能力。
音频组件被配置为输出和/或输入音频信号。例如，音频组件包括一个麦克风(MIC)，当电子设备处于操作模式，如呼叫模式、记录模式和语音识别模式时，麦克风被配置为接收外部音频信号。所接收的音频信号可以被进一步存储在存储器或经由通信组件发送。在一些实施例中，音频组件还包括一个扬声器，用于输出音频信号。 I/O接口为处理组件和外围接口模块之间提供接口，上述外围接口模块可以是键盘，点击轮，按钮等。这些按钮可包括但不限于：主页按钮、音量按钮、启动按钮和锁定按钮。
传感器组件包括一个或多个传感器，用于为电子设备提供各个方面的状态评估。例如，传感器组件可以检测到电子设备的打开/关闭状态，组件的相对定位，例如所述组件为电子设备的显示器和小键盘，传感器组件还可以检测电子设备或电子设备一个组件的位置改变，用户与电子设备接触的存在或不存在，电子设备方位或加速/减速和电子设备的温度变化。传感器组件可以包括接近传感器，被配置用来在没有任何的物理接触时检测附近物体的存在。传感器组件还可以包括光传感器，如CMOS或CCD图像传感器，用于在成像应用中使用。在一些实施例中，该传感器组件还可以包括加速度传感器，陀螺仪传感器，磁传感器，压力传感器或温度传感器。
通信组件被配置为便于电子设备和其他设备之间有线或无线方式的通信。电子设备可以接入基于通信标准的无线网络，如WiFi，2G或3G，或它们的组合。在一个示例性实施例中，通信组件经由广播信道接收来自外部广播管理系统的广播信号或广播相关信息。在一个示例性实施例中，所述通信组件还包括近场通信(NFC)模块，以促进短程通信。例如，在NFC模块可基于射频识别(RFID)技术，红外数据协会(IrDA)技术，超宽带(UWB)技术，蓝牙(BT)技术和其他技术来实现。
在示例性实施例中，电子设备可以被一个或多个应用专用集成电路(ASIC)、数字信号处理器(DSP)、数字信号处理设备(DSPD)、可编程逻辑器件(PLD)、现场可编程门阵列(FPGA)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述方法。
实施例3：
本发明还提供了一种非易失性计算机可读存储介质，计算机可读存储介质中存储有计算机程序或指令，所述计算机程序或指令被处理器执行时，可以实现如实施例1所述的纠正GPS观测得到的K1和K2分潮误差的方法。
具体地，可以提供配有可读存储介质的系统、装置或设备，在该可读存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码，且使该系统、装置或设备的计算机或处理器读出并执行存储在该可读存储介质中的指令。在这种情况下，从可读介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能，因此机器可读代码和存储机器可读代码的可读存储介质构成了本发明的一部分。
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘（如 CD-ROM、CD-R、CD-RW、DVD- 20 ROM、DVD-RAM、DVD-RW、DVD-RW）、磁带等。存储介质可以是通用或专用计算机能够存取的任何可用介质。
应理解，上述处理器可以是中央处理单元(英文：Central Processing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：Digital Signal Processor，简称：DSP)、专用集成电路(英文：Application Specific Integrated Circuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
应理解存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(Application Specific Integrated Circuits，简称：ASIC)中。当然，处理器和存储介质也可以作为分立组件存在于终端或服务器中。
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
上述虽然对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (5)

1.一种纠正GPS观测得到的分潮误差的方法，其特征在于，包括以下步骤：
步骤1，确定分潮对；所述确定分潮对为确定N2对M2、S2对M2、K2对M2、K2对S2、P1对K1、O1对K1、Q1对K1和Q1对O1的分潮对；
步骤2，建立各分潮对的相对关系，获取各个网格点的振幅比
Figure DEST_PATH_IMAGE002A
及迟角差
Figure DEST_PATH_IMAGE004A
；具体为：已知有N个观测点，分别用
Figure DEST_PATH_IMAGE006A
代表，各点的经纬度为
Figure DEST_PATH_IMAGE008A
，其中
Figure DEST_PATH_IMAGE010A
对南纬用负值表示；
Figure DEST_PATH_IMAGE012A
，
Figure DEST_PATH_IMAGE010AA
以度为单位；
Figure DEST_PATH_IMAGE014A
和
Figure DEST_PATH_IMAGE016A
为各验潮站观测点的振幅比和迟角差；
对于观测站
Figure DEST_PATH_IMAGE018A
点，计算它对其它各站点j的距离为
Figure DEST_PATH_IMAGE020A
（8）
求
Figure DEST_PATH_IMAGE018AA
点的加权平均值
Figure DEST_PATH_IMAGE022A
（9）
Figure DEST_PATH_IMAGE024A
（10）
其中，
Figure DEST_PATH_IMAGE026A
代表对所有观测点求和，但不包括
Figure DEST_PATH_IMAGE018AAA
点；
Figure DEST_PATH_IMAGE028A
是修正后的开尔文波波长，其根据
Figure DEST_PATH_IMAGE018AAAA
点附近的水深确定的，方法如下：
设已知海区内及周围各网格点的水深值，即对点
Figure DEST_PATH_IMAGE030A
其经纬度为
Figure DEST_PATH_IMAGE032A
，
Figure DEST_PATH_IMAGE034A
,水深为Dm；对
Figure DEST_PATH_IMAGE018_5A
点，计算网格点
Figure DEST_PATH_IMAGE030AA
与验潮站点
Figure DEST_PATH_IMAGE018_6A
的距离为
Figure DEST_PATH_IMAGE036A
（11）
选取所有
Figure DEST_PATH_IMAGE038A
的点，对其水深值进行平均，记为
Figure DEST_PATH_IMAGE040A
；
计算相应的开尔文波波速
Figure DEST_PATH_IMAGE042A
Figure DEST_PATH_IMAGE044A
（12）
开尔文波波长
Figure DEST_PATH_IMAGE046A
（13）
其中，ui为开尔文波波速，T为潮波周期，即
Figure DEST_PATH_IMAGE048A
（14）
Figure DEST_PATH_IMAGE050A
是潮波角速率；
最后选取
Figure DEST_PATH_IMAGE052A
Figure DEST_PATH_IMAGE054A
（15）
对所有
Figure DEST_PATH_IMAGE056A
各点算出
Figure DEST_PATH_IMAGE058A
和
Figure DEST_PATH_IMAGE060A
，并计算
Figure DEST_PATH_IMAGE062A
（16）
及均方根偏差
Figure DEST_PATH_IMAGE064A
（17）
若
Figure DEST_PATH_IMAGE066A
（18）
则认为第
Figure DEST_PATH_IMAGE018_7A
点数据异常，需要舍去；
为验证去除的异常值后的精确度，由
Figure DEST_PATH_IMAGE058AA
和
Figure DEST_PATH_IMAGE060AA
反算各观测点的振幅比
Figure DEST_PATH_IMAGE014AA
及迟角差
Figure DEST_PATH_IMAGE016AA
Figure DEST_PATH_IMAGE068A
（19）
Figure DEST_PATH_IMAGE070A
Figure DEST_PATH_IMAGE072A
（20）
再由
Figure DEST_PATH_IMAGE074A
（21）
Figure DEST_PATH_IMAGE076A
（22）
计算振幅比与迟角差的均方差
Figure DEST_PATH_IMAGE078A
（23）
Figure DEST_PATH_IMAGE080A
（24）
按公式（18）所记载的方法先舍弃不合理的观测值，按预先给定的分辨率设定各网格点，记各网格点的序号为
Figure DEST_PATH_IMAGE082A
，其对应的经纬度为
Figure DEST_PATH_IMAGE084A
，计算点
Figure DEST_PATH_IMAGE086A
到各观测点
Figure DEST_PATH_IMAGE018_8A
的距离为
Figure DEST_PATH_IMAGE088A
（25）
求
Figure DEST_PATH_IMAGE086AA
点的加权平均值
Figure DEST_PATH_IMAGE090A
（26）
Figure DEST_PATH_IMAGE092A
（27）
其中
Figure 818875DEST_PATH_IMAGE093
计算方法与前述公式（15）中
Figure 549064DEST_PATH_IMAGE094
的计算方法相同，
由
Figure DEST_PATH_IMAGE096A
和
Figure DEST_PATH_IMAGE098A
反算各个网格点的振幅比
Figure DEST_PATH_IMAGE002AA
及迟角差
Figure DEST_PATH_IMAGE004AA
；
Figure DEST_PATH_IMAGE100A
（28）
Figure DEST_PATH_IMAGE102A
Figure DEST_PATH_IMAGE104A
（29）；
步骤3，计算平衡潮调和常数，得到各分潮经地潮校正后的平衡潮振幅
Figure DEST_PATH_IMAGE106A
；
步骤4，计算GPS数据调和常数，获取由GPS数据调和分析得到的各分潮的振幅
Figure DEST_PATH_IMAGE108A
；
步骤5，根据步骤2、步骤3和步骤4的结果并利用相对导纳纠正K1和K2分潮估值误差；具体为：
令
Figure DEST_PATH_IMAGE110A
式中，
Figure DEST_PATH_IMAGE108AA
为由GPS数据调和分析得到的各分潮的振幅，
Figure DEST_PATH_IMAGE106AA
为各分潮平衡潮振幅，
Figure DEST_PATH_IMAGE112A
为模数，将各半日分潮N2、M2、S2和K2的
Figure DEST_PATH_IMAGE114A
值进行曲线拟合，然后将K2分潮的异常
Figure DEST_PATH_IMAGE114AA
值纠正到正常的拟合曲线上，从而完成K2分潮误差的纠正，进而提高K2分潮调和常数的准确度；同理，将各全日分潮K1、P1、O1和Q1的
Figure DEST_PATH_IMAGE114AAA
值进行曲线拟合，然后将K1分潮的异常
Figure DEST_PATH_IMAGE114AAAA
值纠正到正常的拟合曲线上，从而完成K1分潮误差的纠正，进而提高K1分潮调和常数的准确度。
2.如权利要求1所述的一种纠正GPS观测得到的分潮误差的方法，其特征在于，所述步骤1具体过程为：
对某一分潮i，导纳Mi为
Figure DEST_PATH_IMAGE116A
（1）
其中，
Figure DEST_PATH_IMAGE118A
，
Figure DEST_PATH_IMAGE120A
和
Figure DEST_PATH_IMAGE122A
为分潮的调和常数，振幅和迟角，
Figure DEST_PATH_IMAGE124A
为引潮力系数，在同一潮族中，主要分潮n相对于次要分潮m的关系可表示为
Figure DEST_PATH_IMAGE126A
（2）
其中
Figure DEST_PATH_IMAGE128
（3）
Figure DEST_PATH_IMAGE130
（4）
式（3）中
Figure DEST_PATH_IMAGE132
即是传统潮汐学中的振幅比，
Figure DEST_PATH_IMAGE134
为迟角差，
Figure DEST_PATH_IMAGE136
为主要分潮n相对于次要分潮m的振幅关系，
Figure DEST_PATH_IMAGE138
为主要分潮n相对于次要分潮m的迟角关系；如果实际分潮的振幅比与引潮力之比相同，则
Figure DEST_PATH_IMAGE140
等于1；实际海洋中
Figure DEST_PATH_IMAGE140A
不会正好等于1，但对大部分海域它接近于1；如果两个分潮的迟角相等，则
Figure DEST_PATH_IMAGE142
等于零；实际海洋中
Figure DEST_PATH_IMAGE142A
一般不等于零，对大部分海域它大于零或小于零；
确定N2对M2、S2对M2、K2对M2、K2对S2、P1对K1、O1对K1、Q1对K1和Q1对O1分潮对；
Figure DEST_PATH_IMAGE144
（5）
Figure DEST_PATH_IMAGE146
（6）
令振幅比分别为
Figure DEST_PATH_IMAGE148
迟角差分别为
Figure DEST_PATH_IMAGE150
（7）
对某两个分潮，设振幅比为
Figure DEST_PATH_IMAGE152
，迟角差为
Figure DEST_PATH_IMAGE154
,记
Figure DEST_PATH_IMAGE156
，
Figure DEST_PATH_IMAGE158
。
3.如权利要求1所述的一种纠正GPS观测得到的分潮误差的方法，其特征在于，所述步骤3的具体过程为：
对任意观测点
Figure DEST_PATH_IMAGE160
，其对应的平衡潮潮高
Figure DEST_PATH_IMAGE162
可表示为
Figure 12624DEST_PATH_IMAGE163
（30）
式中，
Figure DEST_PATH_IMAGE165A
和
Figure DEST_PATH_IMAGE167A
为分潮的交点因子和交点订正角，
Figure DEST_PATH_IMAGE169A
为平衡潮第
Figure DEST_PATH_IMAGE018_9A
个分潮在
Figure DEST_PATH_IMAGE171A
时刻的位相；
Figure DEST_PATH_IMAGE173A
是分潮角速率；
Figure DEST_PATH_IMAGE175A
为族号，对全日潮
Figure DEST_PATH_IMAGE177
，对半日潮
Figure DEST_PATH_IMAGE179
；
Figure DEST_PATH_IMAGE012AA
为经度；
Figure DEST_PATH_IMAGE181
为北京标准时区号；
Figure DEST_PATH_IMAGE183
为各分潮经地潮校正后的平衡潮振幅，
Figure DEST_PATH_IMAGE185
Figure DEST_PATH_IMAGE187A
Figure DEST_PATH_IMAGE189A
Figure DEST_PATH_IMAGE191A
Figure DEST_PATH_IMAGE193A
Figure DEST_PATH_IMAGE195A
Figure DEST_PATH_IMAGE197A
Figure DEST_PATH_IMAGE199A
（31）
式中，
Figure DEST_PATH_IMAGE010AAA
为纬度。
4.一种纠正GPS观测得到的分潮误差的设备，其特征在于：所述设备包括至少一个处理器和至少一个存储器，所述处理器和存储器相耦合；所述存储器中存储有计算机程序或指令；所述处理器执行所述计算机程序或指令时，可以实现如权利要求1至3任意一项所述的纠正GPS观测得到的分潮误差的方法。
5.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机程序或指令，所述计算机程序或指令被处理器执行时，可以实现如权利要求1至3任意一项所述的纠正GPS观测得到的分潮误差的方法。
CN202210953013.5A 2022-08-10 2022-08-10 一种纠正gps观测得到的分潮误差的方法及设备 Active CN115032670B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211260841.7A CN115616630A (zh) 2022-08-10 2022-08-10 利用gnss研究海潮负荷位移时降低分潮误差的方法
CN202210953013.5A CN115032670B (zh) 2022-08-10 2022-08-10 一种纠正gps观测得到的分潮误差的方法及设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210953013.5A CN115032670B (zh) 2022-08-10 2022-08-10 一种纠正gps观测得到的分潮误差的方法及设备

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN202211260841.7A Division CN115616630A (zh) 2022-08-10 2022-08-10 利用gnss研究海潮负荷位移时降低分潮误差的方法

Publications (2)

Publication Number Publication Date
CN115032670A CN115032670A (zh) 2022-09-09
CN115032670B true CN115032670B (zh) 2022-11-18

Family

ID=83130549

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202210953013.5A Active CN115032670B (zh) 2022-08-10 2022-08-10 一种纠正gps观测得到的分潮误差的方法及设备
CN202211260841.7A Pending CN115616630A (zh) 2022-08-10 2022-08-10 利用gnss研究海潮负荷位移时降低分潮误差的方法

Family Applications After (1)

Application Number Title Priority Date Filing Date
CN202211260841.7A Pending CN115616630A (zh) 2022-08-10 2022-08-10 利用gnss研究海潮负荷位移时降低分潮误差的方法

Country Status (1)

Country Link
CN (2) CN115032670B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117572467B (zh) * 2023-10-18 2024-09-03 澳门大学 一种多海况潮位观测方法、装置、电子设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20170017655A (ko) * 2015-08-06 2017-02-15 대한민국(해양수산부 국립해양조사원장) 수치조류도 소프트웨어 기반 조류도 제작 방법
CN106649991A (zh) * 2016-11-15 2017-05-10 国家海洋局第海洋研究所 一种基于相对导纳计算次要分潮经验同潮图的推算方法
CN113761808A (zh) * 2021-11-08 2021-12-07 长沙理工大学 一种基于gps和经验潮汐模型的地表潮汐位移获取方法及应用方法和系统
CN114444330A (zh) * 2022-04-02 2022-05-06 自然资源部第一海洋研究所 一种顾及平均效应的gnss海潮负荷位移的预测方法
CN114662348A (zh) * 2022-05-25 2022-06-24 山东省科学院海洋仪器仪表研究所 一种基于gnss技术建立时变海潮负载位移模型的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109614711B (zh) * 2018-12-12 2022-10-21 国家海洋局第一海洋研究所 一种基于潮汐潮流数值模拟提高多分潮调和常数计算精度的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20170017655A (ko) * 2015-08-06 2017-02-15 대한민국(해양수산부 국립해양조사원장) 수치조류도 소프트웨어 기반 조류도 제작 방법
CN106649991A (zh) * 2016-11-15 2017-05-10 国家海洋局第海洋研究所 一种基于相对导纳计算次要分潮经验同潮图的推算方法
CN113761808A (zh) * 2021-11-08 2021-12-07 长沙理工大学 一种基于gps和经验潮汐模型的地表潮汐位移获取方法及应用方法和系统
CN114444330A (zh) * 2022-04-02 2022-05-06 自然资源部第一海洋研究所 一种顾及平均效应的gnss海潮负荷位移的预测方法
CN114662348A (zh) * 2022-05-25 2022-06-24 山东省科学院海洋仪器仪表研究所 一种基于gnss技术建立时变海潮负载位移模型的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Determination of Tidal Components and Tidal Types Using Harmonic Analysis in the West Sumatera Waters;Sorja Koesuma等;《2020 IEEE Asia-Pacific Conference on Geoscience, Electronics and Remote Sensing Technology (AGERS)》;20210622;全文 *
潮汐调和分析方法的探讨;张凤烨等;《海洋科学》;20110630;第35卷(第6期);全文 *

Also Published As

Publication number Publication date
CN115032670A (zh) 2022-09-09
CN115616630A (zh) 2023-01-17

Similar Documents

Publication Publication Date Title
CN111983635B (zh) 位姿确定方法及装置、电子设备和存储介质
US20210109228A1 (en) Identifying gnss navigation data as potentially manipulated or as trustworthy at least partially based on an estimated deviation of a second estimate of a satellite state from a first estimate of the satellite state
JP2018163159A (ja) 地点観測を使用して天気予報を精密化する方法及びシステム
CN111127563A (zh) 联合标定方法、装置、电子设备及存储介质
US20130027427A1 (en) Associating Digital Images with Waypoints
CN111624634B (zh) 基于深度卷积神经网络的卫星定位误差评估方法和系统
JPWO2016147569A1 (ja) 衛星測位システム、電子機器及び測位方法
CN115032670B (zh) 一种纠正gps观测得到的分潮误差的方法及设备
CN109492531B (zh) 人脸图像关键点提取方法、装置、存储介质及电子设备
EP3805797A1 (en) Detection of gnss navigation data manipulation based on gnss system time differences
US10587704B2 (en) Location accurate mobile events and social content
CN107272038B (zh) 一种高精度定位的方法及设备
CN106383355A (zh) 基于高精度时间的gnss时差监测方法及系统
TW201205105A (en) Method and apparatus for updating transformation information parameters used in global navigation satellite system
CN105635570A (zh) 拍摄预览方法和系统
CN110168589B (zh) 用于跟踪、捕获和跨多个设备同步活动数据的方法和装置
CN113050030A (zh) 基于到达角度测距的定位方法及装置
Cortés et al. Performance evaluation of adaptive tracking techniques with direct-state Kalman filter
US20210109226A1 (en) Identifying second gnss navigation data as potentially manipulated or as trustworthy at least partially based on first gnss navigation data
US8913993B2 (en) Terminal apparatus with control unit to change accuracy of position information
CN112799110B (zh) 一种顾及Doppler的北斗修正伪距单点定位方法、系统与设备
CN112033544B (zh) 一种热成像仪、测温方法、装置、处理器及存储介质
US20210109231A1 (en) Detection of gnss data manipulation based on correctional gnss data
CN116125495B (zh) 电离层改正数的确定方法、设备、存储介质及程序产品
CN109524025B (zh) 一种唱歌评分方法、装置、电子设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant