JP3933399B2 - プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 - Google Patents
プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 Download PDFInfo
- Publication number
- JP3933399B2 JP3933399B2 JP2001023874A JP2001023874A JP3933399B2 JP 3933399 B2 JP3933399 B2 JP 3933399B2 JP 2001023874 A JP2001023874 A JP 2001023874A JP 2001023874 A JP2001023874 A JP 2001023874A JP 3933399 B2 JP3933399 B2 JP 3933399B2
- Authority
- JP
- Japan
- Prior art keywords
- win
- preamble
- frequency
- components
- power
- 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.)
- Expired - Fee Related
Links
Images
Description
【０００１】
【発明の属する技術分野】
本発明は、受信機に関わり、特に受信信号に既知の信号パターンが含まれているか否かを識別する方法に関するものである。
【０００２】
【従来の技術】
デジタル無線通信システムにおいて、送信機は、伝送する情報（２値信号のビット列）を一定のビット数で分割し、分割されたビット列を含む所定の構成がなされたフレームとし、そのフレーム毎に、例えば16QAM（Quadrature Amplitude Modulation）方式のようなデジタル変調を行い、デジタル変調された信号として送信する。
送信をしていない状態から、送信機が初めて送信を開始するフレームを第１フレームという。以下第１フレームに続いて、第２フレーム、第３フレーム、‥‥‥と送信されていく。
【０００３】
フレームには、受信機が受信し、受信したフレームを用いて、送信される信号に対する同期を行うために、プリアンブルと呼ばれる規定の固定ビットパターンが、定められた位置に、定められた長さで付加されている。
受信機は、送信機が送信した第１フレームの信号を受信する際に、第１フレーム内のプリアンブルを用いて、送信された信号に対して高速同期を行う。ここで、高速同期とは、全く同期が確立されていない状態から同期を行い、同期した情報を現フレームの復調（検波）動作に反映させる処理のことをいう。これに対し、通常同期とは、１つ前のフレームですでに同期が確立されており、現フレームと１つ前のフレームとの間での変動分に対する同期を低速で行い、同期した情報を次のフレームの復調動作に反映する処理のことをいう。
第１フレームは、同期バーストと呼ばれるフレームであり、このフレームにプリアンブルが含まれている。受信機は、例えば、この第１フレームの同期バーストを用いて高速同期を行い、第２フレーム以降では通常同期を行う必要がある。そこで、プリアンブルが含まれている同期バーストを受信しているか否かの判断を行うために、プリアンブルのパターンを識別する方法が考えられている。
このように、同期バーストのフレームは、普通、同期バーストは送信開始後、１フレーム〜数フレームが連続して送信される。例えば、後述する図１６は、第１フレームと第２フレームとが同期バーストのフレームである。
【０００４】
図１４は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に基づいた送信パターンである。SB0とSB1は同期バーストフレーム、TCH0，TCH1，TCH2，‥‥‥，TCHMは主データを載せた通信チャネルフレームである。ここで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に入力されるベースバンド信号は、１シンボル当たりNov回オーバーサンプリングされた複素数（実数部を同相成分(I)、虚数部を直交成分(Q)とする）のデジタル信号である。このオーバーサンプル数Novは２以上の整数である。
入力端子201より入力したベースバンド信号は、相関回路202に入力し、プリアンブルパターンとの相関r(n)を演算し、演算結果を振幅自乗値演算回路203に与える。振幅自乗値演算回路203は、相関r(n)の振幅の自乗値 |r(n)|2を演算し、短時間平均演算回路204に与える。短時間平均演算回路204は相関の自乗値|r(n)|2の短時間平均rave(n)を演算し判定器205に与える。判定器205は、短時間平均rave(n)を所定のしきい値と比較し、短時間平均rave(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の出力である短時間平均rave(n)の時間波形を示す図である。(a)は受信信号、(b)はフェージングなしの場合の短時間平均rave(n)、(c)はフェージング時での短時間平均rave(n)を表す。また、一点鎖線はしきい値である。
図６(a)に示すように、入力端子201に区間601でプリアンブル以外の信号が入力され、区間602でプリアンブルが入力され、区間603で再びプリアンブル以外の信号が入力される場合、短時間平均演算回路204の出力rave(n)は、フェージングなしの場合に図６(b)のようになり、フェージング時には図６(c)のようになる。
【００１３】
短時間平均rave(n)は、相関回路202の段数Nと短時間平均演算回路204の時定数により遅延する。フェージングがない場合には、図６(b)に示すように、高い相関を示すので、識別漏れの確率と誤識別の確率を共に“0”にするしきい値の設定を容易に行うことができる。
しかし、フェージング時には、図６(c)に示すように、受信信号の振幅の落ち込みにより、プリアンブル入力時の相関が落ち込むため、識別漏れの確率を“0”に近づけようとしきい値を低く設定すると誤識別の確率が大きくなり、反対に誤識別の確率を“0”に近づけようとしきい値を高く設定すると識別漏れの確率が大きくなり、両方の確率を“0”にするしきい値を設定することが困難になる。
【００１４】
また、プリアンブルの信号が持つ周波数成分は、シンボルレート（シンボル間隔の逆数）をfbとすると、変調方式がQPSKまたはQAMの場合には、-fb/2と+fb/2の２つの周波数成分であり、変調方式がπ/4シフトQPSKの場合には、-3fb/8と+fb/8の２つの周波数成分である。送信機と受信機との間には、周波数偏差Δfがあり、この周波数偏差Δfが大きくて、プリアンブルの片方の周波数成分が受信機のフィルタの帯域外へ外れる場合には、相関回路202に入力されるプリアンブルの片方の周波数成分が削られるため、周波数偏差Δfによる位相回転の補正を行ったとしても、相関による識別が難しくなる。
【００１５】
【発明が解決しようとする課題】
前述の従来技術には、相関演算結果が受信信号の振幅変動に影響され、識別漏れの確率と誤識別の確率とが相反する傾向にあり、両方の確率を0に近づけるようなしきい値の設定が困難であり、しきい値による識別が難しい欠点があった。
また、送信機と受信機との間で、周波数偏差が大きい場合、プリアンブルの片方の周波数成分が受信機のフィルタの帯域外へ外れ、プリアンブルの２つの周波数成分の片方の成分が削られた場合の識別が困難であるという欠点があった。
本発明の目的は、上記のような欠点を除去し、受信信号の振幅変動に影響されず、送信機と受信機との間で、周波数偏差が大きい場合においても、プリアンブルの識別を容易に行うことができるプリアンブルパターン識別方法を提供することにある。
【００１６】
【課題を解決するための手段】
上記の目的を達成するため、本発明のプリアンブルパターン識別方法は、
(1) 受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2を求め、求めたそれぞれの周波数成分の電力和p1+p2と全周波数成分の電力pとの電力比(p1+p2)/pが所定のしきい値を超えていれば、受信したベースバンド信号がプリアンブルであると判断する。尚、以下の説明で、周波数のずれていない送信ベースバンドでのプリアンブルの周波数成分を、プリアンブル送信時の周波数成分と呼ぶ。
(2) 電力比(p1+p2)/pが所定のしきい値を超え、かつ２つの成分の周波数f1とf2の差の絶対値|f2-f1|が所定の値の範囲内であれば、受信したベースバンド信号がプリアンブルであると判断する。
(3) 抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、２つの周波数f1とf2の成分が残っていると判断した場合で、電力比(p1+p2)/pが所定のしきい値を超え、かつ周波数差の絶対値|f2-f1|が所定の値の範囲内であれば、受信したベースバンド信号がプリアンブルであると判断する。ただし、r1，r2は正の定数で、r1<1，r2>1とする。
【００１７】
(4) また、本発明の周波数偏差検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出して、低い周波数成分の周波数をf1、高い周波数成分の周波数をf2とし、プリアンブル送信時の２成分の周波数のうち低い周波数成分の周波数をf10、高い周波数成分の周波数をf20とし、各成分の周波数偏差f1-f10, f2-f20の平均[(f1-f10)+(f2-f20)]/2を、送信機と受信機との周波数偏差とする。
(5) 更に、本発明の周波数偏差検出方法は、抽出した２成分のうちの一方の成分の周波数をf1'、プリアンブル送信時の２成分の周波数のうち低い周波数成分の周波数をf10、高い周波数成分の周波数をf20とし、f1'-f10またはf1'-f20を、送信機と受信機との周波数偏差とする。
(6) また、本発明の周波数偏差検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断して、２つの周波数f1とf2の成分が残っていると判断した場合には、各成分の周波数偏差f1-f10とf2-f20の平均[(f1-f10)+(f2-f20)]/2を、送信機と受信機との周波数偏差とし、電力比p1/p2が規定値r1未満の場合または規定値r2を超えている場合で、電力の大きい方の成分の周波数をf1'とし、f1'≧0であれば、２成分のうち周波数の低い成分のみが残っているものとし、f1'-f10を送信機と受信機との周波数偏差とし、電力比p1/p2が規定値r1未満の場合または規定値r2を超えている場合で、電力の大きい方の成分の周波数をf1'とし、f1'＜0であれば、２成分のうち周波数の高い成分のみが残っているものとし、f1'-f20を送信機と受信機との周波数偏差とする。
【００１８】
(7) 更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの成分の位相φ1とφ2の差φ1-φ2よりシンボルタイミングを求める。
(8) また更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、抽出したそれぞれの成分の位相φ1とφ2の位相差φ1-φ2よりシンボルタイミングを求める。
(9) また本発明のシンボルタイミング検出方法は、変調方式がQPSKまたはQAMであるベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから周波数成分を抽出し、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1とし、高い周波数の成分の位相をφ2とし、シンボルレートをfbとし、t0=(φ1-φ2)/(2πfb)をシンボルタイミングt0とする。
(10) 更に本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、２つの周波数f1とf2の成分が残っていると判断した場合、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1、高い周波数の成分の位相をφ2、シンボルレートをfbとし、t0=(φ1-φ2)/(2πfb)をシンボルタイミングとする。
(11) また、本発明のシンボルタイミング検出方法は、変調方式がπ/4シフトQPSKであるベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1、高い周波数の成分の位相をφ2、シンボルレートをfbとし、t0=(φ1-φ2-π/2)/(πfb)をシンボルタイミングとする。
(12) 更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、２つの周波数f1とf2の成分が残っていると判断した場合、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1、高い周波数の成分の位相をφ2、シンボルレートをfbとし、t0=(φ1-φ2-π/2)/(πfb)をシンボルタイミングとする。
【００１９】
(13) また更に、本発明のプリアンブル成分抽出方法は、受信信号を１シンボル当たりNov回（Novは2以上の整数）オーバーサンプリングしたベースバンド信号から、任意の連続したNwinサンプル（Nwinは２以上の整数）のベースバンド信号を抽出し、抽出したベースバンド信号に任意の窓関数を乗じ、任意の窓関数を乗じたベースバンド信号x(n) （n=0，1，‥‥‥，(Nwin-1)）のDFT（離散フーリエ変換: Discrete Fourier Transform）X(k)（k=0，1，‥‥‥，(Nwin-1)）を演算し（X(k)は複素数）、DFT出力X(k)（k=0，1，‥‥‥，(Nwin-1)）の全電力p=|X(0)|2+|X(1)|2+…+|X(Nwin-1)|2を算出し、DFT出力X(k)のうち、振幅の自乗値|X(k)|2が最大となるX(k)を検索してそのkをk11とし、X(k11)に隣接する成分から、|X(k11)|2の次に|X(k)|2が大きくなるX(k)を検索してそのkをk12とし、X(k11)とX(k12)の電力和p1=|X(k11)|2+|X(k12)|2を算出し、X(k11)の周辺を補間し、補間した部分のピークにおける周波数をf1、位相をφ1とし、変調方式により定められる自然数の定数をkdとし、X(k)（k=(k11-kd-1) mod Nwin，(k11-kd) mod Nwin，(k11-kd+1) mod Nwin，(k11+kd-1) mod Nwin，(k11+kd) mod Nwin，(k11+kd+1) mod Nwin）から|X(k)|2が最大となるX(k)を検索してそのkをk21とし、X(k21)に隣接する成分のうち、|X(k21)|2の次に|X(k)|2が大きくなるX(k)を検索してそのkをk22とし、X(k21)とX(k22)の電力和p2=|X(k21)|2+|X(k22)|2を算出し、X(k21)の周辺を補間し、補間した部分のピークにおける周波数をf2、位相をφ2とし、算出したf1, f2を比較し、f1>f2であれば、f1，φ1，p1を改めてf2，φ2，p2とし、f2，φ2，p2を改めてf1，φ1，p1とする。
【００２０】
(14) 更にまた、本発明のプリアンブル成分抽出方法は、FFT（高速フーリエ変換：Fast Fourier Transform）アルゴリズムを用いて、任意の窓関数を乗じたベースバンド信号x(n)（n=0，1，‥‥‥，(Nwin-1)）のDFT（離散フーリエ変換：Discrete Fourier Transform）X(k)（k=0，1，‥‥‥，(Nwin-1)）を演算する。
以上において、また、周波数スペクトルから抽出する２つの周波数f1とf2の成分は、振幅が高い順に２つ抽出する
【００２１】
【発明の実施の形態】
本発明の実施形態の説明に入る前に、プリアンブル信号の周波数スペクトルについて説明する。
最初に変調方式が、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はシンボルレートをfbとしωb=2πfbである。
【数３】
【００２２】
送信機が送信する信号は、式(4)に示す伝達関数RR(ω)を持つルートロールオフフィルタで帯域制限され、受信機にも同じRR(ω)の特性のルートロールオフフィルタを持ち、このフィルタの出力信号を復調する。プリアンブルパターン識別では処理遅延を削減するため、受信機のルートロールオフフィルタ前段の信号を用いる。
【数４】
ここで、αをロールオフ率と呼び、0≦α≦1である。
【００２３】
図７(b)は、図７(a)の信号を、式(4)に示す伝達関数RR(ω)を持つルートロールオフフィルタで帯域制限した周波数スペクトルであり、RR(ωb/2) = RR(-ωb/2) = 1/√2であるから、図７(b)の各成分は、式(5)のようになり、ω=±ωb/2の2つの周波数成分を持つ。
【数５】
従って、時間波形は、式(6)のようになる。
【数６】
【００２４】
受信ベースバンド信号がシンボルタイミングより時間t0だけ遅れていて、送信機と受信機との周波数偏差をΔωとすると、式(7)のようになる。また、各成分の角周波数はω1=-ωb/2+Δω、ω2=ωb/2+Δω、ω1とω2各成分の位相φ1，φ2は式(8)のようになる。ここで、θ0は定数である。
【数７】
【数８】
【００２５】
受信ベースバンド信号がプリアンブルの信号である場合の周波数スペクトルは、上記ω1とω2の成分のみであるから、周波数スペクトルの全電力をp、ω1成分の電力をp1、ω2成分の電力をp2とすると、２成分の電力和(p1+p2)が全電力pに近い値となることから、電力比(p1+p2)/pによりプリアンブルパターンの識別ができる。
また、式(8)より(φ1-φ2)を求めてθ0を消去すると、各成分の位相φ1とφ2からシンボルタイミングt0を求めることができる。
【数９】
【００２６】
次に、変調方式がπ/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.9242) / (1+0.9242+0.3832) = 0.927］を占めるため、-3ωb/8の成分の電力をp1とし、ωb/8の成分の電力をp2とし、周波数スペクトルの全電力pに対する比(p1+p2)/pによりプリアンブルパターン識別ができる。
【００２９】
図８(b)及び図８(c)よりプリアンブルの信号の時間波形は、ロールオフ率αが0≦α≦1/4の場合には、式(13)のようになり、ロールオフ率αが1/4≦α≦3/4の場合には、式(14)のようになる。
【数１３】
【数１４】
更に、QPSKまたはQAMの場合と同様に、受信ベースバンド信号がシンボルタイミングより時間t0だけ遅れていて、送信機と受信機との間にΔωの周波数偏差がある場合の時間波形を求めると、ロールオフ率αが0≦α≦1/4の場合には式(15)のようになり、ロールオフ率αが1/4≦α≦3/4の場合には、式(16)のようになる。
【数１５】
【数１６】
【００３０】
式(15)と(16)に示す-3ωb/8+Δω、ωb/8+Δωの各成分の位相φ1、φ2は、式(17)のようになる。ただし、θ0は定数である。
【数１７】
従って、シンボルタイミングt0は数式(18)より求められる。
【数１８】
【００３１】
プリアンブル以外の信号、例えば、通信チャネルTchが受信された場合には、複数の周波数スペクトルが抽出されるが、その場合には、スペクトルの振幅が大きい順に２つ選択し、それらをp1とp2とする。そして、２成分の電力和(p1+p2)を求めると、全電力pに比べて小さな値となることから、電力比(p1+p2)/pが小さくなることによりプリアンブルパターンではないと判定することができる。
【００３２】
次に、本発明の実施形態を、図１を用いて説明する。図１は、本発明のプリアンブルパターン識別、周波数偏差検出、及びシンボルタイミング検出を行なう回路構成の一実施例を示すブロック図である。
図１において、受信されたベースバンド信号が入力端子101に入力される。このベースバンド信号は、従来技術の場合と同様に、１シンボル当たりNov回オーバーサンプリングされた複素数のデジタル信号であり、Novは２以上の正整数である。
【００３３】
入力端子101を介して入力したベースバンド信号は、サンプル毎にシフトレジスタ114に入力される。シフトレジスタ114の段数は、Nwinで、シフトレジスタ114内部のNwin個のレジスタをs(0)，s(1)，‥‥‥，s(Nwin-2)，s(Nwin-1)とすると、入力端子101を介してベースバンド信号が入力されるたびに、s(Nwin-1)，s(Nwin-2)，‥‥‥，s(1)の内容が、それぞれ、s(Nwin-2)，s(Nwin-3)，‥‥‥，s(0)にシフトされ、入力端子101を介して入力されるベースバンド信号がs(Nwin-1)に入力される。シフトレジスタ114の内容s(0)，s(1)，‥‥‥，s(Nwin-2)，s(Nwin-1)は、サンプル毎に、スイッチ113に与えられ、一定の時間間隔Nstep（サンプル）毎に１回、s(0)，s(1)，‥‥‥，s(Nwin-2)，s(Nwin-1)が窓掛け回路102に入力される。（例えば、オーバーサンプル数Nov=2で、Nstep = 8、Nwin = 16である）。このように、スイッチ113は、Nstepサンプル周期に１回閉じていて、Nwinサンプルの信号をNstep間隔で取り出す。
スイッチ113から入力されたNwinサンプルの信号は、窓掛け回路102によって、窓（例えば、ハミング窓（Hamming Window））を掛けられ、FFT演算回路103に入力される。
【００３４】
FFT演算回路103は、窓掛け回路102から入力されるNwinサンプルの信号のDFT（離散フーリエ変換：Discrete Fourier Transform）を、FFT（高速フーリエ変換：Fast Fourier Transform）アルゴリズムにより演算する。
図９は、Nov = 2、Nwin = 16の場合のFFT演算回路103の出力の一実施例で、送信機と受信機との間に周波数偏差がある場合のプリアンブル信号のDFT出力を示す図である。
図９は、角周波数ω1とω2の２つの周波数成分を持つプリアンブル信号が入力された場合のDFT出力であり、-3ωb/8の周波数成分であるX(13)と、5ωb/8の周波数成分であるX(5)がピークとなっている。しかし、DFT出力は、周波数が離散的であるため、実際のプリアンブルの成分ω1（ピーク903）とω2（ピーク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，‥‥‥，(Nwin-1)）から、|X(k)|2が最大となるものを検索し、検索された|X(k)|2が最大となるkをk11（k = k11）としてステップ1002に進む。ステップ1002では、m1とm2を算出し（m1 = (k11 - 1) mod Nwin，m2 = (k11 + 1) mod Nwin）、ステップ1003に進む。
ステップ1003では、算出されたm1とm2から、|X(m1)|2と|X(m2)|2を比較し、|X(m2)|2＞|X(m1)|2であれば、ステップ1004に進み k12 = m2としステップ1006に進む。また、|X(m2)|2≦|X(m1)|2であれば、ステップ1005に進み k12 = m1とし、ステップ1006に進む。ステップ1006では、X(k11)とX(k12)の電力和p1（p1 = |X(k11)|2 + |X(k12)|2）を算出して、ステップ1007に進む。
受信したプリアンブルの実際のピークは、X(k11)とX(k12)の間に存在するため、ステップ1007では、X(m1)、X(k11)、X(m2)の３点でポイント間の補間を行い、補間を行った部分での実際のピークを検索し、その実際のピークでの周波数f1と周波数スペクトルの位相φ1を算出する。
【００３７】
次にもう１つの成分を抽出する。kdはプリアンブルの２成分の周波数差|ω2-ω1|に相当するkの値である。変調方式がQPSKまたはQAMの場合には、|ω2-ω1| = ωbであるから、kd=Nwin/Nov、π/4シフトQPSKの場合には、|ω2-ω1|=ωb/2であるから、kd=Nwin/(2Nov)である。
ステップ1008では、m1=(k11-kd) mod Nwinとm2=(k11+kd) mod Nwinとを算出し、ステップ1009に進む。ステップ1009では、|X(k)|2 (k=(m1-1) mod Nwin，m1，(m1+1) mod Nwin)の最大値を検索し、第１の最大値をx1とし、その時のkをm1（k = m1）としステップ1010に進む。
ステップ1010では、また、|X(k)|2 (k=(m2-1) mod Nwin，m2，(m2+1) mod Nwin)の最大値を検索し、第２の最大値をx2とし、その時のkをm2（k = m2）とする。
【００３８】
ステップ1011では、第１の最大値x1と第２の最大値x2とを比較し、第１の最大値x1が大きければ、ステップ1012で k21=m1を算出しステップ1014に進む。また、そうでなければ、ステップ1013で k21=m2を算出しステップ1014に進む。これにより、k21がもう１つの成分の番号となる。
ステップ1014では、m1=(k21-1) mod Nwinとm2=(k21+1) mod Nwinを算出し、ステップ1015に進む。ステップ1015では、|X(m1)|2と|X(m2)|2を比較する。|X(m1)|2より|X(m2)|2が大きい場合には、ステップ1016でk22=m2としてステップ1018に進む。また、そうでない場合には、ステップ1017でk22=m1としてステップ1018に進む。
ステップ1018では、X(k21)とX(k22)の電力和p2を算出し、ステップ1019に進む。
【００３９】
受信したプリアンブルの実際のピークは、X(k21)とX(k22)の間に存在するため、ステップ1019において、X(m1)、X(k21)、X(m2)の３点でポイント間の補間を行い、補間を行った部分でのピークを検索し、そのピークでの周波数f2と周波数スペクトルの位相φ2を算出し、ステップ1020に進む。
最後に、ステップ1020では、FFT演算回路103の出力の全電力p = |X(0)|2+|X(1)|2+…+|X(Nwin-1)|2を算出する。
【００４０】
ステップ1001〜1020の処理は、図９に示したFFT回路103の出力の場合には、ステップ1001でX(5)が選択され、k11=5、m1=4，m2=6で|X(4)|2と|X(6)|2を比較し（ステップ1003）、|X(4)|2が大きいためk12=4となり（ステップ1005）、p1=|X(5)|2+|X(4)|2が演算される（ステップ1006）。
受信したプリアンブルの実際のピークは，ω2の位置に存在するため、X(4)、X(5)、X(6)の３点の間で補間を行い（破線902）、補間した部分でのピークがω2で検出されるので、その周波数をf1=ω2/(2π)とする。また、ω2での補間値の位相をφ1とする（ステップ1007）。
【００４１】
図９は、変調方式がQPSKまたはQAMの場合の一例であり、kd=Nwin/Nov=8であるから、m1=13, m2=13と同じ値であるため（ステップ1008）、ステップ1009とステップ1010において、|X(k)|2 (k=12，13，14)の最大値を検索し、x1=x2=|X(13)|2，m1=m2=13となり（ステップ1009，1010）、k21=m2=13がもう一つのピーク成分の番号となる（ステップ1013）。m1=12，m2=14として|X(12)|2と|X(14)|2を比較し（ステップ1015）、|X(12)|2が大きいためk22=m1=12となり（ステップ1017）、p2=|X(13)|2+|X(12)|2が演算される（ステップ1018）。
受信したプリアンブルの実際のピークは、ω1の位置に存在するため、X(12)，X(13)，X(14)の３点の間で補間を行い（破線901）、補間した部分でのピークがω1で検出されるので、その周波数をf2=ω1/(2π)とする。また、ω1での補間値の位相をφ2とする（ステップ1019）。
【００４２】
図１のプリアンブルパターン識別回路107、周波数偏差検出回路108、及び、シンボルタイミング検出回路109は、図１０のフローチャートの手順で算出した、DFT出力の全電力p、２つの周波数成分の周波数f1とf2、位相φ1とφ2、電力p1とp2とを用いて、プリアンブルパターン識別、周波数偏差検出、シンボルタイミング検出を行う。
【００４３】
プリアンブルパターン識別回路107の処理は、全電力pに対する２成分の電力和 p1+p2の比γ=(p1+p2)/pを用い、更にプリアンブルの２成分が受信機のフィルタにより削られていいない場合は、２成分の周波数差 f12=|f2-f1|も用いる。入力端子101にプリアンブル信号が入力されていて、２成分ともに削られない場合は、全電力pは、２成分電力和 p1+p2と雑音電力との和であり（ただし、変調方式がπ/4シフトQPSKでロールオフ率αが1/4≦α≦1/2の場合はp1，p2と雑音以外に最大7.3%の成分の電力がpに含まれる）、２成分の内の片方の成分が削られる場合は、p1またはp2のいずれかが雑音電力となる。
従って、片方の成分が削られているか否かの判断は、p1とp2の比較により行う。例えば、p1とp2を比較し、片方が6dB以上低い値（比が4倍以上）ならば片方の成分が削られていると判断する。
【００４４】
変調方式がQPSKまたはQAMの場合で、２成分とも削られていないと判断した場合には、γがしきい値rthを超え、かつf12がfb±f0（fbはシンボルレート、f0は規定値）の範囲内ならば受信信号がプリアンブル信号であると識別する。片方の成分が削られていると判断した場合は、γとしきい値rthの比較のみでプリアンブルパターン識別を行う。
変調方式がπ/4シフトQPSKの場合で、２成分とも削られていないと判断した場合には、γがしきい値rthを超え、かつf12がfb/2±f0の範囲内ならば受信信号がプリアンブル信号であると識別する。片方の成分が削られていると判断した場合は、γとしきい値rthの比較のみでプリアンブルパターン識別を行う。
しきい値rthと規定値f0は、例えばrth=0.8、f0=fb/8に設定する。
受信信号がプリアンブルであると識別した場合には、プリアンブルパターン識別出力端子110を介して“1”を出力し、識別できない場合には“0”を出力する。
【００４５】
周波数偏差検出回路108の処理は、２成分の各々の周波数f1，f2と、プリアンブルの信号の送信ベースバンドでの周波数成分の各成分の周波数f10，f20（低い周波数の成分をf10、高い周波数の成分をf20とする）を用いて行う。２成分とも削られていない場合は両方の周波数を用いて式(19)により算出する。
【数１９】
また、片方の成分が削られている場合は残っている成分の周波数をf1'とし、このf1'を用いて式(20)により算出する。ここで、f1'≧0の場合は周波数の高い成分が削られているものとして、f10との差をΔfとし、f1'<0の場合は周波数の低い成分が削られているものとして、f20との差をΔfとする。
【数２０】
検出した周波数偏差Δfは、周波数偏差出力端子111を介して出力される。
【００４６】
シンボルタイミング検出回路109の処理は、プリアンブルの２成分とも削られていないと判断した場合にのみ行い、２成分の位相φ1とφ2を用いて行う。
変調方式がQPSKまたはQAMの場合、式(9)よりシンボルタイミングfbt0は、式(21)のようになる。
【数２１】
また、変調方式がπ/4シフトQPSKの場合には、式(18)よりシンボルタイミングfbt0は、式(22)のようになる。
【数２２】
検出したシンボルタイミングfbt0は、シンボルタイミング出力端子112を介して出力される。
【００４７】
以上の説明による本発明の実施形態では、入力端子101にプリアンブル以外の信号が入力された場合、全電力pは、２つの周波数スペクトル成分の電力p1とp2以外に信号成分の電力を含んでいるため、電力比γの値は小さくなる。
また、入力端子101にプリアンブルが入力された場合は、振幅が変動したとしても、p1+p2は振幅変動の影響を受けたプリアンブル成分の電力で、全電力pは、振幅変動の影響を受けたプリアンブル成分の電力と雑音電力の和であり、振幅の落ち込みによりp1+p2が小さくなる場合は、pも小さくなる。このため、雑音の影響を受けるものの、振幅変動は直接影響しない。従って、振幅変動の影響を受けないプリアンブルパターン識別が可能である。
【００４８】
更に、送信機と受信機との周波数偏差により、プリアンブルの２つの周波数成分のうち片方の成分が、受信機のフィルタの帯域外へ外れたとしても、p1とp2の一方がプリアンブルの１成分の電力で、他が雑音電力であり、全電力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：ベースバンド信号出力端子。
【発明の属する技術分野】
本発明は、受信機に関わり、特に受信信号に既知の信号パターンが含まれているか否かを識別する方法に関するものである。
【０００２】
【従来の技術】
デジタル無線通信システムにおいて、送信機は、伝送する情報（２値信号のビット列）を一定のビット数で分割し、分割されたビット列を含む所定の構成がなされたフレームとし、そのフレーム毎に、例えば16QAM（Quadrature Amplitude Modulation）方式のようなデジタル変調を行い、デジタル変調された信号として送信する。
送信をしていない状態から、送信機が初めて送信を開始するフレームを第１フレームという。以下第１フレームに続いて、第２フレーム、第３フレーム、‥‥‥と送信されていく。
【０００３】
フレームには、受信機が受信し、受信したフレームを用いて、送信される信号に対する同期を行うために、プリアンブルと呼ばれる規定の固定ビットパターンが、定められた位置に、定められた長さで付加されている。
受信機は、送信機が送信した第１フレームの信号を受信する際に、第１フレーム内のプリアンブルを用いて、送信された信号に対して高速同期を行う。ここで、高速同期とは、全く同期が確立されていない状態から同期を行い、同期した情報を現フレームの復調（検波）動作に反映させる処理のことをいう。これに対し、通常同期とは、１つ前のフレームですでに同期が確立されており、現フレームと１つ前のフレームとの間での変動分に対する同期を低速で行い、同期した情報を次のフレームの復調動作に反映する処理のことをいう。
第１フレームは、同期バーストと呼ばれるフレームであり、このフレームにプリアンブルが含まれている。受信機は、例えば、この第１フレームの同期バーストを用いて高速同期を行い、第２フレーム以降では通常同期を行う必要がある。そこで、プリアンブルが含まれている同期バーストを受信しているか否かの判断を行うために、プリアンブルのパターンを識別する方法が考えられている。
このように、同期バーストのフレームは、普通、同期バーストは送信開始後、１フレーム〜数フレームが連続して送信される。例えば、後述する図１６は、第１フレームと第２フレームとが同期バーストのフレームである。
【０００４】
図１４は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に基づいた送信パターンである。SB0とSB1は同期バーストフレーム、TCH0，TCH1，TCH2，‥‥‥，TCHMは主データを載せた通信チャネルフレームである。ここで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に入力されるベースバンド信号は、１シンボル当たりNov回オーバーサンプリングされた複素数（実数部を同相成分(I)、虚数部を直交成分(Q)とする）のデジタル信号である。このオーバーサンプル数Novは２以上の整数である。
入力端子201より入力したベースバンド信号は、相関回路202に入力し、プリアンブルパターンとの相関r(n)を演算し、演算結果を振幅自乗値演算回路203に与える。振幅自乗値演算回路203は、相関r(n)の振幅の自乗値 |r(n)|2を演算し、短時間平均演算回路204に与える。短時間平均演算回路204は相関の自乗値|r(n)|2の短時間平均rave(n)を演算し判定器205に与える。判定器205は、短時間平均rave(n)を所定のしきい値と比較し、短時間平均rave(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の出力である短時間平均rave(n)の時間波形を示す図である。(a)は受信信号、(b)はフェージングなしの場合の短時間平均rave(n)、(c)はフェージング時での短時間平均rave(n)を表す。また、一点鎖線はしきい値である。
図６(a)に示すように、入力端子201に区間601でプリアンブル以外の信号が入力され、区間602でプリアンブルが入力され、区間603で再びプリアンブル以外の信号が入力される場合、短時間平均演算回路204の出力rave(n)は、フェージングなしの場合に図６(b)のようになり、フェージング時には図６(c)のようになる。
【００１３】
短時間平均rave(n)は、相関回路202の段数Nと短時間平均演算回路204の時定数により遅延する。フェージングがない場合には、図６(b)に示すように、高い相関を示すので、識別漏れの確率と誤識別の確率を共に“0”にするしきい値の設定を容易に行うことができる。
しかし、フェージング時には、図６(c)に示すように、受信信号の振幅の落ち込みにより、プリアンブル入力時の相関が落ち込むため、識別漏れの確率を“0”に近づけようとしきい値を低く設定すると誤識別の確率が大きくなり、反対に誤識別の確率を“0”に近づけようとしきい値を高く設定すると識別漏れの確率が大きくなり、両方の確率を“0”にするしきい値を設定することが困難になる。
【００１４】
また、プリアンブルの信号が持つ周波数成分は、シンボルレート（シンボル間隔の逆数）をfbとすると、変調方式がQPSKまたはQAMの場合には、-fb/2と+fb/2の２つの周波数成分であり、変調方式がπ/4シフトQPSKの場合には、-3fb/8と+fb/8の２つの周波数成分である。送信機と受信機との間には、周波数偏差Δfがあり、この周波数偏差Δfが大きくて、プリアンブルの片方の周波数成分が受信機のフィルタの帯域外へ外れる場合には、相関回路202に入力されるプリアンブルの片方の周波数成分が削られるため、周波数偏差Δfによる位相回転の補正を行ったとしても、相関による識別が難しくなる。
【００１５】
【発明が解決しようとする課題】
前述の従来技術には、相関演算結果が受信信号の振幅変動に影響され、識別漏れの確率と誤識別の確率とが相反する傾向にあり、両方の確率を0に近づけるようなしきい値の設定が困難であり、しきい値による識別が難しい欠点があった。
また、送信機と受信機との間で、周波数偏差が大きい場合、プリアンブルの片方の周波数成分が受信機のフィルタの帯域外へ外れ、プリアンブルの２つの周波数成分の片方の成分が削られた場合の識別が困難であるという欠点があった。
本発明の目的は、上記のような欠点を除去し、受信信号の振幅変動に影響されず、送信機と受信機との間で、周波数偏差が大きい場合においても、プリアンブルの識別を容易に行うことができるプリアンブルパターン識別方法を提供することにある。
【００１６】
【課題を解決するための手段】
上記の目的を達成するため、本発明のプリアンブルパターン識別方法は、
(1) 受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2を求め、求めたそれぞれの周波数成分の電力和p1+p2と全周波数成分の電力pとの電力比(p1+p2)/pが所定のしきい値を超えていれば、受信したベースバンド信号がプリアンブルであると判断する。尚、以下の説明で、周波数のずれていない送信ベースバンドでのプリアンブルの周波数成分を、プリアンブル送信時の周波数成分と呼ぶ。
(2) 電力比(p1+p2)/pが所定のしきい値を超え、かつ２つの成分の周波数f1とf2の差の絶対値|f2-f1|が所定の値の範囲内であれば、受信したベースバンド信号がプリアンブルであると判断する。
(3) 抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、２つの周波数f1とf2の成分が残っていると判断した場合で、電力比(p1+p2)/pが所定のしきい値を超え、かつ周波数差の絶対値|f2-f1|が所定の値の範囲内であれば、受信したベースバンド信号がプリアンブルであると判断する。ただし、r1，r2は正の定数で、r1<1，r2>1とする。
【００１７】
(4) また、本発明の周波数偏差検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出して、低い周波数成分の周波数をf1、高い周波数成分の周波数をf2とし、プリアンブル送信時の２成分の周波数のうち低い周波数成分の周波数をf10、高い周波数成分の周波数をf20とし、各成分の周波数偏差f1-f10, f2-f20の平均[(f1-f10)+(f2-f20)]/2を、送信機と受信機との周波数偏差とする。
(5) 更に、本発明の周波数偏差検出方法は、抽出した２成分のうちの一方の成分の周波数をf1'、プリアンブル送信時の２成分の周波数のうち低い周波数成分の周波数をf10、高い周波数成分の周波数をf20とし、f1'-f10またはf1'-f20を、送信機と受信機との周波数偏差とする。
(6) また、本発明の周波数偏差検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断して、２つの周波数f1とf2の成分が残っていると判断した場合には、各成分の周波数偏差f1-f10とf2-f20の平均[(f1-f10)+(f2-f20)]/2を、送信機と受信機との周波数偏差とし、電力比p1/p2が規定値r1未満の場合または規定値r2を超えている場合で、電力の大きい方の成分の周波数をf1'とし、f1'≧0であれば、２成分のうち周波数の低い成分のみが残っているものとし、f1'-f10を送信機と受信機との周波数偏差とし、電力比p1/p2が規定値r1未満の場合または規定値r2を超えている場合で、電力の大きい方の成分の周波数をf1'とし、f1'＜0であれば、２成分のうち周波数の高い成分のみが残っているものとし、f1'-f20を送信機と受信機との周波数偏差とする。
【００１８】
(7) 更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの成分の位相φ1とφ2の差φ1-φ2よりシンボルタイミングを求める。
(8) また更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、抽出したそれぞれの成分の位相φ1とφ2の位相差φ1-φ2よりシンボルタイミングを求める。
(9) また本発明のシンボルタイミング検出方法は、変調方式がQPSKまたはQAMであるベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから周波数成分を抽出し、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1とし、高い周波数の成分の位相をφ2とし、シンボルレートをfbとし、t0=(φ1-φ2)/(2πfb)をシンボルタイミングt0とする。
(10) 更に本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、２つの周波数f1とf2の成分が残っていると判断した場合、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1、高い周波数の成分の位相をφ2、シンボルレートをfbとし、t0=(φ1-φ2)/(2πfb)をシンボルタイミングとする。
(11) また、本発明のシンボルタイミング検出方法は、変調方式がπ/4シフトQPSKであるベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1、高い周波数の成分の位相をφ2、シンボルレートをfbとし、t0=(φ1-φ2-π/2)/(πfb)をシンボルタイミングとする。
(12) 更に、本発明のシンボルタイミング検出方法は、受信したベースバンド信号の周波数スペクトルを求め、求めた周波数スペクトルから２つの周波数f1とf2の成分を抽出し、抽出したそれぞれの周波数成分の電力p1とp2の電力比p1/p2を算出し、電力比p1/p2が規定値r1以上かつ規定値r2以下ならば２つの周波数f1とf2の成分が残っていると判断し、２つの周波数f1とf2の成分が残っていると判断した場合、抽出した２つの周波数f1とf2の成分のうち、低い周波数の成分の位相をφ1、高い周波数の成分の位相をφ2、シンボルレートをfbとし、t0=(φ1-φ2-π/2)/(πfb)をシンボルタイミングとする。
【００１９】
(13) また更に、本発明のプリアンブル成分抽出方法は、受信信号を１シンボル当たりNov回（Novは2以上の整数）オーバーサンプリングしたベースバンド信号から、任意の連続したNwinサンプル（Nwinは２以上の整数）のベースバンド信号を抽出し、抽出したベースバンド信号に任意の窓関数を乗じ、任意の窓関数を乗じたベースバンド信号x(n) （n=0，1，‥‥‥，(Nwin-1)）のDFT（離散フーリエ変換: Discrete Fourier Transform）X(k)（k=0，1，‥‥‥，(Nwin-1)）を演算し（X(k)は複素数）、DFT出力X(k)（k=0，1，‥‥‥，(Nwin-1)）の全電力p=|X(0)|2+|X(1)|2+…+|X(Nwin-1)|2を算出し、DFT出力X(k)のうち、振幅の自乗値|X(k)|2が最大となるX(k)を検索してそのkをk11とし、X(k11)に隣接する成分から、|X(k11)|2の次に|X(k)|2が大きくなるX(k)を検索してそのkをk12とし、X(k11)とX(k12)の電力和p1=|X(k11)|2+|X(k12)|2を算出し、X(k11)の周辺を補間し、補間した部分のピークにおける周波数をf1、位相をφ1とし、変調方式により定められる自然数の定数をkdとし、X(k)（k=(k11-kd-1) mod Nwin，(k11-kd) mod Nwin，(k11-kd+1) mod Nwin，(k11+kd-1) mod Nwin，(k11+kd) mod Nwin，(k11+kd+1) mod Nwin）から|X(k)|2が最大となるX(k)を検索してそのkをk21とし、X(k21)に隣接する成分のうち、|X(k21)|2の次に|X(k)|2が大きくなるX(k)を検索してそのkをk22とし、X(k21)とX(k22)の電力和p2=|X(k21)|2+|X(k22)|2を算出し、X(k21)の周辺を補間し、補間した部分のピークにおける周波数をf2、位相をφ2とし、算出したf1, f2を比較し、f1>f2であれば、f1，φ1，p1を改めてf2，φ2，p2とし、f2，φ2，p2を改めてf1，φ1，p1とする。
【００２０】
(14) 更にまた、本発明のプリアンブル成分抽出方法は、FFT（高速フーリエ変換：Fast Fourier Transform）アルゴリズムを用いて、任意の窓関数を乗じたベースバンド信号x(n)（n=0，1，‥‥‥，(Nwin-1)）のDFT（離散フーリエ変換：Discrete Fourier Transform）X(k)（k=0，1，‥‥‥，(Nwin-1)）を演算する。
以上において、また、周波数スペクトルから抽出する２つの周波数f1とf2の成分は、振幅が高い順に２つ抽出する
【００２１】
【発明の実施の形態】
本発明の実施形態の説明に入る前に、プリアンブル信号の周波数スペクトルについて説明する。
最初に変調方式が、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はシンボルレートをfbとしωb=2πfbである。
【数３】
送信機が送信する信号は、式(4)に示す伝達関数RR(ω)を持つルートロールオフフィルタで帯域制限され、受信機にも同じRR(ω)の特性のルートロールオフフィルタを持ち、このフィルタの出力信号を復調する。プリアンブルパターン識別では処理遅延を削減するため、受信機のルートロールオフフィルタ前段の信号を用いる。
【数４】
【００２３】
図７(b)は、図７(a)の信号を、式(4)に示す伝達関数RR(ω)を持つルートロールオフフィルタで帯域制限した周波数スペクトルであり、RR(ωb/2) = RR(-ωb/2) = 1/√2であるから、図７(b)の各成分は、式(5)のようになり、ω=±ωb/2の2つの周波数成分を持つ。
【数５】
【数６】
受信ベースバンド信号がシンボルタイミングより時間t0だけ遅れていて、送信機と受信機との周波数偏差をΔωとすると、式(7)のようになる。また、各成分の角周波数はω1=-ωb/2+Δω、ω2=ωb/2+Δω、ω1とω2各成分の位相φ1，φ2は式(8)のようになる。ここで、θ0は定数である。
【数７】
受信ベースバンド信号がプリアンブルの信号である場合の周波数スペクトルは、上記ω1とω2の成分のみであるから、周波数スペクトルの全電力をp、ω1成分の電力をp1、ω2成分の電力をp2とすると、２成分の電力和(p1+p2)が全電力pに近い値となることから、電力比(p1+p2)/pによりプリアンブルパターンの識別ができる。
また、式(8)より(φ1-φ2)を求めてθ0を消去すると、各成分の位相φ1とφ2からシンボルタイミングt0を求めることができる。
【数９】
次に、変調方式がπ/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)となる。
【数１１】
【数１２】
【００２９】
図８(b)及び図８(c)よりプリアンブルの信号の時間波形は、ロールオフ率αが0≦α≦1/4の場合には、式(13)のようになり、ロールオフ率αが1/4≦α≦3/4の場合には、式(14)のようになる。
【数１３】
【数１５】
式(15)と(16)に示す-3ωb/8+Δω、ωb/8+Δωの各成分の位相φ1、φ2は、式(17)のようになる。ただし、θ0は定数である。
【数１７】
【数１８】
プリアンブル以外の信号、例えば、通信チャネルTchが受信された場合には、複数の周波数スペクトルが抽出されるが、その場合には、スペクトルの振幅が大きい順に２つ選択し、それらをp1とp2とする。そして、２成分の電力和(p1+p2)を求めると、全電力pに比べて小さな値となることから、電力比(p1+p2)/pが小さくなることによりプリアンブルパターンではないと判定することができる。
【００３２】
次に、本発明の実施形態を、図１を用いて説明する。図１は、本発明のプリアンブルパターン識別、周波数偏差検出、及びシンボルタイミング検出を行なう回路構成の一実施例を示すブロック図である。
図１において、受信されたベースバンド信号が入力端子101に入力される。このベースバンド信号は、従来技術の場合と同様に、１シンボル当たりNov回オーバーサンプリングされた複素数のデジタル信号であり、Novは２以上の正整数である。
【００３３】
入力端子101を介して入力したベースバンド信号は、サンプル毎にシフトレジスタ114に入力される。シフトレジスタ114の段数は、Nwinで、シフトレジスタ114内部のNwin個のレジスタをs(0)，s(1)，‥‥‥，s(Nwin-2)，s(Nwin-1)とすると、入力端子101を介してベースバンド信号が入力されるたびに、s(Nwin-1)，s(Nwin-2)，‥‥‥，s(1)の内容が、それぞれ、s(Nwin-2)，s(Nwin-3)，‥‥‥，s(0)にシフトされ、入力端子101を介して入力されるベースバンド信号がs(Nwin-1)に入力される。シフトレジスタ114の内容s(0)，s(1)，‥‥‥，s(Nwin-2)，s(Nwin-1)は、サンプル毎に、スイッチ113に与えられ、一定の時間間隔Nstep（サンプル）毎に１回、s(0)，s(1)，‥‥‥，s(Nwin-2)，s(Nwin-1)が窓掛け回路102に入力される。（例えば、オーバーサンプル数Nov=2で、Nstep = 8、Nwin = 16である）。このように、スイッチ113は、Nstepサンプル周期に１回閉じていて、Nwinサンプルの信号をNstep間隔で取り出す。
スイッチ113から入力されたNwinサンプルの信号は、窓掛け回路102によって、窓（例えば、ハミング窓（Hamming Window））を掛けられ、FFT演算回路103に入力される。
【００３４】
FFT演算回路103は、窓掛け回路102から入力されるNwinサンプルの信号のDFT（離散フーリエ変換：Discrete Fourier Transform）を、FFT（高速フーリエ変換：Fast Fourier Transform）アルゴリズムにより演算する。
図９は、Nov = 2、Nwin = 16の場合のFFT演算回路103の出力の一実施例で、送信機と受信機との間に周波数偏差がある場合のプリアンブル信号のDFT出力を示す図である。
図９は、角周波数ω1とω2の２つの周波数成分を持つプリアンブル信号が入力された場合のDFT出力であり、-3ωb/8の周波数成分であるX(13)と、5ωb/8の周波数成分であるX(5)がピークとなっている。しかし、DFT出力は、周波数が離散的であるため、実際のプリアンブルの成分ω1（ピーク903）とω2（ピーク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，‥‥‥，(Nwin-1)）から、|X(k)|2が最大となるものを検索し、検索された|X(k)|2が最大となるkをk11（k = k11）としてステップ1002に進む。ステップ1002では、m1とm2を算出し（m1 = (k11 - 1) mod Nwin，m2 = (k11 + 1) mod Nwin）、ステップ1003に進む。
ステップ1003では、算出されたm1とm2から、|X(m1)|2と|X(m2)|2を比較し、|X(m2)|2＞|X(m1)|2であれば、ステップ1004に進み k12 = m2としステップ1006に進む。また、|X(m2)|2≦|X(m1)|2であれば、ステップ1005に進み k12 = m1とし、ステップ1006に進む。ステップ1006では、X(k11)とX(k12)の電力和p1（p1 = |X(k11)|2 + |X(k12)|2）を算出して、ステップ1007に進む。
受信したプリアンブルの実際のピークは、X(k11)とX(k12)の間に存在するため、ステップ1007では、X(m1)、X(k11)、X(m2)の３点でポイント間の補間を行い、補間を行った部分での実際のピークを検索し、その実際のピークでの周波数f1と周波数スペクトルの位相φ1を算出する。
【００３７】
次にもう１つの成分を抽出する。kdはプリアンブルの２成分の周波数差|ω2-ω1|に相当するkの値である。変調方式がQPSKまたはQAMの場合には、|ω2-ω1| = ωbであるから、kd=Nwin/Nov、π/4シフトQPSKの場合には、|ω2-ω1|=ωb/2であるから、kd=Nwin/(2Nov)である。
ステップ1008では、m1=(k11-kd) mod Nwinとm2=(k11+kd) mod Nwinとを算出し、ステップ1009に進む。ステップ1009では、|X(k)|2 (k=(m1-1) mod Nwin，m1，(m1+1) mod Nwin)の最大値を検索し、第１の最大値をx1とし、その時のkをm1（k = m1）としステップ1010に進む。
ステップ1010では、また、|X(k)|2 (k=(m2-1) mod Nwin，m2，(m2+1) mod Nwin)の最大値を検索し、第２の最大値をx2とし、その時のkをm2（k = m2）とする。
【００３８】
ステップ1011では、第１の最大値x1と第２の最大値x2とを比較し、第１の最大値x1が大きければ、ステップ1012で k21=m1を算出しステップ1014に進む。また、そうでなければ、ステップ1013で k21=m2を算出しステップ1014に進む。これにより、k21がもう１つの成分の番号となる。
ステップ1014では、m1=(k21-1) mod Nwinとm2=(k21+1) mod Nwinを算出し、ステップ1015に進む。ステップ1015では、|X(m1)|2と|X(m2)|2を比較する。|X(m1)|2より|X(m2)|2が大きい場合には、ステップ1016でk22=m2としてステップ1018に進む。また、そうでない場合には、ステップ1017でk22=m1としてステップ1018に進む。
ステップ1018では、X(k21)とX(k22)の電力和p2を算出し、ステップ1019に進む。
【００３９】
受信したプリアンブルの実際のピークは、X(k21)とX(k22)の間に存在するため、ステップ1019において、X(m1)、X(k21)、X(m2)の３点でポイント間の補間を行い、補間を行った部分でのピークを検索し、そのピークでの周波数f2と周波数スペクトルの位相φ2を算出し、ステップ1020に進む。
最後に、ステップ1020では、FFT演算回路103の出力の全電力p = |X(0)|2+|X(1)|2+…+|X(Nwin-1)|2を算出する。
【００４０】
ステップ1001〜1020の処理は、図９に示したFFT回路103の出力の場合には、ステップ1001でX(5)が選択され、k11=5、m1=4，m2=6で|X(4)|2と|X(6)|2を比較し（ステップ1003）、|X(4)|2が大きいためk12=4となり（ステップ1005）、p1=|X(5)|2+|X(4)|2が演算される（ステップ1006）。
受信したプリアンブルの実際のピークは，ω2の位置に存在するため、X(4)、X(5)、X(6)の３点の間で補間を行い（破線902）、補間した部分でのピークがω2で検出されるので、その周波数をf1=ω2/(2π)とする。また、ω2での補間値の位相をφ1とする（ステップ1007）。
【００４１】
図９は、変調方式がQPSKまたはQAMの場合の一例であり、kd=Nwin/Nov=8であるから、m1=13, m2=13と同じ値であるため（ステップ1008）、ステップ1009とステップ1010において、|X(k)|2 (k=12，13，14)の最大値を検索し、x1=x2=|X(13)|2，m1=m2=13となり（ステップ1009，1010）、k21=m2=13がもう一つのピーク成分の番号となる（ステップ1013）。m1=12，m2=14として|X(12)|2と|X(14)|2を比較し（ステップ1015）、|X(12)|2が大きいためk22=m1=12となり（ステップ1017）、p2=|X(13)|2+|X(12)|2が演算される（ステップ1018）。
受信したプリアンブルの実際のピークは、ω1の位置に存在するため、X(12)，X(13)，X(14)の３点の間で補間を行い（破線901）、補間した部分でのピークがω1で検出されるので、その周波数をf2=ω1/(2π)とする。また、ω1での補間値の位相をφ2とする（ステップ1019）。
【００４２】
図１のプリアンブルパターン識別回路107、周波数偏差検出回路108、及び、シンボルタイミング検出回路109は、図１０のフローチャートの手順で算出した、DFT出力の全電力p、２つの周波数成分の周波数f1とf2、位相φ1とφ2、電力p1とp2とを用いて、プリアンブルパターン識別、周波数偏差検出、シンボルタイミング検出を行う。
【００４３】
プリアンブルパターン識別回路107の処理は、全電力pに対する２成分の電力和 p1+p2の比γ=(p1+p2)/pを用い、更にプリアンブルの２成分が受信機のフィルタにより削られていいない場合は、２成分の周波数差 f12=|f2-f1|も用いる。入力端子101にプリアンブル信号が入力されていて、２成分ともに削られない場合は、全電力pは、２成分電力和 p1+p2と雑音電力との和であり（ただし、変調方式がπ/4シフトQPSKでロールオフ率αが1/4≦α≦1/2の場合はp1，p2と雑音以外に最大7.3%の成分の電力がpに含まれる）、２成分の内の片方の成分が削られる場合は、p1またはp2のいずれかが雑音電力となる。
従って、片方の成分が削られているか否かの判断は、p1とp2の比較により行う。例えば、p1とp2を比較し、片方が6dB以上低い値（比が4倍以上）ならば片方の成分が削られていると判断する。
【００４４】
変調方式がQPSKまたはQAMの場合で、２成分とも削られていないと判断した場合には、γがしきい値rthを超え、かつf12がfb±f0（fbはシンボルレート、f0は規定値）の範囲内ならば受信信号がプリアンブル信号であると識別する。片方の成分が削られていると判断した場合は、γとしきい値rthの比較のみでプリアンブルパターン識別を行う。
変調方式がπ/4シフトQPSKの場合で、２成分とも削られていないと判断した場合には、γがしきい値rthを超え、かつf12がfb/2±f0の範囲内ならば受信信号がプリアンブル信号であると識別する。片方の成分が削られていると判断した場合は、γとしきい値rthの比較のみでプリアンブルパターン識別を行う。
しきい値rthと規定値f0は、例えばrth=0.8、f0=fb/8に設定する。
受信信号がプリアンブルであると識別した場合には、プリアンブルパターン識別出力端子110を介して“1”を出力し、識別できない場合には“0”を出力する。
【００４５】
周波数偏差検出回路108の処理は、２成分の各々の周波数f1，f2と、プリアンブルの信号の送信ベースバンドでの周波数成分の各成分の周波数f10，f20（低い周波数の成分をf10、高い周波数の成分をf20とする）を用いて行う。２成分とも削られていない場合は両方の周波数を用いて式(19)により算出する。
【数１９】
【数２０】
【００４６】
シンボルタイミング検出回路109の処理は、プリアンブルの２成分とも削られていないと判断した場合にのみ行い、２成分の位相φ1とφ2を用いて行う。
変調方式がQPSKまたはQAMの場合、式(9)よりシンボルタイミングfbt0は、式(21)のようになる。
【数２１】
【数２２】
【００４７】
以上の説明による本発明の実施形態では、入力端子101にプリアンブル以外の信号が入力された場合、全電力pは、２つの周波数スペクトル成分の電力p1とp2以外に信号成分の電力を含んでいるため、電力比γの値は小さくなる。
また、入力端子101にプリアンブルが入力された場合は、振幅が変動したとしても、p1+p2は振幅変動の影響を受けたプリアンブル成分の電力で、全電力pは、振幅変動の影響を受けたプリアンブル成分の電力と雑音電力の和であり、振幅の落ち込みによりp1+p2が小さくなる場合は、pも小さくなる。このため、雑音の影響を受けるものの、振幅変動は直接影響しない。従って、振幅変動の影響を受けないプリアンブルパターン識別が可能である。
【００４８】
更に、送信機と受信機との周波数偏差により、プリアンブルの２つの周波数成分のうち片方の成分が、受信機のフィルタの帯域外へ外れたとしても、p1とp2の一方がプリアンブルの１成分の電力で、他が雑音電力であり、全電力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 (4)
- 受信したベースバンド信号の周波数スペクトルを求め、該周波数スペクトルから２つの周波数f1とf2の成分を抽出し、該抽出されたそれぞれの周波数f1とf2の成分のそれぞれの電力p1とp2と全周波数成分の電力pを算出し、該それぞれの電力p1とp2の和と全周波数成分の電力pとの電力比(p1+p2)/pを算出し、算出された該電力比(p1+p2)/pが所定のしきい値を超えていれば、受信したベースバンド信号がプリアンブルの信号であると判断することを特徴とするプリアンブルパターン識別方法。
- 請求項１記載のプリアンブルパターン識別方法において、
前記周波数スペクトルから抽出する２つの周波数f1とf2の成分は、振幅が高い順に２つ抽出することを特徴とするプリアンブルパターン識別方法。 - 請求項１または請求項２のいずれかに記載のプリアンブルパターン識別方法において、
前記電力比が所定のしきい値を超えている場合に、前記２つの周波数f1とf2のの差の絶対値|f2-f1|を求め、
該差の絶対値|f2-f1|が所定の第１の値の範囲内にあれば、受信したベースバンド信号がプリアンブルの信号であると判断することを特徴とするプリアンブルパターン識別方法。 - 受信信号を１シンボル当たりNov回（Novは２以上の正整数）オーバーサンプリングしたベースバンド信号から、任意の連続したNwinサンプル（Nwinは２以上の正整数）のベースバンド信号を抽出し、該抽出したベースバンド信号に任意の窓関数を乗じ、該任意の窓関数を乗じたベースバンド信号x(n)（n=0，1，‥‥‥，(Nwin-1)）のDFT（離散フーリエ変換：Discrete Fourier Transform）X(k) （k=0，1，‥‥‥、(Nwin-1)）を演算し（X(k)は複素数）、該DFT出力X(k) （k=0，1，‥‥‥，(Nwin-1)）の全電力p=|X(0)|2+|X(1)|2+…+|X(Nwin-1)|2を算出し、
該DFT出力X(k)のうち、振幅の自乗値|X(k)|2が最大となるX(k)を検索してそのkをk11とし、X(k11)に隣接する成分から、|X(k11)|2の次に|X(k)|2が大きくなるX(k)を検索してそのkをk12とし、X(k11)とX(k12)の電力和p1=|X(k11)|2+|X(k12)|2を算出し、X(k11)の周辺を補間し、該補間した部分のピークにおける周波数をf1、位相をφ1とし、
変調方式により定められる自然数の定数をkdとし、X(k) (k=(k11-kd-1) mod Nwin，(k11-kd) mod Nwin，(k11-kd+1) mod Nwin，(k11+kd-1) mod Nwin，(k11+kd) mod Nwin，(k11+kd+1) mod Nwin)から、|X(k)|2が最大となるX(k)を検索してそのkをk21とし、X(k21)に隣接する成分のうち、|X(k21)|2の次に|X(k)|2が大きくなるX(k)を検索してそのkをk22とし、X(k21)とX(k22)の電力和p2=|X(k21)|2+|X(k22)|2を算出し、X(k21)の周辺を補間し、該補間した部分のピークにおける周波数をf2、位相をφ2とし、
該算出したf1とf2を比較し、f1>f2であれば、f1，φ1，p1を改めてf2，φ2，p2とし、f2，φ2，p2を改めてf1，φ1，p1とすることを特徴とするプリアンブル成分抽出方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001023874A JP3933399B2 (ja) | 2001-01-31 | 2001-01-31 | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001023874A JP3933399B2 (ja) | 2001-01-31 | 2001-01-31 | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 |
Related Child Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007035202A Division JP4493669B2 (ja) | 2007-02-15 | 2007-02-15 | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 |
| JP2007035203A Division JP4540682B2 (ja) | 2007-02-15 | 2007-02-15 | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2002232503A JP2002232503A (ja) | 2002-08-16 |
| JP2002232503A5 JP2002232503A5 (ja) | 2005-09-02 |
| JP3933399B2 true JP3933399B2 (ja) | 2007-06-20 |
Family
ID=18889099
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001023874A Expired - Fee Related JP3933399B2 (ja) | 2001-01-31 | 2001-01-31 | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3933399B2 (ja) |
-
2001
- 2001-01-31 JP JP2001023874A patent/JP3933399B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002232503A (ja) | 2002-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9991930B2 (en) | Configurable correlator for joint timing and frequency synchronization and demodulation | |
| US5787123A (en) | Receiver for orthogonal frequency division multiplexed signals | |
| JP3074103B2 (ja) | Ｏｆｄｍ同期復調回路 | |
| JP2526931B2 (ja) | Ｐｓｋ信号復調装置 | |
| JPH0746218A (ja) | ディジタル復調装置 | |
| JP2002511711A (ja) | 多重搬送波復調システムにおいて精細な周波数同期を行うための方法および装置 | |
| EP2234277A1 (en) | Phase noise correction circuit, transmission device, reception device, radio device, radio communication system, and phase noise correction method | |
| GB2300790A (en) | Symbol clock recovery from received signal having wide frequency error or offset possibilities | |
| JPH10308717A (ja) | 受信装置および受信方法 | |
| US7336723B2 (en) | Systems and methods for high-efficiency transmission of information through narrowband channels | |
| EP1045561B1 (en) | Frequency correction in multicarrier receivers | |
| JP3544147B2 (ja) | Ｏｆｄｍ信号受信装置、ｏｆｄｍ信号通信システム及びその通信制御方法 | |
| JP4493669B2 (ja) | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 | |
| JP3933399B2 (ja) | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 | |
| JP4540682B2 (ja) | プリアンブルパターン識別方法及び周波数偏差検出方法並びにシンボルタイミング検出方法 | |
| EP1392030A2 (en) | Modulators and demodulators which can use a plurality of modulation techniques | |
| KR100429837B1 (ko) | 직교주파수 분할 다중화 신호의 동기화 방법 및 그 장치 | |
| JP3793198B2 (ja) | Ｏｆｄｍ信号通信システム及びｏｆｄｍ信号送信機 | |
| EP0785647A1 (en) | Method and apparatus for the synchronization of an FSK modulated signal in a receiver | |
| JP3852574B2 (ja) | 信号復調装置及び信号復調方法 | |
| JPH118659A (ja) | クロックタイミング再生方法および回路 | |
| JP2003218969A (ja) | 復調装置 | |
| Zicari et al. | A programmable carrier phase independent symbol timing recovery circuit for QPSK/OQPSK signals | |
| JP3735619B2 (ja) | プリアンブルパターン識別方法 | |
| EP1267534A1 (en) | Digital modulation system, radio communication system, radio communication device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050308 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050308 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20061129 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061219 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070215 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20070313 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070313 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100330 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110330 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120330 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120330 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130330 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130330 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140330 Year of fee payment: 7 |
|
| LAPS | Cancellation because of no payment of annual fees |