JP6527745B2 - 航空通信用周波数推定器 - Google Patents

Description

本発明は、例えばＤＶＢ−Ｓ２標準またはその均等な様式に基づき航空条件下でデータをリアルタイムに送受信する周波数推定器に関する。本発明は、低信号対雑音比状況下かつ高ドップラー加速度で、ドップラー効果に起因する振動数シフトを充分に正確な方法で推定するために使用され得る。

本発明は、特に低信号対雑音比と当該分野で高いと考えられる値に達し得るドップラー加速度とを示す空中状況下における周波数同期の問題に関する。ドップラー効果は衛星に対する航空機の移動により引き起こされ、航空チャネルはドップラー効果と呼ばれる受信搬送波における周波数シフトにより乱される。

図１に、地上ステーション１と衛星２と空中システム３または航空機とを有する送信システムを概略的に示す。送り経路Ｌ_ａは、地上ステーション１から航空機などの空中システムへ行く経路として定義され、戻り経路Ｌ_ｒは衛星２を介し航空機３から地上ステーション１へ行く経路として定義される。ドップラー効果は衛星２に対する航空機３の移動により引き起こされる。ドップラー効果およびドップラー効果変動を航空機の速度と加速度に関連づけることが可能である。

送信システムの周波数をｆ_ｅｍ、送信システムの速度をｖ_ｅｍ、受信器の周波数をｆ_ｒｅｃ、受信器の速度をｖ_ｒｅｃとすると、２つの周波数間の関係は光の速度ｃの関数として次式のように表され得る。

衛星は航空機に比較して動かないと考えられるので、送信機と受信器間の周波数差Δｆに対し次式が得られる。

最大ドップラー効果Ｄ_ｍａｘとドップラー効果の最大変化量Ｖ_ｍａｘと、航空機の最大速度Ｓ_ｍａｘと航空機の最大加速度Ａ_ｍａｘとを以下のように関連づけることが可能である。

ここで考えられる「最悪ケース」航空条件は以下のものである。
・旅客機の巡航速度に対応する航空機の最大速度０．９７Ｍａｃｈ、
・次の期間中に旅客機が達し得る最大加速度に対応する航空機の最大加速度２ｇ、
○衛星の方向の離陸と着陸、
○衛星の方向の半径方向加速度での急カーブ、
○航空機に対し衛星が頂点をなすエアポケット。

図２Ａと図２Ｂにおける２つのグラフは、空中ポケットの存在下の加速度Ａおよび高度の変動Ａｌｔと、ドップラー効果Ｄおよびドップラー効果の変動（導出曲線）Ｖとをそれぞれ示す。

無線周波数ＲＦ伝送特性は次の通りである。
・通信衛星において使用される最高周波数範囲Ｋａ帯の典型的な３０ＧＨｚに固定された最大搬送波周波数、
・Ｋａ帯における衛星通信に使用される最小シンボルレートに対応する１Ｍｂａｕｄに固定された最小シンボルレート。

最大ドップラー効果Ｄ_ｍａｘとドップラー効果の最大変化量Ｖ_ｍａｘに対し次の値が定められる。
Ｄ_ｍａｘ＝３３ｋＨｚ、およびＶ_ｍａｘ＝１９６２（Ｈｚ／ｓ）

フレームのシンボルレートＲ_ｓに関するこれらの値を正規化することにより、正規化最大ドップラー

と正規化ドップラー変動

が得られる。

図３に、９０シンボルのヘッダ３００、１４４０データシンボルの第１のブロック３０１_１、続いて３６シンボルのパイロットの第１のブロック３０２_１、次いで１４４０データシンボルの第２のブロック３０１_２、続いて３６シンボルのパイロットの第２のブロック３０２_２等からなるＤＶＢ−Ｓ２フレームの構造を示す。

航空状況下では、Ｋａ帯において、（非特許文献１）に推奨された周波数同期機構に関連付けられたＤＶＢ−Ｓ２フレームの標準化構造は、５ｄＢ（当該分野の典型的値）未満の信号対雑音比の低シンボルレート（５Ｍｂａｕｄ未満）の運用を許可しない。

このアプリケーションの例では、ドップラー推定後の最大残留ドップラーは、位相曖昧性を防止するために２つのパイロットブロック間にπの最大位相シフトを引き起こさねばならない。これは、１４４０データシンボルと３６パイロットシンボルのＤＶＢ−Ｓ２フレームでは次式に等しい最高周波数解像度

に関与する。

図４に、受信器における従来技術による周波数推定の例を概略的に示す。信号の周波数同期は、２段階、すなわちループ構造Ｉにより行われる第１の粗同期と開ループまたは「フィードフォワード」制御モードＩＩにおいて行われる第２の密同期で発生する。受信器により受信された入力信号は、装置の動作中の周波数の補正を可能にする推定周波数値もまた受信するミキサ４００に入る。信号は、ナイキストフィルタブロック４１０と、レート同期モジュール４２０と、フレーム同期モジュール４３０と、第１の周波数同期モジュール４４０とを含むアウトフィットへ通され、第１の周波数同期モジュール４４０の出力は、第１に、位相取得モジュール４７０が続く密同期モジュール４６０へ、第２に、ループフィルタ４６０へ接続される。

ループ同期機構Ｉは周波数同期モジュール４４０とループフィルタのモジュール４５０とを含む。周波数同期モジュールは次式に従って信号

の周波数の推定を行う。

ここで、

であり、ここで、ｒ_ｋは時点ｋに受信されたデータサンプルｒ、

はこの同時点ｋにおける基準シンボルの共役、ａｒｇは複素数の引き数である。雑音に対し非常に敏感であるこの推定は、一次ループフィルタによりフィルタ処理された後、ナイキスト受信フィルタの入力補正として導入される。

ループフィルタのループ帯域幅は第１の周波数同期処理における決定パラメータであり、
− 同期の収束速度に比例し、
− 同様に、雑音に対する推定の感度に比例する。

要約すれば、この推定は、ループ帯域が広くなればなるほど、より急速に周波数の変位へ収束するが、同時に、雑音レベルにより敏感となる。

次の密同期ブロック４６０は、監視モードで動作するアルゴリズムを使用することにより信号の周波数の第２の推定を行う。その推定を行うためにＤＶＢ−Ｓ２フレーム（ヘッダとパイロット）の基準フィールドを使用する。

念のために記しておくと、次式Ｒ_ｌ（ｍ）により、一連の大きさＮにわたる指標ｍの信号ｘの自己相関Ｒ（ｍ）を推定ことが可能である。

当業者に知られたＬｕｉｓｅ＆Ｒｅｇｇｉａｎｉｎｉアルゴリズムからの監視アルゴリズムは、これらの相関をパイロットブロックの数Ｌにわたって平均化する。相関は、パイロットブロックの長さの半分

にわたって行われる、ここで、Ｌ_０はパイロットブロックの長さである。

ここで、
Ｔ_ｓ：シンボル時間
Ｌ_０：パイロットブロックの長さ
ｉ：自己相関の推定の指標
Ｎ：Ｌ_０／２：パイロットブロックにより生成される自己相関の数
ｌ：パイロットブロックの指標
Ｌ：周波数推定が行われるパイロットブロックの数である。

この周波数補正は、相関を平均化するパイロットブロックの数Ｌの値に関する第２の妥協案を生じる。第２の妥協案は推定の精度に比例し、推定の反応度に反比例する。

しかし、このようなスキームでは、特に航空チャネルの周波数変動に追随できるようにするとともに低信号対雑音比０ｄＢでも十分に正確な周波数推定を可能にする周波数同期の調整を見い出すことができない。推定周波数

と実際の周波数ｆ_ｒｅａｌとの周波数差が現れ得、通信中にドロップアウトと従ってフレーム喪失との源となり得る。

公知の従来技術の推定器は通常、ＤＶＢ−Ｓ２フレームのヘッダとパイロットブロックなどの既知の情報を使用する監視アルゴリズムに基づくが、同監視アルゴリズムは、ＤＶＢ−Ｓ２標準に従って動作するのに十分に正確なドップラーの推定と、航空機の移動に起因する著しいドップラー変動に耐えるための十分な反応度との両立を可能にしない。

（非特許文献２）では、長期にわたって平均化されるループＤＡ推定器について説明している。

（特許文献１）では、パイロットを使用してブラインド周波数を推定するアルゴリズムについて説明している。

米国特許出願公開第２００８／０２１１７１９号明細書

ＥＴＳＩ ｄｉｒｅｃｔｉｖｅｓ ＴＲ １０２ ３７６ Ｖ１．１．１，"Ｔｅｃｈｎｉｃａｌ Ｒｅｐｏｒｔ，ＤＶＢ，Ｕｓｅｒ Ｇｕｉｄｅｌｉｎｅｓ ｆｏｒ ｔｈｅ ｓｅｃｏｎｄ ｇｅｎｅｒａｔｉｏｎ ｓｙｓｔｅｍ ｆｏｒ Ｂｒｏａｄｃａｓｔｉｎｇ，Ｉｎｔｅｒａｃｔｉｖｅ Ｓｅｒｖｉｃｅｓ，Ｎｅｗｓ Ｇａｔｈｅｒｉｎｇ ａｎｄ ｏｔｈｅｒ Ｂｒｏａｄｂａｎｄ Ｓａｔｅｌｌｉｔｅ Ａｐｐｌｉｃａｔｉｏｎｓ ＤＶＢ−Ｓ２" ＲＹＵ ｅｔ ａｌ． ｅｎｔｉｔｌｅｄ"Ｈａｒｄｗａｒｅ ｅｆｆｉｃｉｅｎｔ ｆｒｅｑｕｅｎｃｙ ｅｓｔｉｍａｔｏｒ ｂａｓｅｄ ｏｎ ｄａｔａ−ａｉｄｅｄ ａｌｇｏｒｉｔｈｍ ｆｏｒ ｄｉｇｉｔａｌ ｖｉｄｅｏ ｂｒｏａｄｃａｓｔｉｎｇ ｓｙｓｔｅｍ"

本発明の主題である装置は、特に第１に、

の精度を有する周波数の概算見積り

を導出するために、既知データのヘッダに対し監視周波数推定器を使用する。第２に、本装置は、第１の工程中に判断された推定周波数

を中心とする範囲に対応する探索範囲にわたって未知のデータのすべてに対しブラインドアルゴリズムを実行することにより周波数を正確に推定することになる。これは、監視推定の頑強性とブラインド推定の精度との組み合わせを可能にし、極めて正確な周波数の推定が極短い推定期間にわたって得られるようにし、したがって、航空機の移動に起因する非常に大きな変動に追随できるようにするので有利である。

表現「ブラインドアルゴリズム」は、本明細書では、本来通信フレームで送信される未知のデータを使用するアルゴリズムを示す。監視アルゴリズムは既知のデータに作用することになる。

本発明は、第１のステーションと第１のステーションに対して移動する空中システムとの間の航空通信の周波数値を推定する方法に関わり、データは一連のフレームの形式で送信され、フレームは既知のデータの少なくとも１つの第１のヘッダフィールドとそれに続く未知のデータの１つまたは複数のフィールドとを含み、本方法は少なくとも、
ヘッダの既知データに対する監視相関と、

の精度を有する推定周波数

を中心とした第１の周波数範囲の推定とを行う工程１と、
工程１において判断された周波数

を中心とした範囲にわたって相関ピークを求めることにより、および相関ピークに対応する周波数を保持することにより、データフィールドの未知のデータの少なくともすべてに対するブラインド相関を生成する工程２とを含む。

変形実施形態によると、フレームは、データフィールドの前に配置されるパイロットフィールドを含む位相変化による変調により変調されるＤＶＢ−Ｓ２タイプのフレームであり、本変形実施形態の方法は少なくとも、

の精度を有する推定周波数

を中心とした範囲を得るためにＤＶＢ−Ｓ２ヘッダのフーリエ変換により監視相関を行う工程１と、
変調フレームのシンボルを使用し、工程１で定義された周波数

を中心とした範囲に対応する短フレームの変調シンボルのフーリエ変換の４乗を使用し、フーリエ変換のノルムの最大値に対応する周波数を判断しそこからドップラー効果の値を推定する工程２とを含む。

本方法は、
周波数

を推定する工程１であって、

ここで、
− ｃｏｎｊ（）：複素数の共役関数
− ＤＦＴ：フーリエ変換
− ｋ：シンボルの指標
− ｔ：所与の時点
− Δ_ｆ：推定すべきドップラー周波数
− ｒｅｃｔ_９０（ｔ）：長さ９０の直交関数
− ｓｉｎｃ_９０＊π（ｆ）：幅９０^＊πの基本正弦関数
− Ｒｅｆ（ｋ）：ヘッダの受信シンボルＲ（ｋ）の基準シンボル、ここで、１≦ｋ≦９０
− ＢＡＢＧ：付加白色ガウス雑音
である工程１と、
次式からドップラーの値を推定する工程２であって、

ここで、
− ＮＤＡＬｅｎｇｔｈ：ブラインド周波数推定値を生成する長さ、
− Ｒｅｃｔ_{ＮＤＡＬｅｎｇｔｈ}（ｔ）：長さＮＤＡＬｅｎｇｔｈの直交関数、
− Ｓｉｎｃ_{ＮＤＡＬｅｎｇｔｈ＊π}（ｆ）：幅ＮＤＡＬｅｎｇｔｈ^＊πの基本正弦関数、
− ＢＡＢＧＥｑ（ｆ）：雑音のある信号ＳｙｍｂｓＱＰＳＫ＋ＢＡＢＧの４乗から発生する等価付加白色ガウス雑音、
− ＭａｘＩｎｄｅｘＮＤＡ：フーリエ変換のノルムが最大となる周波数範囲

に属する周波数、
である工程２と
を含む。

一例として、フレームはＢＰＳＫまたはＱＰＳＫ変調により変調され、ＮＤＡＬｅｎｇｔｈの値は０ｄＢのＳＮＲに対して８３７０に固定される。

本方法の工程は

に等しいＮＤＡＬｅｎｇｔｈへ適用される。ここで、ＳＮＲＬｉｎは線型で表された信号対雑音比ＳＮＲに対応する。

一実施形態によると、周波数推定は通信を構成するＤＶＢ−Ｓ２フレーム毎に判断される。

本発明はまた、少なくとも第１のステーションと第１のステーションに対して移動する空中システムとを含む航空通信システムにおけるドップラーを推定するための装置に関わり、データは一連のフレームの形式で送信され、フレームは既知のデータの少なくとも１つの第１のヘッダフィールドとそれに続く未知のデータの１つまたは複数のフィールドとを含む。本装置は少なくとも、
・ヘッダの既知のデータに対する監視相関を生成するとともに、

の精度を有する推定周波数

中心とした第１の周波数範囲の推定を行う第１のモジュールと、
・工程１において判断された周波数範囲

にわたって相関ピークを求めることにより、および相関ピークに対応する周波数を保持することにより、データフィールドの未知のデータの少なくともすべてに対するブラインド相関を生成する第２のモジュールとを含む。

データはＱＰＳＫまたはＢＰＳＫ変調ＤＶＢ−Ｓ２データである。

第１の推定モジュールと第２のモジュールは例えばＦＰＧＡ技術を利用して作製される。

本発明の他の特徴と利点は、例示として提供され決して制限するものではない、図面が添付された以下の説明を読むとより良く理解される。

送信システムの例を示す。 空中ポケットの存在下の加速度Ａと高度の変動Ａｌｔを例示する。 ドップラー効果Ｄとドップラー効果の変動すなわち導出曲線Ｖを例示する。 ＤＶＢ−Ｓ２フレームを示す。 従来技術による同期構造の例を示す。 ＤＶＢ−Ｓ２受信器に配置された本発明による周波数推定器の構造の例を示す。

以下に続く例では、本発明の主題についてより良く説明するために、そして本発明の主題の範囲をＤＶＢ−Ｓ２フレームの構造に限定すること無しに、バースト的周波数推定またはそうでなければＱＰＳＫ（４位相シフト変調：Ｑｕａｄｒａｔｕｒｅ Ｐｈａｓｅ Ｓｈｉｆｔ Ｋｅｙｉｎｇ）位相変化による変調を有するＤＶＢ−Ｓ２フレームの場合の周波数推定について記載する。本発明はまた、使用されるフレームが、既知データのヘッダまたはフィールドと、これに続いて複数の未知データのフィールドとを有する通信に適用され得る。一例として、未知データはＢＰＳＫ（バイナリ位相変調方式：Ｂｉｎａｒｙ Ｐｈａｓｅ Ｓｈｉｆｔ Ｋｅｙｉｎｇ）変調またはＱＰＳＫ変調される。以下の詳細な例では、本システムは、既知Ａｒｉｎｃ７９１規準に準拠するために例えば０〜６ｄＢの信号対雑音比で動作することになる。

本発明による方法と周波数推定器は第１の監視周波数推定器と第２のブラインド周波数推定器とを有する。本発明による周波数推定器は、一例として衛星ステーションまたは地上ステーションに配置され例えばプログラマブル回路またはＦＰＧＡ（フィールドプログラマブルゲートアレイ）技術を使用するＤＶＢ−Ｓ２受信器内に実装される。

第１の周波数推定モジュールまたは第１の推定器は、９０シンボルのＤＶＢ−Ｓ２ヘッダのフーリエ変換を利用して監視相関を生成するアルゴリズムを実施する。

本アルゴリズムは、ＤＶＢ−Ｓ２ヘッダから受信されこの同じヘッダの基準シンボル（Ｒｅｆ（ｋ）、１≦ｋ≦９０）に相関付けられた９０シンボル（Ｒ（ｋ）、１≦ｋ≦９０）のフーリエ変換のノルムの最大値に対応する周波数を計算することにより、周波数を推定する。

受信シンボルＲ（ｋ）は次式の通りである。
− Ｒ（ｋ）＝（Ｒｅｆ（ｋ）＋ＢＡＢＧ（ｋ））^＊ｅｘｐ（ｉ^＊２^＊πΔ_ｆ ^＊ｋ）
ここで、
○Ｒｅｆ（ｋ）：シンボルＲ（ｋ）の基準シンボル
○ＢＡＢＧ：付加白色ガウス雑音
○Δ_ｆ：推定すべきドップラー周波数

第１の推定モジュールまたは推定器における周波数推定

は次式で表される。

ここで、
− ｃｏｎｊ（）：複素数の共役関数
− ＤＦＴ：フーリエ変換
− ｋ：シンボルの指標
− ｔ：所与の時点
− ｆ：所与の周波数
− Δ_ｆ：推定すべきドップラー周波数
− ｒｅｃｔ_９０（ｔ）：長さ９０の直交関数
− ｓｉｎｃ_９０＊π（ｆ）：幅９０^＊πの基本正弦関数

推定器は監視推定器であるので動作範囲は＋／−０．５^＊Ｒｓ以下である。

雑音に対する高い耐性を呈示する第１の推定モジュールの結果、周波数値

が得られる。

第２の推定モジュールまたは第２の推定器は、次の周波数範囲に作用する極高い周波数分解能を有するブラインド周波数推定器である。

最大ピークを求めるこの周波数範囲は監視周波数推定器の精度により判断され、したがって、正規化モードでは

を中心に±０．００２５Ｒ_ｓ（１ＭｂａｕｄのＲ_ｓに対して２．５ｋＨｚ）であり、これは０ｄＢにおける監視推定器の最悪の場合の精度に対応する。したがって、０ｄＢでは、

ＳＮＲが０ｄＢである場合、ブラインドアルゴリズムはパイロット（すなわち８３７０シンボル）を有する短ＱＰＳＫフレームのシンボルをすべて使用する。本アルゴリズムは、このフーリエ変換のノルムの最大値に対応する周波数ｆを計算することによりドップラーを推定するために、０ｄＢにおいて、８３７０ＱＰＳＫシンボルの４乗のフーリエ変換を使用することを含む。ドップラーは、４で除したこの推定周波数に対応する。受信ＱＰＳＫシンボルＲ（ｋ）を考慮すると、ここで、１≦ｋ≦ＮＤＡＬｅｎｇｔｈ、
− Ｒ（ｋ）＝（ＳｙｍｂｓＱＰＳＫ（ｋ）＋ＢＡＢＧ（ｋ））^＊ｅｘｐ（ｉ^＊２^＊πΔｆ^＊ｋ）、ここで、
○ＳｙｍｂｓＱＰＳＫ＝ｅｘｐ（ｉ^＊（π／４＋ｎ^＊π／２））、ここで、ｎ＝［０，１，２，３］、
○ＮＤＡＬｅｎｇｔｈ：ブラインド周波数推定値を生成する長さであり、ＮＤＡＬｅｎｇｔｈは、ＳＮＲが０ｄＢである場合８３７０に固定される、
○ＢＡＢＧ：付加白色ガウス雑音
○Δ_ｆ：推定すべきドップラー周波数

第２の推定モジュールまたは推定器における周波数推定

は次式で表される。

ここで、
− Ｒｅｃｔ_{ＮＤＡＬｅｎｇｔｈ}（ｔ）：長さＮＤＡＬｅｎｇｔｈの直交関数、
− ＳｉｎｃＮ_{ＤＡＬｅｎｇｔｈ＊π}（ｆ）：幅ＮＤＡＬｅｎｇｔｈ^＊πの基本正弦関数、
− ＢＡＢＧＥｑ（ｆ）：雑音のある信号ＳｙｍｂｓＱＰＳＫ＋ＢＡＢＧの４乗から発生する等価付加白色ガウス雑音、
− ＭａｘＩｎｄｅｘＮＤＡ：フーリエ変換のノルムが最大となる周波数範囲ＦｒｅｑＲａｎｇｅＮＤＡに属する周波数。

推定器はその入力として信号の４乗を取るブラインド推定器であるので、最大動作範囲は＋／−０．１２５^＊Ｒ_ｓである。

図５は、一例として次の工程を行う２つのモジュールを有する地上ステーションのＤＶＢ−Ｓ２受信器に実装される本発明による周波数推定器の構造を概略的に示す。
・第１に、周波数の粗監視推定は９０シンボルのＤＶＢ−Ｓ２ヘッダに対して行われ、
＋／−２．５ｋＨｚ（１ＭｓｐｓのＲｓ）の精度を有する第１の周波数推定が得られ、
・第２に、ブラインド相関は、粗周波数推定から＋／−２．５ｋＨｚの範囲にわたる短いＱＰＳＫフレームの８３７０シンボルに対して行われる。次に、周波数誤差２０Ｈｚの精度を有する推定が得られる。

換算周波数範囲にわたるブラインド推定器の使用は、ブラインド推定器に起因するＦＥＲを、推定器の全動作範囲とこの換算周波数範囲との比で除算できるようにする。したがって、この場合、推定器の動作範囲は＋／−１２５ｋＨｚであるので、ＦＥＲは５０（１２５／２．５＝５０）で除される。これにより、０ｄＢにおいて１０^−５を越えるＦＥＲから１０^−６近傍のＦＥＲへの変化を可能にする。

本発明に従って実現される周波数推定器は、悪航空条件下の０ｄＢにおいて１０^−５未満のＦＥＲを達成できるようにする。０ｄＢを越える信号対雑音比値ＳＮＲに対し、ブラインド推定の大きさを低減する一方で同じレベルの性能を得ることが可能である。

上に説明された本方法の工程は、周波数が推定される長さが

に対応する場合に適用され、（ここで、ＳＮＲＬｉｎは線型で表現された信号対雑音比ＳＮＲに対応する、性能（最悪の場合の航空ドップラー条件下で１０^−５未満のＦＥＲ（周波数誤差レート））はＤＶＢ−Ｓ２ ＱＰＳＫの典型的なＳＮＲ（０〜７ｄＢ）と同じである。

本発明は特に、極短い推定期間にわたって極めて正確に周波数の推定を提供する利点を有し、したがって、その後の航空機の移動に起因する周波数の非常に大きな変動を許容する。したがって、本発明は完全に「バースト的」な方法で使用され得る、すなわち、ＤＶＢ−Ｓ２フレーム毎に周波数推定が提供され、これは、前フレームにわたって得られる推定とは独立している。したがって、フレームにわたる周波数推定に誤りがあっても、他のフレームに影響を与えなく、その損出は誤った推定のフレームに限定される。

１ 地上ステーション
２ 衛星
３ 空中システム
３００ ヘッダ
３０１_１ １４４０データシンボルの第１のブロック
３０２_１ ３６シンボルのパイロットの第１のブロック
３０１_２ １４４０データシンボルの第２のブロック
３０２_２ ３６シンボルのパイロットの第２のブロック
４００ ミキサ
４１０ ナイキストフィルタ
４２０ レート同期モジュール
４３０ フレーム同期モジュール
４４０ 周波数同期モジュール
４５０ ループフィルタ
４６０ 密同期モジュール
４７０ 位相取得モジュール
５１０ ９０シンボルのヘッダに対する監視相関
５２０ ８３７０シンボルに関するブラインド相関

Claims (11)

1. 第１のステーション（１）と前記第１のステーションに対して移動する空中システム（３）との間の航空通信の周波数値を推定する方法であって、データが一連のフレームの形式で送信され、フレームが既知のデータの少なくとも１つの第１のヘッダフィールド（３００）とそれに続く未知のデータの１つまたは複数のフィールド（３０１_１）とを含み、前記第１のステーションと前記空中システムとのうちの少なくとも１つが、無線通信受信器を備える方法において、前記方法は、前記無線通信受信器によって実行される少なくとも以下の工程、
   前記無線通信受信器による、前記ヘッダの前記既知のデータ（３００）に対する監視相関と、前記無線通信受信器による、
   

   

   の精度を有する推定周波数
   

   

   を中心とした第１の周波数範囲の推定とを行う工程１と、
   前記工程１において判断された前記周波数
   

   

   を中心とした範囲にわたって相関ピークを求めることにより、および前記相関ピークに対応する周波数を保持することにより、前記無線通信受信器によって、少なくとも前記データフィールド（３０１_１）のすべての未知のデータに対するブラインド相関を生成する工程２とを含み、前記保持される周波数は、前記第１のステーションと前記空中システムとの間でデータを送信するために使用される、方法。
2. 前記フレームは、前記データフィールド（３０１_１）の前に配置されるパイロットフィールド（３０２_１）を含む位相変化による変調により変調されるＤＶＢ−Ｓ２タイプのフレームであることと、前記方法は、少なくとも、
   

   

   の精度を有する推定周波数
   

   

   を中心とした範囲を得るために、前記ＤＶＢ−Ｓ２ヘッダのフーリエ変換により監視相関を行う工程１と、
   前記変調フレームのシンボルを使用し、前記工程１で定義された前記周波数
   

   

   を中心とした範囲に対応する短フレームの変調シンボルの４乗からのフーリエ変換を使用し、前記フーリエ変換のノルムの最大値に対応する周波数を判断しそこからドップラー効果の値を推定する工程２とを含むこととを特徴とする、請求項１に記載の方法。
3. 次式の周波数
   

   

   を推定する工程１であって、
   

   

   ここで、
   − ｃｏｎｊ（）：複素数の共役関数
   − ＤＦＴ：フーリエ変換
   − ｋ：シンボルの指標
   − ｔ：所与の時点
   − Δ_ｆ：推定すべきドップラー周波数
   − ｒｅｃｔ_９０（ｔ）：長さ９０の直交関数
   − ｓｉｎｃ_９０＊π（ｆ）：幅９０^＊πの基本正弦関数
   − Ｒｅｆ（ｋ）：ヘッダの受信シンボルＲ（ｋ）の基準シンボル、ここで、１≦ｋ≦９０
   − ＢＡＢＧ：付加白色ガウス雑音
   である工程１と、
   次式からドップラーの値を推定する工程２であって、
   

   

   ここで、
   − ＮＤＡＬｅｎｇｔｈ：ブラインド周波数推定値を生成する長さ、
   − Ｒｅｃｔ_{ＮＤＡＬｅｎｇｔｈ}（ｔ）：長さＮＤＡＬｅｎｇｔｈの直交関数、
   − ｓｉｎｃ_{ＮＤＡＬｅｎｇｔｈ＊π}（ｆ）：幅ＮＤＡＬｅｎｇｔｈ^＊πの基本正弦関数、
   − ＢＡＢＧＥｑ（ｆ）：雑音のある信号ＳｙｍｂｓＱＰＳＫ＋ＢＡＢＧの４乗から発生する等価付加白色ガウス雑音、
   − ＭａｘＩｎｄｅｘＮＤＡ：フーリエ変換のノルムが最大となる周波数範囲ＦｒｅｑＲａｎｇｅＮＤＡに属する周波数
   である工程２と
   を少なくとも含むことを特徴とする、請求項２に記載の方法。
4. 前記フレームはＢＰＳＫまたはＱＰＳＫ変調により変調されることを特徴とする、請求項２または３に記載の方法。
5. ＮＤＡＬｅｎｇｔｈの値は０ｄＢのＳＮＲに対して８３７０に固定されることを特徴とする、請求項４に記載の方法。
6. 前記ＮＤＡＬｅｎｇｔｈの値は
   

   

   に等しく、ここで、ＳＮＲＬｉｎは線型に表された信号対雑音比ＳＮＲに対応することを特徴とする、請求項５に記載の方法。
7. 周波数推定値は前記通信を構成するＤＶＢ−Ｓ２フレーム毎に判断されることを特徴とする、請求項２乃至６のいずれか一項に記載の方法。
8. 少なくとも第１のステーション（１）と前記第１のステーションに対して移動する空中システム（３）とを含む航空通信システムにおけるドップラーを推定するための装置において、データが一連のフレームの形式で送信され、フレームが既知のデータの少なくとも１つの第１のヘッダフィールド（３００）とそれに続く未知のデータの１つまたは複数のフィールド（３０１_１）とを含み、前記装置は少なくとも、
   ・前記ヘッダの既知のデータに対する監視相関を生成するとともに、
   

   

   の精度を有する推定周波数
   

   

   を中心とした第１の周波数範囲の推定を行うよう構成された第１のモジュール（５１０）が実装された無線通信受信器を備え、
   ・前記無線通信受信器はさらに、前記推定周波数
   

   

   にわたって相関ピークを求めることにより、および前記相関ピークに対応する周波数を保持することにより、少なくともデータフィールド（３０１_１）のすべての未知のデータに対するブラインド相関を生成するよう適合された第２のモジュール（５２０）が実装され、前記保持される周波数は、前記第１のステーションと前記空中システムとの間でデータを送信するために使用される、装置。
9. 前記フレームはＱＰＳＫまたはＢＰＳＫ変調ＤＶＢ−Ｓ２データであることを特徴とする請求項８に記載の装置。
10. 前記第１の推定モジュールと前記第２のモジュールはＦＰＧＡ技術を使用して生成されることを特徴とする請求項８または９に記載の装置。
11. 少なくとも第１のステーションと前記第１のステーションに対して移動する空中システムとを含む航空通信システムにおけるドップラーを推定するための装置において、データが一連のフレームの形式で送信され、フレームが既知のデータの少なくとも１つの第１のヘッダフィールドとそれに続く未知のデータの１つまたは複数のフィールドとを含み、前記装置は少なくとも、
    前記ヘッダの既知のデータに対する監視相関装置が実装されるとともに、
    

    

    の精度を有する推定周波数
    

    

    を中心とした第１の周波数範囲の推定をする無線通信受信器を備え、
    前記無線通信受信器はさらに、前記推定周波数
    

    

    にわたって相関ピークを求めることにより、および前記相関ピークに対応する周波数を保持することにより、少なくともデータフィールドのすべての未知のデータに対するブラインド相関を行う装置が実装され、前記保持される周波数は、前記第１のステーションと前記空中システムとの間でデータを送信するために使用される、装置。

Images