(57)【特許請求の範囲】
【請求項１】
受信信号をキャリア周波数の信号で周波数変換し、
前記周波数変換した前記受信信号について、長さＬの拡散系列を長さΔＬの部分系列に分割し、
各部分系列毎に逆拡散ベクトルを算出し、
各部分系列毎に算出された前記逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させて、２本の枝に分岐し、
前記拡散系列の最後まで前記２本の枝に分岐を続けて２分木を生成し、
前記２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された前記逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、
前記複数のパスのうち、前記合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出し、
前記キャリア周波数をF[Hz]とし、送信側での周波数安定度をAt[ppm]とし、受信側での周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、前記拡散系列の長さをL[chip]とし、前記逆拡散ベクトルの算出過程で許容される位相回転をα[rad]としたときに、
前記拡散系列を等間隔に分割する周波数ピッチΔfが次式であらわされ、
【数１】
0005745293.tif 000023
これにより前記受信信号の長さＬの拡散系列が、次式でｎ分割され、
【数２】
0005745293.tif 000024
前記部分系列の長さΔＬが、次式に示す値以下になることを特徴とするスペクトル拡散通信システムの受信方法。
【数３】
0005745293.tif 000025
【請求項２】
前記各部分系列毎に前記逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、前記逆拡散ベクトルの位相回転の累積による時間軸変化が、逆拡散ベクトルの算出に際し最小スライド間隔に達したときに、
次式で定義される時間ずれΔtを考慮して、逆拡散する前に入力する前記部分系列をずらすことを特徴とする請求項１記載のスペクトル拡散通信システムの受信方法。
【数４】
0005745293.tif 000026
【請求項３】
キャリア周波数を生成するキャリア周波数生成手段と、
前記キャリア周波数によって周波数変換された受信信号について、長さＬの拡散系列を長さΔＬの部分系列に分割する拡散系列分割手段と、
前記拡散系列分割手段が分割した各部分系列毎に、前記受信信号を逆拡散して逆拡散ベクトルを算出する逆拡散手段であるマッチドフィルタと、
前記逆拡散ベクトルの位相を正負に一定角度回転させて、２本の枝に分岐し、前記拡散系列の最後まで前記２本の枝に分岐を続ける２分木生成手段と、
前記２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された前記逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、前記複数のパスのうち、前記合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出する最大値検出手段と、
前記最大値検出手段の逆拡散データに基づいて、受信信号を復調する復調手段と、を備え、
キャリア周波数をF[Hz]とし、送信装置の周波数安定度をAt[ppm]とし、受信装置の周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、前記拡散系列の長さをL[chip]とし、前記逆拡散ベクトルの算出過程で許容される位相回転α[rad]としたときに、
前記拡散系列を等間隔に分割する周波数ピッチΔfが次式であらわされ、
【数５】
0005745293.tif 000027
これにより前記受信信号の長さＬの拡散系列が、次式でｎ分割され、
【数６】
0005745293.tif 000028
前記部分系列の長さΔＬが、次式に示す値以下になるように前記拡散系列分割手段を構成したことを特徴とするスペクトル拡散通信システムの受信装置。
【数７】
0005745293.tif 000029
【請求項４】
前記各部分系列毎に前記逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、前記逆拡散ベクトルの位相回転の累積による時間軸変化が、前記マッチドフィルタの最小スライド間隔に達したときに、次式で定義される時間ずれΔtを考慮して、前記マッチドフィルタに入力する前記部分系列をずらす時間ずれ補償手段を備えたことを特徴とする請求項３記載のスペクトル拡散通信システムの受信装置。
【数８】
0005745293.tif 000030
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、スペクトル拡散通信システムの受信方法および受信装置に関し、特に長距離無線通信に好適なスペクトル拡散通信システムの受信方法および受信装置に関するものである。
【背景技術】
【０００２】
ユビキタス情報化社会の到来とともに実用化の期待されているセンサネットワークでは、通信速度は、遅くても良いが通信距離は長くしたいとの要求がある。例えば、無線通信を応用した倉庫の在庫管理や農場管理および鳥獣監視などがある。
【０００３】
一方、センサに搭載する通信装置として有望視されている現在のRFID(Radio Frequency IDentification)の通信距離は、数十ｍしかない。
【０００４】
したがって、RFIDのメリットを生かしたまま通信距離を拡大させることができれば、ユビキタス情報化社会の実現は一層促進されると思われる。
【０００５】
シャノンの定理では、通信容量Cの理論限界は次式で規定される。
【０００６】
【数１】
0005745293.tif 000002
【０００７】
ここで、Wは通信路の帯域幅，Sは信号電力，Nは雑音電力である。いま、受信電力（信号電力S）が通信距離Dのα乗に比例して減衰すると仮定すれば、次式を導出できる。
【０００８】
【数２】
0005745293.tif 000003
【０００９】
ここで、aは定数である。したがって、通信距離Dは、通信容量Cを小さく（通信速度を低く）することで拡大できることが分かる。
【００１０】
しかし、実際の通信では、周波数ずれと時間ずれによる性能劣化によって通信距離が制限される。以下、周波数ずれと時間ずれの性能劣化について説明する。
【００１１】
（１）周波数ずれによる性能劣化
図１は、一般的なスペクトル拡散通信のシステム構成図である。スペクトル拡散通信システムは、送信装置１と受信装置２を備えて構成される。
【００１２】
送信装置１は、送信信号３を拡散符号４によりスペクトル拡散する拡散手段５と、キャリア周波数生成手段６と、前記キャリア周波数生成手段６で生成される信号と前記拡散手段５の出力信号を乗算して、前記出力信号を周波数変換するための乗算器７と、乗算器７の出力信号を伝送信号として無線で送信する送信用アンテナ８と、によって構成される。
【００１３】
受信装置２は、送信装置１からの伝送信号を受信する受信用アンテナ９と、キャリア周波数生成手段10と、前記キャリア周波数生成手段10で生成される信号と前記受信用アンテナ９の出力信号を乗算して、前記出力信号を周波数変換するための乗算器11と、当該乗算器11の出力信号を拡散符号12によりスペクトル逆拡散する逆拡散手段13と、によって構成される。そして、その逆拡散手段13の出力が受信信号14となる。
【００１４】
数２に示すシャノンの変形定理を帯域幅Wのスペクトル拡散通信に当てはめると、狭帯域化、すなわち拡散系列長を長くする（通信速度を低速にする）ことで、通信距離Dを拡大できることが分かる。
【００１５】
しかし、実際の通信システムでは、送信装置１と受信装置２との間には僅かなキャリア周波数のずれがあり、これが原因で逆拡散期間に受信信号の位相が回転する。その結果、正しい処理利得が得られなくなり、通信距離Dが制限される。
【００１６】
スペクトル拡散の正規化処理利得Gpは、逆拡散期間に発生する位相回転をθとするとき次式で与えられる。
【００１７】
【数３】
0005745293.tif 000004
【００１８】
図２にキャリア周波数を40kHz，チップ速度を1kc/s，送信装置１及び受信装置２の周波数安定度をそれぞれ10ppmとした場合において、拡散系列長Lと正規化処理利得Gpの関係を示す。この図では、拡散系列長Lが大きくなると、正規化処理利得Gpが１から急激に落ち込んでいることがわかる。また図３に、拡散系列長Lと処理利得Gp（正規化してない）の関係を示す。この図では、拡散系列長Lが大きくなると、処理利得Gpが頭打ちになり、拡散比が制限されることがわかる。図中の横軸から垂直に伸びた点線は、拡散比の制限が始まる位置を示している。これらの図から、拡散系列長Lを大きくすると十分な処理利得Gpが得られないことがわかる。
【００１９】
（２）時間ずれによる性能劣化
キャリア周波数以外にも、送信装置１と受信装置２との間には僅かなデータクロック周波数のずれが生じる。その結果、逆拡散期間に拡散系列のチップずれが起こり、同様に正しい処理利得が得られなくなり、通信距離Dが制限される。
【００２０】
図４にキャリア周波数を40kHz，チップ速度を1kc/s，送信装置１及び受信装置２のキャリア周波数の周波数安定度をそれぞれ10ppmとした場合において、送信装置１と受信装置２との間において、逆拡散期間に生ずるチップずれの影響による拡散系列長Lと処理利得の関係を示す。この図でも、拡散系列長Lが大きくなると、正規化処理利得Gpが１から急激に落ち込んでいることがわかる。
【００２１】
実際の通信では、周波数ずれと時間ずれの両方の影響を同時に受け、通信距離Dは、周波数ずれや時間ずれの性能劣化により制限を受ける。
【００２２】
このような問題に対し、従来技術では、受信側で、予め保持している期待信号と、受信信号搬送波との相関演算によって、受信信号搬送波の中に存在するトグル点の候補を検出し、その検出結果に基づいてシフトさせた拡散符号を受信信号に乗算することにより、受信信号を復調する方法が開示されている（特許文献１）。
【００２３】
また、遅延検波によって周波数ずれを伴う信号を復調する方法（特許文献２）や、搬送波周波数オフセット補正方法（特許文献３）などが開示されている。
【先行技術文献】
【特許文献】
【００２４】
【特許文献１】特許３６３９８３９号公報
【特許文献２】特開平９−９３２１７号公報
【特許文献３】特開２００７−２０２０８８号
【発明の概要】
【発明が解決しようとする課題】
【００２５】
しかしながら、上記特許文献１〜３は、いずれも送受信機のキャリア周波数ずれの影響を受けない範囲にスペクトル拡散系列長を設定しており、S／Nの良い受信信号から送信機のキャリア周波数を検出する方法等を採用している。そして、長いスペクトル拡散系列長を使用する必要のある場合は、例えば温度補償機能付発振器のような高精度の発振器を使用することで対応している。
【００２６】
したがって、送受信機間のキャリア周波数ずれの影響が大きい範囲にスペクトル拡散系列長を設定し、送信信号の狭帯域化すなわち通信速度を低速化することにより、送信電力を上げることなく無線長距離伝送するものではなかった。
【００２７】
例えば図５では、スペクトル拡散通信システムで送信する信号Sと雑音Nの関係について、ｎ個に分割された部分系列の送信信号Sと雑音Nがあった場合、部分系列の送信信号Sを加算し、その振幅をｎ倍にすれば、信号電力はｎ２倍、雑音電力はｎ倍となることを示している。ここで、送信信号Sの出力を増加させれば、S／Nは良くなるが、高出力電波は電波法の規制等により上限が決まっている。
【００２８】
図６は、スペクトル拡散通信システムで受信するｎ個に分割された部分系列の信号Sと雑音Nの関係を示しているが、受信側では送信側との間にキャリア周波数ずれ（各部分系列の逆拡散ベクトルの位相のずれ）があるため、部分系列の受信信号Sを加算したとしてもS／Nが悪く、受信側では、S／Nが悪い条件での復調が必要となる。
【００２９】
そこで、本発明では、信号の狭帯域化の制限を無くし、無線長距離通信が可能なスペクトル拡散通信システムの受信方法および受信装置を提供することを目的とする。
【課題を解決するための手段】
【００３０】
本発明の請求項１に係る発明は、受信信号をキャリア周波数の信号で周波数変換し、前記周波数変換した前記受信信号について、長さＬの拡散系列を長さΔＬの部分系列に分割し、各部分系列毎に逆拡散ベクトルを算出し、各部分系列毎に算出された前記逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させて、２本の枝に分岐し、前記拡散系列の最後まで前記２本の枝に分岐を続けて２分木を生成し、前記２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された前記逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、前記複数のパスのうち、前記合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出し、
前記キャリア周波数をF[Hz]とし、送信側でのキャリア周波数の周波数安定度をAt[ppm]とし、受信側でのキャリア周波数の周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、前記拡散系列の長さをL[chip]とし、前記逆拡散ベクトルの算出過程で許容される位相回転をα[rad]としたときに、
前記拡散系列を等間隔に分割する周波数ピッチΔfが次式であらわされ、
【００３１】
【数４】
0005745293.tif 000005
【００３２】
これにより前記受信信号の長さＬの拡散系列が、次式でｎ分割され、
【００３３】
【数５】
0005745293.tif 000006
【００３４】
前記部分系列の長さΔＬが、次式に示す値以下になる
【００３５】
【数６】
0005745293.tif 000007
【００３６】
ことを特徴とする。
【００３７】
本発明の請求項２に係る発明は、前記各部分系列毎に前記逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、前記逆拡散ベクトルの位相回転の累積による時間軸変化が、逆拡散ベクトルの算出に際し最小スライド間隔に達したときに、次式で定義される時間ずれΔtを考慮して、逆拡散する前に入力する前記部分系列をずらす
【００３８】
【数７】
0005745293.tif 000008
【００３９】
ことを特徴とする。
【００４０】
本発明の請求項３に係る発明は、キャリア周波数を生成するキャリア周波数生成手段と、前記キャリア周波数によって周波数変換された受信信号について、長さＬの拡散系列を長さΔＬの部分系列に分割する拡散系列分割手段と、前記拡散系列分割手段が分割した各部分系列毎に、前記受信信号を逆拡散して逆拡散ベクトルを算出する逆拡散手段であるマッチドフィルタと、前記逆拡散ベクトルの位相を正負に一定角度回転させて、２本の枝に分岐し、前記拡散系列の最後まで前記２本の枝に分岐を続ける２分木生成手段と、前記２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された前記逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、前記複数のパスのうち、前記合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出する最大値検出手段と、前記最大値検出手段の逆拡散データに基づいて、受信信号を復調する復調手段と、を備え、
前記受信信号の長さＬの拡散系列を、キャリア周波数をF[Hz]とし、送信装置の周波数安定度をAt[ppm]とし、受信装置の周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、前記拡散系列の長さをL[chip]とし、前記逆拡散ベクトルの算出過程で許容される位相回転α[rad]としたときに、前記拡散系列を等間隔に分割する周波数ピッチΔfが次式であらわされ、
【００４１】
【数８】
0005745293.tif 000009
【００４２】
これにより前記受信信号の長さＬの拡散系列が、次式でｎ分割され、
【００４３】
【数９】
0005745293.tif 000010
【００４４】
前記部分系列の長さΔＬが、次式に示す値以下になるように前記拡散系列分割手段を構成したことを特徴とする。
【００４５】
【数１０】
0005745293.tif 000011
【００４６】
本発明の請求項４に係る発明は、前記各部分系列毎に前記逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、前記逆拡散ベクトルの位相回転の累積による時間軸変化が、前記マッチドフィルタの最小スライド間隔に達したときに、次式で定義される時間ずれΔtを考慮して、前記マッチドフィルタに入力する前記部分系列をずらす時間ずれ補償手段を備えたことを特徴とする。
【００４７】
【数１１】
0005745293.tif 000012
【発明の効果】
【００４８】
本発明の請求項１によれば、２分木の最初から最後までの複数のパスのうち、合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出することにより、結果的に各部分系列毎に算出された逆拡散ベクトルの位相ずれ、すなわち周波数ずれが最小となる信号を検出することとなり、キャリア周波数ずれを最小化した受信信号の復調が可能となるため、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【００４９】
本発明の請求項２によれば、各部分系列毎に逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、逆拡散ベクトルの位相回転の累積による時間軸変化が逆拡散ベクトルの算出に際し最小スライド間隔に達したとき、時間ずれΔtを考慮して、逆拡散する前に入力する部分系列をずらすことにより、時間ずれを補償することができ、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【００５０】
本発明の請求項３によれば、受信信号の拡散系列を分割して部分系列毎に逆拡散ベクトルを算出し、前記逆拡散ベクトルの位相を正負に一定角度回転させ、２本の枝に分岐し、前記拡散系列の最後まで前記２本の分岐を続ける２分木生成手段と、前記２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された前記逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、前記複数のパスのうち、前記合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出する最大値検出手段と、を備えたことにより、キャリア周波数ずれを最小化した受信信号の復調が可能となるため、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【００５１】
本発明の請求項４によれば、各部分系列毎に逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、逆拡散ベクトルの位相回転の累積による時間軸変化がマッチドフィルタの最小スライド間隔に達したときに、時間ずれΔt分だけ、マッチドフィルタの入力信号をずらす時間ずれ補償手段を備えたことにより、時間ずれを補償することができ、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【図面の簡単な説明】
【００５２】
【図１】従来のスペクトル拡散通信のシステム構成図である。
【図２】拡散系列長Lと正規化処理利得Gpの関係を示す図である。
【図３】拡散系列長Lと処理利得Gpの関係を示す図である。
【図４】送受信機間の逆拡散期間に生ずるチップずれと処理利得の関係を示す図である。
【図５】スペクトル拡散通信システムの送信信号と雑音の関係を示す図である。
【図６】スペクトル拡散通信システムの受信信号と雑音の関係を示す図である。
【図７】比較例の周波数軸アプローチの概略構成図である。
【図８】本発明の一実施例における時間軸アプローチの概略構成図である。
【図９】同上、スペクトル拡散系列の分割を示す概念図である。
【図１０】部分系列毎の位相回転がない場合の合成逆拡散ベクトルを示す図である。
【図１１】部分系列毎の位相回転がある場合の合成逆拡散ベクトルを示す図である。
【図１２】同上、ｍ個の部分系列の合成逆拡散ベクトルを示す図である。
【図１３】本発明の実施例の２分木の構成図である。
【図１４】同上、時間軸アプローチによる二分木生成による探索動作を説明する概念図である。
【図１５】同上、二分木生成による探索結果から合成逆拡散ベクトルのノルムが最大となるパスを逆拡散データとして検出する動作を説明する概念図である。
【図１６】比較例の周波数軸アプローチ特性の計算結果を示す図である。
【図１７】本発明の実施例の時間軸アプローチ特性の計算結果を示す図である。
【図１８】比較例の周波数軸アプローチ特性のシミュレーション結果を示す図である。
【図１９】本発明の実施例の時間軸アプローチ特性のシミュレーション結果を示す図である。
【図２０】有線実験（α＝π/2）において比較例の周波数軸アプローチと本発明の実施例の時間軸アプローチを適用した場合のS/N特性を示す図である。
【図２１】有線実験（α＝π/4）において比較例の周波数軸アプローチと本発明の実施例の時間軸アプローチを適用した場合のS/N特性を示す図である。
【図２２】比較例の周波数軸アプローチと本発明の実施例の時間軸アプローチの演算量の比較を示す図である。
【発明を実施するための形態】
【００５３】
以下、図面を参照して本発明の実施形態について詳細に説明する。なお、本発明は、時間軸アプローチによって、送信装置に対する受信装置のキャリア周波数ずれを補償するものであるが、わかりやすく説明するために、従来技術である周波数軸アプローチによるキャリア周波数ずれの補償方法を、比較例として最初に説明する。
【００５４】
１．周波数軸および時間軸アプローチの動作原理
（１）周波数軸アプローチ
まず、比較例として、キャリア周波数ずれを周波数軸から見てずれ補償する周波数軸アプローチについて説明する。
【００５５】
図７は、周波数軸アプローチを実現する受信装置20の概略構成図である。具体的には、無線で伝送されてきた信号を受信する受信用アンテナ21と、受信側のキャリア周波数Fの上限fHから下限fLまでのずれ周波数を周波数ピッチΔf間隔にｎ分割した周波数に対応したｎ個のキャリア周波数生成手段22と、それらのｎ個のキャリア周波数生成手段21で生成される信号と受信用アンテナ21で受信した受信信号を乗算して、当該受信信号を周波数変換するための乗算器23と、乗算器23の各出力信号を逆拡散する逆拡散手段であるマッチドフィルタ24と、逆拡散された各信号を比較し、キャリア周波数ずれの少ない最大の信号を受信データとして選択する比較器25と、その選択された信号を復調する復調手段26と、によって構成される。
【００５６】
図示しない送信装置と受信装置20との間の周波数差は、キャリア周波数生成手段２に組み込まれた発振器（図示せず）の製造バラツキや温度等の使用環境が原因で変動する。しかし、逆拡散期間は十分短く、ここでは逆拡散期間内の周波数変動はないと仮定する。
【００５７】
周波数軸アプローチでは、キャリア周波数ずれを周波数軸から補償する。いま、キャリア周波数をF[Hz]，送信装置のキャリア周波数の周波数安定度をAt[ppm]とし、受信装置のキャリア周波数の周波数安定度をAr[ppm]，チップ速度をRC[cps]，拡散系列長をL[chip]，逆拡散過程で許容される位相回転をα[rad]とすると、キャリア周波数Fの上限fHおよび下限fLは次式で表わされる。
【００５８】
【数１２】
0005745293.tif 000013
【００５９】
【数１３】
0005745293.tif 000014
【００６０】
図７に示す受信装置20は、上限周波数fHから下限周波数fLまでのずれ周波数を周波数ピッチΔf間隔に分割し、分割した各周波数で受信信号のキャリア復調及び逆拡散をｎ個のキャリア周波数生成手段22，乗算器23，マッチドフィルタ24で行い、各マッチドフィルタ24からの出力値を比較器25で比較することにより、その最大値を受信データとして検出する。
【００６１】
周波数軸アプローチの正規化処理利得Gpsは、逆拡散過程で位相がα回転する場合、数３のθにαを代入することで与えられる。
【００６２】
つぎに、周波数軸アプローチの演算量を求める。
【００６３】
拡散系列長Lは、前記周波数ピッチΔfを用いると、
【００６４】
【数１４】
0005745293.tif 000015
【００６５】
で与えられる。これから周波数ピッチΔfは、
【００６６】
【数１５】
0005745293.tif 000016
【００６７】
となる。帯域2(At+Ar)Fのずれ周波数は周波数ピッチΔfによりｎ分割される。
【００６８】
【数１６】
0005745293.tif 000017
【００６９】
したがって、逆拡散にマッチドフィルタ24を使用するとき、周波数軸アプローチの演算量NSは次式で与えられる。
【００７０】
【数１７】
0005745293.tif 000018
【００７１】
周波数軸アプローチの演算量は、Lの3乗に比例して増加する。 
【００７２】
（２）時間軸アプローチ
つぎに、本発明の実施例として、キャリア周波数ずれを時間軸から見てずれ補償する時間軸アプローチについて説明する。
【００７３】
図８は、本発明の受信方法（時間軸アプローチ）を実現する受信装置２の概略構成図である。具体的には、無線で伝送されてきた信号を受信する受信用アンテナ31と、受信側で所定のキャリア周波数の信号を生成するキャリア周波数生成手段32と、受信された受信信号とキャリア周波数生成手段32の出力信号を乗算して、当該受信信号を周波数変換する乗算器33と、乗算器33の出力を逆拡散する逆拡散手段であるマッチドフィルタ34と、送信装置（図示しない）と受信装置２のキャリア周波数のずれ範囲を探索する二分木生成手段35と、最大値検出手段36を備えた探索手段37と、前記探索手段37の探索結果に基づいて、受信信号を復調する復調手段38と、を備えて構成される。
【００７４】
上記受信装置２が行なう時間軸アプローチの動作手順は、以下の通りである。
【００７５】
（Ｉ）拡散系列を部分系列に分割
長さLの拡散系列を、位相回転がα[rad]となる長さΔLの部分系列に分割する。ここで、αは処理利得Gpの劣化が許容できる位相回転である。数３で与えられる正規化処理利得Gpは、逆拡散過程に発生する位相回転θが小さいとき、ほとんど劣化しない。例えば、θがπ/4のときのGpは約0.974である。時間軸アプローチはこの特徴を利用する。図９に拡散系列の分割を示すが、これは乗算器33の前段または後段に設けた拡散系列分割手段（図示せず）で行なわれる。
【００７６】
（ＩＩ）部分系列の逆拡散
受信用アンテナ31で受信信号に対して、前記拡散系列分割手段が分割した各部分系列ごとにI,Q直交軸で逆拡散を行う。ここでいうI軸は直交変調による波形の同相成分を意味し、Q軸は直交位相成分を意味する。マッチドフィルタ34による逆拡散結果は、（逆拡散ベクトル）として扱う。
【００７７】
例えば、図１０は、各部分系列毎に算出されたｍ個の逆拡散ベクトルΔL１、ΔL２・・・ΔLｍに位相回転が無い場合において、その逆拡散ベクトルΔL１、ΔL２・・・ΔLｍを合成した合成逆拡散ベクトルを示している。このような理想的な場合、正規化処理利得GpはロスがないLとなりI軸と同軸上になる。一方、図１１のように逆拡散ベクトルΔL１、ΔL２・・・ΔLｍにそれぞれθずつ位相回転がある場合、ｍ個分では、図１２のように、位相のずれが累積していき、合成逆拡散ベクトルは小さくなる。例えば、ｍが20で、θがπ/12の場合、正規化処理利得Gpは0.19Lとなる。なお、これら図１０〜図１２は逆拡散ベクトルについて説明するための図であり、本発明を適用したものではない。
【００７８】
（ＩＩＩ）二分木の構成
マッチドフィルタ34で得られた部分系列の逆拡散ベクトルは、次の二分木生成手段35で位相を±α回転した２本の枝に分岐される。二分木生成手段35は、この操作を拡散系列の最後まで続け二分木を完成させる。
【００７９】
図１３に二分木の構成を示す。なお、この図１３は、位相回転の操作を示すものであるため、実際の動作上のキャリア周波数ずれに伴う位相ずれは、ここでは含めていない。したがって、例えばΔL1から＋α位相が回転（ΔL2）した後、さらに−α位相が回転した場合（ΔL3）、位相回転の累積はゼロとなる。また、ΔL3の部分系列の逆拡散ベクトルの先端では、2本の枝が位相回転がゼロとなる点で交差しており、ΔL4の部分系列の逆拡散ベクトルの先端では、位相回転が＋α、−αの点で2本の枝が交差する。そのため、実際の二分木生成手段35による探索では、その交点以後で重複する計算を省略して、一本の枝として計算することができる。
【００８０】
（ＩＶ）ベクトル和の最大値の検出
最大値検出手段36は、二分木の初期状態から最終枝までのn個の部分系列からなるパスについて、パス毎に合成逆拡散ベクトルのノルムを求める。その後、ノルムが最大となるパスを復調すべき逆拡散データとして検出し、復調手段38に送出する。図１４は、時間軸アプローチによる二分木生成手段35による探索動作を説明する概念図である。この図は、部分系列毎に位相ずれを伴いながら二分木探索する様子を示している。ここで、図１３と図１４を比較すると、図１３では、ΔL4の部分系列の逆拡散ベクトルが６本記載されているのに対し、図１４では、８本記載されている。しかし、位相ずれを伴う図１４においても、８本のうち、２本の逆拡散ベクトルが、他の２本の逆拡散ベクトルと大きさも向きも同じベクトルであるために計算が省略され、実際の計算は６本分である。これは、図１３において交点以後で重複する計算を省略することと同様に説明できるが、より詳細な説明をすると、各部分系列ΔL中において、累積された位相の回転角度が同じ逆拡散ベクトルどうしは、その算出が１回の計算で済む、という意味である。
【００８１】
図１５は、最大値検出手段36によって、部分系列毎の逆拡散ベクトルを合成した逆拡散ベクトルのノルムが最大となるパスを検出する動作を示す概念図である。この図のように、合成逆拡散ベクトルのノルムが最大となるパスを選択し、その正規化処理利得をGpとすると、結果的に部分系列毎の逆拡散ベクトルの位相ずれを最小化するようにパスを選択したこととなり、キャリア周波数のずれも最小化する。
【００８２】
時間軸アプローチの正規化処理利得Gptは、部分系列ごとに位相を＋α，−α交互に回転する場合に得られ、数３で求めたのと同様に次式で与えられる。
【００８３】
【数１８】
0005745293.tif 000019
【００８４】
例えば、αがπ/４の場合には、Gptは約0.9Lとなる。
【００８５】
なお、最大値検出手段36は、図示しない記憶手段であるメモリ40に各部分系列毎の逆拡散ベクトルやパス毎の合成逆拡散ベクトル等のデータを保存したり、読み込んだりして、パス毎に合成逆拡散ベクトルのノルムを求めたり、ノルムが最大となるパスを復調すべき逆拡散データとして検出したりする。
【００８６】
つぎに、図８の受信装置２で行なわれる時間軸アプローチの演算量を求める。
【００８７】
各部分系列の長さΔLは、拡散系列長Lを数９で規定するｎ個に分割したものであり、次式であらわせる。なお、上記（１）周波数軸アプローチの場合と同様に、送信装置でのキャリア周波数の周波数安定度をAt[ppm]とし、受信装置でのキャリア周波数の周波数安定度をAr[ppm]とする。
【００８８】
【数１９】
0005745293.tif 000020
【００８９】
時間軸アプローチの演算量NTは、各部分相関のマッチドフィルタ演算，二分木の構成，ベクトルノルムの最大値の検出から求まる。ここで、各部分相関のマッチドフィルタの演算量の総和は、拡散系列長Lのマッチドフィルタ34の演算量に等しく、残りの演算はこれに比べると十分小さい。したがって、時間軸アプローチの演算量NTは、次式で近似できる。
【００９０】
【数２０】
0005745293.tif 000021
【００９１】
時間軸アプローチの演算量NTは、拡散系列長Lの2乗となり、拡散系列長Lの3乗に比例する周波数軸アプローチに比べ、管理センターとなる受信装置２側の信号処理負荷は小さい。つまり、図７の構成ではなく、図８の構成を採用することで、受信用アンテナ31が受信信号を受信してから、キャリア周波数のずれを最小化するような逆拡散データ（受信データ）を出力するまでの処理速度を、拡散系列長Lの3乗の比例から2乗に減らすことができる。
【００９２】
つぎに、前述した時間ずれの補償について検討する。時間ずれを補償する手段として、図７に示すマッチドフィルタ24の前段や、図８に示すマッチドフィルタ34の前段には、それぞれ時間ずれ補償手段29、39が配設される。
【００９３】
送信装置１や受信装置２において、キャリア周波数とデータクロック周波数は、それぞれ別個の発振器で生成する場合もあるが、ここでは同じ発振器として、上記キャリア周波数生成手段22，32で生成する場合を考える。
【００９４】
このとき、双方のデータクロック周波数のずれがΔfであれば、送信装置１と受信装置２との間の拡散系列の時間軸は、送信キャリア周波数をfS，受信キャリア周波数fRとするとき、(fR−Δf)/fR倍に伸縮する。したがって、図７に示す周波数軸アプローチを実現する受信装置２では、時間ずれ補償手段29により周波数のずれに応じて時間軸を伸縮させたのち逆拡散することで、周波数ずれ及び時間ずれの補償が可能となる。
【００９５】
一方、図８に示す時間軸アプローチを実現する受信装置２では、部分系列毎にする位相回転αに伴う時間ずれΔtが、次式のようにあらわせる。
【００９６】
【数２１】
0005745293.tif 000022
【００９７】
したがって、時間軸アプローチでは、位相回転の累積による時間軸変化がマッチドフィルタ34の最小スライド間隔に達したとき、上記数１１の時間ずれΔt分だけ、部分系列のマッチドフィルタ34への入力をずらせば、周波数ずれ及び時間ずれの補償が可能となる。
【００９８】
２．周波数軸および時間軸アプローチの性能比較
（１）周波数軸及び時間軸アプローチ特性
キャリア周波数を40kHz，チップ速度を1kc/s，送信装置１及び受信装置２のキャリア周波数の周波数安定度をそれぞれ10ppmとした場合の、周波数軸アプローチ特性の計算結果を、α＝π/8，π/4，π/2の３つの場合について図１６に示す。また、時間軸アプローチ特性の計算結果を同様に図１７に示す。また、周波数軸アプローチ特性のシミュレーション結果を図１８に示し、時間軸アプローチ特性のシミュレーション結果を図１９に示す。周波数軸アプローチでは、前記図７に示す受信装置２を用い、時間軸アプローチでは、前記図８に示す受信装置２を用いている。
【００９９】
これらの図１６〜図１９から、周波数軸アプローチ及び時間軸アプローチともに、拡散系列長を大きくしたときの性能劣化はなく、キャリア周波数ずれ及び時間ずれを補償出来ていることが分かる。また、性能の劣化量は許容位相回転αにより調整できる。αを大きくすれば、許容誤差も大きくなるため、劣化量が大きくなるが、αを小さくすれば、劣化量が小さくなるため、正規化処理利得が１に近づく。
【０１００】
比較例および本発明の実施例は、無線に限らず、有線でも適用できるため、有線に適用した場合のシミュレーション結果を図２０、図２１に示す。
【０１０１】
図２０及び図２１は、αをπ/2、π/4とし、S/Nが良い場合（60dB）と悪い場合（-30dB）について、シミュレーション実験した結果であり、周波数軸アプローチ（比較例）、時間軸アプローチ（本発明の実施例）、理論計算値の特性をそれぞれ示している。いずれの結果も、拡散系列長Lが103以上になると、理論計算値と同様の結果となる。なお、S/Nが悪い場合（-30dB）においては、拡散系列長Lが小さいと、理論計算値とずれがあるが、これは実験における雑音の発生状況によるものと考えられる。
【０１０２】
（２）演算量の比較
図２２に、それぞれ数１７と数２０に示す周波数軸アプローチと時間軸アプローチの演算量を示す。
【０１０３】
拡散系列長Lを大きくするとき、周波数軸アプローチ式はLの3乗に比例するのに対し、時間軸アプローチはLの2乗となる。また、時間軸アプローチの演算量は、許容位相回転αに依存しない。
【０１０４】
したがって、拡散系列長を大きくした長距離通信では、時間軸アプローチに必要な演算量は周波数軸アプローチの数桁分の１になり、管理センター側の信号処理負荷は小さく済むことが分かる。
【０１０５】
（３）まとめ
シャノンの定理から通信距離と通信速度の関係を導出し、通信速度を低くすることで通信距離が拡大できることを示した。しかし、この変形シャノンの定理をスペクトル拡散通信に応用すると、送受信機間に生ずる周波数ずれと時間ずれが問題となることが分かった。
【０１０６】
そこで、スペクトル拡散を応用した長距離無線通信として、周波数軸アプローチと時間軸アプローチによる。周波数ずれと時間ずれの補償方法について検討した。これら補償方法は、拡散系列を長くしても性能劣化はほとんどなく、劣化量は調整できることが分かった。
【０１０７】
また、ずれ補償に必要な演算量は、周波数軸アプローチは拡散系列長の3乗に比例するのに対し、時間軸アプローチは拡散系列長の2乗で済むことが分かった。技術革新の著しいCPU性能の向上を考慮すると、特に時間軸アプローチは十分に実用可能と思われる。また、本発明の受信方法は、センサに搭載する無線機をそのまま使用できるため有利である。
【０１０８】
以上のように、本実施例の受信方法では、受信信号をキャリア周波数の信号で周波数変換し、周波数変換した前記受信信号について、長さＬの拡散系列を長さΔＬの部分系列に分割し、各部分系列毎に逆拡散ベクトルを算出し、各部分系列毎に算出された逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させて、２本の枝に分岐し、拡散系列の最後まで前記２本の枝に分岐を続けて２分木を生成し、２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、複数のパスのうち、合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出している。
【０１０９】
また、キャリア周波数をF[Hz]とし、送信側での周波数安定度をAt[ppm]とし、受信側での周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、拡散系列の長さをL[chip]とし、逆拡散ベクトルの算出過程で許容される位相回転をα[rad]としたときに、拡散系列を等間隔に分割する周波数ピッチΔfが上記数４であらわされ、これにより受信信号の長さＬの拡散系列が、上記数５でｎ分割され、部分系列の長さΔＬが、上記数６に示されるように分割している。
【０１１０】
この場合、２分木の最初から最後までの複数のパスのうち、合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出することにより、結果的に各部分系列毎に算出された逆拡散ベクトルの位相ずれ、すなわち周波数ずれが最小となる信号を検出することとなり、キャリア周波数ずれを最小化した受信信号の復調が可能となるため、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【０１１１】
また、本実施例の受信方法では、キャリア周波数をF[Hz]とし、送信側での周波数安定度をAt[ppm]とし、受信側での周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、拡散系列長をL[chip]とし、逆拡散過程で許容される位相回転をα[rad]としたときに、各部分系列毎に逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、逆拡散ベクトルの位相回転の累積による時間軸変化が、逆拡散ベクトルの算出に際し最小スライド間隔に達したときに、上記数２１で定義される時間ずれΔtを考慮して、逆拡散する前に入力する前記部分系列をずらしている。
【０１１２】
この場合、各部分系列毎に逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、逆拡散ベクトルの位相回転の累積による時間軸変化が逆拡散ベクトルの算出に際し最小スライド間隔に達したとき、時間ずれΔtを考慮して、逆拡散する前に入力する部分系列をずらすことにより、時間ずれを補償することができ、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【０１１３】
また、本実施例の受信装置２では、キャリア周波数を生成するキャリア周波数生成手段32と、キャリア周波数によって周波数変換された受信信号を逆拡散する逆拡散手段であるマッチドフィルタ34と、受信信号の拡散系列を分割して部分系列毎に逆拡散ベクトルを算出し、逆拡散ベクトルの位相を正負に一定角度回転させて、２本の枝に分岐し、拡散系列の最後まで分岐を続ける２分木生成手段35と、２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、複数のパスのうち、合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出する最大値検出手段36と、最大値検出手段36の逆拡散データに基づいて、受信信号を復調する復調手段38と、を備えて構成されている。
【０１１４】
また、前記受信信号の長さＬの拡散系列を、キャリア周波数をF[Hz]とし、送信装置１の周波数安定度をAt[ppm]とし、受信装置２の周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、拡散系列長をL[chip]とし、逆拡散過程で許容される位相回転をα[rad]としたときに、拡散系列を等間隔に分割する周波数ピッチΔfが上記数４であらわされ、これにより受信信号の長さLの拡散系列が、上記数５でｎ分割され、前記部分系列の長さΔＬが、上記数６に示される長さとする拡散系列分割手段を備えて構成されている。
【０１１５】
この場合、受信信号の拡散系列を分割して部分系列毎に逆拡散ベクトルを算出し、前記逆拡散ベクトルの位相を正負に一定角度回転させ、２本の枝に分岐し、前記拡散系列の最後まで前記２本の分岐を続ける２分木生成手段35と、前記２分木の最初から最後までの複数のパスについて、各部分系列毎に算出された前記逆拡散ベクトルを合成した合成逆拡散ベクトルのノルムをそれぞれ算出し、前記複数のパスのうち、前記合成逆拡散ベクトルのノルムが最大となるパスを復調すべき逆拡散データとして検出する最大値検出手段36と、を備えたことにより、キャリア周波数ずれを最小化した受信信号の復調が可能となるため、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【０１１６】
また、本実施例の受信装置２は、キャリア周波数をF[Hz]とし、送信装置１の周波数安定度をAt[ppm]とし、受信装置２の周波数安定度をAr[ppm]とし、チップ速度をRC[cps]とし、拡散系列長をL[chip]とし、逆拡散過程で許容される位相回転をα[rad]としたときに、各部分系列毎に逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、前記逆拡散ベクトルの位相回転の累積による時間軸変化が前記マッチドフィルタ34の最小スライド間隔に達したときに、上記数７で定義される時間ずれΔtを考慮して、前記マッチドフィルタ34に入力する部分系列をずらす時間ずれ補償手段39を備えて構成されている。
【０１１７】
この場合、各部分系列毎に逆拡散ベクトルの位相をそれぞれ正負に一定角度回転させるのに伴い、逆拡散ベクトルの位相回転の累積による時間軸変化がマッチドフィルタ34の最小スライド間隔に達したときに、時間ずれΔt分だけ、マッチドフィルタ34の入力信号をずらす時間ずれ補償手段39を備えたことにより、時間ずれを補償することができ、信号の狭帯域化の制限を無くし、無線長距離通信が可能となる。
【０１１８】
なお、本発明は上記実施例に限定されるものではなく、本発明の趣旨を逸脱しない範囲で変更可能である。例えば、上記実施例の二分木生成による探索では、すべての部分系列の分木させた枝に対して逆拡散ベクトル演算している。つまり、すべての枝について、最後の深さまで演算している。しかし、二分木生成の途中過程で、明らかに位相ずれが大きい枝については演算を打ち切って、演算量を減らしても良い。
【符号の説明】
【０１１９】
２  受信装置
32  周波数変動手段
34  マッチドフィルタ
35  二分木生成手段
36  最大値検出手段
37  探索手段
38  復調手段
39  時間ずれ補償手段
【図１】
0005745293.tif 000031
【図２】
0005745293.tif 000032
【図３】
0005745293.tif 000033
【図４】
0005745293.tif 000034
【図５】
0005745293.tif 000035
【図６】
0005745293.tif 000036
【図７】
0005745293.tif 000037
【図８】
0005745293.tif 000038
【図９】
0005745293.tif 000039
【図１０】
0005745293.tif 000040
【図１１】
0005745293.tif 000041
【図１２】
0005745293.tif 000042
【図１３】
0005745293.tif 000043
【図１４】
0005745293.tif 000044
【図１５】
0005745293.tif 000045
【図１６】
0005745293.tif 000046
【図１７】
0005745293.tif 000047
【図１８】
0005745293.tif 000048
【図１９】
0005745293.tif 000049
【図２０】
0005745293.tif 000050
【図２１】
0005745293.tif 000051
【図２２】
0005745293.tif 000052