利用GNSS研究海潮负荷位移时降低分潮误差的方法与流程

利用gnss研究海潮负荷位移时降低分潮误差的方法
技术领域
1.本发明属于海洋监测技术领域，尤其涉及一种利用gnss研究海潮负荷位移时降低分潮误差的方法。


背景技术：

2.随着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研究海潮负荷位移的难点和热点问题。


技术实现要素：

3.本发明利用s2对m2，n2对m2，k2对s2，o1对k1，p1对k1和q1对o1分潮的相对导纳纠正gps观测得到k1和k2分潮估值误差过大问题。
4.本发明第一方面提供了一种利用gnss研究海潮负荷位移时降低分潮误差的方法，其特征在于，包括以下步骤：
5.步骤1，确定分潮对；所述确定分潮对为确定n2对m2、s2对m2、k2对m2、k2对s2、p1对k1、o1对k1、q1对k1和q1对o1的分潮对；
6.步骤2，建立各分潮对的相对关系，获取各个网格点的振幅比及迟角差
7.步骤3，计算平衡潮调和常数，得到各分潮经地潮校正后的平衡潮振幅
8.步骤4，计算gps数据调和常数，获取由gps数据调和分析得到的各分潮的振幅h
gpsi
；根据潮汐理论，假设ζ为海潮负荷对观测站中某一个方向的瞬时位移，则主要分潮的调和常数可表示为，
[0009][0010]
式中，s0为位移的平均位置；h和g为各分潮的调和常数；fi和ui为分潮的交点因子和交点订正角，v0为平衡潮第i个分潮在t＝0时刻的位相；ωi是分潮角速率；这里f、ωt及ν0+u与时间有关，与地点无关；而h和g与地点有关，与时间无关；故上式可写作
[0011][0012]
其中，
[0013][0014][0015]
按最小二乘法原理，解算出xi和yi，再按照下式计算出调和常数h
gps
和g
gps
，
[0016][0017]
步骤5，根据步骤2、步骤3和步骤4的结果并利用相对导纳纠正k1和k2分潮估值误差。
[0018]
在一种可能的设计中，所述步骤1具体过程为：
[0019]
对某一分潮i，导纳mi为
[0020][0021]
其中，hi和gi为分潮的调和常数，振幅和迟角，ci为引潮力系数，在同一潮族中，主要分潮n相对于次要分潮m的关系可表示为
[0022][0023]
其中
[0024][0025]am/n
＝g
m-gnꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0026]
式(3)中hm/hn即是传统潮汐学中的振幅比，g
m-gn为迟角差，a
*m/n
为主要分潮n相对于次要分潮m的振幅关系，a
m/n
为主要分潮n相对于次要分潮m的迟角关系；如果实际分潮的振幅比与引潮力之比相同，则等于1；实际海洋中不会正好等于1，但对大部分海域它接近于1；如果两个分潮的迟角相等，则a
m/n
等于零；实际海洋中a
m/n
一般不等于零，对大部分海域它大于零或小于零；
[0027]
确定n2对m2、s2对m2、k2对m2、k2对s2、p1对k1、o1对k1、q1对k1和q1对o1分潮对；
[0028][0029]gm
＝a
m/n
+gnꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
令振幅比分别为迟角差分别为
[0030]
对某两个分潮，设振幅比为h，迟角差为g,记a＝hcosg，b＝hsing；
[0031]
所述步骤2的具体过程为：
[0032]
已知有n个观测点，分别用i(i＝1,2...,n)代表，各点的经纬度为其中对南纬用负值表示；λ，以度为单位；和为各验潮站观测点的振幅比和迟角差；
[0033]
对于观测站i点，计算它对其它各站点j的距离为
[0034][0035]
求i点的加权平均值
[0036][0037][0038]
其中，代表对所有观测点求和，但不包括i点；ri是修正后的开尔文波波长，其根据i点附近的水深确定的，方法如下：
[0039]
设已知海区内及周围各网格点的水深值，即对点m其经纬度为设已知海区内及周围各网格点的水深值，即对点m其经纬度为水深为dm；对i点，计算网格点m与验潮站点i的距离为
[0040][0041]
选取所有r
im
＜150
×
103m的点，对其水深值进行平均，记为
[0042]
计算相应的开尔文波波速
[0043][0044]
开尔文波波长
[0045]
li＝uit
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0046]
其中，ui为开尔文波波速，t为潮波周期，即
[0047]
t＝360
°
/ω
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0048]
ω是潮波角速率；
[0049]
最后选取
[0050]ri
＝κli(κ＝1/4)
ꢀꢀꢀꢀꢀ
(15)
[0051]
对所有i＝1,2......,n各点算出和并计算
[0052][0053]
及均方根偏差
[0054][0055]
若
[0056]
δi≥mσ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0057]
则认为第i点数据异常，需要舍去；
[0058]
为验证去除的异常值后的精确度，由和反算各观测点的振幅比及迟角差
[0059][0060][0061]
再由
[0062][0063][0064]
计算振幅比与迟角差的均方差
[0065][0066][0067]
按公式(18)所记载的方法先舍弃不合理的观测值，按预先给定的分辨率设定各网
格点，记各网格点的序号为k(k＝1,2......)，其对应的经纬度为计算点k到各观测点i的距离为
[0068][0069]
求k点的加权平均值
[0070][0071][0072]
其中rk计算方法与前述公式(15)中ri的计算方法相同，
[0073]
由和反算各个网格点的振幅比及迟角差
[0074][0075][0076]
在一种可能的设计中，所述步骤3的具体过程为：
[0077]
对任意观测点其对应的平衡潮潮高ζ
eq
可表示为
[0078][0079]
式中，fi和ui为分潮的交点因子和交点订正角，v0为平衡潮第i个分潮在t＝0时刻的位相；ωi是分潮角速率；p为族号，对全日潮p＝1，对半日潮p＝2；λ为经度；s＝-8h为北京标准时区号；为各分潮经地潮校正后的平衡潮振幅，
[0080][0081][0082][0083][0084]
式中，为纬度。
[0085]
在一种可能的设计中，所述步骤5的具体过程为：
[0086]
令
[0087]
式中，h
gpsi
为由gps数据调和分析得到的各分潮的振幅，为各分潮平衡潮振幅，ψi为模数，将各半日分潮n2、m2、s2和k2的ψ值进行曲线拟合，然后将k2分潮的异常ψ值纠正到正常的拟合曲线上，从而完成k2分潮误差的纠正，进而提高k2分潮调和常数的准确度；同理，将各全日分潮k1、p1、o1和q1的ψ值进行曲线拟合，然后将k1分潮的异常ψ值纠正到正常的拟合曲线上，从而完成k1分潮误差的纠正，进而提高k1分潮调和常数的准确度。
[0088]
本发明第二方面还提供了一种利用gnss研究海潮负荷位移时降低分潮误差的设备，所述设备包括至少一个处理器和至少一个存储器，所述处理器和存储器相耦合；所述存储器中存储有计算机程序或指令；所述处理器执行所述计算机程序或指令时，可以实现如第一方面所述的利用gnss研究海潮负荷位移时降低分潮误差的方法。
[0089]
本发明第三方面还提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序或指令，所述计算机程序或指令被处理器执行时，可以实现如第一方面所述的利用gnss研究海潮负荷位移时降低分潮误差的方法。
[0090]
与现有技术相比，本发明提供了一种利用gnss研究海潮负荷位移时降低分潮误差的方法，并产生如下的有益效果：
[0091]
1.针对gps轨道周期与k2分潮周期一致、gps星座的重复周期与k1分潮的周期一致从而产生耦合作用，导致k2分潮和k1分潮估值误差过大的难题。本发明提出了一种基于相对导纳纠正gps观测得到的k1和k2分潮误差的方法，能够有效改善gps观测得到的k1和k2分潮的误差并提高k1和k2分潮的准确度。
[0092]
2.通过gnss准确地获取观测站位由海潮负荷引起的位移形变，并进而得到海潮负荷位移调和常数，本发明为提高研究海潮负荷效应的准确度提供一种新的方法。
[0093]
3.为提高研究海潮负荷效应的准确度，本发明采用建立8个分潮对相对关系的方法，能有效得到各潮族的相对导纳关系。结合本发明采用的gps数据分析调和常数与相应平衡潮调和常数的比值，能有效的纠正k1和k2分潮的误差并提高k1和k2分潮的准确度。
附图说明
[0094]
图1为本发明纠正gps观测得到的k1和k2分潮误差方法的流程框图。
[0095]
图2为本发明纠正k2分潮误差的示意图。
[0096]
图3为本发明纠正gps观测得到的k1和k2分潮误差的设备结构简易示意图。
具体实施方式
[0097]
下面结合具体实施例对发明进行进一步说明。
[0098]
实施例1：
[0099]
为实现基于相对导纳纠正gps观测得到的k1和k2分潮误差的推算技术，其包括相对导纳计算和纠正gps观测得到的k1和k2分潮误差的计算步骤，大体流程如图1所示。
[0100]
(一)确定分潮对
[0101]
对某一分潮i，导纳mi为
[0102][0103]
其中hi和gi为分潮的调和常数，振幅和迟角；ci为引潮力系数。在同
一潮族中，主要分潮n相对于次要分潮m的关系可表示为
[0104][0105]
其中
[0106][0107]am/n
＝g
m-gnꢀꢀꢀ
(4)式(3)中hm/hn即是传统潮汐学中的振幅比，g
m-gn为迟角差。如果实际分潮的振幅比与引潮力之比相同，则等于1。实际海洋中不会正好等于1，但对大部分海域它接近于1。如果两个分潮的迟角相等，则a
m/n
等于零。实际海洋中a
m/n
一般不等于零，对大部分海域它大于零或小于零。
[0108]
确定n2对m2，s2对m2，k2对m2，k2对s2，p1对k1，o1对k1，q1对k1和q1对o1分潮对。
[0109][0110]gm
＝a
m/n
+gnꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0111]
令振幅比分别为迟角差分别为
[0112]
对某两个分潮，设振幅比为h，迟角差为g，记a＝hcosg，b＝hsing。
[0113]
(二)建立各分潮对的相对关系
[0114]
已知有n个观测点，分别用i(i＝1,2...,n)代表，各点的经纬度为其中对南纬用负值表示。λ，以度为单位。hi和gi为各验潮站观测点的振幅比和迟角差。
[0115]
对于观测站i点，计算它对其它各站点j的距离为
[0116][0117]
求i点的加权平均值
[0118][0119][0120]
其中，代表对所有观测点求和，但不包括i点。ri是根据i点附近的水深确定的，方
法如下：
[0121]
设已知海区内及周围各网格点的水深值，即对点m其经纬度为设已知海区内及周围各网格点的水深值，即对点m其经纬度为水深为dm。对i点，计算网格点m与验潮站点i的距离为
[0122][0123]
选取所有r
im
＜150
×
103m的点，对这些点的水深值进行平均，记为
[0124]
计算相应的开尔文波波速
[0125][0126]
开尔文波波长
[0127]
li＝uit
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
[0128]
其中，t为潮波周期，即
[0129]
t＝360
°
/ω
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(14)
[0130]
这里ω是潮波角速率。
[0131]
最后选取
[0132]ri
＝κli(κ＝1/4)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(15)
[0133]
对所有i＝1,2......,n各点算出和并计算
[0134][0135]
及均方根偏差
[0136][0137]
若
[0138]
δi≥mσ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0139]
则认为第i点数据可能异常，考虑是否舍去，这里m可取2。
[0140]
为验证去除的异常值后的精确度，由和反算各观测点的振幅比及迟角差
[0141]
[0142][0143]
再由
[0144][0145][0146]
计算振幅比与迟角差的均方差
[0147][0148][0149]
按上述方法先舍弃不合理的观测值，按预先给定的分辨率设定各网格点，记各网格点的序号为k(k＝1,2......)，其对应的经纬度为计算点k到各观测点i的距离为
[0150][0151]
求k点的加权平均值
[0152][0153][0154]rk
计算方法与前述方法类似。
[0155]
由和反算各个网格点的振幅比及迟角差
[0156]
[0157][0158]
(三)计算平衡潮调和常数
[0159]
对任意观测点其对应的平衡潮潮高ζ
eq
可表示为，
[0160][0161]
式中，fi和ui为分潮的交点因子和交点订正角，v0为平衡潮第i个分潮在t＝0时刻的位相；ωi是分潮角速率；p为族号，对全日潮p＝1，对半日潮p＝2；λ为经度；s＝-8h为北京标准时区号；为各分潮经地潮校正后的平衡潮振幅，
[0162][0163][0164][0165][0166]
式中，为纬度。
[0167]
(四)计算gps数据调和常数
[0168]
分析gps观测数据，根据潮汐理论，假设ζ为海潮负荷对观测站东向、北向和径向中某一个方向的瞬时位移，则主要分潮的调和常数可表示为，
[0169][0170]
式中，s0为位移的平均位置；h和g为各分潮的调和常数；fi和ui为分潮的交点因子和交点订正角，v0为平衡潮第i个分潮在t＝0时刻的位相；ωi是分潮角速率；这里f、ωt及ν0+u与时间有关，与地点无关；而h和g与地点有关，与时间无关；故上式可写作
[0171][0172]
其中，
[0173][0174][0175]
按最小二乘法原理，解算出xi和yi，再按照下式计算出调和常数h
gps
和g
gps
，
[0176][0177]
在实际计算过程中，由于gps卫星轨道的周期(约为11.967小时)、gps星座的重复周期(约为23.9319小时)与部分分潮的周期(k2周期约为11.9672小时，k1周期约为23.9345小时)非常相近，容易产生耦合现象。gps卫星轨道的误差极易被这些分潮吸收，导致k1/k2分潮负荷位移的估算结果误差较大。
[0178]
(五)利用相对导纳纠正k1和k2分潮估值误差
[0179]
对于同一个潮族中，都有主要分潮和次要分潮，主要分潮占比大，次要分潮占比小，但同一个潮族中分潮的潮汐特征基本一致。
[0180]
令
[0181][0182]
式中，h
gpsi
为由gps数据调和分析得到的各分潮的振幅；为各分潮平衡潮振幅，ψi为模数。
[0183]
对同一个潮族，各分潮的相对导纳特征具有类似特征，对半日潮族，s2对m2，n2对m2，k2对s2相对导纳特征相似，所以m2、n2、s2和k2分潮对应的r值也具有相似性，将各半日分潮的ψ值进行曲线拟合，将k2分潮的异常ψ值纠正到正常的拟合曲线上，这样就解决了k2分潮负荷位移的估算结果误差较大的问题，如图2所示。同理，对全日潮族，o1对k1，p1对k1和q1对o1分潮的相对导纳特征相似，所以k1、p1、o1和q1分潮对应的r值也具有相似性，将各半日分潮的ψ值进行曲线拟合，将k1分潮的异常ψ值纠正到正常的拟合曲线上，这样就解决了k1分潮负荷位移的估算结果误差较大的问题。
[0184]
实施例2：
[0185]
如图3所示，本发明同时提供了一种利用gnss研究海潮负荷位移时降低分潮误差的设备，设备包括至少一个处理器和至少一个存储器，同时还包括通信接口和内部总线；所述存储器中存储有计算机程序或指令；所述处理器执行所述计算机程序或指令时，可以实现如实施例1所述的利用gnss研究海潮负荷位移时降低分潮误差的方法。其中内部总线可以是工业标准体系结构(industry standard architecture，isa)总线、外部设备互连(peripheral component，pci)总线或扩展工业标准体系结构(.xtendedindustry standard architecture，eisa)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，本技术附图中的总线并不限定仅有一根总线或一种类型的总线。其中存储器可能包含高速ram存储器，也可能还包括非易失性存储nvm，例如至少一个磁盘存储器，还可以为u盘、移动硬盘、只读存储器、磁盘或光盘等。
[0186]
设备可以被提供为终端、服务器或其它形态的设备。
[0187]
图3是为示例性示出的一种设备的框图。设备可以包括以下一个或多个组件：处理组件，存储器，电源组件，多媒体组件，音频组件，输入/输出(i/o)的接口，传感器组件，以及通信组件。处理组件通常控制电子设备的整体操作，诸如与显示，电话呼叫，数据通信，相机操作和记录操作相关联的操作。处理组件可以包括一个或多个处理器来执行指令，以完成上述的方法的全部或部分步骤。此外，处理组件可以包括一个或多个模块，便于处理组件和
其他组件之间的交互。例如，处理组件可以包括多媒体模块，以方便多媒体组件和处理组件之间的交互。
[0188]
存储器被配置为存储各种类型的数据以支持在电子设备的操作。这些数据的示例包括用于在电子设备上操作的任何应用程序或方法的指令，联系人数据，电话簿数据，消息，图片，视频等。存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(sram)，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，可编程只读存储器(prom)，只读存储器(rom)，磁存储器，快闪存储器，磁盘或光盘。
[0189]
电源组件为电子设备的各种组件提供电力。电源组件可以包括电源管理系统，一个或多个电源，及其他与为电子设备生成、管理和分配电力相关联的组件。多媒体组件包括在所述电子设备和用户之间的提供一个输出接口的屏幕。在一些实施例中，屏幕可以包括液晶显示器(lcd)和触摸面板(tp)。如果屏幕包括触摸面板，屏幕可以被实现为触摸屏，以接收来自用户的输入信号。触摸面板包括一个或多个触摸传感器以感测触摸、滑动和触摸面板上的手势。所述触摸传感器可以不仅感测触摸或滑动动作的边界，而且还检测与所述触摸或滑动操作相关的持续时间和压力。在一些实施例中，多媒体组件包括一个前置摄像头和/或后置摄像头。当电子设备处于操作模式，如拍摄模式或视频模式时，前置摄像头和/或后置摄像头可以接收外部的多媒体数据。每个前置摄像头和后置摄像头可以是一个固定的光学透镜系统或具有焦距和光学变焦能力。
[0190]
音频组件被配置为输出和/或输入音频信号。例如，音频组件包括一个麦克风(mic)，当电子设备处于操作模式，如呼叫模式、记录模式和语音识别模式时，麦克风被配置为接收外部音频信号。所接收的音频信号可以被进一步存储在存储器或经由通信组件发送。在一些实施例中，音频组件还包括一个扬声器，用于输出音频信号。i/o接口为处理组件和外围接口模块之间提供接口，上述外围接口模块可以是键盘，点击轮，按钮等。这些按钮可包括但不限于：主页按钮、音量按钮、启动按钮和锁定按钮。
[0191]
传感器组件包括一个或多个传感器，用于为电子设备提供各个方面的状态评估。例如，传感器组件可以检测到电子设备的打开/关闭状态，组件的相对定位，例如所述组件为电子设备的显示器和小键盘，传感器组件还可以检测电子设备或电子设备一个组件的位置改变，用户与电子设备接触的存在或不存在，电子设备方位或加速/减速和电子设备的温度变化。传感器组件可以包括接近传感器，被配置用来在没有任何的物理接触时检测附近物体的存在。传感器组件还可以包括光传感器，如cmos或ccd图像传感器，用于在成像应用中使用。在一些实施例中，该传感器组件还可以包括加速度传感器，陀螺仪传感器，磁传感器，压力传感器或温度传感器。
[0192]
通信组件被配置为便于电子设备和其他设备之间有线或无线方式的通信。电子设备可以接入基于通信标准的无线网络，如wifi，2g或3g，或它们的组合。在一个示例性实施例中，通信组件经由广播信道接收来自外部广播管理系统的广播信号或广播相关信息。在一个示例性实施例中，所述通信组件还包括近场通信(nfc)模块，以促进短程通信。例如，在nfc模块可基于射频识别(rfid)技术，红外数据协会(irda)技术，超宽带(uwb)技术，蓝牙(bt)技术和其他技术来实现。
[0193]
在示例性实施例中，电子设备可以被一个或多个应用专用集成电路(asic)、数字
信号处理器(dsp)、数字信号处理设备(dspd)、可编程逻辑器件(pld)、现场可编程门阵列(fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述方法。
[0194]
实施例3：
[0195]
本发明还提供了一种非易失性计算机可读存储介质，计算机可读存储介质中存储有计算机程序或指令，所述计算机程序或指令被处理器执行时，可以实现如实施例1所述的利用gnss研究海潮负荷位移时降低分潮误差的方法。
[0196]
具体地，可以提供配有可读存储介质的系统、装置或设备，在该可读存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码，且使该系统、装置或设备的计算机或处理器读出并执行存储在该可读存储介质中的指令。在这种情况下，从可读介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能，因此机器可读代码和存储机器可读代码的可读存储介质构成了本发明的一部分。
[0197]
上述存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(sram)，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，可编程只读存储器(prom)，只读存储器(rom)，磁存储器，快闪存储器，磁盘或光盘(如cd-rom、cd-r、cd-rw、dvd-20rom、dvd-ram、dvd-rw、dvd-rw)、磁带等。存储介质可以是通用或专用计算机能够存取的任何可用介质。
[0198]
应理解，上述处理器可以是中央处理单元(英文：central processing unit，简称：cpu)，还可以是其他通用处理器、数字信号处理器(英文：digital signal processor，简称：dsp)、专用集成电路(英文：application specific integrated circuit，简称：asic)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合发明所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
[0199]
应理解存储介质耦合至处理器，从而使处理器能够从该存储介质读取信息，且可向该存储介质写入信息。当然，存储介质也可以是处理器的组成部分。处理器和存储介质可以位于专用集成电路(application specific integrated circuits，简称：asic)中。当然，处理器和存储介质也可以作为分立组件存在于终端或服务器中。
[0200]
这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。
[0201]
用于执行本公开操作的计算机程序指令可以是汇编指令、指令集架构(isa)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如smalltalk、c++等，以及常规的过程式编程语言—诸如“c”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包
括局域网(lan)或广域网(wan)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(fpga)或可编程逻辑阵列(pla)，该电子电路可以执行计算机可读程序指令，从而实现本公开的各个方面。
[0202]
以上所述仅为本技术的优选实施例而已，并不用于限制本技术，对于本领域的技术人员来说，本技术可以有各种更改和变化。凡在本技术的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本技术的保护范围之内。
[0203]
上述虽然对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。