(57)【要約】【課題】  無線通信において干渉信号の情報がない状態でも、また干渉信号と希望信号の信号帯域幅が異なっていても干渉信号が抑圧された希望信号の抽出を可能とする。【解決手段】  時間相関抽出部１５で受信信号の時間相関が演算され、受信信号の時間相関波形から、希望信号の時間相関成分以外が、干渉信号の時間相関波形として抽出され、信号波形再生部１７でその干渉信号の時間相関波形が参照信号で割算されて干渉信号が出力され、判定切替部１８で信号波形再生部１７の出力信号が所定レベル以上になったか否かが判定され、判定レベルになっていなければ受信信号が参照信号とされ、所定レベル以上であれば再生部１７の出力信号が参照信号とされる。再生部１７から再生干渉信号が減算部２３で受信信号から差し引かれて干渉が抑圧された希望信号が得られる。  【選択図】  図２
2005151461.tif 000002
【特許請求の範囲】
【請求項１】
無線受信信号（以下単に受信信号という）及び抽出対象情報が入力され、上記受信信号の時間相関が演算され、抽出信号の時間相関波形を抽出する時間相関抽出部と、
上記時間相関波形と参照信号が入力され、その時間相関波形を参照信号で割算して上記抽出信号を出力する信号波形再生部と、
上記信号波形再生部の出力信号が入力され、その出力信号が上記抽出信号であるか否かを判定し、抽出信号でなければ上記受信信号を上記参照信号とし、抽出信号であれば上記出力信号を上記参照信号とする判定切替部と
を具備する信号抽出装置。
【請求項２】
請求項１記載の装置において、
上記抽出信号は干渉信号であり、
上記受信信号から上記信号再生部からの抽出信号を減算して希望信号として出力する減算部を備えることを特徴とする歪除去信号抽出装置。
【請求項３】
請求項２記載の装置において、
上記受信信号に希望信号が含まれていない状態を検出する希望信号有無検出部と、
希望信号が含まれていない受信信号の上記時間相関抽出部で得られる相関値に所定レベル以上のものがあれば上記抽出対象情報を干渉信号とし、上記減算部の出力を、所定レベル以上のものがなければ上記抽出対象情報を希望信号とし上記信号波形再生部の出力を希望信号として出力する干渉有無判定部と
を備える信号抽出装置。
【請求項４】
請求項２記載の装置において、
上記時間相関抽出部は希望信号の時間相関波形を抽出する希望抽出部と希望信号以外の信号の時間相関波形を抽出する干渉抽出部とを備え、
上記信号波形再生部として、上記希望抽出部の出力が入力される希望波形再生部と、上記干渉抽出部の出力が入力され、出力を上記減算部へ出力する干渉波形再生部とを備え、
上記判定切替部として、上記希望波形再生部に対し希望判定切替部が、上記干渉波形抽出部に対し干渉判定切替部が設けられ、
上記受信信号に希望信号が含まれていない状態を検出する希望信号有無検出部と、
上記希望信号が含まれていない受信信号の上記時間相関抽出部で得られる相関値に所定レベル以上のものがあれば上記抽出対象情報を干渉信号とし、上記減算部の出力を希望信号として出力し、所定レベル以上のものがなければ上記抽出対象情報を希望信号とし、上記希望波形再生部の出力を希望信号として出力する干渉有無判定部とを備える信号抽出装置。
【請求項５】
請求項１〜４のいずれかに記載の装置において、
上記時間相関抽出部は上記受信信号の周期自己相関値を計算する周期自己相関計算部と、
上記周期自己相関値から上記抽出対象情報に応じて抽出信号の成分を抽出する相関ピーク抽出部と、
上記相関ピーク抽出部から抽出された相関値を逆フーリエ変換して抽出信号の時間相関波形を出力する逆フーリエ変換部とを具備することを特徴とする信号抽出装置。
【請求項６】
請求項１〜４のいずれかに記載の装置において、
上記時間相関抽出部は上記受信信号の互いに異なる遅延時間についての時間相関値を計算してそれぞれ周期自己相関値として出力する複数時間相関計算部と、
上記各時間相関計算部よりの周期自己相関値がそれぞれ入力され、上記抽出対象情報に応じて上記抽出信号の成分を抽出して時間相関波形をそれぞれ出力する複数の帯域通過フィルタ部とを具備することを特徴とする信号抽出装置。
【請求項７】
請求項１〜４のいずれかに記載の装置において、
上記時間相関抽出部は上記受信信号の時間相関値を計算する時間相関計算部と、
上記時間相関計算部よりの時間相関値が入力され、上記抽出対象情報に応じて上記抽出信号の成分を時間相関波形として出力する帯域通過フィルタ部とを具備することを特徴とする信号抽出装置。
【請求項８】
受信信号の周期自己相関値を計算する周期自己相関計算部と、
通信要求信号を受信し、通信要求の受付けを表わす制御信号により、上記周期自己相関計算部より希望信号が含まれない受信信号の周期自己相関値（以下事前周期自己相関値という）を出力させる通信要求受付部と、
上記事前周期自己相関値が入力され、その事前周期自己相関値が相対的に小さい成分の近くに時間相関値の成分をもつ希望信号の変調方式、信号帯域幅又はパルスレート及び搬送波周波数を送信パラメータとして決定する送信パラメータ決定部と、
上記送信パラメータを送信する送信部と
を具備する受信側装置。
【請求項９】
請求項８に記載の装置において、
上記送信パラメータ決定部は上記事前周期自己相関値が所定値以下の遅延時間−周期周波数領域を使用可能領域として検出する使用可能領域検出部と、
使用可能変調方式と使用可能信号帯域幅又はパルスレートと使用可能搬送波周波数との組み合せの複数のパラメータ組のそれぞれについて、その組み合せをもつ信号の周期自己相関値が所定値以上となる遅延時間と周期周波数との組をそれぞれ求める周期自己相関出現領域計算部と、
上記遅延時間と周期周波数の組中の上記使用可能領域に属するものと対応するパラメータ組を抽出する使用可能パラメータ抽出部と、
上記通信要求信号の受信電力と上記使用可能領域中の周期自己相関値から受信品質を推定する品質推定部と、
上記受信品質と上記抽出された各パラメータ組とから使用可能な符号化率及び拡散率をそれぞれ抽出し、そのパラメータ組に組み込みパラメータ組候補とする符号化率・拡散率決定部と、
上記パラメータ組候補についてそれぞれのビットレートを計算するビットレート計算部と、
上記計算されたビットレート中の最大値を検出し、その最大値となるパラメータ組候補を上記送信パラメータとして選択する送信パラメータ選択部とを備えることを特徴とする受信側装置。
【請求項１０】
請求項９に記載の装置において、
上記送信パラメータ決定部は上記事前周期自己相関値が所定値以下の遅延時間−周期周波数領域を使用可能領域として検出する使用可能領域検出部と、
使用可能変調方式と使用可能信号帯域幅又はパルスレートと使用可能搬送波周波数との組み合せの複数のパラメータ組のそれぞれについて、その組み合せをもつ信号の周期自己相関値が所定値以上である遅延時間と周期周波数との組をそれぞれ求める周期自己相関出現領域計算部と、
上記遅延時間と周期周波数の組中の上記使用可能領域に属するものと対応するパラメータ組を抽出する使用可能パラメータ抽出部と、
上記通信要求信号の受信電力と上記使用可能領域中の周期自己相関値から受信品質を推定する品質推定部と、
上記受信品質と上記抽出された各パラメータ組とから使用可能な符号化率及び拡散率をそれぞれ抽出し、そのパラメータ組に組み込みパラメータ組候補とする符号化率・拡散率決定部と、
上記パラメータ組候補についてそれぞれのビットレートを計算するビットレート計算部と、
上記通信要求信号に含まれた要求品質情報を分離する要求品質分離部と、
上記計算されたビットレート中の上記要求品質情報が示す要求品質を満足する範囲で、最低のビットレートを得ることができるパラメータ組候補を上記送信パラメータとして選択する送信パラメータ選択部とを備えることを特徴とする受信側装置。
【請求項１１】
請求項８〜１０の装置において、
請求項１〜６のいずれかに記載した信号抽出装置が信号処理部として設けられ、その信号処理部における時間相関抽出部中の周期自己相関値を求める方法が上記周期自己相関計算部と兼用されていることを特徴とする受信側装置。
【請求項１２】
請求項１〜１１のいずれかに記載の装置としてコンピュータを機能させるためのプログラム。
【請求項１３】
送信側装置から通信要求信号を受信側装置へ送信し、
受信側装置は上記通信要求信号を受信すると、受信信号の周期自己相関値を計算し、
上記計算された周期自己相関値が小さい成分の近くに希望信号の周期自己相関値の成分が現われる、その希望信号の変調方式、信号帯域幅又はパルスレート、及び搬送波周波数を送信パラメータ組として決定し、この送信パラメータ組を含む送信許可通知を上記送信側装置へ送信し、
上記送信側装置は上記送信許可通知を受信すると、これに含まれる送信パラメータ組を分離し、
送信データにより変調され、上記送信パラメータ組の各パラメータを満す希望信号を生成して上記受信側装置へ送信することを特徴とする無線通信方法。
【請求項１４】
請求項１３記載の方法において、
上記周期自己相関値が所定値以下の遅延時間・周期周波数領域を使用可能領域として検出し、
上記通信要求信号の受信電力と上記使用可能領域の周期自己相関値から受信品質を推定し、
上記推定した受信品質に基づき、上記使用可能領域に、周期自己相関値の成分が現われる信号の変調方式、信号帯域幅、搬送波周波数、拡散率、符号化率のパラメータ組候補を求め、
上記各パラメータ組候補について送信データのビットレートを計算し、
これらビットレート中の最大を与えるパラメータ組候補を上記送信パラメータとすることを特徴とする無線通信方法。
【請求項１５】
請求項１３記載の方法において、
上記送信側装置は上記通信要求信号に要求品質情報を付加して送信し、
上記周期自己相関値が所定値以下の遅延時間・周期周波数領域を使用可能領域として検出し、
上記通信要求信号の受信電力と上記使用可能領域の周期自己相関値から受信品質を推定し、
上記推定した受信品質に基づき、上記使用可能領域に、周期自己相関値の成分が現われる信号の変調方式、信号帯域幅、搬送波周波数、拡散率、符号化率のパラメータ組候補を求め、
上記各パラメータ組候補について送信データのビットレートを計算し、
これらビットレート中で上記通信要求信号中の要求品質情報以上で最低のビットレートを与えるパラメータ組候補を上記送信パラメータとすることを特徴とする無線通信方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、移動通信方式、その他の各種無線通信方式に適用され、無線受信信号から抽出したい希望信号あるいは干渉信号などを抽出する信号抽出装置、特に歪を抑圧した信号として抽出する装置、及びそのプログラム、またその装置を用いた無線通信方法、その装置、及びそのプログラムに関する。
【背景技術】
【０００２】
無線通信システムにおいては、限られた周波数資源を効率的に使用するために、受信信号に対する干渉量を抑える技術が重要となる。従来の干渉量抑圧技術として、図１に示すようなレプリカ生成型干渉キャンセラがある。
受信信号ｚ（ｔ）は希望信号ｘ（ｔ）と干渉信号ｗ（ｔ）と雑音ｎ（ｔ）とが重畳されたものであり、通信路推定部１で希望信号ｘ（ｔ）及び干渉信号ｗ（ｔ）の各伝搬路特性を推定し、これらに希望信号ｘ（ｔ）及び干渉信号ｗ（ｔ）の取りうる全てのシンボル系列候補を希望レプリカ生成部２及び干渉レプリカ生成部３でそれぞれ乗算することにより、全てのシンボル系列候補に対して希望信号レプリカと干渉信号レプリカを生成し、これらを加算部４で加算して受信信号レプリカとし、減算部５でこの受信信号レプリカを受信信号ｚ（ｔ）から差し引き、その誤差信号を通信路推定部１に入力すると共に、自乗部６で自乗して最尤系列推定部７に入力する。最尤系列推定部７で各シンボル系列候補を生成して各レプリカ生成部２及び３に入力し、実際の受信信号ｚ（ｔ）と最も近い受信信号レプリカを与える希望信号ｘ（ｔ）及び干渉信号ｗ（ｔ）のシンボル系列候補を最尤系列推定処理により判定し、希望信号ｘ（ｔ）のシンボル系列候補を受信信号ｚ（ｔ）の判定結果として出力する。このように受信信号ｚ（ｔ）から適応的に干渉信号ｗ（ｔ）を除去することにより、異なる複数の信号が同一時刻に同一周波数を同一空間で使用することが可能となり、周波数利用効率を向上させることが可能となる（特許文献１  参照）。
【特許文献１】特開平１０−１７３５７３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
上述したレプリカ生成型干渉キャンセラは、希望信号及び干渉信号の各レプリカを生成することにより、実効的に干渉を除去することが可能であるが、レプリカを生成するときに希望信号だけでなく干渉信号の変調方式等の情報を把握している必要がある。このため、干渉信号の情報を把握している場合は高度な干渉除去効果を得ることが可能であるが、この情報を有さない場合は干渉除去効果を得ることができない。
さらに、従来の干渉キャンセラは希望信号と干渉信号の信号帯域幅が等しい場合を想定しているため、広帯域を要求する希望信号と、狭帯域を要求する信号を対象とする干渉信号が合成されて受信される場合、十分な干渉抑圧効果を得ることができない。このため、干渉キャンセラによる周波数利用効率改善効果を得るためには、広帯域信号を使用している周波数帯域では、同時に狭帯域信号を利用することができない。
【０００４】
この発明の目的は干渉信号の情報を有さない状態であっても、干渉信号と希望信号の信号帯域幅が異なっていても干渉信号が抑圧された希望信号の抽出を可能とする、または希望信号を含まない干渉信号の抽出を可能とする信号抽出装置、そのプログラム及び無線通信方法を提供することにある。
【課題を解決するための手段】
【０００５】
この発明の信号抽出装置によれば、時間相関抽出部で受信信号の時間相関が演算され、抽出対象情報が表わす抽出信号の時間相関波形が抽出され、信号波形再生部で上記時間相関波形が参照信号で割算されて上記抽出信号が出力され、判定切替部で上記信号波形再生部の出力信号が上記抽出信号であるか否かが判定され、抽出信号でなければ上記受信信号が上記参照信号とされ、抽出信号であればその抽出信号が上記参照信号とされる。
この発明の方法によれば送信側装置から通信要求信号を受信側装置へ送信し、受信側装置は上記通信要求信号を受信すると、受信信号の周期自己相関値を計算し、上記計算された周期自己相関値が小さい成分の近くに希望信号の周期自己相関値の成分が現われる、その希望信号の変調方式、信号帯域幅又はパルスビットレート、及び搬送波周波数を送信パラメータ組として決定し、この送信パラメータ組を含む送信許可通知を上記送信側装置へ送信し、上記送信側装置は上記送信許可通知を受信すると、これに含まれる送信パラメータ組を分離し、送信データにより変調され、上記送信パラメータ組の各パラメータを満す希望信号を生成して上記受信側装置へ送信する。
【発明の効果】
【０００６】
この発明の信号抽出装置では抽出信号の時間相関波形のみを利用するため、例えば抽出信号が希望信号の場合、干渉信号の情報を知らない場合でも、希望信号を抽出することができる。また希望信号と干渉信号との信号帯域幅が異なっていても希望信号を抽出することができる。抽出信号が干渉信号の場合は、その干渉信号を同様に抽出することができるため、この抽出した干渉信号を受信信号から差し引くことにより希望信号を抽出することができる。干渉信号の抽出は時間相関抽出部で演算した時間相関値から希望信号の成分以外を干渉信号の時間相関波形として抽出することにより得られる。
【０００７】
この発明の方法によれば、同一時刻に同一周波数を多数の送受信機で共用することが可能となり、周波数利用効率が改善される。
【発明を実施するための最良の形態】
【０００８】
まずこの発明による信号抽出装置を適用した無線受信機の実施形態を説明する。
［実施形態１］
原理
図２を参照して実施形態１の原理的構成を説明する。アンテナ１１により受信された無線受信信号は必要に応じて周波数変換器１２により中間周波数の受信信号に変換され、この受信信号は一般にはＡ／Ｄ変換器１３により周期的にサンプリングされ、その各サンプルがディジタル値に変換され、ディジタル信号系列の受信信号に変換される。前記中間周波数ｆｉは無線受信信号中のベースバンド信号周波数より高ければよいが、処理のし易さからなるべく低い周波数がよい。Ａ／Ｄ変換器１３でのサンプリング周波数ｆｓは２ｆｉ 以上とされる。以下前記ディジタル信号系列の受信信号を単に受信信号という。この受信信号がこの発明の信号抽出装置の入力端子１４に入力される。
【０００９】
入力端子１４よりの受信信号は時間相関抽出部１５へ供給される。時間相関抽出部１５では受信信号ｚ（ｎ）の時間相関値を演算し、抽出対象指定部１６からの抽出信号識別・情報に基づき、抽出信号の時間相関成分が抽出される。ｎは受信信号系列における各サンプル時刻、つまり離散時刻を表わす。抽出信号は受信信号ｚ（ｎ）中から抽出したい信号であって、一般には希望信号ｘ（ｎ）または干渉信号ｗ（ｎ）である。希望信号ｘ（ｎ）の時間相関波形はｘ（ｎ）ｘ*(ｎ−τ)であり、干渉信号ｗ（ｚ）の時間相関波形はｗ（ｎ）ｘ*(ｎ−τ)である。τは受信信号ｚ（ｎ）のサンプル周期（１／ｆｓ）の整数倍であり、「* 」はｘ（ｎ−τ）の複素共役を表わす。時間相関抽出部１５の具体的構成は後で説明する。
【００１０】
時間相関抽出部１５で抽出された時間相関波形は信号波形再生部１７へ入力され、参照信号ｕ（ｎ−τ）により割算され抽出信号が信号波形再生部１７から出力される。抽出信号が希望信号ｘ（ｎ）の場合、信号波形再生部１７で下記の計算が行われる。
ｘ＾（ｎ）＝ｘ（ｎ）ｘ*(ｎ−τ)／ｕ*(ｎ−τ)                （１）
参照信号として過去の再生信号波形ｘ（ｎ−τ）を用いることができれば出力される再生信号波形は次のように希望信号波形が再現されたものになり、つまり受信信号ｚ（ｎ）から抽出したい信号が抽出される。
【００１１】
ｘ＾（ｎ）＝ｘ（ｎ）ｘ*(ｎ−τ)／ｘ*(ｎ−τ)
＝ｘ（ｎ）
抽出信号が干渉信号ｗ（ｎ）の場合も、同様に干渉信号が再現され、受信信号ｚ（ｎ）から抽出したい信号が抽出される。
抽出したい信号波形は既知ではなく、無線信号の受信を開始した時には過去の再生信号を参照信号ｕ（ｎ−τ）として用いることができない。この場合は、受信信号ｚ（ｎ）を用いる。図示例では判定切替部１８内の切替スイッチ１８ａは接点ｎｃ側に接続され、入力端子１４からの受信信号ｚ（ｎ）が切替スイッチ１８ａを通じて参照信号ｕ（ｎ−τ）として信号波形再生部１７へ入力される。この場合、必要に応じて時間相関抽出部１５から出力される時間相関波形の時刻と、タイミングが一致するように、遅延部１９により受信信号ｚ（ｎ）が遅延される。信号波形再生部１７の出力信号が判定切替部１８内の判定部１８ｂに入力され、判定部１８ｂはその入力された出力信号のレベルが所定値より大になったか否かを監視し、所定値より大になれば、信号波形再生部１７より目的とする抽出信号が再生出力されたと判断して、切替スイッチ１８ａを接点ｎｏ側に切替える。これによって信号波形再生部１７よりの出力再生信号が切替スイッチ１８ａを通じて参照信号ｕ（ｎ−τ）として信号波形再生部１７に入力されるようになる。
【００１２】
このようにして、受信信号ｚ（ｎ）から抽出したい信号のみを抽出することができる。図２に示した例では特に希望信号ｘ（ｎ）を抽出して、復調器２１へ供給するようにした場合である。抽出対象指定部１６において希望信号が指定されると、抽出対象指定部１６により切替部２２が制御され、切替部２２内の切替スイッチ２２ａ及び２２ｂが共に接点２２ｓ側に切替えられ、信号波形再生部１７よりの再生希望信号波形ｘ＾（ｎ）は切替スイッチ２２ａ及び２２ｂの各接点２２ｓ側を通じて復調器２１へ供給される。抽出対象指定部１６で干渉信号が指定されると、抽出対象指定部１６により切替部２２が制御され、その切替スイッチ２２ａ及び２２ｂは共に接点ｉ側に切替えられる。信号波形再生部１７よりの再生干渉信号ｗ＾（ｎ）は切替スイッチ２２ａの接点ｉ側を通じて減算部２３へ入力され、減算部２３で入力端子１４からの受信信号ｚ（ｎ）から再生干渉信号ｗ＾（ｎ）が差し引かれて、この受信信号ｚ（ｎ）から干渉信号ｗ＾（ｎ）が除去されたｚ（ｎ）−ｗ＾（ｎ）が切替スイッチ２２ｂの接点ｉを通じて復調器２１へ供給される。このようにして復調器２１では干渉に影響されずに信号を復調することができる。
【００１３】
この場合必要に応じて、再生干渉信号ｗ＾（ｎ）と受信信号ｚ（ｎ）のタイミングを合わせるため受信信号ｚ（ｎ）を遅延部２４で遅延させて減算部２３へ供給する。遅延部１９や２４としては遅延回路を用いる、あるいは受信信号を一括してメモリに蓄積し、遅延部２４としては信号再生部１７から出力される再生干渉信号の時刻と合せて、前記メモリから受信信号を読み出して減算部２３へ供給し、遅延部１９としては時間相関抽出部１５よりの時間相関波形の時刻と合せてメモリから受信信号を読み出して切替スイッチ１８ａへ供給してもよい。
以上のように信号の時間相関波形の周波数特性のみを利用することにより、干渉信号の情報を把握していない場合であっても、希望信号を抽出することが可能となる。また、干渉信号を再生する際には、受信している干渉信号を一括して取り扱うため、干渉信号の波数に関わらず回路規模を一定に保つことが可能となる。
【００１４】
［実施形態２］
実施形態２は実施形態１における時間相関抽出部１５として周期自己相関を利用するものである。
実施例１
周期自己相関を利用する時間相関抽出部１５の具体例また抽出対象指定部１６及び信号再生部１７の各具体例を図３に示す。実施例１における時間相関抽出部１５、抽出対象指定部１６及び信号再生部１７以外の構成は図２に示した構成と同様であるから特に図としては示していない。
【００１５】
時間相関抽出部１５は周期自己相関計算部１５ａと相関ピーク抽出部１５ｂと逆フーリエ変換部１５ｃとより構成される。信号ｖ（ｔ）の時間差τでの周期自己相関（Cyclic Autocorrelation）は次式で定義される。
【数１】
2005151461.tif 000003
α：信号ｖ（ｔ）の周期周波数、
干渉信号ｗk(ｔ)がＫ（ｋ＝１，２，…，Ｋ）波存在する場合の受信信号ｚ（ｔ）を
ｚ（ｔ）＝ｘ（ｔ）＋Σk=1Kｗk(ｔ)＋ｎ（ｔ）
とすると、希望信号と干渉信号が無相関であれば、受信信号ｚ（ｔ）の周期自己相関値Ｒz (τ，α)は以下のようにそれぞれの周期自己相関値の和となる。
Ｒz(τ，α)＝Ｒx(τ，α)＋Σk=1KＲwk（τ，α）
実際の演算においてはディジタルに処理する。この場合の周期自己相関Ｒz （τ,α）は次式の離散フーリエ変換ＤＦＴにより演算される。
【００１６】
Ｒz(τ，α)＝ＤＦＴ（ｚ（ｎ）ｚ*（ｎ−τ））
受信信号ｚ（ｎ）のサンプリング周期をＴs ＝１／ｆｓ、ａを１以上の整数とするとτ＝…，−２ａＴs ，−ａＴs ，０，ａＴs ，２ａＴs ，…である。τがその受信信号ｚ（ｎ）におけるパルス長（ほぼ信号帯域幅の逆数）までのｚ（ｎ）とｚ（ｎ−τ）との相関はかなり大きいが、τがこれより大きくなると相関は比較的急に小さくなる。従ってτは最大でパルス長の２〜３倍程度であればよく、τの刻み間隔ａＴs は小さい方がよいが、受信信号に対するサンプリングは一般にオーバーサンプリングとされているからパルスレートの４〜８倍の逆数が現実的である。なおパルスレートはＣＤＭＡのような符号多重通信方式ではチップレートであり、拡散率ＳＦ＝１の場合やＴＤＭＡの場合などはパルスレートはシンボルレートである。
【００１７】
シンボルレートが異なる２つのＢＰＳＫ信号とＡＭ信号とが混在する受信信号ｚ（ｎ）の周期自己相関値を計算した例を、τ−α平面に対し垂直に相関値を表わすと、例えば図４に示すようになる。相関値の遅延時間τに対する特性はτ＝０の軸に対し対称に表われ、ＢＰＳＫ信号（１）は周期周波数α１で−τ１とτ１をそれぞれ中心としたピークが現われ、ＢＰＳＫ信号（１）よりもシンボルレートが速い、ＢＰＳＫ信号（２）は周期周波数α２で−τ２とτ２をそれぞれ中心としたピークが現われる。ＡＭ信号は周期周波数α３に、τ＝０を中心にピークが漸次小さくなる小山が連続的に現われる。
【００１８】
希望信号及び干渉信号の信号帯域幅や搬送波周波数が異なると、希望信号と干渉信号の周期自己相関値のピークが現れる周期周波数αと時間差τの範囲も異なるため、相関ピーク抽出部１５ｂでは、ある時間差τ及び希望信号また干渉信号のピークが出現する周期周波数の範囲についてのみ周期自己相関値Ｒx(τ，α)またはＲw(τ，α)を取り出すことが可能となる。以下の説明では干渉信号（干渉波）は１波の場合とするが、複数の干渉信号（干渉波）の場合も同様に対応できることは容易に理解されよう。例えば受信したい希望信号が予め決められた複数種のいずれかを選択的に抽出したい場合は図３中の抽出対象指定部１６内に、これら各希望信号ｘ１，ｘ２，…についてその中心のτの値と、各周期周波数成分をそれぞれ記憶部１６ａに記憶しておき、抽出したい希望信号ｘを抽出対象指定部１６に入力すると、そのｘで記憶部１６ａの対応するτと周期周波数α…とを読み出し、そのパルス長と対応するτ−Δτ〜τ＋Δτと周期周波数α…を相関ピーク抽出部１５ｂに入力する。相関ピーク抽出部１５ｂには、周期自己相関計算部１５ａ中の各τ１，…，τＰについての演算部１５ａ1 ，…，１５ａP よりの各周期自己相関値の各周波数成分が入力されている。演算部１５ａ1 ，…，１５ａPと対応してその周波数成分ごとのゲートを備えるゲート群１５ｂ1 ，…、１５ｂP が相関ピーク抽出部１５ｂに設けられ、これらゲート群１５ｂの各ゲート中の入力されたτ−Δτ〜τ＋Δτとα，…とそれぞれ対応するものが開とされ、他のゲートは全て閉とされ、周期自己相関波形Ｒz (τ１，α)，…，Ｒz (τＰ，α)中の入力されたτ及び周期周波数αと対応する成分が通過され、希望信号ｘの周期自己相関値Ｒx(τ，α)が抽出される。
【００１９】
このようにして相関ピーク抽出部１５ｂで抽出された希望信号の周期自己相関成分は逆フーリエ変換部１５ｃにおける、ゲート群１５ｂ1 ，…，１５ｂP とそれぞれ対応するτ１変換部１５ｃ1 ，…,τｐ変換部１５ｃP でそれぞれ逆フーリエ変換されて、希望信号ｘの時間相関波形ｘ(ｎ)ｘ*(ｎ−（τ＋Δτ）)，…，ｘ(ｎ)ｘ*(ｎ−（τ−Δτ）)が出力される。
干渉信号を抽出する場合は、希望信号ｘに対する干渉信号として抽出対象指定部１６に指定入力され、抽出対象指定部１６は希望信号ｘと対応するτとα，…を記憶部１６ａから読み出し、相関ピーク抽出部１５ｂに対し、各τ−Δτ〜τ＋Δτ、各α，…と対応するゲートを閉とし、他の全てのゲートを開とする。このようにして希望信号ｘに対する干渉信号ｗ（ｎ）の各τについての周期自己相関値Ｒw(τ，α)の各周波数成分が相関ピーク抽出部１５ｂにより抽出され、これらが逆フーリエ変換部５ｃで逆フーリエ変換されて干渉信号の時間相関波形ｗ（ｎ）ｗ*（ｎ−τ１），…，ｗ（ｎ）ｗ*（ｎ−τＰ）が得られる。
時間相関波形が得られた後の信号再生部１７での処理その他は実施形態１と同様である。
【００２０】
信号波形再生部
信号波形再生部１７では時間相関抽出部１５の出力のあるτｐについての例えばｘ（ｎ）ｘ*（ｎ−τｐ）を、既に再生された同一τｐの信号を用い、ｘ＾*（ｎ−τｐ）で割算すればよいが、再生信号波形ｘ＾（ｘ）のゼロクロス点付近の時刻に基づくｘ＾*（ｎ−τ）で割算する状態が生じることがある。この場合この割算を十分大きなビット数で処理することができればよいが、取り扱う処理ビット数に制限されて、割算結果が不定となるおそれがある。このような問題を避けるには信号再生部１７を例えば図３中に示すように構成すればよい。この例ではτの刻み間隔（例えばτ２−τ１）ａＴs はａ＝１／２とした場合である。つまり１つのτについて、２つのサンプリング時刻が生じる場合である。
【００２１】
時間相関抽出部１５よりの各τ１，…，τＰについての時間相関波形が選択部１７ａに入力され、選択部１７ａにより選択された１つの時間相関値が除算部１７ｂにおいて参照信号検出部１７ｃよりの参照信号ｕ*(ｎ−τｐ)により割算されて信号再生部１７の出力とされると共に、バッファ１７ｄ内の各τ１，…，τＰ領域１７ｄ１，…，１７ｄＰに格納されている割算結果（再生信号）を１サンプル時刻分過去の領域へ、図３で左側に移し、かつ領域１７ｄ１，…，１７ｄＰ内の、各τｐごとの対応サンプル時刻（この例では０又は１）の部分に、その割算結果を切替スイッチ１８ａを通じて格納する。
【００２２】
次の処理時刻（サンプル時刻）では参照信号検出部１７ｃでバッファ１７ｄに格納されている過去の割算結果（再生信号）中の除算部１７ｂでの安定な割算を行うことができるもの、例えばレベルが最も大きいものが検出され、参照信号ｕ*(ｎ−τｐ)として出力され、またその検出した過去の割算結果（再生信号）のτｐと対応した、時間相関抽出部１５の出力時間相関波形を選択部１７ａで選択して除算部１７ｂへ出力する。以下同様の処理をする。
なお実施形態１で説明したようにこの信号抽出装置の動作開始時においてはバッファ１７ｄになんら信号が格納されていない。そこで入力端子１４からの受信信号ｚ（ｎ）が、少くとも数サンプル（処理時刻）分、切替スイッチ１８ａを通じてバッファ１７ｄに順次格納され、その後、参照信号検出部１７ｃで前述したようにバッファ１７ｄ内に格納された信号から参照信号を検出し、また選択部１７ａでの選択を行って除算部１７での割算を行うようにする。
【００２３】
実施例２
実施例２は時間相関抽出部１５において、周期自己相関を計算する手法を用いる具体的構成としている為実施例１と異なるだけであり、その他、抽出対象指定部１６、信号再生部１７などは実施例１又は実施形態１と同様に構成される。よってその時間相関抽出部１５の構成を図５に示してこの部分についてのみ説明する。
入力された受信信号ｚ（ｎ）は時間相関計算部１５ｄ１，…，１５ｄＰでその各遅延部１５ｅによりそれぞれτ１，…，τＰ遅延されたものと遅延されないものとが各乗算部１５ｆでそれぞれ乗算されて、時間差τ１，…，τＰについての時間相関波形ｚ（ｎ）ｚ*(ｎ−τ１),…，ｚ（ｎ）ｚ*(ｎ−τＰ)が計算される。これら時間相関計算部１５ｄ１，…，１５ｄＰで計算された時間相関波形はそれぞれ帯域通過フィルタ部１５ｇ１，…，１５ｇＰへ供給され、これら各帯域通過フィルタ部は通過周波数特性が可変の小帯域部の複数よりなるフィルタバンク１５ｈが設けられ、入力されたα，…に応じて小帯域部の周波数特性が制御され、フィルタバンク１５ｈを通過した周波数成分はそれぞれ合成部１５ｉで合成されて通過帯域フィルタ部１５ｇ１，…，１５ｇＰより時間相関波形ｘ（ｎ）ｘ*(ｎ−τ１)，…，ｘ（ｎ）ｘ*(ｎ−τＰ)又はｗ（ｎ）ｗ*(ｎ−τ１)，…，ｗ（ｎ）ｗ*(ｎ−τＰ)がそれぞれ出力される。
【００２４】
帯域通過フィルタ部１５ｇ１，…，１５ｇＰには、例えば図３を参照して説明した抽出対象指定部１６からの希望信号ｘに対するτｐ−Δτ〜τｐ＋Δτとα，…が入力され、そのτｐ−Δτ〜τｐ＋Δτと対応する帯域通過フィルタ部内のフィルタバンク１５ｈにおいて、入力されたα，…と対応する周波数成分のみが通過される。希望信号ｘに対する干渉信号が指定入力された場合は、前記希望信号ｘ指定時に通過されたτ、αの成分のみが遮断され、他のτ，α成分は通過されるように帯域通過フィルタ部１５ｇ１，…，１５ｇＰが制御される。これら帯域通過フィルタ部１５ｇ１，…，１５ｇＰの各出力は、それぞれ合成部１５ｉ１，…，１５ｉＰで合成されて希望信号ｘ（ｎ）の時間相関波形又は干渉信号ｗ（ｎ）の時間相関波形として出力される。
【００２５】
帯域通過フィルタ部１５ｇ１，…，１５ｇＰとして例えば可変特性ディジタルフィルタが用いられ、その各ディジタルフィルタのフィルタ係数を、通過又は遮断τ，αの組に応じて設定してもよい。この場合、これらフィルタ係数を、抽出対象指定部１６内の記憶部に、希望信号ｘ１，ｘ２，…ごとに、これを抽出する場合と遮断する場合とをそれぞれ予め計算して格納しておき、希望信号の抽出か、希望信号に対する干渉信号の抽出かに応じて対応するフィルタ係数を読み出し、帯域通過フィルタ部１５ｇ１，…，１５ｇＰの対応するものに設定してもよい。
上述した実施例１及び２のいずれも希望信号と干渉信号が混ざり合った信号から、時間相関計算に用いる遅延時間及び時間相関波形の周波数特性を利用して、高精度に希望信号（干渉信号）を分離・抽出することが可能となる。
【００２６】
実施例３
この実施例３は１つの時間相関計算部と１つの帯域通過フィルタ部とのみで時間相関抽出部１５を構成する場合である。例えば図５中に２点鎖線枠で示すように時間相関計算部１５ｄ２と帯域通過フィルタ部１５ｇ２とにより時間相関抽出部１５が構成される。この場合の遅延部１５ｅの遅延量τは例えば希望信号のパルス長の２分の１程度とされる。つまり、時間相関抽出部１５では予め決めた特定のτについてのみ時間相関波形を抽出する。
【００２７】
この時、得られる時間相関波形は希望信号と干渉信号との間で変調方式、信号帯域幅、又は搬送波周波数のいずれかが異なれば、例えば図６Ａに示すように受信信号ｚ（ｎ）の時間相関波形ｚ（ｎ）ｚ*(ｎ−τ)中の希望信号ｘ（ｎ）の時間相関成分ｘ（ｎ）ｘ*(ｎ−τ)の周波数成分と、干渉信号ｗ（ｎ）の時間相関成分ｗ（ｎ）ｗ*(ｎ−τ)の周波数成分とが分離して現われる。よって抽出対象指定部１６から指定された抽出対象が希望信号であれば、図６Ｂに示すように希望信号成分のみを通過させるように帯域通過フィルタ部１５ｇを制御することにより、希望信号ｘ（ｎ）の時間相関波形が抽出される。また指定された抽出対象が干渉信号であれば、帯域通過フィルタ部１５ｇを、希望信号成分のみを遮断し、他の全ての成分を通過させるように制御することにより、干渉信号ｗ（ｎ）の時間相関波形が抽出される。なお図６Ａは図４に示す時間相関波形特性においてτ軸中の特定の値における周期周波数特性を表わしている。
【００２８】
抽出したい信号の時間相関波形を抽出後の処理は、前述した各種実施例と同様である。ただし信号波形再生部１７においては、時間相関抽出部１５で用いられるτが１つの固定値であるから、図３に示したように過去の複数の再生信号をバッファ１７ｄに格納して、除算部１７ｂから正しい割算結果が得られるものを参照信号として選択することはできない。しかし、除算部１７ｂでの演算桁数（ビット数）がある程度大きければ、演算が不定になる程参照信号のレベルが著しく小さくなることは例外的になり、それ程問題はない。それよりも実施例３では演算処理量が実施例１や２よりも著しく少なくて済む利点がある。つまりある遅延時間に対する時間相関波形の周波数特性を利用して、容易に希望信号と干渉信号を分離・抽出することが可能となる。
【００２９】
［実施形態３］
この実施形態３は希望信号の復調に影響を与えるような干渉信号が存在するか否かを予め検出するものである。
実施例４
この実施例４は先に述べた各実施例において抽出対象指定部１６により抽出対象を希望信号にするか干渉信号にするかを、自動的に指定するものである。通常の無線受信機は常時は待ち受け状態になっており、相手送信局から通信要求信号を受信した後に、希望信号が受信され、その通信要求信号は例えば制御チャネルにより送信されるなど希望信号と分離が可能とされている。
入力端子１４よりの受信信号は信号分離部２７に入力され、信号分離部２７で通信要求信号が分離されて抽出対象指定部１６内の通信要求受付部１６ｂに入力される。通信要求信号の分離はＡ／Ｄ変換器１３（図２）より前段で行ってもよい。通信要求受付部１６ｂで通信要求信号より通信要求を検出すると、干渉有無判定部１６ｃに制御信号を送る。干渉有無判定部１６ｃはこの制御信号を受けると、時間相関抽出部１５における希望信号が受信されていない状態での演算された時間相関波形レベルが所定値より大きいか否かを判定する。例えば周期自己相関計算部１５ａから、希望信号送信前の干渉信号と雑音のみからなる受信信号
ｚ（ｎ）＝Σk=1Kｗk(ｎ)＋ｎ（ｎ）
についての周期自己相関値Ｒz(τ１，α)〜Ｒz(τＰ，α)を取り込み、その各ピーク値を測定し（図８、Ｓ１）、この測定したピーク値中に所定値より大きいものがあるか否か調べ（図８、Ｓ２）、所定値以上のものがあれば、復調に影響を与える干渉信号があると判定して抽出対象を干渉信号とする指定を抽出対象決定部１６ｄへ入力し、ピーク値中に所定値より大きいものがなければ、復調に影響を与える干渉信号がないと判定して抽出対象を希望信号とする指定を抽出対象決定部１６ｄへ入力する（図８、Ｓ２）。
【００３０】
抽出対象決定部１６ｄは指定入力が干渉信号であれば干渉信号の時間相関波形を抽出するように時間相関抽出部１５を設定する。例えば先に説明したように相関ピーク抽出部１５ｂにおいて希望信号の各相関値成分のみを阻止し、他の成分は通過させるように相関ピーク抽出部１５ｂを設定する。また切替部２２を接点ｉ側に設定する（図８、Ｓ３）。抽出対象決定部１６ｄは指定入力が希望信号であれば、希望信号の各相関値成分のみを通過させ、他の成分を阻止するように相関ピーク抽出部１５ｂを設定し、切替部２２を接点ｓ側に設定する（図８、Ｓ４）。抽出対象決定部１６ｄには相関ピーク値が所定より大のものがあるか否かのみを入力してもよく、必要に応じて時間相関抽出部１５に対する設定のための起動指令を入力する。切替部２２に対する設定も干渉有無判定部１６ｃの出力によってもよい。時間相関抽出部１５に対する設定のためには、例えば抽出対象決定部１６ｄに、図３中に示した記憶部１６ａが設けられ、希望信号に対するτ，α，…についての情報が保持されている。
【００３１】
受信信号ｚ（ｎ）中の希望信号ｘ（ｎ）の形式は予め決められた１つの場合は、希望信号が抽出対象とされると、抽出対象決定部１６ｄでその形式に対するτ，α，…を出力すればよい。希望信号ｘ（ｎ）の形式が決められていない場合は、無線受信機からどのような形式の希望信号ｘ（ｎ）を送信するべきかを送信局に通知する場合、逆に送信局から通信要求信号中に希望信号の形式情報を付加するなどにより抽出対象決定部１６ｄでは希望信号の形式がわかるようにされる。時間相関抽出部１５としては図５に示したものあるいは実施例３で説明したものでもよく、これらの場合は時間相関抽出部１５で抽出可能な全ての成分が出力されるように時間相関抽出部１５内の帯域通過フィルタ部１５ｇを全帯域通過特性に設定してその出力を図７中に破線で示すように干渉有無判定部１６ｃに入力し、その出力中のピーク値が所定値を越えるものがあるか否かを判定すればよい。いずれにしても、この判定は、例えば１０００〜２０００サンプル程度の受信信号がｚ（ｎ）が時間相関抽出部１５に入力されれば行うことができる。
【００３２】
以上のように時間相関抽出部１５、切替部２２に対する設定が行われた後、希望信号が受信され、干渉信号があれば、これを除去して復調器２１へ供給することができる。このようにすることにより、干渉信号の受信レベルが低く、希望信号送信開始後の受信信号がほとんど希望信号と雑音からなる受信信号ｚ（ｎ）＝ｘ（ｎ）＋ｎ（ｎ）となる場合においても、信号波形再生部１７で、この受信信号を参照信号とすることにより、希望信号波形を高精度に再生することが可能となる。しかも、このようにして希望信号の波形を再生するようにしている通信中において干渉信号が新たに発生する場合であっても、干渉信号の影響を受けずに希望信号成分のみを抽出することが可能となる。一方、強く受信される干渉信号が存在する場合は、希望信号送信開始前の干渉信号と雑音からなる受信信号ｚ（ｎ）＝ｗ（ｎ）＋ｎ（ｎ）を参照信号とすることによって、干渉信号波形の再生を開始することにより、干渉信号を高精度に再生することが可能となりここで、信号波形再生部１７が再生信号波形を出力し始めると、切替スイッチ１８ａを接点ｎｏ側に切り替えて再生した波形を参照信号として帰還することにより、干渉信号波形を高精度に再生することができ、送信を開始した後も、その影響を受けずに干渉信号波形を再生し、これを用いて、受信信号から干渉信号波形を除去することにより希望信号波形を高精度に再生することができる。希望信号を送信する時刻が予めわかっている時など、希望信号の送信開始がわかる場合は、その送信開始になるべく直前に受信される干渉信号の状態を検出して、時間相関抽出部１５，切替部２２などに対する設定を行えばよく、要は希望信号が送信されていない状態に干渉信号の有無を検出すればよく通信要求信号は必ずしも必要としない。
実施例５
この実施例５は時間相関抽出部１５で希望信号の時間相関波形と干渉信号の時間相関波形とを同時に抽出可能とされ、かつ信号波形再生部１７として、希望信号の希望波形再生部１７ｓと干渉信号の干渉波形再生部１７ｉとが設けられる。時間相関抽出部１５としては先に示したいずれのものを使用してもよいが、この図９では図３に示した例を用いる。
【００３３】
実施例４と同様に、通信要求が通信要求受付部１６ｂで検出され、その出力により干渉有無判定部１６ｃが制御され、干渉有無判定部１６ｃは時間相関抽出部１５で生成された希望信号送信前の干渉信号と雑音のみからなる受信信号についての時間相関波形、この例では周期自己相関値Ｒz(τ，α）のピーク値が予め定めた所定値より大であるか否かを判定して、所定値より大であれば干渉信号がありと、所定値より大でなければ干渉信号なしと判定して抽出対象決定部１６ｄへ入力する。
抽出対象決定部１６ｄにその判定結果が入力されると抽出対象決定部１６ｄは相関ピーク抽出部１５ｂを制御して希望信号の周期自己相関値成分を希望用逆フーリエ変換部１５ｃｓへ入力し、それら以外の周期自己相関値成分を干渉用フーリエ変換部１５ｃｉへ入力するように受信信号ｚ（ｎ）の周期自己相関値成分を振り分ける。この振り分け制御は、先に述べた希望信号の周期自己相関値成分の選出、干渉信号の周期自己相関値成分の選出と同様に行うことができる。希望用逆フーリエ変換部１５ｃｓから希望信号時間相関波形が、干渉用逆フーリエ変換部１５ｃｉから干渉信号時間相関波形がそれぞれ出力され、前者は希望波形再生部１７ｓへ、後者は干渉波形再生部１７ｉへそれぞれ供給される。なお希望信号送信前では希望信号時間相関波形は抽出されない。
【００３４】
抽出対象指定部１６で抽出対象が希望信号と指定（判定）されていれば、抽出対象指定部１６は判定切替部１８の切替スイッチ１８ａｓを接点ｎｃ側に、切替スイッチ１８ａｉを接点ｎｏ側に切替えて、受信信号ｚ（ｎ）が切替スイッチ１８ａｓの接点ｎｃ側を通じて希望波形再生部１７ｓへ参照信号として供給され、また切替部２２の切替スイッチ２２ｂは接点ｓ側に接続され、抽出対象が干渉信号と指定（判定）されていれば切替スイッチ１８ａｓを接点ｎｏ側に、切替スイッチ１８ａｉを接点ｎｃ側に切替えて、受信信号ｚ（ｎ）が切替スイッチ１８ａｉの接点ｎｃ側を通じて干渉波形再生部１７ｉに参照信号として供給され、また切替スイッチ２２ｂは接点ｉ側に接続される。
【００３５】
判定切替部１８内の判定部１８ｂｓ及び１８ｂｉには希望波形再生部１７ｓ及び干渉波形再生部１７ｉの各出力信号がそれぞれ入力され、各判定部１８ｂｓ及び１８ｂｉはその出力信号のレベルが所定値を超えると、それぞれ再生信号を出力している状態と判定して切替スイッチ１８ａｓ及び１８ａｉを接点ｎｏ側に切替える。希望波形再生部１７ｓでは切替スイッチ２２ｂの出力が切替スイッチ１８ａｓの接点ｎｏ側を通じて参照信号として入力され、干渉波形再生部１７ｉではその出力信号が切替スイッチ１８ａｉの接点ｎｏ側を通じて参照信号として入力される。切替スイッチ２２ｂの出力信号は判定部２２ｃにも入力され、判定部２２ｃはこの入力信号を監視し、そのレベルが所定値を超えると希望信号の再生信号が入力されたと判定して、切替スイッチ２２ｂを接点ｓ側に切替える。これにより希望波形再生部１７ｓよりの再生信号波形が切替スイッチ２２ｂの接点ｓ側を通じて復調器２１へ供給されることになる。切替スイッチ１８ａｉ，１８ａｓ，２２ｂの各状態における設定状態を図１０に示す。
【００３６】
このように構成されているから、図１１に示すように希望信号送信前に干渉有無判定部１６ｃで干渉信号が有る、つまり抽出対象が希望信号でないと判定されると（Ｓ１）、切替スイッチ１８ａｓ，１８ａｉ及び２２ｂはそれぞれ接点ｎｏ側、ｎｃ側及びｉ側に設定され（Ｓ２）、干渉波形再生部１７ｉに受信信号ｚ（ｎ）が参照信号として入力され、干渉信号波形が精度よく再生され、受信信号ｚ（ｎ）に希望信号が重畳された状態になると、この受信信号ｚ（ｎ）から再生された干渉信号が減算部２３で差し引かれ、干渉信号が抑圧された希望信号が切替スイッチ２２ｂを通じて復調器２１へ供給される。このような状態になるとその干渉信号が抑圧された希望信号は判定部２２ｃに入力されると共に切替スイッチ１８ａｓを通じて希望波形再生部１７ｓに参照信号として入力され、またこの希望波形再生部１７ｓに時間相関抽出部１５で抽出された希望信号の時間相関波形が入力され、かつ判定部２２ｃはその入力信号が所定レベルを超えたことを検出して（Ｓ３）、切替スイッチ２２ｂを接点ｓ側に切替える（Ｓ４）。従って希望波形再生部１７ｓで精度よく希望信号が再生されるようになり、この再生希望信号が切替スイッチ２２ｂを通じて復調器２１へ供給される。
【００３７】
つまり希望信号送信開始時の抽出対象が干渉信号であっても、希望信号を高精度に再生することができ、切替スイッチ２２ｂが接点ｓ側に切替えられると、希望信号のみに着目して信号を再生する状態になり、希望信号に新たに別の干渉信号が重畳されるようになるなど伝搬路環境が大きく変化しても、安定かつ高精度に希望信号を再生することができる。
一方、希望信号送信前に干渉有無判定部１６ｃで干渉信号がないと判定され、つまり抽出対象が希望信号であれば（Ｓ１）、切替スイッチ１８ａｓ，１８ａｉ及び２２ｂはそれぞれ接点ｎｃ側、ｎｏ又はｎｃ側及びａｓ側に設定され（Ｓ５）、図７に示した実施例におけるその場合と同様に希望波形再生部１７ｓには受信信号ｚ（ｎ）が参照信号として入力され、受信信号ｚ（ｎ）に希望信号が重畳された状態になると希望波形再生部１７ｓから希望信号が精度よく再生されるようになり、希望波形再生部１７ｓの出力レベルが所定値以上となり（Ｓ６）、切替スイッチ１８ａｓが接点ｎｏ側に切替わり（Ｓ７）、希望波形再生部１７ｓの出力が切替スイッチ２２ｂを通じて復調器２１へ供給される。
【００３８】
時間相関抽出部１５として図５に示した例を用いる場合は希望信号時間相関波形抽出用の帯域通過フィルタ部１５ｇ１〜１５ｇＰと干渉信号時間相関波形抽出用の帯域通過フィルタ部１５ｇ１〜１５ｇＰとを設ければよく、実施例３で説明したものを用いる場合は希望信号成分のみを通過させる帯域通過フィルタ部と、希望信号成分を遮断し、他の成分のみを通過させる帯域通過フィルタ部と受信信号の時間相関波形を供給するようにすればよい。また抽出対象指定部１６において、図２、図３を参照して説明したように通信要求受付部１６ｂ及び干渉有無判定部１６ｃを設けることなく、抽出対象を決定する場合にもこの実施例５を適用することもできる。
【００３９】
［実施形態４］
この実施形態はこの発明の無線通信方法の実施形態であり、これを適用したシステム構成例を図１２に、その処理の流れを図１３にそれぞれ示す。送信側装置１００と受信側装置２００との間で送受信が行われる。
送信側装置１００は、情報信号の送信を始める前に通信要求生成部１１０で、通信要求信号ＳD を生成して、合成部１１１，アンテナ１１２を通じて受信側装置２００へ送信する。この送信は例えば制御チャネルを通じて行う。受信側装置２００では、アンテナ２１０で受信し、信号分離部２１１で分離した通信要求信号を通信要求受付部２１２で受け付ける。通信要求受付部２１２は周期自己相関計算部２１３へ通信要求が到来したことを示す制御信号を送る。制御信号を受けた周期自己相関計算部２１３では、希望信号の送信を開始する前の、アンテナ２１０で受信され、信号分離部２１１で分離された干渉信号と雑音のみからなる受信信号について周期自己相関値を計算し、その計算結果をパラメータ決定部２１４へ入力する。パラメータ決定部２１４は、入力された周期自己相関値から、好ましいつまり干渉信号の周期自己相関値がゼロ又は最も小さい遅延時間τと周期周波数αに希望信号の周期自己相関のピークをもつ希望信号を探し、その希望信号の送信に使用する送信信号パラメータ（変調方式Ｍｏｄ、信号帯域幅Ｂ、シンボルレートＲｓ、符号化率ｒ及び搬送波周波数ｆｃ）Ｐｓを決定し、送信許可と共に送信部２１５、アンテナ２１６を通じて送信側装置１００例えば制御チャネルを通じて送信する。また、同時に決定した送信信号パラメータＰｓを信号処理部２１７及び復調器２１８へ入力し、同期検波の基準周波数や帯域制限フィルタの通過帯域等を設定させる。この送信信号パラメータＰｓの決定のやりかたは後で説明する。
【００４０】
送信側装置１００では、アンテナ１１３、受信部１１４を通じて受信した送信信号パラメータＰｓを分配部１１５により変調方式制御部１１６、発振周波数制御部１１７、信号帯域幅・シンボルレート制御部１１８及び誤り訂正符号化部１１９へ分配入力する。誤り訂正符号化部１１９は、受信したパラメータＰｓ中の符号化率に関する情報に基づき、送信データを誤り訂正符号化してベースバンド変調器１２１へ入力する。誤り訂正符号化方法は予め決めたものを用い、その符号化率がパラメータＰｓ中の符号化率ｒになるようにする。変調方式制御部１１６及び信号帯域幅・シンボルレート制御部１１８は、ベースバンド変調器１２１を、受信側装置２００で定められた変調方式Ｍｏｄ、信号帯域幅Ｂ及びシンボルレートの信号を出力するように制御する制御信号を出力する。ベースバンド変調器１２１は、誤り訂正符号化部１１９から入力された誤り訂正符号化された送信データを、変調方式制御部１１６及び信号帯域幅・シンボルレート制御部１１８からの制御信号に従って変調する。このとき、ベースバンド変調器１２１は、信号帯域幅・シンボルレート制御部１１８から入力された信号帯域幅制御信号に応じてパルスレートを制御するだけでなく、変調信号にかける帯域通過フィルタの通過領域・遮断領域を制御する。発振周波数制御部１１７は、周波数変換部１２２から出力される信号の搬送波周波数ｆｃが受信側装置２００から指定された搬送波周波数ｆｃとなるように、局部発振器１２３を制御する。このようにして、変調方式制御部１１６、発振周波数制御部１１７、信号帯域幅・シンボルレート制御部１１８及び誤り訂正符号化部１１９は、ベースバンド変調器１２１及び局部発振器１２３を制御することにより、送信信号を受信側装置２００において決定されたパラメータＰｓを有する信号にして送信を開始する。
【００４１】
ＣＤＭＡのような符号分割多重を行う場合の各制御部１１７〜１１８、誤り訂正符号化部１１９、ベースバンド変調器１２１についての具体的動作例を、図１４を参照して説明する。分配部１１５は、通知された送信信号パラメータ（変調方式Ｍｏｄ、信号帯域幅Ｂ又はシンボルレートＲｓ、搬送波周波数ｆｃ、拡散率ＳＦ、符号化率ｒ）を、Ｍｏｄを変調方式制御部１１６へ、ＢとＳＦをシンボルレート制御部１１８へ、ｒを誤り訂正符号化部１１９へ、ｆｃを周波数制御部１１７へそれぞれ分配入力する。送信データは、まず誤り訂正符号化部１１９へ入力される。誤り訂正符号化部１１９では、複数の符号化率ｒ＝１／２，３／３，９／１６，…を実現する符号化器１１９ａ１，１１９ａ２，…を用意しておき、スイッチ制御部１１９ｂにより指定されたパラメータｒの符号化率を実現する符号化器に送信データが入力されるように、スイッチ１１９ｃの切り替えを行う。これにより、送信データは符号化率ｒで符号化されたデータとなり、ベースバンド変調器１２１へ入力される。変調方式制御部１１６は、入力されたパラメータＭｏｄで指定された変調方式で変調された信号が出力されるように、ベースバンド変調器１２１内の変調部１２１ａの変調方式を設定し、また、シンボルレート制御部１１８は、パラメータＢからチップレートＲｃを計算し、一般には便宜的にＢ＝Ｒｃとして求めたチップレートＲｃとＳＦからシンボルレートＲｓ＝Ｒｃ／ＳＦを求め、そのシンボルレートＲｓに応じて、変調部１２１ａが出力する信号のシンボルレートを制御する誤り訂正符号化された送信データは変調部１２１ａでいわゆるシンボルへのマッピングが行われる。変調方式制御部１１６とシンボルレート制御部１１８は、変調方式Ｍｏｄと拡散率ＳＦを実現する拡散符号を生成するように拡散符号を生成し、この拡散符号が変調部１２１ａから出力される信号と乗算部１２１ｃで乗算される。このようにすることで、符号化されたデータは、ベースバンド変調器１２１において、変調方式Ｍｏｄで変調され、帯域幅Ｂを持つ拡散率ＳＦのベースバンド信号に変換される。続いて、このベースバンド信号へ入力される。周波数制御部１１７は、入力されたパラメータｆｃを実現するために、局部発振器１２３の発振周波数を制御し、この局部発振器１２３の発振信号が周波数変換部１２２へ入力され、周波数変換部１２２から出力された信号の、中心周波数がｆｃに変換される。なお拡散率ＳＦを１とすることにより、他の多重化方式ＴＤＭＡ，ＦＤＭＡなどにも上記手法は適用できる。この場合はＲｃ＝Ｒｓ＝Ｒｐ(Ｒｐはパルスレート)となる。
送信信号パラメータ決定１
図１２中のパラメータ決定部２１４の機能構成例を図１５に示す。送信側で使用可能なパラメータ（変調方式、帯域幅、及び搬送波周波数）をＭｏｄ格納部２１４ａ、ｆｃ格納部２１４ｂ、Ｂ格納部２１４ｃにデータとして保有する。周期自己相関出現領域計算部２１４ｄは、保有しているパラメータＭｏｄ，ｆｃ，Ｂの全組み合わせに対し、周期自己相関値Ｒ′x(τ，α)が出現する（τ′，α′）を計算により求め、使用可能パラメータ抽出部２１４ｅにこれら各出現τ′，α′とパラメータ（Ｍｏｄ，ｆｃ，Ｂ）との組を入力する。つまり例えば前記パラメータＭｏｄ，ｆｃ，Ｂの全組合せのそれぞれについて、そのパラメータ（Ｍｏｄ，ｆｃ，Ｂ）から無線信号の式ｘ（ｎ）を作り、このｘ（ｎ）について前述したように周期自己相関Ｒ′x(τ，α)を計算すればよい。あるいは変調方式Ｍｏｄが決るとピークが立つ周期周波数αは搬送波周波数ｆｃと帯域幅Ｂを用い、又は帯域幅Ｂのみを用いる簡単な計算で求まることが比較的多くの変調方式について知られている。従って、使用可能な各変調方式Ｍｏｄについて、各使用可能なｆｃとＢ又はＢを用いて周期周波数αをそれぞれ求めればよい。これらのようにして求めた（τ，α）を、予め記憶部に記憶しておき、（τ，α）と（Ｍｏｄ，ｆｃ，Ｂ）との一方をアドレスとして他方を読み出すことができるようにしておいてもよい。なお、このような記憶部が図３中の抽出対象指定部１６内の記憶部１６ａとして用いられる。
【００４２】
図１２中の送信パラメータ決定部２１４における送信信号パラメータＰｓの決定方法の例を説明する。周期自己相関計算部２１３において計算された干渉信号ｗ（ｎ）と雑音ｎ（ｎ）からなる受信信号ｚ（ｎ）の周期自己相関値Ｒz(τ，α)をτ−α空間上に表わした例を図１６に示す。つまり図１６はτ−α空間上における干渉信号のピークが立っている状況を表している。送信パラメータ決定部２１４では図１６中の点線で囲った部分のように周期自己相関値Ｒz(τ，α)が十分に小さい領域を使用可能領域検出部２１４ｆで探索する。具体的には、干渉信号と雑音のみからなる受信信号の周期自己相関値を、予め定めておいた所定値と比較し、所定値以下となる領域（以下使用可能領域という）を探索し、使用可能パラメータ抽出部２１４ｅへ入力する。
【００４３】
使用可能パラメータ抽出部２１４ｅは入力されたτ′，α′の組中の前記検出された使用可能領域に属するものの１つと対応する（Ｍｏｄ，ｆｃ，Ｂ）を希望信号の送信パラメータＰｓとして抽出し、送信側装置１００、信号処理部２１７及び復調器２１８へ入力する。拡散率ＳＦや符号化率ｒは一般に使用される値がいくつかに限られているから、受信側装置２００で適当なものを決めて送信側装置１００へ通知してもよく、あるいは送信側装置１００で適当なものを選択してもよい。以上のようにして実施形態１〜３で説明した信号抽出装置において希望信号または干渉信号の時間相関波形を容易に抽出することが可能となる。
【００４４】
使用可能領域の決定に用いる前記所定値は希望信号の変調方式に応じて、変調多値数が多く、復調にＣＩＮＲ（Carrier to Interference plus background Noise power Ratio：搬送波電力／（干渉波電力＋背景雑音電力））等の受信品質が高品質であることが必要な条件では所定値を低く、逆に変調多値数が少なく、復調に必要なＣＩＮＲ等の受信品質に対する要求が高くない条件では所定値を高く設定するなど、適切な値に設定するのが望ましい。図１６に示した相関値−τ−α特性を求めるには、例えば図３中に示した周期自己相関値の計算方法によるか、あるいは図５に示した時間相関抽出部１５において各帯域通過フィルタ部１５ｇ１〜１５ｇＰを、全相関値周波数成分をそれぞれ分離して出力するように構成してこれらの出力により求めることもできる。
【００４５】
送信信号パラメータ決定２
パラメータ決定部２１４で符号化率ｒも決定する場合であり、まず可能な限り最大のビットレートの通信を可能とする例を説明する。
図１７にその機能構成の要部を示し、図１８にその処理手順を示す。図１２中の通信要求受付部２１２において、送信側装置１００からの通信要求を受け付けると（Ｓ１）、周期自己相関計算部２１３は、希望信号送信前の干渉信号と雑音からなる受信信号について、周期自己相関値を計算し、パラメータ決定部２１４にその計算結果を入力する（Ｓ２）。パラメータ決定部２１４は、前述したようにして、使用可能領域に希望信号の周期自己相関値を出現させることのできるパラメータ組（変調方式Ｍｏｄ、信号帯域幅Ｂ、中心周波数ｆｃ）の全てをパラメータ組候補として抽出する（Ｓ３）。
【００４６】
一方、通信要求受付部２１２に入力された通信要求信号の電力を受信電力検出部２１２ａにより検出し、これを希望信号送信開始後の受信信号電力と推定し、この受信信号電力と、使用可能領域の周期自己相関値で割算して使用可能領域での受信品質（ＣＩＮＲ）を、ＣＩＮＲ推定部２１４ｇで推定する（Ｓ４）。誤り率推定部２１４ｈでは、推定したＣＩＮＲと抽出したパラメータ組の変調方式Ｍｏｄとにより誤り訂正符号化されていない希望信号送信開始後のビット誤り率を推定する（Ｓ５）。このビット誤り率の推定は、変調方式Ｍｏｄに応じた、ＣＩＮＲを変数とする関数計算により、あるいは、予め求まっているビット誤り率とＣＩＮＲ特性曲線を参照して求めることができる。
【００４７】
拡散率・符号化率決定部２１４ｉでは、その推定ビット誤り率が所定値以下となるように、誤り訂正符号の符号化率ｒや、ＣＤＭＡ等の符号分割多重システムを用いる場合は拡散率ＳＦを定め、ビットレート計算部２１４ｋへ入力する（Ｓ６）。拡散率ＳＦは一般にいくつかしかとらないから、その各ＳＦに対する誤り訂正符号化が行われていない信号に対するビット誤り率の改善度が知られており、この拡散率ＳＦによる改善と符号化率ｒによる改善とにより、推定ビット誤り率より、希望信号のビット誤り率が所定値以下になるようにする。ビットレート計算部２１４ｋでは、入力された符号化率ｒ及び拡散率ＳＦと、対応するパラメータ組候補の変調方式Ｍｏｄ及び信号帯域幅Ｂとにより実効的なビットレートを推定する（Ｓ７）。つまり信号帯域幅Ｂから、パルスレートＲｐ又はチップレートＲｃが決り、更に拡散率ＳＦからシンボルレートＲｓが決り、シンボルレートＲｓと変調方式Ｍｏｄから誤り訂正符号化送信データのビットレートが決る。この誤り訂正符号のビットレートは例えばＢＰＳＫの場合シンボルレートＲｓと同じでありＱＰＳＫであれば２Ｒｓとなり、１６ＱＡＭであれば４Ｒｓとなる。この誤り訂正符号のビットレートに対し符号化率ｒを乗算すれば、送信データのビットレートＲｂが求まる。
【００４８】
このようにして推定ビットレートＲｂを使用可能パラメータ抽出部２１４ｅにおいて抽出された全てのパラメータ組候補に対して求め、送信パラメータ選択部２１４ｊは、求めた推定ビットレート中の最大となるパラメータの組［Ｍｏｄ，Ｂ，ｆｃ，ＳＦ，ｒ］を希望信号の送信パラメータＲｓとして決定し（Ｓ８）、決定したパラメータ組を送信側装置１００へ送信し、通信開始状態にする（Ｓ９）。このようにすることにより、希望信号は、常に可能な限り最大のビットレートで通信を行うことが可能となる。
符号分割多重システムを用いる場合、予め決めたいくつかの拡散率ＳＦについて通信要求信号を逆拡散した時の各ＣＩＮＲを推定するがこれは例えば受信電力とＳＦの積により求めることができる。このＣＩＮＲを用いて誤り率推定部２１４ｈで誤り率を推定する。使用される符号化率ｒはいくつかに限られており、その各符号化率について、予め求められているビット誤り率−ＣＩＮＲ特性曲線図を参照して、図１７に破線で示すように、拡散率・符号化率決定部２１４ｉで推定ＣＩＮＲから所定値以下の符号化率ｒを求めてもよい。つまり図１８中に破線で示すようにステップＳ４のＣＩＮＲ推定からステップＳ６の符号化率・拡散率決定に直ちに移ってもよい。
【００４９】
次に送信側装置１００から送信データに要求されるビットレートが通知される場合の送信パラメータＰｓの決定について説明する。この前提として、送信側装置１００においてその通信要求生成部１１０（図１２参照）で通信要求信号を生成する際に所定のビットレートも要求品質を示す情報もその通信要求信号に付加される。この要求品質情報は図１７中に示すように受信側装置２００における通信要求受付部２１２内の要求品質分離部２１２ｂで分離出力される。
受信側装置２００におけるパラメータ決定部２１４では前述したようにして利用可能領域から抽出された全てのパラメータ組候補について、拡散率ＳＦ、符号化率ｒを決定し、実効的なビットレートＲｂを推定する。送信パラメータ選択部２１４ｊは図１７中に破線で示すように、要求品質分離部２１２ｂからの要求品質情報、つまり要求ビットレートが入力され、ビットレート計算部２１４ｋで計算された推定ビットレート中の要求ビットレートを超えるパラメータ組［Ｍｏｄ，Ｂ，ｆｃ，ＳＦ，ｒ］を抽出し、その中で推定ビットレートが最低となるパラメータ組を希望信号の送信パラメータＰｓとして決定し、その決定したパラメータＰｓ送信側装置１００へ送信する。このようにして、要求品質（ビットレート）の低い信号に、相関ピーク抽出部において抽出する周期自己相関値の歪みを許容させることにより、τ−α空間上にできるだけ多くの信号を収容することが可能となる。要求品質を満す送信パラメータＰｓの決定における処理手順は図１８中のステップＳ８において括弧書きで示すように、ステップＳ８で要求品質（要求ビットレート）以上で最低ビットレートを与えるパラメータ組を選択することになる。
【００５０】
図１２に示した受信側装置２００は図２、図３、図５、図９に示した信号抽出装置のいずれかを備えるものとして構成することができる。その例として図２、図７を備えたものを図１９に示し、対応する構成で参照番号が対応して異なるものは一方を括弧書きで示した。
送信側装置１００からの通信要求信号を受信すると、その通信要求信号は信号分離部２１０（２７）から通信要求受付部２１２（１６ｂ）へ送られる。通信要求受付部２１２（１６ｂ）は送信側装置からの通信要求を受けると、干渉有無判定部１６ｃへ制御信号を送る。干渉有無判定部１６ｃは、周期自己相関計算部２１３（１５ａ）において計算された、希望信号送信前の干渉信号と雑音のみからなる受信信号についての周期自己相関値Ｒz(τ，α)について、相関値の最大値（ピーク）が、予め定めておいたしきい値以上であるか否かを判定することにより、強い電力で受信されている干渉信号の有無を判定し、判定結果を抽出対象決定部１６ｄへ入力する。抽出対象指定部１６ｄは入力された判定結果に応じて相関ピーク抽出部１５ｂ及びスイッチ制御部１８ａを制御する。また、干渉有無判定部１６ｃは、干渉信号と雑音のみからなる受信信号の周期自己相関値Ｒz(τ，α)を予め定めておいた所定値と比較し、相関値が所定値以下となるτ−α空間上の領域、つまり使用可能領域を抽出し、その抽出した使用可能領域に関する情報を送信パラメータ決定部２１４へ入力する。送信パラメータ決定部２１４では、その使用可能領域に希望信号の周期自己相関値にピークを立たせることのできる送信信号のパラメータ（変調方式Ｍｏｄ、信号帯域幅Ｂ及び搬送波周波数ｆｃ）を決定する。送信パラメータ決定部２１４において決定された送信信号パラメータＰｓは、送信側装置１００へ通知されると同時に、相関ピーク抽出部１５ｂへ入力される。
【００５１】
送信側装置１００が送信パラメータ決定部２１４より通知された送信信号パラメータＰｓで信号（希望信号）の送信を開始すると、受信側装置２００で希望信号と干渉信号と雑音からなる信号が受信される。その受信信号は信号分離部２１０（２７）を通じて周期自己相関計算部２１３（１５ａ）、減算部２３、切替スイッチ１８ａに入力される。周期自己相関計算部２１３（１５ａ）では、受信信号の周期自己相関値を計算する。抽出対象決定部１６ｄで決定した抽出対象に応じて、相関ピーク抽出部１５ｂと切り替え部２２が同時に制御される。相関ピーク抽出部１５ｂは、送信パラメータ決定部２１４から入力された希望信号のパラメータを用いて、抽出対象決定部１６ｄから入力された抽出対象の情報に従い、希望信号または干渉信号の周期自己相関値のピークを抽出し、取り出した相関ピークＲx(τ，α)またはＲw(τ，α)を逆フーリエ変換部へ入力する。逆フーリエ変換部１５ｃでは、入力されたαの周期自己相関関数値を逆フーリエ変換し、希望信号（又は干渉信号）の時間相関値ｘ（ｎ）ｘ*(ｎ−τ)（又はｗ（ｎ）ｗ*(ｎ−τ)）を出力する。出力された時間相関値は、信号波形再生部１７で切替スイッチ１８ａからの出力信号より割算されて信号波形に再生される。
【００５２】
このような構成の受信側装置を用いることにより、希望信号と干渉信号の周期自己相関値のピークが重ならないように、希望信号のパラメータを設定し、周期自己相関値を用いて、それぞれの信号の相関ピークを抽出して、信号波形を再生させることにより、多くの無線送受信機が、同一時刻に同一周波数を使用することが可能となるため、周波数利用効率を改善することが可能となる。
図２、図３、図５、図９、図１９などに示した装置をコンピュータにより機能させることもできる。この場合はコンピュータにこれら装置のいずれかとして機能させるためのプログラムをＣＤ−ＲＯＭ、磁気ディスクなどの記録媒体からインストールまたは通信回線を介してダウンロードし、そのプログラムを実行させればよい。
【図面の簡単な説明】
【００５３】
【図１】従来のレプリカ生成型干渉キャンセラの機能構成を示すブロック図。
【図２】この発明による信号抽出装置の機能構成例を示すブロック図。
【図３】図２中の時間相関抽出部１５、信号波形再生部１７の具体的機能構成例を示す図。
【図４】周期自己相関関数計算により求めた相関値分布の概念を示す図。
【図５】図２中の時間相関抽出部１５の他の具体的機能構成例を示す図。
【図６】Ａは時間相関波形の周波数特性を、Ｂは抽出したその希望信号成分を示す図である。
【図７】この発明の信号抽出装置の変形例の一部の機能構成を示すブロック図。
【図８】図７に示した装置の処理手順の例を示す流れ図。
【図９】この発明の信号抽出装置の他の機能構成例を示すブロック図。
【図１０】図９中の切替スイッチの各処理状態では切替え状態を示すブロック図。
【図１１】図９中の切替スイッチの各状態における制御手順の例を示す図。
【図１２】この発明方法が適用されるシステムの構成例を示す図。
【図１３】この発明方法手順の例を示す流れ図。
【図１４】図１２中の送信側装置１００の機能構成例を示す図。
【図１５】図１２中の受信側装置２００における送信パラメータ決定部２１４の機能構成例を示すブロック図。
【図１６】図１５中の使用可能領域検出部２１４ｆの動作を説明するための周期自己相関値分布例を示す図。
【図１７】送信パラメータ決定部２１４の他の機能構成例を示すブロック図。
【図１８】受信側装置２００の処理手順の例を示す図。
【図１９】受信側装置２００の具体的機能構成例を示すブロック図。
【図１】
2005151461.tif 000004
【図２】
2005151461.tif 000005
【図３】
2005151461.tif 000006
【図４】
2005151461.tif 000007
【図５】
2005151461.tif 000008
【図６】
2005151461.tif 000009
【図７】
2005151461.tif 000010
【図８】
2005151461.tif 000011
【図９】
2005151461.tif 000012
【図１０】
2005151461.tif 000013
【図１１】
2005151461.tif 000014
【図１２】
2005151461.tif 000015
【図１３】
2005151461.tif 000016
【図１４】
2005151461.tif 000017
【図１５】
2005151461.tif 000018
【図１６】
2005151461.tif 000019
【図１７】
2005151461.tif 000020
【図１８】
2005151461.tif 000021
【図１９】
2005151461.tif 000022