JP4493669B2 - プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 - Google Patents

Description

本発明は、受信機に関わり、特に受信信号に既知の信号パターンが含まれているか否かを識別する方法に関するものである。

デジタル無線通信システムにおいて、送信機は、伝送する情報（２値信号のビット列）を一定のビット数で分割し、分割されたビット列を含む所定の構成がなされたフレームとし、そのフレーム毎に、例えば16QAM（Quadrature AmplitudeModulation）方式のようなデジタル変調を行い、デジタル変調された信号として送信する。送信をしていない状態から、送信機が初めて送信を開始するフレームを第１フレームという。以下第１フレームに続いて、第２フレーム、第３フレーム、‥‥‥と送信されていく。

フレームには、受信機が受信し、受信したフレームを用いて、送信される信号に対する同期を行うために、プリアンブルと呼ばれる規定の固定ビットパターンが、定められた位置に、定められた長さで付加されている。受信機は、送信機が送信した第１フレームの信号を受信する際に、第１フレーム内のプリアンブルを用いて、送信された信号に対して高速同期を行う。ここで、高速同期とは、全く同期が確立されていない状態から同期を行い、同期した情報を現フレームの復調（検波）動作に反映させる処理のことをいう。これに対し、通常同期とは、１つ前のフレームですでに同期が確立されており、現フレームと１つ前のフレームとの間での変動分に対する同期を低速で行い、同期した情報を次のフレームの復調動作に反映する処理のことをいう。第１フレームは、同期バーストと呼ばれるフレームであり、このフレームにプリアンブルが含まれている。受信機は、例えば、この第１フレームの同期バーストを用いて高速同期を行い、第２フレーム以降では通常同期を行う必要がある。そこで、プリアンブルが含まれている同期バーストを受信しているか否かの判断を行うために、プリアンブルのパターンを識別する方法が考えられている。このように、同期バーストのフレームは、普通、同期バーストは送信開始後、１フレーム〜数フレームが連続して送信される。例えば、後述する図１６は、第１フレームと第２フレームとが同期バーストのフレームである。

図１４はSCPC（Single Channel Per Carrier）の標準規格であるARIB STD-T61による同期バーストフレームのフレーム構造を示す図である。更にまた、図１５は通信チャネルフレームのフレーム構造を示す図である。図１４または図１５において、LP+Rはリニアライザプリアンブル・ランプアップ部、Pbはプリアンブル部、RIは通信情報チャネル部、SWは同期ワードパターン部、PIはパラメータ情報チャネル部、Gはガードタイム部、Tchは通信チャネル部、UDは未定義部である。また、それぞれの記号の下に記した数値はそれぞれの領域のビット（bit）数を表す。プリアンブルは、例えば、全て“0”、または、“1，0，0，1”を繰り返すような固定パターンからなる。

また、移動無線通信システムで送信を行う場合の送信パターンの一例を図１６に示す、図１６はARIB STD-T61に基づいた送信パターンである。SB₀とSB₁は同期バーストフレーム、TCH₀，TCH₁，TCH₂，‥‥‥，TCH_Mは主データを載せた通信チャネルフレームである。ここでMは自然数である。ARIB STD-T61では、送信を開始するときには、通信チャネルフレームの送信に先立ち、同期バーストフレームをm回（mは２以上の正整数）送信することになっており（通常は２回以上）、図１６では、まず同期バーストフレームを２個送信し、次から送信チャネルフレームの送信を開始し、送信終了となるまで送信する。即ち、通話は２個の同期バーストとM個の通信チャネルで構成される。ただし、M個の通信チャネルは通話の長さにより変化する。

プリアンブルパターン識別方法の説明に入る前に、図１１を用いて受信機について簡単に説明する。図１１は、従来の受信機の簡単な構成を示すブロック図である。受信信号入力端子1101には、アンテナから受信した信号が入力され、高周波部回路1102に入力される。高周波部回路1102は、受信された信号である無線帯域の高周波の信号を、A/D変換器1103でサンプリング可能な低い周波数へ周波数変換してA/D変換器1103に与える。A/D変換器1103は、周波数変換された受信信号をサンプリング及び量子化してデジタル信号とし、乗算器1104-1と1104-2とに入力する。

正弦波発生回路1106は、角周波数ωの正弦波cosωtを出力し、乗算器1104-1と移相器1105とに入力する。移相器1105は、正弦波cosωtの位相をπ/2（rad）だけ進ませて（-sinωtとし）乗算器1104-2へ入力する。乗算器1104-1は、A/D変換器1103の出力と正弦波発生回路1106の出力との積を演算し、ローパスフィルタ（LPF）1107-1へ入力し、乗算器1107-2はA/D変換器1103の出力と移相器1105の出力との積を演算し、ローパスフィルタ（LPF）1107-2へ入力する。乗算器1104-1および1104-2の出力信号には不要な高周波成分を含んでいる。ローパスフィルタ1107-1は、乗算器1104-1の出力信号から高周波成分を除去し、ルートロールオフフィルタ1108-1へ入力する。同様に、ローパスフィルタ1107-2は、乗算器1104-2の出力信号から高周波成分を除去し、ルートロールオフフィルタ1108-2へ入力する。ルートロールオフフィルタは、ローパスフィルタ1107-1と1107-2のそれぞれの出力信号の帯域制限を行い、その出力をベースバンド信号出力端子1109-1および1109-2を介してそれぞれ出力する。

ベースバンド信号出力端子1109-1を介して出力される信号をベースバンド信号の同相成分と呼び、1109-2を介して出力される信号をベースバンド信号の直交成分と呼ぶ。ベースバンド信号は、実数部を同相成分、虚数部を直交成分とする複素数の信号である。ここでは、ベースバンド信号の同相成分と直交成分をそれぞれ別の端子で表しているが、以下の説明では簡単のため１つの端子で表す。

従来のプリアンブルパターン識別方法を、図２によって説明する。図２は、従来のプリアンブルパターン識別のための回路の構成を示すブロック図である。図２において、受信したベースバンド信号が入力端子201に入力される。入力端子201に入力されるベースバンド信号は、１シンボル当たりN_ov回オーバーサンプリングされた複素数（実数部を同相成分(I)、虚数部を直交成分(Q)とする）のデジタル信号である。このオーバーサンプル数N_ovは２以上の整数である。入力端子201より入力したベースバンド信号は、相関回路202に入力し、プリアンブルパターンとの相関r(n)を演算し、演算結果を振幅自乗値演算回路203に与える。振幅自乗値演算回路203は、相関r(n)の振幅の自乗値 |r(n)|²を演算し、短時間平均演算回路204に与える。短時間平均演算回路204は相関の自乗値|r(n)|²の短時間平均r_ave(n)を演算し判定器205に与える。判定器205は、短時間平均r_ave(n)を所定のしきい値と比較し、短時間平均r_ave(n)がしきい値を上回っていれば入力信号がプリアンブルであると判断し、プリアンブル識別出力端子206を介して“1”を出力する。ここで、nはサンプル番号で、正整数ある。

入力端子201に、プリアンブルが入力される場合のベースバンド信号は、QPSK（Quaternary Phase Shift Keying）やQAMといった変調方式の場合には、図４に示すように、（I，Q）=（A，-A）と（I，Q）=（-A，A）をシンボル間隔で交互に繰り返す波形であり、変調方式がπ/4シフトQPSKの場合、図５に示すような８シンボル周期の波形である（ここで、Aは定数）。図４は、変調方式がQPSKまたはQAMである場合のプリアンブル信号のコンスタレーションおよび時間波形を示す図であり、図５は、変調方式がπ/4シフトQPSKである場合のプリアンブル信号のコンスタレーションおよび時間波形を示す図である。図４と図５において、それぞれ、(a)は、横軸をベースバンド信号の同相成分(I)、縦軸を直交成分(Q)として、I-Q座標平面にプロットしたコンスタレーションである。また(b)は、横軸を時間、縦軸をベースバンド信号の同相成分(I)としてプロットした時間波形である。更に(c)は、横軸を時間、縦軸をベースバンド信号の直交成分(Q)としてプロットした時間波形である。また、図中の黒点はシンボルタイミングでのベースバンド信号の瞬時値である。

図３は、図２で説明した相関回路202の構成を、より詳細に図示したブロック図である。図３の相関回路は、N-1個の１シンボル遅延回路302、N個の複素乗算器304とN-1個の複素加算器303で構成される。ここでNは相関演算を行うシンボル数で、２以上の正整数であり、N≧8が望ましい。pは、プリアンブルパターン（複素数）で、π/4シフトQPSKの場合には、次の式(1)で表され、16QAMの場合には、式(2)で表される。また、p^*はpの共役複素数であり、N個の乗算器304の該当するものに与えられる係数である。

図６は、図２に示した従来のプリアンブルパターン識別のための回路の、短時間平均演算回路204の出力である短時間平均r_ave(n)の時間波形を示す図である。(a)は受信信号、(b)はフェージングなしの場合の短時間平均r_ave(n)、(c)はフェージング時での短時間平均r_ave(n)を表す。また、一点鎖線はしきい値である。図６ (a)に示すように、入力端子201に区間601でプリアンブル以外の信号が入力され、区間602でプリアンブルが入力され、区間603で再びプリアンブル以外の信号が入力される場合、短時間平均演算回路204の出力r_ave(n)は、フェージングなしの場合に図６ (b)のようになり、フェージング時には図６ (c)のようになる。

短時間平均r_ave(n)は、相関回路202の段数Nと短時間平均演算回路204の時定数により遅延する。フェージングがない場合には、図６ (b)に示すように、高い相関を示すので、識別漏れの確率と誤識別の確率を共に“0”にするしきい値の設定を容易に行うことができる。しかし、フェージング時には、図６ (c)に示すように、受信信号の振幅の落ち込みにより、プリアンブル入力時の相関が落ち込むため、識別漏れの確率を“0”に近づけようとしきい値を低く設定すると誤識別の確率が大きくなり、反対に誤識別の確率を“0”に近づけようとしきい値を高く設定すると識別漏れの確率が大きくなり、両方の確率を“0”にするしきい値を設定することが困難になる。

また、プリアンブルの信号が持つ周波数成分は、シンボルレート（シンボル間隔の逆数）をf_bとすると、変調方式がQPSKまたはQAMの場合には、-f_b/2と+f_b/2の２つの周波数成分であり、変調方式がπ/4シフトQPSKの場合には、-3f_b/8と+f_b/8の２つの周波数成分である。送信機と受信機との間には、周波数偏差Δfがあり、この周波数偏差Δfが大きくて、プリアンブルの片方の周波数成分が受信機のフィルタの帯域外へ外れる場合には、相関回路202に入力されるプリアンブルの片方の周波数成分が削られるため、周波数偏差Δfによる位相回転の補正を行ったとしても、相関による識別が難しくなる。

前述の従来技術には、相関演算結果が受信信号の振幅変動に影響され、識別漏れの確率と誤識別の確率とが相反する傾向にあり、両方の確率を0に近づけるようなしきい値の設定が困難であり、しきい値による識別が難しい欠点があった。また、送信機と受信機との間で、周波数偏差が大きい場合、プリアンブルの片方の周波数成分が受信機のフィルタの帯域外へ外れ、プリアンブルの２つの周波数成分の片方の成分が削られた場合の識別が困難であるという欠点があった。本発明の目的は、上記のような欠点を除去し、受信信号の振幅変動に影響されず、送信機と受信機との間で、周波数偏差が大きい場合においても、プリアンブルの識別を容易に行うことができるプリアンブルパターン識別方法を提供することにある。

上記の目的を達成するため、本発明のプリアンブルパターン識別方法は、
(1)受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出したそれぞれの周波数成分の電力p₁とp₂を求め、求めたそれぞれの周波数成分の電力和p₁+p₂と全周波数成分の電力pとの電力比(p₁+p₂)/pが所定のしきい値を超えていれば、受信したベースバンド信号がプリアンブルであると判断する。尚、以下の説明で、周波数のずれていない送信ベースバンドでのプリアンブルの周波数成分を、プリアンブル送信時の周波数成分と呼ぶ。
(2) 電力比(p₁+p₂)/pが所定のしきい値を超え、かつ２つの成分の周波数f₁とf₂の差の絶対値|f₂-f₁|が所定の値の範囲内であれば、受信したベースバンド信号がプリアンブルであると判断する。
(3) 抽出したそれぞれの周波数成分の電力p₁とp₂の電力比p₁/p₂を算出し、電力比p₁/p₂が規定値r₁以上かつ規定値r₂以下ならば２つの周波数f₁とf₂の成分が残っていると判断し、２つの周波数f₁とf₂の成分が残っていると判断した場合で、電力比(p₁+p₂)/pが所定のしきい値を超え、かつ周波数差の絶対値|f₂-f₁|が所定の値の範囲内であれば、受信したベースバンド信号がプリアンブルであると判断する。ただし、r₁，r₂は正の定数で、r₁<1，r₂>1とする。

(4) また、本発明の周波数偏差検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出して、低い周波数成分の周波数をf₁、高い周波数成分の周波数をf₂とし、プリアンブル送信時の２成分の周波数のうち低い周波数成分の周波数をf₁₀、高い周波数成分の周波数をf₂₀とし、各成分の周波数偏差f₁-f₁₀, f₂-f₂₀の平均[(f₁-f₁₀)+(f₂-f₂₀)]/2を、送信機と受信機との周波数偏差とする。
(5) 更に、本発明の周波数偏差検出方法は、抽出した２成分のうちの一方の成分の周波数をf₁'、プリアンブル送信時の２成分の周波数のうち低い周波数成分の周波数をf₁₀、高い周波数成分の周波数をf₂₀とし、f₁'-f₁₀またはf₁'-f₂₀を、送信機と受信機との周波数偏差とする。
(6) また、本発明の周波数偏差検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出したそれぞれの周波数成分の電力p₁とp₂の電力比p₁/p₂を算出し、電力比p₁/p₂が規定値r₁以上かつ規定値r₂以下ならば２つの周波数f₁とf₂の成分が残っていると判断して、２つの周波数f₁とf₂の成分が残っていると判断した場合には、各成分の周波数偏差f₁-f₁₀とf₂-f₂₀の平均[(f₁-f₁₀)+(f₂-f₂₀)]/2を、送信機と受信機との周波数偏差とし、電力比p₁/p₂が規定値r₁未満の場合または規定値r₂を超えている場合で、電力の大きい方の成分の周波数をf₁'とし、f₁'≧0であれば、２成分のうち周波数の低い成分のみが残っているものとし、f₁'-f₁₀を送信機と受信機との周波数偏差とし、電力比p₁/p₂が規定値r₁未満の場合または規定値r₂を超えている場合で、電力の大きい方の成分の周波数をf₁'とし、f₁'＜0であれば、２成分のうち周波数の高い成分のみが残っているものとし、f₁'-f₂₀を送信機と受信機との周波数偏差とする。

(7) 更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出したそれぞれの成分の位相φ₁とφ₂の差φ₁-φ₂よりシンボルタイミングを求める。
(8) また更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出したそれぞれの周波数成分の電力p₁とp₂の電力比p₁/p₂を算出し、電力比p₁/p₂が規定値r₁以上かつ規定値r₂以下ならば２つの周波数f₁とf₂の成分が残っていると判断し、抽出したそれぞれの成分の位相φ₁とφ₂の位相差φ₁-φ₂よりシンボルタイミングを求める。
(9) また本発明のシンボルタイミング検出方法は、変調方式がQPSKまたはQAMであるベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから周波数成分を抽出し、抽出した２つの周波数f₁とf₂の成分のうち、低い周波数の成分の位相をφ₁とし、高い周波数の成分の位相をφ₂とし、シンボルレートをf_bとし、t₀=(φ₁-φ₂)/(2πf_b)をシンボルタイミングt₀とする。
(10) 更に本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出したそれぞれの周波数成分の電力p₁とp₂の電力比p₁/p₂を算出し、電力比p₁/p₂が規定値r₁以上かつ規定値r₂以下ならば２つの周波数f₁とf₂の成分が残っていると判断し、２つの周波数f₁とf₂の成分が残っていると判断した場合、抽出した２つの周波数f₁とf₂の成分のうち、低い周波数の成分の位相をφ₁、高い周波数の成分の位相をφ₂、シンボルレートをf_bとし、t₀=(φ₁-φ₂)/(2πf_b)をシンボルタイミングとする。
(11) また、本発明のシンボルタイミング検出方法は、変調方式がπ/4シフトQPSKであるベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出した２つの周波数f₁とf₂の成分のうち、低い周波数の成分の位相をφ₁、高い周波数の成分の位相をφ₂、シンボルレートをf_bとし、t₀=(φ₁-φ₂-π/2)/(πf_b)をシンボルタイミングとする。
(12) 更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f₁とf₂の成分を抽出し、抽出したそれぞれの周波数成分の電力p₁とp₂の電力比p₁/p₂を算出し、電力比p₁/p₂が規定値r₁以上かつ規定値r₂以下ならば２つの周波数f₁とf₂の成分が残っていると判断し、２つの周波数f₁とf₂の成分が残っていると判断した場合、抽出した２つの周波数f₁とf₂の成分のうち、低い周波数の成分の位相をφ₁、高い周波数の成分の位相をφ₂、シンボルレートをf_bとし、t₀=(φ₁-φ₂-π/2)/(πf_b)をシンボルタイミングとする。

(13) また更に、本発明のプリアンブル成分抽出方法は、受信信号を１シンボル当たりN_ov回（N_ovは2以上の整数）オーバーサンプリングしたベースバンド信号から、任意の連続したN_winサンプル（N_winは２以上の整数）のベースバンド信号を抽出し、抽出したベースバンド信号に任意の窓関数を乗じ、任意の窓関数を乗じたベースバンド信号x(n) （n=0，1，‥‥‥，(N_win-1)）のDFT（離散フーリエ変換: Discrete Fourier Transform）X(k)（k=0，1，‥‥‥，(N_win-1)）を演算し（X(k)は複素数）、DFT出力X(k)（k=0，1，‥‥‥，(N_win-1)）の全電力p=|X(0)|²+|X(1)|²+…+|X(N_win-1)|²を算出し、DFT出力X(k)のうち、振幅の自乗値|X(k)|²が最大となるX(k)を検索してそのkをk₁₁とし、X(k₁₁)に隣接する成分から、|X(k₁₁)|²の次に|X(k)|²が大きくなるX(k)を検索してそのkをk₁₂とし、X(k₁₁)とX(k₁₂)の電力和p₁=|X(k₁₁)|²+|X(k₁₂)|²を算出し、X(k₁₁)の周辺を補間し、補間した部分のピークにおける周波数をf₁、位相をφ₁とし、変調方式により定められる自然数の定数をk_dとし、X(k)（k=(k₁₁-k_d-1) mod N_win，(k₁₁-k_d) mod N_win，(k₁₁-k_d+1) mod N_win，(k₁₁+k_d-1) mod N_win，(k₁₁+k_d) mod N_win，(k₁₁+k_d+1) mod N_win）から|X(k)|²が最大となるX(k)を検索してそのkをk₂₁とし、X(k₂₁)に隣接する成分のうち、|X(k₂₁)|²の次に|X(k)|²が大きくなるX(k)を検索してそのkをk₂₂とし、X(k₂₁)とX(k₂₂)の電力和p₂=|X(k₂₁)|²+|X(k₂₂)|²を算出し、X(k₂₁)の周辺を補間し、補間した部分のピークにおける周波数をf₂、位相をφ₂とし、算出したf₁, f₂を比較し、f₁>f₂であれば、f₁，φ₁，p₁を改めてf₂，φ₂，p₂とし、f₂，φ₂，p₂を改めてf₁，φ₁，p₁とする。

(14) 更にまた、本発明のプリアンブル成分抽出方法は、FFT（高速フーリエ変換：Fast Fourier Transform）アルゴリズムを用いて、任意の窓関数を乗じたベースバンド信号x(n)（n=0，1，‥‥‥，(N_win-1)）のDFT（離散フーリエ変換：Discrete Fourier Transform）X(k)（k=0，1，‥‥‥，(N_win-1)）を演算する。以上において、また、周波数スペクトルから抽出する２つの周波数f₁とf₂の成分は、振幅が高い順に２つ抽出する。

本発明によれば、プリアンブルパターン識別に用いる電力比が、受信信号の振幅変動に影響せず、振幅変動の影響を受けないプリアンブルパターン識別が可能である。また、送信機と受信機との周波数偏差によりプリアンブルの２つの周波数成分のうち、片方の周波数成分が受信機のフィルタの帯域外へ外れて削られる場合においても、プリアンブルパターン識別が可能である。更に、本発明によれば、送信機と受信機との周波数偏差の検出が可能で、プリアンブルの２つの周波数成分が残っている場合には、シンボルタイミング検出を行うこともできる。

本発明の実施形態の説明に入る前に、プリアンブル信号の周波数スペクトルについて説明する。最初に変調方式が、QPSKまたはQAMの場合を、図１２と図７によって説明する。図１２は、変調方式がQPSKまたはQAMである場合のプリアンブル信号のマッピングとゼロ補間（0補間）を説明するための図であり、図７は、変調方式がQPSKまたはQAMである場合のプリアンブル信号のDFT（離散フーリエ変換: Discrete Fourier Transform）出力を示す図である。図１２に示すように、シンボル毎に√2・A exp(-jπ/4)、√2・A exp(j3π/4)を交互にマッピングし、１シンボル当たり２回オーバーサンプリングしてゼロ補間したサンプル値x(0)〜x(15)のDFTを計算すると式(3)のようになり、周波数スペクトルは図７ (a)に示すようになる。ここで、Aは図４と同じ値で、ω_bはシンボルレートをf_bとしω_b=2πf_bである。

送信機が送信する信号は、式(4)に示す伝達関数R_R(ω)を持つルートロールオフフィルタで帯域制限され、受信機にも同じR_R(ω)の特性のルートロールオフフィルタを持ち、このフィルタの出力信号を復調する。プリアンブルパターン識別では処理遅延を削減するため、受信機のルートロールオフフィルタ前段の信号を用いる。

ここで、αをロールオフ率と呼び、0≦α≦1である。

図７ (b)は、図７ (a)の信号を、式(4)に示す伝達関数R_R(ω)を持つルートロールオフフィルタで帯域制限した周波数スペクトルであり、R_R(ω_b/2) = R_R(-ω_b/2) = 1/√2であるから、図７ (b)の各成分は、式(5)のようになり、ω=±ω_b/2の2つの周波数成分を持つ。

従って、時間波形は、式(6)のようになる。

受信ベースバンド信号がシンボルタイミングより時間t₀だけ遅れていて、送信機と受信機との周波数偏差をΔωとすると、式(7)のようになる。また、各成分の角周波数はω₁=-ω_b/2+Δω、ω₂=ω_b/2+Δω、ω₁とω₂各成分の位相φ₁，φ₂は式(8)のようになる。ここで、θ₀は定数である。

受信ベースバンド信号がプリアンブルの信号である場合の周波数スペクトルは、上記ω₁とω₂の成分のみであるから、周波数スペクトルの全電力をp、ω₁成分の電力をp₁、ω₂成分の電力をp₂とすると、２成分の電力和(p₁+p₂)が全電力pに近い値となることから、電力比(p₁+p₂)/pによりプリアンブルパターンの識別ができる。また、式(8)より(φ₁-φ₂)を求めてθ₀を消去すると、各成分の位相φ₁とφ₂からシンボルタイミングt₀を求めることができる。

次に、変調方式がπ/4シフトQPSKの場合のプリアンブル周波数スペクトルについて、図１３と図８によって説明する。図１３は、変調方式がπ/4シフトQPSKである場合のプリアンブルの信号のマッピングとゼロ補間を説明するための図であり、図８は、変調方式がπ/4シフトQPSKである場合のプリアンブルの信号のDFT出力を示す図である。図１３に示すように、シンボル毎にA exp(j0)、A exp(-jπ/4)、A exp(jπ/2)、A exp(jπ/4)、A exp(jπ)、A exp(j3π/4)、A exp(-jπ/2)、A exp(-j3π/4)でマッピングし、１シンボル当たり２回オーバーサンプリングしてゼロ補間したサンプル値x(0)〜x(15)のDFTを計算すると式(10)のようになり、周波数スペクトルは図８ (a)のようになる。ここで、Aは図５と同じ値である。

この場合も、変調方式がQPSKまたはQAMである場合と同様に、式(4)の伝達関数を持つルートロールオフフィルタで帯域制限された周波数スペクトルを考えると、周波数スペクトルは図８ (b)〜(d)のようになる。この場合、ルートロールオフフィルタのロールオフ率αにより周波数スペクトルが変化する。ロールオフ率αが0≦α≦1/4の場合は、図８ (b)に示すように、角周波数が-3ω_b/8とω_b/8の2成分X(13)とX(1)のみが残り、これらの成分の振幅は変化しない。

ロールオフ率αが1/4≦α≦3/4の場合は、図８ (c)に示すように、角周波数がω_b/8であるX(1)はそのまま残り、-3ω_b/8，5ω_b/8の成分はそれぞれ、式(11)のようになり、ロールオフ率α=0.5の場合はX'(13) = 0.924 X(13)、X'(5) = 0.383 X(5)となる。

ロールオフ率αが3/4≦α≦1の場合は同図(d)に示すように、-7ω_b/8、-3ω_b/8、ω_b/8、5ω_b/8の4成分が残り、各成分はそれぞれ式(12)のようになる。

ロールオフ率αは、通常、0.5 以下であることが殆どであるため、図８ (b)及び図８ (c)の場合のみ考えると、ロールオフ率α=0.5の場合プリアンブルの信号入力時は、-3ω_b/8とω_b/8の２成分の電力が全電力の92.7%以上［α=0.5の場合、(1+0.924²) / (1+0.924²+0.383²) = 0.927］を占めるため、-3ω_b/8の成分の電力をp₁とし、ω_b/8の成分の電力をp₂とし、周波数スペクトルの全電力pに対する比(p₁+p₂)/pによりプリアンブルパターン識別ができる。

図８ (b)及び図８ (c)よりプリアンブルの信号の時間波形は、ロールオフ率αが0≦α≦1/4の場合には、式(13)のようになり、ロールオフ率αが1/4≦α≦3/4の場合には、式(14)のようになる。

更に、QPSKまたはQAMの場合と同様に、受信ベースバンド信号がシンボルタイミングより時間t₀だけ遅れていて、送信機と受信機との間にΔωの周波数偏差がある場合の時間波形を求めると、ロールオフ率αが0≦α≦1/4の場合には式(15)のようになり、ロールオフ率αが1/4≦α≦3/4の場合には、式(16)のようになる。

式(15)と(16)に示す-3ω_b/8+Δω、ω_b/8+Δωの各成分の位相φ₁、φ₂は、式(17)のようになる。ただし、θ₀は定数である。

従って、シンボルタイミングt₀は数式(18)より求められる。

プリアンブル以外の信号、例えば、通信チャネルTchが受信された場合には、複数の周波数スペクトルが抽出されるが、その場合には、スペクトルの振幅が大きい順に２つ選択し、それらをp₁とp₂とする。そして、２成分の電力和(p₁+p₂)を求めると、全電力pに比べて小さな値となることから、電力比(p₁+p₂)/pが小さくなることによりプリアンブルパターンではないと判定することができる。

次に、本発明の実施形態を、図１を用いて説明する。図１は、本発明のプリアンブルパターン識別、周波数偏差検出、及びシンボルタイミング検出を行なう回路構成の一実施例を示すブロック図である。図１において、受信されたベースバンド信号が入力端子101に入力される。このベースバンド信号は、従来技術の場合と同様に、１シンボル当たりN_ov回オーバーサンプリングされた複素数のデジタル信号であり、N_ovは２以上の正整数である。

入力端子101を介して入力したベースバンド信号は、サンプル毎にシフトレジスタ114に入力される。シフトレジスタ114の段数は、N_winで、シフトレジスタ114内部のN_win個のレジスタをs(0)，s(1)，‥‥‥，s(N_win-2)，s(N_win-1)とすると、入力端子101を介してベースバンド信号が入力されるたびに、s(N_win-1)，s(N_win-2)，‥‥‥，s(1)の内容が、それぞれ、s(N_win-2)，s(N_win-3)，‥‥‥，s(0)にシフトされ、入力端子101を介して入力されるベースバンド信号がs(N_win-1)に入力される。シフトレジスタ114の内容s(0)，s(1)，‥‥‥，s(N_win-2)，s(N_win-1)は、サンプル毎に、スイッチ113に与えられ、一定の時間間隔N_step（サンプル）毎に１回、s(0)，s(1)，‥‥‥，s(N_win-2)，s(N_win-1)が窓掛け回路102に入力される。（例えば、オーバーサンプル数N_ov=2で、N_step= 8、N_win = 16である）。このように、スイッチ113は、N_stepサンプル周期に１回閉じていて、N_winサンプルの信号をN_step間隔で取り出す。スイッチ113から入力されたN_winサンプルの信号は、窓掛け回路102によって、窓（例えば、ハミング窓（Hamming Window））を掛けられ、FFT演算回路103に入力される。

FFT演算回路103は、窓掛け回路102から入力されるN_winサンプルの信号のDFT（離散フーリエ変換：Discrete Fourier Transform）を、FFT（高速フーリエ変換：Fast Fourier Transform）アルゴリズムにより演算する。図９は、N_ov= 2、N_win = 16の場合のFFT演算回路103の出力の一実施例で、送信機と受信機との間に周波数偏差がある場合のプリアンブル信号のDFT出力を示す図である。図９は、角周波数ω₁とω₂の２つの周波数成分を持つプリアンブル信号が入力された場合のDFT出力であり、-3ω_b/8の周波数成分であるX(13)と、5ω_b/8の周波数成分であるX(5)がピークとなっている。しかし、DFT出力は、周波数が離散的であるため、実際のプリアンブルの成分ω₁（ピーク903）とω₂（ピーク904）は、それぞれ、X(12)とX(13)のポイント間とX(4)とX(5)のポイント間に存在する。従って、実際のピーク903を求めるには、例えば、ポイントX(12)，X(13)，X(14)間でデータ補間し、包絡線901を求めて行う。同様に、実際のピーク904を求めるには、例えば、ポイントX(4)，X(5)，X(6)間でデータ補間し、包絡線902を求めて行う。

FFT演算回路103での演算結果（DFT出力）は、全電力算出回路104、及び、プリアンブル成分抽出回路105と106とに与えられる。全電力算出回路104、及び、プリアンブル成分抽出回路105と106の処理は、図１０に示すフローチャートの手順に従って、全電力算出およびプリアンブル成分の抽出を行う。図１０は、本発明のDFT出力からプリアンブル成分を抽出する手順の一実施例を示すフローチャートである。

図１０において、ステップ1001では、FFT演算回路103の出力X(k)（k = 0，1，‥‥‥，(N_win-1)）から、|X(k)|²が最大となるものを検索し、検索された|X(k)|²が最大となるkをk₁₁（k = k₁₁）としてステップ1002に進む。ステップ1002では、m₁とm₂を算出し（m₁= (k₁₁ - 1) mod N_win，m₂ = (k₁₁ + 1) mod N_win）、ステップ1003に進む。ステップ1003では、算出されたm₁とm₂から、|X(m₁)|²と|X(m₂)|²を比較し、|X(m₂)|²＞|X(m₁)|²であれば、ステップ1004に進み k₁₂ = m₂としステップ1006に進む。また、|X(m₂)|²≦|X(m₁)|²であれば、ステップ1005に進み k₁₂ = m₁とし、ステップ1006に進む。ステップ1006では、X(k₁₁)とX(k₁₂)の電力和p₁（p₁= |X(k₁₁)|² + |X(k₁₂)|²）を算出して、ステップ1007に進む。受信したプリアンブルの実際のピークは、X(k₁₁)とX(k₁₂)の間に存在するため、ステップ1007では、X(m₁)、X(k₁₁)、X(m₂)の３点でポイント間の補間を行い、補間を行った部分での実際のピークを検索し、その実際のピークでの周波数f₁と周波数スペクトルの位相φ₁を算出する。

次にもう１つの成分を抽出する。k_dはプリアンブルの２成分の周波数差|ω₂-ω₁|に相当するkの値である。変調方式がQPSKまたはQAMの場合には、|ω₂-ω₁|= ω_bであるから、k_d=N_win/N_ov、π/4シフトQPSKの場合には、|ω₂-ω₁|=ω_b/2であるから、k_d=N_win/(2N_ov)である。ステップ1008では、m₁=(k₁₁-k_d) mod N_winとm₂=(k₁₁+k_d) mod N_winとを算出し、ステップ1009に進む。ステップ1009では、|X(k)|²(k=(m₁-1) mod N_win，m₁，(m₁+1) mod N_win)の最大値を検索し、第１の最大値をx₁とし、その時のkをm₁（k= m₁）としステップ1010に進む。ステップ1010では、また、|X(k)|² (k=(m₂-1) mod N_win，m₂，(m₂+1) mod N_win)の最大値を検索し、第２の最大値をx₂とし、その時のkをm₂（k = m₂）とする。

ステップ1011では、第１の最大値x₁と第２の最大値x₂とを比較し、第１の最大値x₁が大きければ、ステップ1012で k₂₁=m₁を算出しステップ1014に進む。また、そうでなければ、ステップ1013で k₂₁=m₂を算出しステップ1014に進む。これにより、k₂₁がもう１つの成分の番号となる。ステップ1014では、m₁=(k₂₁-1) mod N_winとm₂=(k₂₁+1) mod N_winを算出し、ステップ1015に進む。ステップ1015では、|X(m₁)|²と|X(m₂)|²を比較する。|X(m₁)|²より|X(m₂)|²が大きい場合には、ステップ1016でk₂₂=m₂としてステップ1018に進む。また、そうでない場合には、ステップ1017でk₂₂=m₁としてステップ1018に進む。ステップ1018では、X(k₂₁)とX(k₂₂)の電力和p₂を算出し、ステップ1019に進む。

受信したプリアンブルの実際のピークは、X(k₂₁)とX(k₂₂)の間に存在するため、ステップ1019において、X(m₁)、X(k₂₁)、X(m₂)の３点でポイント間の補間を行い、補間を行った部分でのピークを検索し、そのピークでの周波数f₂と周波数スペクトルの位相φ₂を算出し、ステップ1020に進む。最後に、ステップ1020では、FFT演算回路103の出力の全電力p = |X(0)|²+|X(1)|²+…+|X(N_win-1)|²を算出する。

ステップ1001〜1020の処理は、図９に示したFFT回路103の出力の場合には、ステップ1001でX(5)が選択され、k₁₁=5、m₁=4，m₂=6で|X(4)|²と|X(6)|²を比較し（ステップ1003）、|X(4)|²が大きいためk₁₂=4となり（ステップ1005）、p₁=|X(5)|²+|X(4)|²が演算される（ステップ1006）。受信したプリアンブルの実際のピークは，ω₂の位置に存在するため、X(4)、X(5)、X(6)の３点の間で補間を行い（破線902）、補間した部分でのピークがω₂で検出されるので、その周波数をf₁=ω₂/(2π)とする。また、ω₂での補間値の位相をφ₁とする（ステップ1007）。

図９は、変調方式がQPSKまたはQAMの場合の一例であり、k_d=N_win/N_ov=8であるから、m₁=13, m₂=13と同じ値であるため（ステップ1008）、ステップ1009とステップ1010において、|X(k)|² (k=12，13，14)の最大値を検索し、x₁=x₂=|X(13)|²，m₁=m₂=13となり（ステップ1009，1010）、k₂₁=m₂=13がもう一つのピーク成分の番号となる（ステップ1013）。m₁=12，m₂=14として|X(12)|²と|X(14)|²を比較し（ステップ1015）、|X(12)|²が大きいためk₂₂=m₁=12となり（ステップ1017）、p₂=|X(13)|²+|X(12)|²が演算される（ステップ1018）。受信したプリアンブルの実際のピークは、ω₁の位置に存在するため、X(12)，X(13)，X(14)の３点の間で補間を行い（破線901）、補間した部分でのピークがω₁で検出されるので、その周波数をf₂=ω₁/(2π)とする。また、ω₁での補間値の位相をφ₂とする（ステップ1019）。

図１のプリアンブルパターン識別回路107、周波数偏差検出回路108、及び、シンボルタイミング検出回路109は、図１０のフローチャートの手順で算出した、DFT出力の全電力p、２つの周波数成分の周波数f₁とf₂、位相φ₁とφ₂、電力p₁とp₂とを用いて、プリアンブルパターン識別、周波数偏差検出、シンボルタイミング検出を行う。

プリアンブルパターン識別回路107の処理は、全電力pに対する２成分の電力和p₁+p₂の比γ=(p₁+p₂)/pを用い、更にプリアンブルの２成分が受信機のフィルタにより削られていいない場合は、２成分の周波数差 f₁₂=|f₂-f₁|も用いる。入力端子101にプリアンブル信号が入力されていて、２成分ともに削られない場合は、全電力pは、２成分電力和 p₁+p₂と雑音電力との和であり（ただし、変調方式がπ/4シフトQPSKでロールオフ率αが1/4≦α≦1/2の場合はp₁，p₂と雑音以外に最大7.3%の成分の電力がpに含まれる）、２成分の内の片方の成分が削られる場合は、p₁またはp₂のいずれかが雑音電力となる。従って、片方の成分が削られているか否かの判断は、p₁とp₂の比較により行う。例えば、p₁とp₂を比較し、片方が6dB以上低い値（比が4倍以上）ならば片方の成分が削られていると判断する。

変調方式がQPSKまたはQAMの場合で、２成分とも削られていないと判断した場合には、γがしきい値r_thを超え、かつf₁₂がf_b±f₀（f_bはシンボルレート、f₀は規定値）の範囲内ならば受信信号がプリアンブル信号であると識別する。片方の成分が削られていると判断した場合は、γとしきい値r_thの比較のみでプリアンブルパターン識別を行う。変調方式がπ/4シフトQPSKの場合で、２成分とも削られていないと判断した場合には、γがしきい値r_thを超え、かつf₁₂がf_b/2±f₀の範囲内ならば受信信号がプリアンブル信号であると識別する。片方の成分が削られていると判断した場合は、γとしきい値r_thの比較のみでプリアンブルパターン識別を行う。しきい値r_thと規定値f₀は、例えばr_th=0.8、f₀=f_b/8に設定する。受信信号がプリアンブルであると識別した場合には、プリアンブルパターン識別出力端子110を介して“1”を出力し、識別できない場合には“0”を出力する。

周波数偏差検出回路108の処理は、２成分の各々の周波数f₁，f₂と、プリアンブルの信号の送信ベースバンドでの周波数成分の各成分の周波数f₁₀，f₂₀（低い周波数の成分をf₁₀、高い周波数の成分をf₂₀とする）を用いて行う。２成分とも削られていない場合は両方の周波数を用いて式(19)により算出する。

また、片方の成分が削られている場合は残っている成分の周波数をf₁'とし、このf₁'を用いて式(20)により算出する。ここで、f₁'≧0の場合は周波数の高い成分が削られているものとして、f₁₀との差をΔfとし、f₁'<0の場合は周波数の低い成分が削られているものとして、f₂₀との差をΔfとする。

検出した周波数偏差Δfは、周波数偏差出力端子111を介して出力される。

シンボルタイミング検出回路109の処理は、プリアンブルの２成分とも削られていないと判断した場合にのみ行い、２成分の位相φ₁とφ₂を用いて行う。変調方式がQPSKまたはQAMの場合、式(9)よりシンボルタイミングf_bt₀は、式(21)のようになる。

また、変調方式がπ/4シフトQPSKの場合には、式(18)よりシンボルタイミングf_bt₀は、式(22)のようになる。

検出したシンボルタイミングf_bt₀は、シンボルタイミング出力端子112を介して出力される。

以上の説明による本発明の実施形態では、入力端子101にプリアンブル以外の信号が入力された場合、全電力pは、２つの周波数スペクトル成分の電力p₁とp₂以外に信号成分の電力を含んでいるため、電力比γの値は小さくなる。また、入力端子101にプリアンブルが入力された場合は、振幅が変動したとしても、p₁+p₂は振幅変動の影響を受けたプリアンブル成分の電力で、全電力pは、振幅変動の影響を受けたプリアンブル成分の電力と雑音電力の和であり、振幅の落ち込みによりp₁+p₂が小さくなる場合は、pも小さくなる。このため、雑音の影響を受けるものの、振幅変動は直接影響しない。従って、振幅変動の影響を受けないプリアンブルパターン識別が可能である。

更に、送信機と受信機との周波数偏差により、プリアンブルの２つの周波数成分のうち片方の成分が、受信機のフィルタの帯域外へ外れたとしても、p₁とp₂の一方がプリアンブルの１成分の電力で、他が雑音電力であり、全電力pはプリアンブルの１成分の電力と雑音電力の和であるため、電力比γを用いてプリアンブルパターン識別を行うことが可能である。また更に、抽出した２つの周波数成分の周波数より、送信機と受信機との周波数偏差を検出することができ、プリアンブルの２つの周波数成分が受信機のフィルタにより削られずに残っている場合は、抽出した２つの成分の位相よりシンボルタイミングを検出することも可能である。

本発明の一実施例のプリアンブルパターン識別および周波数偏差・シンボルタイミング検出のための回路構成を示すブロック図。 従来のプリアンブルパターン識別のための回路構成を示すブロック図。 従来のプリアンブルパターン識別で用いる相関回路の構成を示すブロック図。 変調方式がQPSKまたはQAMである場合のプリアンブル信号のコンスタレーションおよび時間波形を示す図。 変調方式がπ/4シフトQPSKである場合のプリアンブル信号のコンスタレーションおよび時間波形を示す図。 図２の回路での相関出力の自乗値の短時間平均の時間波形を示す図。 変調方式がQPSKまたはQAMである場合のプリアンブル信号のDFT出力を示す図。 変調方式がπ/4シフトQPSKである場合のプリアンブル信号のDFT出力を示す図。 本発明のFFT演算回路の出力の一実施例を示す図。 本発明のDFT出力からプリアンブル成分を抽出する手順の一実施例を示すフローチャート。 従来の受信機の簡単な構成を示すブロック図。 変調方式がQPSKまたはQAMである場合のプリアンブル信号のマッピングとゼロ補間を説明するための図。 変調方式がπ/4シフトQPSKである場合のプリアンブル信号のマッピングとゼロ補間を説明するための図。 従来の同期バーストフレームのフレーム構造を示す図。 従来の通信チャネルフレームのフレーム構造を示す図。 従来のARIB STD-T61に基づく送信パターンを説明する図。

符号の説明

101：入力端子、 102：窓掛け回路、 103：FFT演算回路、 104：全電力算出回路、 105，106：プリアンブル成分抽出回路、 107：プリアンブルパターン識別回路、 108：周波数偏差検出回路、 109：シンボルタイミング検出回路、 110：プリアンブルパターン識別出力端子、 111：周波数偏差出力端子、112：シンボルタイミング出力端子、 113：スイッチ、 114：シフトレジスタ、 201：入力端子、 202：相関回路、 203：振幅自乗値演算回路、 204：短時間平均演算回路、 205：判定器、 206：プリアンブル識別出力端子、 301：入力端子、 302：１シンボル遅延回路、 303：複素加算器、 304：複素乗算器、 305：相関出力端子、 601，603：プリアンブル以外の信号を受信する区間、 602：プリアンブルを受信する区間、 901：DFT出力のX(12)，X(13)，X(14)を補間した包絡線、 902：DFT出力のX(4)，X(5)，X(6)を補間した包絡線、903，904：実際のピーク、 1101：受信信号入力端子、 1102：高周波部回路、1103：A/D変換器、 1104：乗算器、 1105：移相器、 1106：正弦波発生回路、 1107：ローパスフィルタ、 1108：ルートロールオフフィルタ、 1109：ベースバンド信号出力端子。

Claims (1)

1. 受信したベースバンド信号の周波数スペクトルをＦＦＴによって求め、該周波数スペクトルから振幅の高い順に２つの周波数成分を抽出し、
   該抽出した２つの周波数成分のうちの低い方の周波数をf₁、高い方の周波数をf₂とし、
   前記周波数f ₁ と前記周波数f ₂ のそれぞれの電力p ₁ とp ₂ の電力比p ₁ /p ₂ を算出し、
   該電力比p ₁ /p ₂ が規定値r ₁ 以上かつ規定値r ₂ 以下ならば、前記周波数f ₁ からプリアンブルの信号の送信時の周波数成分の２成分のうち低い方の周波数f ₁₀ を減じた周波数偏差f ₁ -f ₁₀ を求めると共に、前記周波数f ₂ からプリアンブルの信号の送信時の周波数成分の２成分のうち高い方の周波数f ₂₀ を減じた周波数偏差f ₂ -f ₂₀ を求め、周波数偏差f ₁ -f ₁₀ とf ₂ -f ₂₀ の平均[(f ₁ -f ₁₀ )+(f ₂ -f ₂₀ )]/2を送信機と受信機との周波数偏差とし、
   電力比p ₁ /p ₂ が規定値r ₁ 未満の場合または規定値r ₂ を超えている場合は、振幅の高い方の成分の周波数をf ₁ ′とし、f ₁ ′≧0であればf ₁ ′-f ₁₀ を送信機と受信機との周波数偏差とする一方、f ₁ ′＜0であればf ₁ ′-f ₂₀ を送信機と受信機との周波数偏差とする、
   ことを特徴とする周波数偏差検出方法。

Images