空中湍流位置系统和方法

专利名称：空中湍流位置系统和方法
技术领域：
本发明通常涉及飞行器安全测量，更具体地涉及用于测量湍流
(turbulence )位置的系统。
背景技术：
晴空湍流(CAT， clear airturbulence)由于大量的不同因素而发生在典型的 喷气式飞机巡航姿态，所述因素包括射流的位置、对流层顶(tropopause) 边界的位置、局部(localized)恶劣天气(包括雷暴)的位置、和山岭地区的 位置。然而，所有晴空湍流现象都具有下面共同的一点所述现象发生在晴 空中，并且飞行器一旦遭遇到湍流区域则除了操纵或改变飞行高度以外，目 前几乎不具有避免这些现象的方式。
湍流(尤其是CAT)是空中旅行中要严重关注的事。 一些航空公司采取 了强大的测量来避免所预测的或先前发现的湍流的区域，有时会飞行许多额 外英里来避免湍流。在大多数国家地区中，航空公司可用于航班改航的选项 的数量由于空域性能或空中交通控制系统约束而受到限制。已经对提高湍气 流的识别进行了许多努力，但是，这些努力仅在一个或多个航班已经飞入乱 流之后有帮助，并且它们没有提供不断扫描所述空域的方式。目前，对于飞 机没有承担得起的可靠方式来及时检测CAT以采取回避行动。
CAT的问题是它发生在缺少特定物质(例如云)的晴空中。许多现有的 检测湍流的方法取决于空气中特定物质(水蒸气、灰尘等)的存在，这不助 于在飞行器(商业、军用、或无人驾驶的空中飞行器(UAV))能够飞行所处 的巡航姿态中检测CAT。例如，在班机飞行的高晴空中，对于依赖于空中的 颗粒的多普勒雷达或激光雷达系统来说通常不能检测到CAT。
早期的工作已经论证了使用GPS信号的闪烁现象(scintillation)来检测 CAT在技术上是切实可行的。然而，这些早期的方法不能够解决关于到GPS 接收机的距离来精确地定位湍气流的问题。现有方法仅能够随着信号的路径 (沿着到卫星的视线)确定某处湍流的存在。其它估计的方法太昂贵或者繁重而而不能在飞行器上安装实践。
因此，期望具有一种成本有效的系统和方法来^r测CAT，以^便在飞行器 贯穿湍气流区域之前通过更精确地测定湍流的距离和高度以增加空中交通系 统性能。进一步期望该系统和方法使得空中交通系统和操作员更精确地计划 避免这些湍气流区域。
结合附图以及前面的技术领域和背景技术内容，从随后详细的描述和所 附权利要求来看，本发明实施例的其它期望特征和特性将变得明显。

发明内容
公开了 一种用于计算从天线平台到空气湍流的距离的系统和方法。所述 方法包括确定与在从第一天线到第二天线的干扰图案的传播中的时移对应 的相关滞后(lag)值，所述第一天线与第二天线相距一段距离。基于所述相关 滞后值和天线间距，所述方法计算从天线平台到空中湍流的距离，并且发出 指示所述距离的警告。
所述方法通过下列步骤确定所述相关滞后值在第一天线处取样第一干 扰图案幅值/相位序列以获得第 一幅值/相位图案；和在第二天线处取样第二干 扰图案幅值/相位序列以获得第二幅值/相位图案。然后，所述方法以时间滞后 序列，计算所述第一干扰图案幅值/相位序列和所述第二干扰图案幅值/相位序 列之间的时间相关系数，以便确定所述幅值/相位图案在每个时间滞后处是否 匹配良好。所述方法选择与计算出的时间相关系数的最高值对应的相关滞后 值，其中在所述最高值所述幅值/相位图案最接近地匹配。具有最高相关性的 时间滞后值是相关滞后值。
所述相关滞后值和从第 一天线到第二天线的天线间距被用来计算到湍气 流区域的距离。可以发出指示到所述湍气流区域的距离的警告。
此处所述的本发明的各个示例实施例计算湍气流区域的距离以便适应下 列方案天线平台是否固定或移动和/或湍气流区域是否固定或移动；天线平 台是否处于卫星的^L线中；或者天线是否水平和/或垂直地分离。
此处所述的技术可以通过更精确地测定湍气流区域的位置(距离和高度) 来提高空中交通系统性能，提供了用于空中交通系统提供商和操作员的一种 方式来更精确地计划避免这些湍气流区域。


通过参考结合下列附图考虑的详细描述和权利要求，可以推导出本发明 的更全面的理解，其中在整个附图中，相同的附图标记指向类似的元件。
图l是用于定位空中湍流的环境的示意表示图2是^t配置成定位空中湍流的系统的示意表示图3是图示用于定位空中湍流的处理的流程图4是图示用于定位空中湍流的另一处理的流程图5图示了在每个天线取样的幅值图案的示例；
图6图示了时间相关系数相对于相关滞后值的示例；和
图7图示了到湍流的各个距离的相关波峰。
具体实施例方式
下面的详细描述实际上仅仅是图示性的，而往往不限制发明或者申请的 实施例以及所述实施例的使用。而且，往往不通过先前的技术领域、背景技 术、发明内容或下面的详细描述中呈现的明确或暗示的理论来进行约束。
此处可以在概念和/或逻辑块组件和各种处理步骤方面来描述本发明的 实施例。应当理解，所述块组件可以通过任意数目的硬件、软件、和/或被配 置成执行特定功能的固件组件来实现。例如，本发明的一个实施例可以利用 各种集成电路组件，例如存储器元件、数字信号处理元件、逻辑元件、查找 表等等，它们可以在一个或多个微处理器或其它控制期间的控制下实施各种 各样的功能。另外，本领域的普通技术人员将会理解，本发明的实施例可以 结合任意数目的通信系统和任意数目的平台来实践，并且此处描述的系统仅 仅是本发明的一个示例实施例。
为了简洁，此处可能不详细描述关于信号处理、数据传输、信令、测量 的传统技术以及系统的其它功能方面(和系统的各个操作组件)。而且，此处 包含的各个附图中示出的连接线意欲表示各个元件之间的示例功能关系和/ 或物理耦合。应当注意，在本发明的实施例中可以呈现许多替换物或附加功 能关系或物理连接。
下面的描述提及元件或节点或特征是"连接"或者"耦接"在一起。如 此处所使用的，除非明确地相反限定，否则"连接"意思是一个元件/节点/ 特征直接结合到(或者与...直接通信)另一个元件/节点/特征，而不必机械地连接。同样地，除非明确地相反限定，否则"耦接"意思是一个元件/节点/ 特征直接或间接地结合到(或者与...直接或间接地通信)另一个元件/节点/
特征，而不必机械地连接。因此，尽管图1中图示的环境和图2中示出的示
意图描述了元件的示例结构，但是附加中间元件、部件、特征或组件可被呈 现在本发明的实施例中(假设电路的功能未被不利地影响)。
在大多数大气湍流中，称作乱流(turbule )的空气单元由于它们不同的 温度而以不同的比率上升和下降。(类似地，认定乱流是在冷空气海洋中上升 的暖空气的气泡)。湍气流区域包含许多各种尺寸的乱流。每个乱流就象镜头 那样(它当波穿过它时，它折射光或者RF波)影响电磁波。这种折射稍微 改变波的相位和行进方向。来自诸如卫星的远距离源的电f兹信号就象平面波 一样动作，直到波阵面到达湍气流区域为止在沿着基本平行的路径行进的 波阵面上的所有点。然而，在湍流区域内，每条信号路径会引起相位和方向 上的变化，所述变化与由相邻路径引起的变化不同。因此，从湍流区域涌现 的信号不再是平面波；相反，波阵面上的每个点以不同相位在稍微不同的方 向上行进。随着它们行进得越远，这些波元素就会彼此形成结构性的和破坏 性的干扰的区域。结构性的和破坏性的干扰的空间图案定义了强或弱信号强 度的位置、以及相关在前或在后信号相位的位置。信号强度和信号相位的空 间干扰图案可以通过在不同的位置处放置接收^L来测量。在乱流在湍流区域 内移动或者在诸如卫星的信号源相对于湍流区域移动时，空间干扰图案随着 时间变化。
为了说明目的，假设湍气流区域被单个大固定镜头代替，并且接收卫星 信号的平台也不移动。从坐在平台观察卫星穿过天空行进的观察员的视角来 看，来自单个卫星的卫星信号源看起来是横向行进，因为卫星本身具有横向 速度。当卫星移动到镜头位于平台和位置之间的位置时，平台不会看到卫星 也不能检测它的信号，除非卫星、镜头和平台排成一行，因此镜头将信号聚 焦在平台位置。(焦点对应于上述的干扰图案中的波峰)。后来，当卫星穿过 天空移动地更远时，镜头将卫星信号聚焦在不同的位置。(这对应于整个干扰 图案中的横向偏移)。
本发明的示例实施例确定与在从第一天线在第一位置到第二天线在第二 位置的干扰图案的接收传播中的时移对应的时间间隔，并且基于所述时间间 隔计算湍气流距离。图1是根据本发明示例实施例的环境100的示意表示图。环境100使用 位于平台114 (其可以是停飞或者飞行的飞机，飘浮或者运动的太空船、静 止或者运动的地面车辆、或者任何适当配置的结构)上的多个卫星接收天线 118/116 (前部/后部)来接收来自卫星102的信号。在这个示例中，干扰图案 108 (如在上面段落中所解释的)受湍气流区域106的影响，在第一时间值 120处穿过后部天线116和在第二时间值122处穿过前部天线118。如下描述 的伪测距(ranging)算法通过分析第一时间值与第二时间值之间的时间间隔 (时移)来估计湍气流区域106的距离126。
图2是根据本发明示例实施例配置的系统200的示意表示图。图2中描 绘的各个块能够以任意数目的在平台或系统200上到处放置的物理组件或模 块来实现。物理系统200可以包括大量电控制单元(ECU)、通信系统、机载 计算机系统、测量结构、网络和图2中所示组件以外的组件。这里将不详细 描述系统200的传统子系统、特征和各方面。
系统200通常包括多个天线(包括第一天线214和第二天线216)、接收 机220、包括相关器模块224的处理逻辑元件222、存储器模块226和通信网 络架构230,所述存储器模块226包括时间相关系数和相应相关滞后值228 的表或者与所述表通信。
第一天线214和第二天线216耦接到接收机220，并且被配置成提供如 下所述的湍流计算所需的数据。所述天线可以沿着平台的长度或宽度相距一 段距离丄。，和/或沿着平台的高度相距一段距离仏。每个天线响应于入射电磁 能而产生接收干扰图案，其中所述电磁能沿着一条须经湍气流区域106的路 径行进。
本发明的示例实施例要求最少两个天线。然而，根据平台设计，平台的 正常操纵可能导致卫星信号被阻挡或者遭受严重的多径衰落。在这点上，可 能需要替换位置的天线。所述平台可以包括但是不限于飞机飞行器、船、地 面车辆或者基于地面的结构。例如，对于飞行器，但是不限于，可以使用三 个天线 一个天线可位于垂直水平尾翼和/或翼尖，以及两个天线可沿着机身 放置以获得从特定卫星输入的多个信号。也可以添加另外的天线来提供附加 的数据源以补偿由于飞行器运动引起的信号的变化(例如，转向、爬升或者 下降)。卫星信号可以但是不限于GPS信号、GPS和GLONASS和伽利略信 号、或者GPS和iGPS信号。接收机220耦接到所述天线并且被配置成测量接收的干扰图案的属性。 通过一个算法来分析天线接收到的信号，所述算法将平台(例如飞行器)的 相对运动与卫星的相对运动进行比较。该比较算法利用了以下事实即，例 如GPS卫星的卫星广播一组代码，该组代码包含关于卫星的位置和速度的信 息、GPS系统时间、时钟误差、和数据的良好状态和精确度。所述信息允许 高精度地预测GPS卫星的相对位置和运动。例如，通常，飞行器机载的GPS 接收机将每个卫星的粗获取(C/A)代码与接收机数据库中包含的相同代码副 本进行匹配。通过转变接收机的代码以匹配卫星信号以及通过将这种转换与 卫星内部时钟相比较，接收机能够估计到卫星的距离。然而，通过使用多个 天线和此处描述的伪测距算法，也能够通过分析GPS信号的相位和时移来估 计湍气流区域106的3巨离。
另外，通过机载飞行管理计算机系统(图2中未示出)高精度地测量飞 行器的位置和相对运动。使用下面详细解释的伪测距算法，该数据连同来自 飞行器的卫星距离和卫星的速度随后被用来计算到湍气流区域106的距离。 接收机220可以而不限于GPS接收机、GPS和GLONASS和伽利略接收机、 或GPS和iGPS接收机。
处理逻辑元件222耦接到接收机220,并且可以包括任意数目的独特的 处理模块、相关器、或者被配置成执行此处更详细描述的任务、处理和操作 的组件。尽管图2中仅示出了一个处理块，但是现实的实现可以利用任意数 目的独特的物理和/和逻辑处理器，它们可以散布在系统200的各处。现实中， 可以使用被设计来执行此处所描述的功能的通用处理器、内容可寻址存储器、 数字信号处理器、特定用途集成电路、现场可编程门阵列、任何合适的可编 程逻辑器件、离散门或者晶体管逻辑、离散硬件组件或者它们的任何组合， 来实现或者执行处理逻辑元件222。处理器可被实现为微处理器、控制器、 微控制器或状态机。处理器也可被实现为计算器件的组合，例如，数字信号 处理器和微处理器、多个微处理器、 一个或多个微处理器结合数字信号处理 器核心、或者任何其它所述配置的组合。
相关器模块224耦接到处理逻辑元件222或者被实现在处理逻辑元件222
中，并且被配置成计算时间相关系数K^^和相应相关滞后^H并且选择相 关滞后值，如下面在图4的上下文中更详细解释的。相关器模块224可以包 括但是不限于并行处理相关器或者串行处理相关器。相关滞后值是在前部天线处的信号变化与后来在后部天线处的信号变化A之间的时间间隔，并且时 间相关系数是幅值/相位图案在每个时间滞后A在每个天线处匹配良好的测
量，如在下面的附图4-7的上下文中更详细解释的。
存储器模块226可以是具有适当数量的存储器的任何适当数据存储区， 被格式化来支持系统200的操作。存储器模块226被配置成存储、维护和提 供所需的数据，以便以下面描述的方式支持系统200的功能性。在实际的实 施例中，存储器模块226可被实现为RAM存储器、闪存、ROM存储器、EPROM 存储器、EEPROM存储器、寄存器、硬盘、可移动盘或者现有技术中已知的 任何其它形式的存储介质。存储器模块226可以耦接到处理逻辑元件222， 并且被配置成存储但不限于对应于下列测量的数据值平台位置和速度、卫 星的位置和速度、干扰图案的幅值/相位序列和取样间隔。另外地，存储器模 块226可以表示动态地更新包含时间相关系数和相应相关滞后228的表的数 据块，用以使用伪测距算法来计算到湍气流区域106的距离，如下面详细解 释的。
通信网络架构230可以耦接到处理逻辑元件222和测量结构(图2中未 示出)。通信网络架构230被配置成通信CAT数据并且报告CAT数据。CAT 数据是可被集成来改善实时报告和预测飞行天气两者的许多天气信息源之 一。经由通信网络架构230将CAT数据从CAT数据服务提供商发送到软件 应用程序。通信网络架构230可以但是不限于地面通信网络或者无线通信网 络，例如卫星通信网络。 一组中间件，例如系统宽信息管理(SWIM)系统， 经由通信网络架构230自动地读取来自许多源的数据，并且通过许多不同的 用户和应用程序重定格式使用所需的数据。SWIM允许集成大气天气简档信 息、飞行器状态配置、实时CAT跟踪和多普勒天气雷达。在这点上，系统200 可以向通信网络架构230的用户发出表示湍气流区域106的警告。
本发明的示例实施例基于下列信息来计算到湍气流区域的距离平台相 对于卫星速度矢量的被测速度矢量；卫星到平台114的距离124;卫星的速 率；和受湍气流区域106影响的干扰图案K.的速率/速度，如下面详细解释的。
参考图1,如果卫星102以速度Ks在距平台124的恒定距离124 (A) 处移动到右边，并且平台114和湍气流区域106静止，则从平台114到湍气 流区域106的距离126 (A)固定，并且唯一移动的物体是卫星102。在这点 上，干扰图案108看起来行进的水平速度F'(假设卫星信号在空间中是一条直线)通过下列关系来描绘:
可以通过重新排列上述等式(1)的项来获得到湍气流区域106的距离
凡
如等式(2)所示，湍气流区域106与平台114之间的距离i ,是已知距离 &、卫星的速度(两者都可以从卫星提供的导航数据中确定)以及由湍气 流区域106计算的卫星图像的未知速度F,的函数。值K.可以被确定为在前部 天线处的信号变化与随后后部天线处的信号改变之间的时间间隔A的函数， 如下面图3的上下文中详细所解释的。这些信号变化可以是信号幅值/相位、 信号相位、或者两者的变化。
图3是图示根据本发明示例实施例的用于定位湍气流区域的空中湍流位 置处理300的流程图。与处理器300相关执行的各个任务可以通过软件、硬 件、固件或者它们的组合来执行。为了图示目的，处理300的下列描述可以 参考结合图1和图2的上述元件。在实践的实施例中，处理300的各部分可 以通过系统200的下列不同元件来执行例如，第一天线214、第二天线216、 接收机220、处理逻辑元件222、相关器模块224、存储器模块226和通信网 络架构230。
处理300可以开始于获取测量数据用以基于相关滞后值来计算所述距 离(任务302)。例如，测量数据可以包括(但不限于)平台的水平速度、 卫星的水平速度、湍气流区域的水平速度、天线的水平间距、平台距卫星的 距离、平台的垂直速度、卫星的垂直速度、湍气流区域的垂直速度、和天线 的垂直间距。
如上所解释的，干扰图案108的未知速度值K如等式(1)中所述，可以 被确定为时间间隔A的函数。因此，随后，处理300确定与在第一天线和第 二天线处的干扰图案的接收传播中的时间偏移对应的时间间隔A (相关滞后 值)(任务304)，如在下面的附图4的上下文中详细描述的，并且继续基于 该时间间隔A计算距离i ,(任务306)。
此处描述的发明的各个示例实施例基于下列方案来计算距湍气流区域的 距离i ,:平台和/或湍气流区域是否固定和/或移动；平台是否处于卫星的近水平视线中；或者天线是否水平和/或垂直地分离，如下面更详细解释的。
在一个示例实施例中，参考图1，后部天线116和前部天线U8被水平分 离，并且平台114和湍气流区域106两者静止。所述天线可以在平台114上 相距距离￡0。在该示例中，干扰图案108在时间^ (用附图标记120表示) 穿过后部天线116,并且在时间间隔A之后，图像在时间^ (用附图标记122 表示)穿过前部天线118。
根据观察到的穿过两个天线的干扰图案之间的时间间隔A来计算干扰图 案速度K。干扰图案正在移动的速度基于下列关系来确定 ^ (3)
将等式(3)的^代入到上面的等式2，会得到一个作为时间间隔At的函 数的、用于计算距湍气流区域的距离/ ,的关系式。
因此，根据第一示例实施例，所述方法基于下列关系式计算距离《
i , ~^_ (4)
其中A是从静止的湍气流区域到静止的平台的距离，凡是静止的平台与
卫星之间的距离，z。是第一天线与第二天线之间的水平间距，J^是卫星的水 平速度，AZ是相关滞后。
在第二示例实施例中，如上所解释的，后部天线和前部天线被水平分离
并且湍气流区域静止，然而，平台是移动的(与上面的第一示例实施例不同)。 因此，需要考虑平台的移动来调整前面的方程式。在这种情况下，干扰图案 在时间^到达后部天线。前部天线距该时刻u)达一段距离丄。。然而，当干 扰图案(在时间G)到达前向天线的初始位置(在^时的位置)，第二前部天线 已经移动(图1中未示出)。图像最后在时间^ (图1中未示出)跟上前部天
线。这种情形下的定时通过下列等式(5)来描绘
在等式(5)中，K是平台的速度。因此干扰图案速度基于下列关系式来 计算
K=~^ (6) 将等式(6)代入到等式(2)得到距湍气流区域106的距离《。因此，根据第二示例实施例，所述方法基于下列关系式计算所述距离
<formula>formula see original document page 16</formula>其中凡是从静止的湍气流区域到移动平台的距离，A、 Z。和A"口上所解
释的，K是平台的水平速度。
第一和第二示例实施例都基于静止的紊气团来计算到湍气流区域的距离
《。然而，因为风可能移动湍气流区域，因此并不总是这种情况。
在第三示例实施例中，后部天线和前部天线如上所解释地水平分离，但 是不同于第一和第二示例实施例，湍气流区域和平台都移动。这种情形数学
地表示为下列等式(8):
<formula>formula see original document page 16</formula>将等式(8 )代入到等式(2 )得到距离A。
因此，根据第三示例实施例，所述方法基于下列关系式来计算距离i^:
<formula>formula see original document page 16</formula>
(9)
其中^是A^移动的湍气流区域到移动的平台的距离，A、 A)、 A 和r。如
上所解释的，K是湍气流区域的水平速度。
等式(9)给出距离A是时间间隔"和风速^的函数。当风速已知时，可 以照现在情形使用等式(9)。这可以是在良好测量的区域(类似北美、欧洲、 澳大利亚和靠近日本的地区)上的情况。如果已知风沿着从平台到卫星的矢 量而粗略地一致，则可以将值F,插入到等式(9)以使_得到距离凡。在风沿着 从平台到卫星的矢量平滑地变化(除湍流区中的风之外)的情况，使用三维 模型，可以使用第一猜想值^来精选可以估计F,处的位置。将该值^插入到 等式(9)以便估计新的值凡。所述处理重复，直到它在单个值A和r,上收敛。
得克萨斯大小的气团的水平运动通常以每秒几十米为单位进行测量。类 似气团的垂直运动通常以每秒几厘米为单位进行测量。在这点上，在几乎所 有情形中，假设湍气流区域具有为零的有效垂直速度是安全的。因此，垂直 分离的卫星天线可被用来测量干扰图案的垂直速度，相信，由于已知的卫星运动、已知的平台运动和从平台到湍流区域的未知距离，垂直速度几乎是完 整的。
在第四示例实施例中，与上面解释的第一、第二和第三示例实施例相反， 第一天线和第二天线被垂直地分离，并且平台和湍气流区域静止。
根据第四示例实施例，所述方法基于下列关系式计算距离i ,:
尺=
仏
(10)
其中A是从静止的湍气流区域106到静止的平台的距离，i ,和A,如上所 //。是第一天线与第二天线之间的垂直间隔，Rz是卫星的垂直速度。 在第五示例实施例中，第一天线和第二天线垂直分离并且湍气流区域静 止，如在上面第四实施例的上下文中所解释的，但是与第四实施例相反，在 这种情况下平台是移动的。
根据第五示例实施例，所述方法基于下列关系式来计算距离《 仏
(11)
<formula>formula see original document page 17</formula>
其中/ ,是从静止的湍气流区域到移动的平台的距离，A、仏、A/和Kz 如上所解释的，Fm是平台的垂直速度。
在第六示例实施例中，第一天线和第二天线垂直分离并平台是移动的， 如在上面第四和第五实施例的上下文中所解释的，但是与第四和第五实施例 相反，在这种情况下平台和湍气流区域都是移动的。
根据第六示例实施例，所述方法基于下列关系式来计算距离<formula>formula see original document page 17</formula>
(12)
其中7 ,是从移动的湍气流区域到移动的平台的距离，仏、△,、 Kz 和^如上所解释的，F,z是湍气流区域的垂直速度。
在第七示例实施例中，飞机到卫星的视线从水平线来看是几度。而且， 在第七示例实施例中，第一天线和第二天线垂直分离并且平台和湍气流区域 是移动的，如在上面第六实施例的上下文中所解释的。然而，与第六示例实 施例相反，参数受由于平台到卫星的视线51起的校正影响。因此，根据第七示例实施例，所述方法基于下列关系式来计算距离<formula>formula see original document page 18</formula>
其中A是从移动的湍气流区域到移动的平台的距离，A针对平台到卫星
的视线而被校正，凡是移动的平台与卫星之间的距离，F,—e/是移动的湍气流
区域投影到仰角单位(elevation unit)矢量的垂直速度，e,是投影到仰角单 位矢量的第一天线与第二天线之间的垂直间距，At是如上所解释的相关滞后，
F。—e/是投影到仰角单位矢量的移动平台的垂直速度，^e,是投影到仰角单位矢
量的卫星的垂直速度，并且e是飞机到卫星的矢量凡的仰角单位矢量E。 E
是在增加的仰角的方向上。当凡水平时，E是垂直的，但是当凡不水平时，
即，当从平台看来卫星不是刚好在水平面上时，E远离垂直地倾斜。E净皮定
义为力口下
<formula>formula see original document page 18</formula>
其中，￡、 ^和A是向量，A是z方向(垂直)上的单位矢量，"x，，表示 向量积，并且垂直条表示矢量幅值。为了计算任意矢量A到E的投影，使用 点积<formula>formula see original document page 18</formula>
例如，值e,在等式(13b)中可被计算为 <formula>formula see original document page 18</formula>注意，在平台到卫星矢量基本非水平的情况下，Kj不能被^暇设为零。 在计算距湍气流区域的距离之后，根据上述实施例，处理300随后发出
指示所述距离的警告(任务308)并且停止。
如上所述，湍流气团移动(shift)卫星信号的相位，并且因为气团具有 "凸凹不平的"折射率，因此相移根据信号的路径而改变。这将有序均匀的
卫星信号转换为干扰图案。通常，干扰图案是当从相同源的稍微不同路径行
进的两个或更多个信号被集合时产生的波峰和波谷的图案。该干扰图案随着 时间和位置而变化。前导天线检测由尾随(trailing)天线第一跟随的千扰图案， 并且确定^r测前导天线处和尾随天线处的干扰之间的时间差，以便计算距湍 气流区域的距离i ,，如下所解释的。图4是图示根据本发明示例实施例的空中湍流位置处理400的流程图。 与处理400相关执行的各个任务可以通过软件、硬件、固件或者它们的任何 组合来执行。为了图示目的，处理400的下列描述可以是指结合图1、图2 和图3的上述元件。在特定的实施例中，处理400的各部分可以通过系统200 的下列不同元件来执行例如第一天线214、第二天线216、 *接收机220、处 理逻辑元件222、相关器模块224、存储器模块226和通信网络架构230。
处理400可以开始于获得测量数据，用以基于如上面图3的上下文中 解释的相关滞后值来计算所述距离(任务402)。该示例实施例使用时间相关 技术来测量前导(前部)天线和尾随(后部)天线之间的时间差。来自前导 天线和尾随天线的信号的相关性在接收前导天线和尾随天线处的干扰图案中 (针对飞行器运动、风等校正了卫星图像速度)，在与由湍气流区域造成的时 移对应的值A 处产生相关波峰。本发明的该示例实施例对于许多值A/计算相 关系数K的，并且选择具有最高值的值(波峰相关值)标记为到最强烈湍流 的距离。该实施例使用在前导天线和尾随天线处的信号幅值测量来计算相关 系数。替换的实施例可以使用信号相位测量或者信号相位和信号幅值的组合 来计算相关系数。
当干扰图案速度^与平台速度几乎相同时，干扰图案从一个天线移动到 另一个天线的A,可以是几秒。最大相关系数KAO可以随着更长"而减小，因 为湍流图案随着时间变化。时间连贯性(coherence)描述了在时间的不同时 刻处观察的信号之间的相关性或者可预测的关系。优选实施例使用湍流的时 间连贯性间隔在短A,处缩小r(AO值，并且在长"处i文大KA,H直，以便在全部 距离处给出湍流强度的没有偏差的比较。
因此，处理400通过取样对应于后部天线处的干扰图案的第 一信号幅值/ 相位序列(任务406)和取样对应于前部天线处的干扰图案的第二信号幅值/ 相位序列(任务408)，继续计算相关滞后(附图标记404)。在该实施例中， 第一信号幅值/相位序列和第二信号幅值/相位序列是以相同的取样间隔取样 的。图5图示了在每个天线处取样的幅值图案的示例。干扰图案108的第一 信号幅值序列502在前部天线118处被取样，干扰图案108的第二信号幅值 序列504在后部天线116处被取样，在每个天线使用256个取样。
该示例中的卫星处于低环地轨道中。该示例中的平台是飞机。卫星的行 进方向远离飞机头部大约125。，即，卫星看起来是相对于飞才几向后运动。卫星的速度^是大约每秒7000米，但是在与从飞机到卫星的视线垂直的方向上 近大约4000 m/s。从飞机到卫星的距离i ,是大约4200 km。飞才几的实际空速r。 是大约200 m/s。前部天线和后部天线相距大约20米。
然后，处理400计算第一信号幅值序列520与第二信号幅值序列504之 间的时间相关系数KA^ (任务410),并且继续计算与该时间相关系数对应的 相关滞后值A"(任务412)。相关滞后值是在前部天线处的信号变化与后在来 后部天线处的信号变化之间的时间间隔Ar (如上面图1的上下文中所解释 的)，并且通过测量信号幅值到达(祐j企测)前部天线的时间和将前面的时间 与信号幅值到达后部天线的时间相减来计算该值。时间相关性描述了在时间 上的不同时刻处观察到的信号之间的相关性或可预测的关系。因此，在各个 天线处测量的信号幅值的统计协议确定在所计算的值A/处的时间相关系翁: (幅值匹配多好)。替换实施例可以使用信号相位测量、或者信号相位和信号 幅值的组合来计算时间相关系数。
处理400存储计算出的时间相关系数和相应的相关滞后值(任务414 )， 并且选择与计算出的时间相关系数的最高值对应的相关滞后值(任务416 )。 例如，在图6中，时间相关系数的最高值在大约4个取样602 (在相关系凄史 的波峰)处发生。这意味着如果在时间^到达前部天线的信号穿过特定位置的 湍气流区域，则相同位置将在时间"被到达后部天线的信号穿过。假定大约 0.05秒的取样间隔，这是大约0.2秒的相关滞后值。假设天线(丄。)之间相距大 约20米，这对应于朝向飞机后部的大约100m/s的相对图像速度^ 。飞机以 大约200 m/s ( F。)运动，因此真实图像速度^是前向大约100m/s。假设卫 星在大约4200 km的距离A处以大约400 m/s(^)向后运动，则大约100 m/s 的图像速度F;表示湍气流区域在距飞机大约105km的位置处被穿过。
图7图示了以千米计的距湍气流区域的各个距离的相关波峰(用附图标 记702-712表示)。如图7所示，波峰相关值可以发生在正值A,和负值A 两者 处。当干扰图案速度K为正并且比平台快时，即，干扰图案(湍气流区域影 响的幅值序列)在前导(前部)天线之前影响尾随(后部)天线，负值A,发 生。这可以使用iGPS来发生，其中在低且快速运动的轨道中的卫星生成GPS 信号。
特别地，如图7所示，相关波峰随着距离的增加(由附图标记708-712 表示)移动到逐渐增加的正时间间隔A Q人前部天线到后部天线)。然而，,爻定倒退的卫星，则在特定距离出现奇点(singularity )。当距离接近该奇点， 在这种情况下大约200 km (图7中未示出)，相关波峰快速地领先，然后弯 曲到负Aa。然后，它接近左边的原点(△/=())(由附图标记702-706表示)， 随着距离增加，渐近地接近原点。
随后处理400继续计算所述距离(任务418)，并且发出指示所述距离的 警告(任务420),如上面图3的上下文中所解释的。
上述的本发明的示例实施例确定了到湍气流区域的中心的距离。然而， 根据本发明的另一示例实施例，湍气流区域的厚度可以通过4企查相关波峰的 宽度来确定。因为所有的气团在相同距离处并且因此计算以相同速度移动的 图像， 一个薄区域具有一个高且窄的相关波峰(在A内较窄)。 一个厚区i或具 有一个更宽的波峰，因为一些气团比其它气团相距更远，因此图像具有许多 以稍微不同的速度移动的分量。对于给定数量的卫星信号闪烁，强度与所测 量的湍气流区域的厚度成反比。在这点上，如果时间相关曲线的波峰较宽(在 A内较宽)，则处理400指示厚湍气流区域具有相应较低的湍流强度，如果时 间相关曲线的波峰较窄，则处理400指示薄湍气流区域具有相应较高的湍流 强度。
此处描述的发明实施例可以通过更精确地测定湍气流区i或(包括CAT) 的位置(距离和高度)来提高空中交通系统性能，提供一种用于空中交通系 统提供商和操作员的方式以便更精确地计划避免湍气流区域以提高乘客在交 通或旅行期间的安全性和舒适性。
尽管在前面的详细描述中已经呈现了至少一个示例性实施例，但是应当 理解存在大量的变型。还应当理解， 一个或多个示例性实施例仅仅是示例， 并且决不用于以任何方式限制发明的范围、应用或者配置。相反，前面的详 细描述将向本领域的普通技术人员提供用于实现所述一个或多个示例性实施 例的便利指导。应当理解，在不背离所附权利要求书及其合法等效物所阐述 的发明范畴的情况下，可以在元件的功能和结构上进行各种变化。
权利要求
1. 一种用于计算从湍气流区域到平台的距离的方法，所述方法包括确定在从第一天线到第二天线的干扰图案的传播中的相关滞后值，其中所述相关滞后值对应于由所述湍气流区域造成的时移，并且其中所述湍气流区域处于从卫星发送到所述平台的信号的路径中；和基于所述相关滞后值而计算所述距离。
2. 如权利要求l的方法，其中确定相关滞后值的步骤进一步包括以一取样间隔在所述第 一天线处取样对应于所述干扰图案的第 一信号幅 值序列；以所述取样间隔在所述第二天线处取样对应于所述干扰图案的第二信号 幅^直序列；计算所述第 一信号幅值序列和所述第二信号幅值序列之间的时间相关系数；计算对应于所述时间相关系数的相关滞后值；和 选择与所述时间相关系数的最高值对应的相关滞后值。
3. 如权利要求1或2的方法，进一步包括发出指示所述距离的警告。
4. 如任意在前权利要求的方法，其中 所述湍气流区域是静止的湍气流区域； 所述平台是静止的平台；并且基于关系式i ,丄n计算距离，其中&是从静止的湍气流区域到所述平台的距离，A是静止的平台与卫星之间的距离，4是第一天线与第 二天线之间的水平间距，^;是卫星的水平速度，并且Af是相关滞后。
5.如权利要求1-3中任一项的方法，其中所述湍气流区域是静止的湍气流区域；所述平台是移动的平台；并且基于关系式^=-&、AZ.计算所述距离，其中&是从静止的湍气流区域到移动的平台的距离，尺是移动的平台与卫星之间的距离，丄。是第一天线与第二天线之间的水平间距，A,是相关滞后，r。是移动的平台的水平 速度，并且^是卫星的水平速度。
6.如权利要求1-3中任一项的方法，其中 所述湍气流区域是移动的湍气流区域； 所述平台是移动的平台；并且基于关系式《=-凡乙->计算所述距离，其中^是从移动的湍气流区域到移动的平台的距离，i ,是移动的平台与卫星之间的距离，^是湍 气流区域的水平速度，￡。是第一天线与第二天线之间的水平间距，A,是相关 滞后，^是卫星的水平速度，并且K。是移动的平台的水平速度。
7.如权利要求1-3中任一项的方法，其中所述湍气流区域是静止的湍气流区域；所述平台是静止的平台；并且基于关系式<formula>formula see original document page 0</formula>.计算所述距离，其中A是从静止的湍气流区域到静止的平台的距离，^是静止的平台与卫星之间的距离，仏是第一天线 与第二天线之间的垂直间隔，^是卫星的垂直速度，并且&是相关滞后。
8.如权利要求1-3中任一项的方法，其中 所述湍气流区域是静止的湍气流区域； 所述平台是移动的平台；并且仏基于关系式i ,.计算所述距离，其中是从静止的湍气流区域到移动的平台的距离，/ ,是移动的平台与卫星之间的距离，仏是第 一天线与第二天线之间的垂直间隔，M是相关滞后，f^是移动的平台的垂直 速度，并且^是卫星的垂直速度。
9.如权利要求1-3中任一项的方法，其中所述湍气流区域是移动的湍气流区域；所述平台是移动的平台；并且基于关系式i , = -A "~^计算所述距离，其中A是从移动的湍气流区域到移动的平台的距离，^是移动的平台与卫星之间的距离，F,z是移动 的湍气流区域的垂直速度，仏是第一天线与第二天线之间的垂直间隔，"是 相关滞后，^是移动的湍气流区域的垂直速度，P^是移动的平台的垂直速度， 并且^是卫星的垂直速度。
10.如权利要求1-3中任一项的方法，其中 所述湍气流区域是移动的湍气流区域； 所述平台是移动的平台；并且基于关系式i ,<formula>formula see original document page 4</formula>计算所述距离，其中/ ,是从移动的湍气流区域到平台的距离，^针对平台到卫星的视线而被校正，^是移动的 平台与卫星之间的距离，e/是移动的湍气流区域投影到仰角单位矢量的垂 直速度，e/是投影到仰角单位矢量的第一天线与第二天线之间的垂直间距， At是相关滞后，Fa—e/是投影到仰角单位矢量的移动的平台的垂直速度，并且e/是投影到仰角单位矢量的卫星的垂直速度。
11. 如任意在前权利要求的方法，进一步包括基于所述相关滞后值获取 用于计算所述距离的数据，所述数据包括所述平台的水平速度；卫星的水平速度；所述湍气流区域的水平速度；所述第 一天线和所述第二天线之间的水平间距；所述平台距所述卫星的距离；所述平台的垂直速度；所述卫星的垂直速度；所述湍气流区域的垂直速度；和所述第一天线和所述第二天线之间的垂直间距。
12. 如权利要求2-11中任一项的方法，进一步包括存储 所述第 一信号幅值序列和所述第二信号幅值序列；计算出的时间相关系数； 所述相关滞后值；和所述取才羊间隔。
13. 如权利要求2-12中任一项的方法，其中计算出的时间相关系数的 最高值处于时间相关曲线的波峰处。
14. 如权利要求13的方法，进一步包括基于所述时间相关曲线的波峰确 定所述湍气流区域的厚度。
15. 如权利要求14的方法，进一步包括如果所述时间相关曲线的波峰相对较宽，则指示厚湍气流区域；和 如果所述时间相关曲线的波峰相对较窄，则指示薄湍气流区域。
16. —种用于计算从湍气流区域到平台的距离的系统，所述系统包括 多个天线，每个天线响应于入射电磁能而产生接收干扰图案，其中所述电磁能沿着一条须经过所述湍气流区域的路径行进；耦接到所述天线的接收机，其被配置成测量所述接收干扰图案的属性；和耦接到所述接收机的处理逻辑元件，其被配置成确定在从第 一天线到第二天线的接收干扰图案的传播中的相关滞后 值，其中所述相关滞后值对应于由所述湍气流区域造成的时移；和 基于所述相关滞后值而计算所述距离。
17. 如权利要求16的系统，其中所述处理逻辑元件被进一步配置成以 一取样间隔在第 一天线处取样对应于所述接收干扰图案的第 一信号幅 值序列；和以所述取样间隔在第二天线处取样对应于所述接收干扰图案的第二信号 幅^直序列。
18. 如权利要求17的系统，进一步包括耦接到所述处理逻辑元件的相关 器模块，该相关器模块被配置成计算所述第 一信号幅值序列和所述第二信号幅值序列之间的时间相关系数；计算与计算出的时间相关系数对应的相关滞后值；和 选择与计算出的时间相关系数的最高值对应的相关滞后值。
19. 如权利要求18的系统，进一步包括耦接到所述处理逻辑元件的存储器模块，该存储器模块被配置成存储所述第 一信号幅值序列和所述第二信号幅值序列；计算出的时间相关系数；所述相关滞后值；和所述取样间隔。
20. 如权利要求16-19中任一项的系统，其中所述接收机包括卫星接收机
21. 如权利要求16-20中任一项的系统，进一步包括耦接到所述处理逻 辑元件的通信网络，所述通信网络被配置成发出指示所述距离的警告。
22. 如权利要求16-21中任一项的系统，其中所述处理逻辑元件被进一 步配置成以 一取样间隔在第 一天线处取样对应于所述接收干扰图案的第 一信号相 位序列；和以所述取样间隔在第二天线处取样对应于所述接收干扰图案的第二信号 相位序列。
23. 如权利要求18-22中任一项的系统，其中所述相关器模块被进一步 配置成计算所述第 一信号相位序列和所述第二信号相位序列之间的时间相关系数；计算与计算出的时间相关系数对应的相关滞后值；和 选择与计算出的时间相关系数的最高值对应的相关滞后值。
全文摘要
公开了一种用于计算从天线平台到湍气流区域的距离的系统和方法。所述方法，确定在响应于入射电磁能的天线处沿着经湍气流区域的路径行进的干扰图案的接收传播中的、与时移相对应的相关滞后值。基于所述相关滞后值，所述方法计算从天线平台到湍气流区域的距离，并且发出指示该距离的警告。
文档编号G01W1/00GK101535846SQ200780040788
公开日2009年9月16日 申请日期2007年10月23日 优先权日2006年10月31日
发明者布赖恩·蒂洛森, 马里萨·辛格尔顿 申请人:波音公司