(57)【要約】【課題】一次傾斜の最適な補正を容易に行う。【解決手段】複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、入力信号の複素乗算を行って出力する、直列に接続された乗算器３１０−１〜３１０−ｎと、乗算器３１０−１〜３１０−ｎのうちの最後段の乗算器３１０−ｎからの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する判定器３２０と、乗算器３１０−ｎからの出力信号と誤差信号とに基づいて、複数の第２のタップ係数を算出する相関器３３０と、複数の第２のタップ係数の直交成分の総和を算出し、その総和が最小となる場合の傾斜補正値に属する複数の第１のタップ係数を乗算器３１０−１〜３１０−ｎに設定する制御部３４０とを有し、乗算器３１０−１〜３１０−ｎは、制御部３４０が設定した傾斜補正値に属する複数の第１のタップ係数を用いて入力信号の複素乗算を行う。【選択図】図１３
2018098597.tif 000002
【特許請求の範囲】
【請求項１】
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、入力信号の複素乗算を行って出力する、直列に接続された複数の乗算器と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する判定器と、
前記乗算器からの出力信号と、前記判定器が生成した誤差信号とに基づいて、複数の第２のタップ係数を算出する相関器と、
前記相関器が算出した複数の第２のタップ係数の直交成分の総和を算出し、該算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する制御部とを有し、
前記乗算器は、前記制御部が設定した傾斜補正値に属する複数の第１のタップ係数を用いて、前記入力信号の複素乗算を行う受信信号調整回路。
【請求項２】
請求項１に記載の受信信号調整回路において、
前記複数の乗算器は、サンプリング周波数を用いて動作する遅延素子を介して互いに直列に接続されている受信信号調整回路。
【請求項３】
請求項１または請求項２に記載の受信信号調整回路において、
前記制御部は、あらかじめ設定された複数のタイミングで相関器が算出した複数の第２のタップ係数の直交成分の平均値を算出し、該算出した平均値の総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する受信信号調整回路。
【請求項４】
請求項１から３のいずれか１項に記載の受信信号調整回路において、
前記判定器は、前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点との差分が最も少ない格子点を判定値として算出し、該判定値と前記乗算器からの出力信号との差分値を前記誤差信号として生成する受信信号調整回路。
【請求項５】
請求項１から４のいずれか１項に記載の受信信号調整回路において、
前記相関器は、前記乗算器からの出力信号を、直列に接続された複数の遅延素子を用いて遅延させたそれぞれの信号と、前記判定器が生成した誤差信号とに基づいて、前記複数の第２のタップ係数を算出する受信信号調整回路。
【請求項６】
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器が入力信号の複素乗算を行う処理と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する処理と、
前記乗算器からの出力信号と、前記誤差信号とに基づいて、複数の第２のタップ係数を算出する処理と、
前記複数の第２のタップ係数の直交成分の総和を算出する処理と、
前記算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する処理とを行う受信信号調整方法。
【請求項７】
コンピュータに、
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器で入力信号の複素乗算を行う手順と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する手順と、
前記乗算器からの出力信号と、前記誤差信号とに基づいて、複数の第２のタップ係数を算出する手順と、
前記複数の第２のタップ係数の直交成分の総和を算出する手順と、
前記算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する手順とを実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、受信信号調整回路、受信信号調整方法およびプログラムに関する。
【背景技術】
【０００２】
図１５は、一般的な屋外受信装置と屋内受信装置とから構成された通信システムの一例を示す図である。図１５に示した通信システムは、アンテナ３０００が受信した信号を屋外で受信する屋外受信装置２０００と、屋内に設置され、同軸ケーブル４０００を介して屋外受信装置２０００と接続された屋内受信装置１０００とから構成される。同軸ケーブル４０００は、屋外受信装置２０００から屋内受信装置１０００へＩＦ（Ｉｎｔｅｒｍｅｄｉａｔｅ  Ｆｒｅｑｕｅｎｃｙ）信号を受け渡す。
屋外受信装置２０００は、ＢＰＦ２１００と、受信周波数変換部２２００とを有する。ＢＰＦ２１００は、アンテナ３０００から受信したＲＦ（Ｒａｄｉｏ  Ｆｒｅｑｕｅｎｃｙ）信号を受けるフィルタである。受信周波数変換部２２００は、ＢＰＦ２１００を通過したＲＦ信号をＩＦ帯の信号へダウンコンバートする。
屋内受信装置１０００は、直交復調部１１００と、Ａ／Ｄ変換部１２００と、線形自動等化器１４００とを有する。直交復調部１１００は、同軸ケーブル４０００を介して受信したＩＦ信号を直交復調し、Ｉｃｈ成分とＱｃｈ成分とに分離する。Ａ／Ｄ変換部１２００は、直交復調部１１００が直交復調した信号をデジタル信号へ変換する。線形自動等化器１４００は、無線伝搬路で発生するフェージングに起因する符号間干渉を取り除く。線形自動等化器１４００は既知の技術であり、線形フィルタから構成され、受信した信号の判定値と誤差情報との相関を取り、そのフィルタタップ係数を適応的に制御する。
【０００３】
図１６は、図１５に示した直交復調部１１００に入力されるＩＦ帯のＱＡＭ（Ｑｕａｄｒａｔｕｒｅ  Ａｍｐｌｉｔｕｄｅ  Ｍｏｄｕｌａｔｉｏｎ）信号における周波数特性の一例を示す図である。図１６に示すように、ＱＡＭ信号の先頭部分が、同軸ケーブル４０００固有の一次傾斜特性により、右肩下がりになっている。この同軸ケーブル４０００に起因する一次傾斜は、線形自動等化器１４００でも等化することが可能である。
【０００４】
また、屋外装置からケーブルを介して送信されてきた信号を受信した屋内装置が、その信号の受信レベルを検出して屋外装置へフィードバックし、屋外装置がフィードバックされた情報に基づいて信号の増幅率を決定するシステムが開示されている（例えば、特許文献１参照。）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】国際公開第２００７／０６９４５０号
【発明の概要】
【発明が解決しようとする課題】
【０００６】
図１５に示したシステムにおいて、線形自動等化器１４００のリソースを一次傾斜の補正の処理に割いてしまうと、本来補償すべきフェージングに起因する符号間干渉を十分に除去することができないという問題が生じてしまう。これは、ＱＡＭ信号が多値になる程深刻になる。
また、特許文献１に記載されたシステムにおいては、屋外装置と屋内装置との間で制御信号を送受信する必要があり、その仕組みを構築しなければならないという問題点がある。
【０００７】
本発明の目的は、上記課題を解決する受信信号調整回路、受信信号調整方法およびプログラムを提供することにある。
【課題を解決するための手段】
【０００８】
本発明の受信信号調整回路は、
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、入力信号の複素乗算を行って出力する、直列に接続された複数の乗算器と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する判定器と、
前記乗算器からの出力信号と、前記判定器が生成した誤差信号とに基づいて、複数の第２のタップ係数を算出する相関器と、
前記相関器が算出した複数の第２のタップ係数の直交成分の総和を算出し、該算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する制御部とを有し、
前記乗算器は、前記制御部が設定した傾斜補正値に属する複数の第１のタップ係数を用いて、前記入力信号の複素乗算を行う。
また、本発明の受信信号調整方法は、
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器が入力信号の複素乗算を行う処理と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する処理と、
前記乗算器からの出力信号と、前記誤差信号とに基づいて、複数の第２のタップ係数を算出する処理と、
前記複数の第２のタップ係数の直交成分の総和を算出する処理と、
前記算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する処理とを行う。
また、本発明のプログラムは、
コンピュータに実行させるためのプログラムであって、
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器で入力信号の複素乗算を行う手順と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する手順と、
前記乗算器からの出力信号と、前記誤差信号とに基づいて、複数の第２のタップ係数を算出する手順と、
前記複数の第２のタップ係数の直交成分の総和を算出する手順と、
前記算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する手順とを実行させる。
【発明の効果】
【０００９】
以上説明したように、本発明においては、一次傾斜の最適な補正を容易に行うことができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の受信信号調整回路を適用した通信システムの第１の実施の形態を示す図である。
【図２】図１に示した複素線形フィルタにおける周波数ｆに対する振幅｜Ｈ（ｆ）｜および位相ψＨ（ｆ）の特性の一例を示した図である。
【図３】図１に示した複素線形フィルタの内部構成の一例を示す図である。
【図４】図１に示したメモリにあらかじめ格納された第１のタップ係数の一例を示す図である。
【図５】図１に示した複素線形フィルタにおいて、傾斜補正値を「９」とした場合の周波数ｆに対する振幅｜Ｈ（ｆ）｜および位相ψＨ（ｆ）の特性の一例を示した図である。
【図６】図１に示した線形自動等化器の内部構成の一例を示す図である。
【図７】図６に示した適応フィルタの内部構成の一例を示す図である。
【図８】図１に示した線形自動等化器における、一次傾斜が無い場合のタップ係数の一例を示す図である。
【図９】図１に示した線形自動等化器における、一次傾斜がある場合のタップ係数の一例を示す図である。
【図１０】図１に示した通信システムにおける受信信号調整方法を説明するためのフローチャートである。
【図１１】本発明の受信信号調整回路を適用した通信システムの第２の実施の形態を示す図である。
【図１２】本発明の受信信号調整回路を適用した通信システムの第３の実施の形態を示す図である。
【図１３】本発明の受信信号調整回路の第４の実施の形態を示す図である。
【図１４】図１３に示した受信信号調整回路における受信信号調整方法を説明するためのフローチャートである。
【図１５】一般的な屋外受信装置と屋内受信装置とから構成された通信システムの一例を示す図である。
【図１６】図１５に示した直交復調部に入力されるＩＦ帯のＱＡＭ信号における周波数特性の一例を示す図である。
【発明を実施するための形態】
【００１１】
本発明は、受信部に一次傾斜補正に特化した複素線形フィルタを設けた構成をとる。この複素線形フィルタが、屋外装置と屋内装置とを接続する同軸ケーブルに起因した一次傾斜を補正し、線形自動等化器が、無線伝搬路で生じるフェージングに起因する符号間干渉を補償することで、適切なリソース配分が可能となる。
（第１の実施の形態）
【００１２】
本形態における一次傾斜の自動補正は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）における制御で行う。ＣＰＵと接続された不揮発性の外部メモリに、複素線形フィルタのタップ係数候補をテーブルとしてあらかじめ格納しておく。ＣＰＵが外部メモリから、タップ係数を逐次読み出して複素線形フィルタに設定することで一次傾斜の補正を行いつつ、一次傾斜の大きさを推定する。そして、ＣＰＵが、その推定した値が最も小さくなるようなタップ係数を、外部メモリに格納されたタップ係数候補から選び出すことで最適な補正を実現する。一次傾斜量の推定は、同軸ケーブルに依存する一次傾斜の大きさに応じて、受信部の線形自動等化器におけるタップ係数の直交成分が成長する性質を利用して行う。
この一次傾斜の自動調整は、装置立ち上げ直後に所定のトレーニング期間を設けて実施する。このため、トレーニング期間中に限り、線形自動等化器を一次傾斜の推定に充てる。一度調整が完了すれば、同軸ケーブルを取り換えない限り、装置を立ち上げる毎に再調整する必要性は無い。トレーニング完了後は、装置の動作モードを実運用モードに切り替え、複素線形フィルタが、トレーニング時に選出した最適なタップ係数を用いて傾斜補正を行う。これにより、実運用時は、線形自動等化器が同軸ケーブルの一次傾斜補正のために無駄にタップ係数を成長させることがなくなる。そのため、線形自動等化器がフェージングに起因する符号間干渉の補償に十分機能を発揮することができる。
【００１３】
図１は、本発明の受信信号調整回路を適用した通信システムの第１の実施の形態を示す図である。本形態は図１に示すように、アンテナ３００が受信した信号を屋外で受信する屋外受信装置２００と、屋内に設置され、同軸ケーブル４００を介して屋外受信装置２００と接続された屋内受信装置１００とから構成される。同軸ケーブル４００は、屋外受信装置２００から屋内受信装置１００へＩＦ信号を受け渡す。
屋外受信装置２００は、ＢＰＦ２１０と、受信周波数変換部２２０とを有する。ＢＰＦ２１０は、アンテナ３００から受信したＲＦ信号を受けるフィルタである。受信周波数変換部２２０は、ＢＰＦ２１０を通過したＲＦ信号をＩＦ帯の信号へダウンコンバートする。
屋内受信装置１００は、直交復調部１１０と、Ａ／Ｄ変換部１２０と、複素線形フィルタ１３０と、線形自動等化器１４０と、メモリ１５０と、ＣＰＵ１６０とを有する。この複素線形フィルタ１３０と、線形自動等化器１４０と、メモリ１５０と、ＣＰＵ１６０とで本発明の受信信号調整回路を構成する。
直交復調部１１０は、同軸ケーブル４００を介して受信したＩＦ信号を直交復調し、Ｉｃｈ成分とＱｃｈ成分とに分離する。Ａ／Ｄ変換部１２０は、直交復調部１１０が直交復調した信号をデジタル信号へ変換する。複素線形フィルタ１３０は、ＣＰＵ１６０が設定したタップ係数を用いて、同軸ケーブル４００の一次傾斜を補正する。詳細は後述する。線形自動等化器１４０は、無線伝搬路で発生するフェージングに起因する符号間干渉を取り除く。詳細は後述する。メモリ１５０は、第１のタップ係数の候補をテーブルとしてあらかじめ格納しておく不揮発性の外部メモリである。ＣＰＵ１６０は、傾斜補正値を自律で変更し、その補正値に応じたタップ係数をメモリ１５０から読み出して複素線形フィルタ１３０に設定する機能を有する制御部である。さらに、ＣＰＵ１６０は、線形自動等化器１４０からタップ係数を読み出して一次傾斜の大きさを推定する機能および推定値に基づいてタップ係数の最適値を選択する機能を有する。
【００１４】
図２は、図１に示した複素線形フィルタ１３０における周波数ｆに対する振幅｜Ｈ（ｆ）｜および位相ψＨ（ｆ）の特性の一例を示した図である。ここで、複素線形フィルタ１３０は、サンプリング周波数Ｆｓで動作するデジタルフィルタを想定する。Ｈ（ｆ）は複素関数であり、｜Ｈ（ｆ）｜は関数Ｈ（ｆ）の絶対値、ψＨ（ｆ）は、Ｈ（ｆ）の実数部をＲｅ｛Ｈ（ｆ）｝、虚数部をＩｍ｛Ｈ（ｆ）｝とすると、
ψＨ（ｆ）＝ｔａｎ−１Ｉｍ｛Ｈ（ｆ）｝／Ｒｅ｛Ｈ（ｆ）｝
で表現される。図２に示すように、｜Ｈ（ｆ）｜は、複素線形フィルタ１３０で表現可能な周波数範囲内（−Ｆｓ／４〜Ｆｓ／４）において右上がり（周波数増加に伴い振幅が増加）の一次傾斜特性を有する。同軸ケーブル４００の一次傾斜特性は右下がり（周波数増加に伴い振幅が減少）を示すため、ケーブル固有の傾斜特性を打ち消す特性となる。Ｋ［ｄＢｐｐ］は、上述した表現領域（−Ｆｓ／４〜Ｆｓ／４）における振幅の振れ幅（以下、傾斜補正値と称する）を示す。ψＨ（ｆ）は、位相の補正は実施されないため、上述した表現領域（−Ｆｓ／４〜Ｆｓ／４）において周波数に依存せず一定の値となる。複素線形フィルタ１３０のタップ係数は、複素関数Ｈ（ｆ）の逆フーリエ変換により生成される。生成されたタップ係数も複素数である。
【００１５】
図３は、図１に示した複素線形フィルタ１３０の内部構成の一例を示す図である。図１に示した複素線形フィルタ１３０は図３に示すように、サンプリング周波数を用いて動作する複数のＤフリップフロップ（遅延素子）１３２−１〜１３２−１３をそれぞれ介して直列に接続された複数の複素乗算器１３１−１〜１３１−１３を有するデジタルＦＩＲフィルタである。複素乗算器１３１−１〜１３１−１３は、センタータップに加え、その前後６タップから成り、その動作は、サンプリング周波数Ｆｓでの動作を想定している。図３に示したＲｅ＿ｊ（Ｋ）はタップ係数の実数部を表し、Ｉｍ＿ｊ（Ｋ）はタップ係数の虚数部を表す。また、Ｋは傾斜補正値を表し、ｊはセンタータップ位置を基準（ｊ＝０）としたときのそれぞれのフィルタタップ位置を表す。図３に示した構成では、ｊの取り得る値は−６〜＋６の整数であり、正の値はセンタータップに対して時間的に遅れた位置を示し、負の値はセンタータップに対して時間的に進んだ位置を示す。
【００１６】
図４は、図１に示したメモリ１５０にあらかじめ格納された第１のタップ係数の一例を示す図である。図４に示すように、図１に示したメモリ１５０には、０〜９［ｄＢｐｐ］の１ｄＢ刻みの傾斜補正値に対応した複素線形フィルタ１３０のタップ係数がテーブルとして格納されている。
【００１７】
図５は、図１に示した複素線形フィルタ１３０において、傾斜補正値を「９」とした場合の周波数ｆに対する振幅｜Ｈ（ｆ）｜および位相ψＨ（ｆ）の特性の一例を示した図である。図４に示したタップ係数において、例えば、傾斜補正値「９」に対応したタップ係数はＲｅ＿ｊ（９）およびＩｍ＿ｊ（９）であり、その係数は図５に示した複素関数Ｈ（ｆ）の逆フーリエ変換で導出可能である。
傾斜補正値が「０」である場合は、Ｈ（ｆ）の振幅および位相特性が表現周波数領域において一定となる。このときのＨ（ｆ）を逆フーリエ変換した値であるタップ係数は、実数部のセンタータップＲｅ＿０（０）＝１を除く全てのタップ係数が値「０」となり、複素線形フィルタ１３０をスルーした場合と等価になる。
【００１８】
図６は、図１に示した線形自動等化器１４０の内部構成の一例を示す図である。図１に示した線形自動等化器１４０は図６に示すように、適応フィルタ１４１ａ，１４１ｂと、判定器１４４ａ，１４４ｂと、相関器１４５ａ，１４５ｂとを有する。
適応フィルタ１４１ａは、複素線形フィルタ１３０から出力された信号のＩｃｈの信号およびＱｃｈの信号を入力として、符号間干渉の影響を除去したＩｃｈ信号を出力する。適応フィルタ１４１ｂは、複素線形フィルタ１３０から出力された信号のＩｃｈの信号およびＱｃｈの信号を入力として、符号間干渉の影響を除去したＱｃｈ信号を出力する。
判定器１４４ａは、適応フィルタ１４１ａから出力されるＩｃｈ信号を入力とし、このＩｃｈ信号と、ＱＡＭ（直角位相振幅変調）の理想格子点（例えば、４０９６ＱＡＭの場合、４０９６通りの格子点が存在する）上のＩｃｈ信号との差分をそれぞれ算出し、算出した差分が最も小さな格子点を判定値として検出する。また、判定器１４４ａは、Ｉｃｈ信号と判定値との差分をＩｃｈ誤差信号として生成し、相関器１４５ｂへ出力する。また、判定器１４４ｂは、適応フィルタ１４１ｂから出力されるＱｃｈ信号を入力とし、判定器１４４ａと同様の処理でＱｃｈ誤差信号を生成し、相関器１４５ａへ出力する。
相関器１４５ｂは、適応フィルタ１４１ａの各タップから得られるＩｃｈ信号の極性信号（Ｉｃｈ信号の正負をバイナリで表現したもの）とＩｃｈ誤差信号とで相関をとることにより、第２のタップ係数の同相成分を生成し、適応フィルタ１４１ａへ出力する。また、相関器１４５ａは、適応フィルタ１４１ａの各タップから得られるＱｃｈ信号の極性信号（Ｑｃｈ信号の正負をバイナリで表現したもの）とＱｃｈ誤差信号とで相関をとることにより、第２のタップ係数の直交成分を生成し、適応フィルタ１４１ａへ出力する。
【００１９】
ここで、Ｉｃｈ信号の極性信号ＤｉおよびＱｃｈ信号の極性信号Ｄｑの定義を、信号値の極性が正の場合は「１」、信号値の極性が負の場合は「０」のバイナリ表現とする。また、Ｉｃｈ誤差信号ＥｉおよびＱｃｈ誤差信号Ｅｑの定義を、適応フィルタ出力が判定値よりも大きな場合、「１」、小さな場合は「０」のバイナリ表現とする。すると、以下のように相関値が表現される。
相関値＝ｎｏｔ（Ｄｉ・Ｅｉ）
相関値＝ｎｏｔ（Ｄｑ・Ｅｑ）
線形自動等化器１４０は、相関値が「０」である場合、第２のタップ係数を減少させ、相関値が「１」である場合は、第２のタップ係数を増加させることにより、符号間干渉を打ち消す方向にタップ係数を成長させていく。
【００２０】
図７は、図６に示した適応フィルタ１４１ａの内部構成の一例を示す図である。図６に示した適応フィルタ１４１ａは図７に示すように、サンプリング周期Ｔｓで動作するＤフリップフロップ（遅延素子）１４１１ａ〜１４１１ｌと、各Ｄフリップフロップ出力にタップ係数Ｃｒｋ，Ｃｉｋ（図７に示した構成においては、ｋ＝−６〜＋６の整数）を乗じるための乗算器１４２１ａ〜１４２１ｈと、各乗算器１４２１ａ〜１４２１ｈの出力の総和を求めるための加算器１４３１，１４３１ａ，１４３１ｂとを有する。ここで、タップ係数Ｃｒｋは、線形自動等化器係数Ｃｋ（図７に示した構成においては、ｋ＝−６〜＋６の整数）の実数成分Ｃｒｋ＝Ｒｅ｛Ｃｋ｝に相当し、本係数をタップ係数の同相成分と表現する。また、タップ係数Ｃｉｋは、線形自動等化器係数Ｃｋの虚数成分Ｃｉｋ＝Ｉｍ｛Ｃｋ｝に相当し、本係数をタップ係数の直交成分と表現する。なお、図７には、Ｄフリップフロップの一部および乗算器の一部を省略して示しているが、ｋの値に応じた数のＤフリップフロップおよび乗算器が設けられている。
【００２１】
図７に示した形態において、Ｄフリップフロップ１４１１ｃの出力信号がセンタータップのＩｃｈ信号に相当する。このＩｃｈ信号は加算器１４３１ａに直接接続されるため、本形態ではセンタータップのタップ係数Ｃｒ０は常時「１」となる。また、Ｄフリップフロップ１４１１ｉの出力信号がセンタータップのＱｃｈ信号に相当する。このＱｃｈ信号は未使用であるため、本形態ではセンタータップのタップ係数Ｃｉ０は常時「０」となる。
【００２２】
また、図７に示した形態において、Ｄフリップフロップ１４１１ｂ，１４１１ｈの出力信号はそれぞれ、センタータップに対してＴｓ遅れたＩｃｈ信号、Ｑｃｈ信号を表し、それぞれ相関器１４５ｂ，１４５ａに入力される。相関器１４５ｂは、入力したＩｃｈ信号とＩｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｒ−１の値を初期値「０」から増減することで自動調整を行う。相関器１４５ａは、Ｑｃｈ信号とＱｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｉ−１の値を初期値「０」から増減することで自動調整を行う。
【００２３】
また、図７に示した形態において、Ｄフリップフロップ１４１１ａ，１４１１ｇの出力信号はそれぞれ、センタータップに対して６Ｔｓ遅れたＩｃｈ信号、Ｑｃｈ信号を表し、それぞれ相関器１４５ｂ，１４５ａに入力される。相関器１４５ｂは、入力したＩｃｈ信号とＩｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｒ−６の値を初期値「０」から増減することで自動調整を行う。相関器１４５ａは、Ｑｃｈ信号とＱｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｉ−６の値を初期値「０」から増減することで自動調整を行う。
【００２４】
また、図７に示した形態において、Ｄフリップフロップ１４１１ｄ，１４１１ｊの出力信号はそれぞれ、センタータップに対してＴｓ進んだＩｃｈ信号、Ｑｃｈ信号を表し、それぞれ相関器１４５ｂ，１４５ａに入力される。相関器１４５ｂは、入力したＩｃｈ信号とＩｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｒ−１の値を初期値「０」から増減することで自動調整を行う。相関器１４５ａは、Ｑｃｈ信号とＱｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｉ−１の値を初期値「０」から増減することで自動調整を行う。
【００２５】
また、図７に示した形態において、Ｄフリップフロップ１４１１ｆ，１４１１ｌの出力信号はそれぞれ、センタータップに対して６Ｔｓ進んだＩｃｈ信号、Ｑｃｈ信号を表し、それぞれ相関器１４５ｂ，１４５ａに入力される。相関器１４５ｂは、入力したＩｃｈ信号とＩｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｒ−６の値を初期値「０」から増減することで自動調整を行う。相関器１４５ａは、Ｑｃｈ信号とＱｃｈ誤差信号とを用いて相関値を算出し、算出した相関値に応じてタップ係数Ｃｉ−６の値を初期値「０」から増減することで自動調整を行う。
【００２６】
適応フィルタ１４１ａの出力は、適応フィルタ１４１ａに入力されるＩｃｈ信号から、この信号に対して時間的に±６Ｔｓ、±５Ｔｓ、±４Ｔｓ、±３Ｔｓ、±２Ｔｓ、±Ｔｓ前後したＩｃｈ信号、Ｑｃｈ信号が及ぼす符号間干渉を除去したものに相当する。
【００２７】
図８は、図１に示した線形自動等化器１４０における、一次傾斜が無い場合のタップ係数の一例を示す図である。図９は、図１に示した線形自動等化器１４０における、一次傾斜がある場合のタップ係数の一例を示す図である。図８および図９においては、タップ係数Ｃｋのｋの値を−５〜＋５の整数としている。
図８および図９に示すように、対象の等化器タップ係数Ｃｋは、センタータップの係数Ｃ０を中心に前後５タップ分存在する。Ｃｋは複素数であり、実数部がタップ係数の同相成分、虚数部が直交成分を表すものとする。図８に示すように、同軸ケーブル等に起因する一次傾斜が無く、無線伝搬路が良好な場合、等化器タップは成長しない（タップ係数のデフォルト値は、センタータップＣ０の同相成分に限り「１」であり、その他は「０」である）。また、図９に示すように、同軸ケーブル等に起因する一次傾斜がある場合、線形自動等化器１４０は傾斜補正すべくタップ係数を成長させる。図１６を用いて説明したように、一次傾斜は信号の中心周波数に対して左右非対称の周波数特性であるため、タップ係数の同相成分ではなく直交成分が成長する。また、本タップ係数値は、一次傾斜の度合い（傾きの急峻さ）に比例して増加する。これらの性質を利用すると、タップ係数の直交成分の総和を計算することにより、一次傾斜の影響度合いを推定することが可能となる。ただし、センタータップの値に限り、一次傾斜の影響と無相関であるため、この計算対象からは外す。
【００２８】
以下に、図１に示した通信システムにおける受信信号調整方法について説明する。図１０は、図１に示した通信システムにおける受信信号調整方法を説明するためのフローチャートである。図１に示した屋内受信装置１００は、同軸ケーブル４００の受信信号調整（一次傾斜補正）を所定のトレーニング期間において実施する。送信側である図１に示した屋外受信装置２００における一次傾斜補正は理想的に実施されているものとする。
【００２９】
まず、ＣＰＵ１６０が、メモリ１５０にあらかじめ格納されているテーブルから、傾斜補正値「０」に対応付けられた第１のタップ係数Ｒｅ＿ｊ（０）、Ｉｍ＿ｊ（０）を読み出す。ＣＰＵ１６０は、メモリ１５０から読み出した各タップ係数Ｒｅ＿ｊ（０）、Ｉｍ＿ｊ（０）を複素線形フィルタ１３０の複素乗算器１３１−１〜１３１−１３に設定する。具体的には、図３に示すように、ＣＰＵ１６０は、タップ係数Ｒｅ＿−６（０）、Ｉｍ＿−６（０）を図３に示した複素乗算器１３１−１に設定し、タップ係数Ｒｅ＿−５（０）、Ｉｍ＿−５（０）を図３に示した複素乗算器１３１−２に設定するといった、それぞれのタップ係数をそれに応じた複素乗算器に設定する。タップ係数が設定された複素乗算器１３１−１〜１３１−１３が入力信号Ｉｃｈ、Ｑｃｈの複素乗算を行う（ステップＳ１）。
続いて、線形自動等化器１４０の適応フィルタ１４１ａが、複素線形フィルタ１３０から出力されてきたＩｃｈ信号およびＱｃｈ信号から符号化干渉の影響を除去したＩｃｈ信号を出力する。すると、判定器１４４ａが、適応フィルタ１４１ａから出力されたＩｃｈ信号と、ＱＡＭの理想格子点上のＩｃｈ信号との差分をそれぞれ算出する（ステップＳ２）。判定器１４４ａは、差分が最小となる格子点を判定値として検出する（ステップＳ３）。そして、判定器１４４ａは、Ｉｃｈ信号と判定値との差分をＩｃｈ誤差信号として算出し（ステップＳ４）、相関器１４５ｂへ出力する。
同様に、線形自動等化器１４０の適応フィルタ１４１ｂが、複素線形フィルタ１３０から出力されてきたＱｃｈ信号およびＩｃｈ信号から符号化干渉の影響を除去したＱｃｈ信号を出力する。すると、判定器１４４ｂが、適応フィルタ１４１ｂから出力されたＱｃｈ信号と、ＱＡＭの理想格子点上のＱｃｈ信号との差分をそれぞれ算出する。判定器１４４ｂは、差分が最小となる格子点を判定値として検出する。そして、判定器１４４ｂは、Ｑｃｈ信号と判定値との差分をＱｃｈ誤差信号として算出し、相関器１４５ａへ出力する。
相関器１４５ａは、適応フィルタ１４１ａから出力されたＱｃｈの極性信号と、判定器１４４ｂから出力されたＱｃｈ誤差信号とに基づいて、複数の第２のタップ係数を算出する（ステップＳ５）。具体的には、図７を用いて説明した通りである。
【００３０】
ＣＰＵ１６０は、線形自動等化器１４０から第２のタップ係数Ｃｋ（ｋは−６〜＋６の０を除く整数）の直交成分Ｉｍ｛Ｃｋ｝を取り出し、所定の演算回数Ｎ（Ｎは自然数）において、タップ係数毎の移動平均値を算出する（ステップＳ６）。このとき、ＣＰＵ１６０は、センタータップＩｍ｛Ｃ０｝の情報に限り、一次傾斜の推定に用いることは出来ず、ノイズ源になるだけなので、演算対象に含めない。
時刻ｔにおける演算結果Ｉｍ｛ａｖｅ＿Ｃｋ｝ｔは、以下に示す式で表現される。
Ｉｍ｛ａｖｅ＿Ｃｋ｝ｔ＝｛Ｉｍ｛Ｃｋ｝ｔ＋Ｉｍ｛Ｃｋ｝ｔ−１＋・・・＋Ｉｍ｛Ｃｋ｝ｔ−Ｎ＋１｝／Ｎ
この各タップ係数の直交成分平均化は、ノイズ除去を目的として行う。ＣＰＵ１６０は、平均値の演算処理が完了した後、平均値の各タップ係数の総和を算出する（ステップＳ７）。算出式を以下に示す。
γＫ＝｛Ｉｍ｛ａｖｅ＿Ｃ−６｝＋Ｉｍ｛ａｖｅ＿Ｃ−５｝＋・・・＋Ｉｍ｛ａｖｅ＿Ｃ＋５｝＋Ｉｍ｛ａｖｅ＿Ｃ＋６｝
ここで、「γＫ（Ｋは０〜９の整数）」は、傾斜補正値Ｋ［ｄＢｐｐ］に対応するタップ係数Ｒｅ＿ｊ（Ｋ）、Ｉｍ＿ｊ（Ｋ）を用いて、複素線形フィルタ１３０で補正したときの、一次傾斜推定値（補正後）に相当する。
現時点では、ＣＰＵ１６０が傾斜補正値「０」を自律で設定しているため、演算結果は「γ０」と表現される。この「γ０」は、複素線形フィルタ１３０が傾斜補正値「０」で補正を行った時の一次傾斜推定値に相当する。
【００３１】
ＣＰＵ１６０は、傾斜補正値「０」と、それに対応する傾斜推定値「γ０」を自身のレジスタに格納する。その後、ＣＰＵ１６０は、傾斜補正値を１ｄＢインクリメントさせ、メモリ１５０に格納されたテーブルから、傾斜補正値「１」に対応したタップ係数Ｒｅ＿ｊ（１）、Ｉｍ＿ｊ（１）を読み出す。そして、ＣＰＵ１６０は、タップ係数Ｒｅ＿ｊ（１）、Ｉｍ＿ｊ（１）を複素線形フィルタ１３０に設定し、上述した補正を行う。
以下同様の手順で、ＣＰＵ１６０は、傾斜補正値を０［ｄＢｐｐ］から９［ｄＢｐｐ］まで、１ｄＢずつ増やしながら、傾斜補正値Ｋ［ｄＢｐｐ］に応じた傾斜推定量「γＫ」［真値］を自身のレジスタ領域に格納していく。
その後、以下に示す式で定義される「γＫ」の最小値「γＫｍｉｎ」を与える傾斜補正値「Ｋｍｉｎ」を検出する。
γＫｍｉｎ＝ｍｉｎｉｍｕｍ（γ０，・・・，γＫ，・・・，γ９）
以上の処理を経て、同軸ケーブル４００に起因する一次傾斜を自動補正するためのトレーニング過程は完了し、屋内受信装置１００は動作モードを実運用の動作モードに切り替える。実運用時は、ＣＰＵ１６０が、メモリ１５０に格納されたテーブルから、傾斜補正量「Ｋｍｉｎ」に対応したタップ係数Ｒｅ＿ｊ（Ｋｍｉｎ）、Ｉｍ＿ｊ（Ｋｍｉｎ）を読み出し、読み出したタップ係数を複素線形フィルタ１３０に設定し（ステップＳ８）、傾斜補正を実施する。
【００３２】
本形態においては、−６〜＋６の整数をｊやｋの値の一例として説明したが、他の整数を用いるものであっても良い。ただし、タップ係数の数が少なくなると、そのタップ係数を保有するデジタルフィルタが表現できる周波数特性の精度が落ちるため、あまり小さな整数を用いることは好ましくない。
【００３３】
以上説明したように、複素線形フィルタ１３０に複数の第１のタップ係数を設定し、複素線形フィルタ１３０の出力信号と直角位相振幅変調の理想格子点とに基づいた第２のタップ係数の直交成分の総和の最小値を複素線形フィルタ１３０に設定する第１のタップ係数とする。これにより、屋内受信装置１００と屋外受信装置２００とを接続する同軸ケーブル４００の材質、または、長さによりケーブル固有の一次傾斜特性が異なった場合であっても、一次傾斜を最適に自動補正することが可能となる。また、フェージングの補償に用いる為の自動線形等化器の性能を劣化させることなく、一次傾斜の補正が可能となる。
（第２の実施の形態）
【００３４】
図１１は、本発明の受信信号調整回路を適用した通信システムの第２の実施の形態を示す図である。本形態においては、図１１に示すように、屋内受信装置１０１にて、第１の実施の形態におけるＣＰＵ１６０とメモリ１５０とを周辺回路（不揮発性メモリを含む）を内蔵したカスタムＩＣ１７１で実現する。これにより、１個のデバイス（ＩＣ：Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）でＣＰＵ１６０とメモリ１５０との機能を実現することも可能である。
（第３の実施の形態）
【００３５】
図１２は、本発明の受信信号調整回路を適用した通信システムの第３の実施の形態を示す図である。本形態においては、図１２に示すように、屋内受信装置１０２がＣ／Ｎ推定器１７２を有する。Ｃ／Ｎ推定器１７２は、線形自動等化器１４０から出力される誤差信号を利用して、Ｃ／Ｎ（信号電力対雑音電力の比）推定値を算出し、算出した結果をＣＰＵ１６２へ出力する。ＣＰＵ１６２は、Ｃ／Ｎ推定器１７２から出力されてきたＣ／Ｎ推定結果が所定の閾値よりも高い場合は一次傾斜の推定を行い、所定の閾値よりも低い場合は一次傾斜の推定を取りやめる制御を行う。これにより、無線伝搬路が良好な状態に限り一次傾斜の推定が可能となり、一次傾斜の推定精度を高めることができる効果がある。
（第４の実施の形態）
【００３６】
図１３は、本発明の受信信号調整回路の第４の実施の形態を示す図である。本形態における受信信号調整回路は図１３に示すように、直列に接続された複数の乗算器３１０−１〜３１０−ｎ（ｎは自然数）と、判定器３２０と、相関器３３０と、制御部３４０とを有する。
乗算器３１０−１〜３１０−ｎは、あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、入力信号の複素乗算を行って出力する。乗算器３１０−１〜３１０−ｎは、制御部３４０が設定した傾斜補正値に属する複数の第１のタップ係数を用いて、入力信号の複素乗算を行う。
判定器３２０は、乗算器３１０−１〜３１０−ｎのうちの最後段の乗算器３１０−ｎからの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する。
相関器３３０は、乗算器３１０−ｎからの出力信号と、判定器３２０が生成した誤差信号とに基づいて、複数の第２のタップ係数を算出する。
制御部３４０は、相関器３３０が算出した複数の第２のタップ係数の直交成分の総和を算出する。また、制御部３４０は、算出した総和が最小となる場合の傾斜補正値に属する複数の第１のタップ係数を乗算器３１０−１〜３１０−ｎに設定する。
【００３７】
以下に、図１３に示した受信信号調整回路における受信信号調整方法について説明する。図１４は、図１３に示した受信信号調整回路における受信信号調整方法を説明するためのフローチャートである。
【００３８】
乗算器３１０−１〜３１０−ｎが、あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器で入力信号の複素乗算を行う（ステップＳ１１）。続いて、判定器３２０が、乗算器３１０−１〜３１０−ｎのうちの最後段の乗算器３１０−ｎからの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する（ステップＳ１２）。すると、相関器３３０が、乗算器３１０−ｎからの出力信号と、判定器３２０が生成した誤差信号とに基づいて、複数の第２のタップ係数を算出する（ステップＳ１３）。制御部３４０が、相関器３３０が算出した複数の第２のタップ係数の直交成分の総和を算出する（ステップＳ１４）。また、制御部３４０は、算出した総和が最小となる場合の傾斜補正値に属する複数の第１のタップ係数を乗算器３１０−１〜３１０−ｎに設定する（ステップＳ１５）。
【００３９】
このように、乗算器３１０−１〜３１０−ｎに複数の第１のタップ係数を設定し、乗算器３１０−ｎの出力信号と直角位相振幅変調の理想格子点とに基づいた第２のタップ係数の直交成分の総和の最小値を乗算器３１０−１〜３１０−ｎに設定する第１のタップ係数とする。これにより、一次傾斜の最適な補正を容易に行うことができる。
【００４０】
上述した受信信号調整回路に設けられた各構成要素が行う処理は、目的に応じてそれぞれ作製された論理回路で行うようにしても良い。また、処理内容を手順として記述したコンピュータプログラム（以下、プログラムと称する）を受信信号調整回路にて読取可能な記録媒体に記録し、この記録媒体に記録されたプログラムを受信信号調整回路に読み込ませ、実行するものであっても良い。受信信号調整回路にて読取可能な記録媒体とは、フロッピー（登録商標）ディスク、光磁気ディスク、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｃ）、ＣＤ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ）、Ｂｌｕ−ｒａｙ（登録商標）  Ｄｉｓｃなどの移設可能な記録媒体の他、受信信号調整回路に内蔵されたＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）等のメモリやＨＤＤ（Ｈａｒｄ  Ｄｉｓｃ  Ｄｒｉｖｅ）等を指す。この記録媒体に記録されたプログラムは、受信信号調整回路に設けられたＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）にて読み込まれ、ＣＰＵの制御によって、上述したものと同様の処理が行われる。ここで、ＣＰＵは、プログラムが記録された記録媒体から読み込まれたプログラムを実行するコンピュータとして動作するものである。
【００４１】
上記の実施の形態の一部または全部は、以下の付記のようにも記載され得るが、以下には限られない。
（付記１）あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、入力信号の複素乗算を行って出力する、直列に接続された複数の乗算器と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する判定器と、
前記乗算器からの出力信号と、前記判定器が生成した誤差信号とに基づいて、複数の第２のタップ係数を算出する相関器と、
前記相関器が算出した複数の第２のタップ係数の直交成分の総和を算出し、該算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する制御部とを有し、
前記乗算器は、前記制御部が設定した傾斜補正値に属する複数の第１のタップ係数を用いて、前記入力信号の複素乗算を行う受信信号調整回路。
（付記２）前記複数の乗算器は、サンプリング周波数を用いて動作する遅延素子を介して互いに直列に接続されている、付記１に記載の受信信号調整回路。
（付記３）前記制御部は、あらかじめ設定された複数のタイミングで相関器が算出した複数の第２のタップ係数の直交成分の平均値を算出し、該算出した平均値の総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する、付記１または付記２に記載の受信信号調整回路。
（付記４）前記判定器は、前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点との差分が最も少ない格子点を判定値として算出し、該判定値と前記乗算器からの出力信号との差分値を前記誤差信号として生成する、付記１から３のいずれか１項に記載の受信信号調整回路。
（付記５）前記相関器は、前記乗算器からの出力信号を、直列に接続された複数の遅延素子を用いて遅延させたそれぞれの信号と、前記判定器が生成した誤差信号とに基づいて、前記複数の第２のタップ係数を算出する、付記１から４のいずれか１項に記載の受信信号調整回路。
（付記６）あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器が入力信号の複素乗算を行う処理と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する処理と、
前記乗算器からの出力信号と、前記誤差信号とに基づいて、複数の第２のタップ係数を算出する処理と、
前記複数の第２のタップ係数の直交成分の総和を算出する処理と、
前記算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する処理とを行う受信信号調整方法。
（付記７）コンピュータに、
あらかじめ設定された複数の傾斜補正値に属する複数の第１のタップ係数をそれぞれ用いて、直列に接続された複数の乗算器で入力信号の複素乗算を行う手順と、
前記複数の乗算器のうちの最後段の乗算器からの出力信号と直角位相振幅変調の理想格子点とに基づいて、誤差信号を生成する手順と、
前記乗算器からの出力信号と、前記誤差信号とに基づいて、複数の第２のタップ係数を算出する手順と、
前記複数の第２のタップ係数の直交成分の総和を算出する手順と、
前記算出した総和が最小となる場合の前記傾斜補正値に属する複数の第１のタップ係数を前記乗算器に設定する手順とを実行させるためのプログラム。
【符号の説明】
【００４２】
１００〜１０２    屋内受信装置
１１０    直交復調部
１２０    Ａ／Ｄ変換部
１３０    複素線形フィルタ
１３１−１〜１３１−１３    複素乗算器
１３２−１〜１３２−１３，１４１１ａ〜１４１１ｌ    Ｄフリップフロップ
１４０    線形自動等化器
１４１ａ，１４１ｂ    適応フィルタ
１４４ａ，１４４ｂ，３２０  判定器
１４５ａ，１４５ｂ，３３０    相関器
１４２１ａ〜１４２１ｈ，３１０−１〜３１０−ｎ    乗算器
１４３１，１４３１ａ，１４３１ｂ    加算器
１５０    メモリ
１６０，１６２    ＣＰＵ
１７１    カスタムＩＣ
１７２    Ｃ／Ｎ推定器
２００    屋外受信装置
２１０    ＢＰＦ
２２０    受信周波数変換部
３００    アンテナ
３４０    制御部
４００    同軸ケーブル
【図１】
2018098597.tif 000003
【図２】
2018098597.tif 000004
【図３】
2018098597.tif 000005
【図４】
2018098597.tif 000006
【図５】
2018098597.tif 000007
【図６】
2018098597.tif 000008
【図７】
2018098597.tif 000009
【図８】
2018098597.tif 000010
【図９】
2018098597.tif 000011
【図１０】
2018098597.tif 000012
【図１１】
2018098597.tif 000013
【図１２】
2018098597.tif 000014
【図１３】
2018098597.tif 000015
【図１４】
2018098597.tif 000016
【図１５】
2018098597.tif 000017
【図１６】
2018098597.tif 000018