(57)【要約】  重複多重に基づく一種類の変調復調方法と装置、変調方法では、プリコーディング構成を採用し、発送側は入力情報シーケンスに対しまずパリティチェックプロダクトコードエンコードを行い、エンコード後の結果でエンコード規則により因子グラフを生成し、それから重複多重変調エンコードを行い、エンコード後の信号をアンテナで発射する。復調方法では、信号はチャンネルにより伝送され、受信側はアンテナにより信号を受信すると、まずデジタル信号処理を行う。それは同期、均一化等のプロセスを含む。その後、それに対し重複多重復調デコードを行い、最後にデコード後の結果に対し因子グラフ確率伝搬法デコードを行い、最終的にデコード後のシーケンスを得る。本発明では、プロダクトコードデコード方法を採用し、パリティチェックコードをサブコードとし、因子グラフ確率伝搬法の考えをデコード側に使う。それにより、パリティチェックプロダクトコードを採用し、その構成が簡単で、且つデコードプロセスで因子グラフを採用する方法を実現したので、運算複雑度を下げている。【選択図】図１１
2019535196.tif 000002
【特許請求の範囲】
【請求項１】
重複多重に基づく一種類の変調方法であって、下記を含む変調方法：
入力情報を取得する；
入力情報に対しパリティチェックプロダクトコードエンコードを行い、因子グラフを生成する；
重複多重変調エンコードを行う；
エンコード後の信号を発射する。
【請求項２】
請求項1記載の方法において、入力情報に対しパリティチェックプロダクトコードエンコードを行い、下記のことを含む方法：
入力情報をエンコード構成の情報ビットに記入する；
情報ビットの中の情報に対しラインコーディングを行う；
情報ビットの中の情報に対しコラムコーディングを行う；
エンコード後の結果に対しエンコード規則により因子グラフを生成する。
【請求項３】
請求項2記載の方法において、上記エンコード構成は対角エンコード構成、二次元エンコード構成、三次元エンコード構成又は四次元エンコード構成である方法。
【請求項４】
請求項2記載の方法において、入力情報に対しパリティチェックプロダクトコードエンコードを行い、下記のことを含む方法：
入力情報シーケンスを
2019535196.tif 000092
の二次元構成で対応する情報ビットに記入する。その中、入力情報長さ
2019535196.tif 000093
；
2019535196.tif 000094
は行数、
2019535196.tif 000095
は列数である。
各行の第
2019535196.tif 000096
ビットの情報を現在行の前
2019535196.tif 000097
列のモジューロツーサム結果としラインコーディングする；
各列の第
2019535196.tif 000098
ビットの情報を現在列の前
2019535196.tif 000099
行モジューロツーサム結果としコラムコーディングする；
エンコード後の結果に対しエンコード規則により因子グラフを生成する。
【請求項５】
重複多重に基づく一種類の復調方法であって、下記を含む方法：
入力信号を取得する；
入力信号に対し重複多重復調デコードを行う；
因子グラフ確率伝搬法デコードを行う；
デコード結果を出力する。
【請求項６】
請求項5記載の方法において、因子グラフ確率伝搬法デコードが下記を含む方法：
初期対数尤度比を計算する；
反復の最大回数を設定する；
検証ノードを計算し、検証情報を更新する；
変数ノードを計算し、情報メッセージを更新する；
すべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算する；
判決をする；
一定のプリセット条件を満足してから、デコード結果を出力する。
【請求項７】
請求項6記載の方法において、公式
2019535196.tif 000100
又は
2019535196.tif 000101
を採用し検証ノードを計算し、検証情報を更新する；その中、
2019535196.tif 000102
は検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000103
は情報メッセージであり、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000104
は検証ノードの制約する局部コードエレメント情報セットである；
2019535196.tif 000105
は
2019535196.tif 000106
が第j個変数ノードのサブセットを含まないことを示す；
2019535196.tif 000107
は連乗運算である；
及び/又は、上記プリセット条件は最大反復回数となることである方法。
【請求項８】
請求項6記載の方法において、下記のことを含む方法：
公式
2019535196.tif 000108
を採用し変数ノードを計算し、情報メッセージを更新する；
公式
2019535196.tif 000109
を採用しすべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算する；
その中、
2019535196.tif 000110
は発信機の発送する信号の中の発送コードワードである；
2019535196.tif 000111
は受信機の受信する入力信号の中の受信コードワードである；
2019535196.tif 000112
は変数ノードの参加する検証セットで、
2019535196.tif 000113
は
2019535196.tif 000114
が第i個検証ノードのサブセットを含まないことを示す；
2019535196.tif 000115
は検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000116
は受信機の初期受信のチャンネル情報の対数尤度比の表現形式である；
2019535196.tif 000117
は情報メッセージであり、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000118
はすべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を示す。
【請求項９】
重複多重に基づく一種類の変調装置であって、下記のことを含む変調装置：
入力情報取得モジュールは、入力情報を取得することに使われる；
パリティチェックプロダクトコードエンコードモジュールは、入力情報に対しパリティチェックプロダクトコードエンコードを行うことに使われる；
重複多重変調エンコードモジュールは、重複多重変調エンコードを行うことに使われる；
信号発射モジュールは、エンコード後の信号を発射することに使われる。
【請求項１０】
重複多重に基づく一種類の復調装置であって、下記のことを含む復調装置：
入力信号取得モジュールは、入力信号を取得することに使われる；
重複多重復調デコードモジュールは、入力信号に対し重複多重復調デコードを行うことに使われる；
因子グラフ確率伝搬法デコードモジュールは、因子グラフ確率伝搬法デコードを行うことに使われる；
デコード結果出力モジュールは、デコード結果を出力することに使われる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は通信分野に関し、具体的には重複多重に基づく一種類の変調復調方法と装置に関する。
【背景技術】
【０００２】
重複多重（OvXDM：Overlapped X Division Multiplexing）に基づく変調復調技術は多種類の具体的な実現案を含む。例えば、重複時分割多重（OvTDM：Overlapped Time Division Multiplexing）に基づく変調復調技術、重複周波数分割多重（OvFDM：Overlapped Frequency Division Multiplexing）に基づく変調復調技術、重複コード分割多重（OvCDM：Overlapped Code Division Multiplexing）に基づく変調復調技術、重複スペース分割多重（OvSDM：Overlapped Space Division Multiplexing）に基づく変調復調技術、重複混合多重（OvHDM：Overlapped Hybrid Division Multiplexing）に基づく変調復調技術などである。
【０００３】
説明する必要のあるのは、本発明で言及されたOvXDMでは、Xは任意領域を示す。例えば時間T、スペースS、周波数F、コード分割C、混合H等。
【０００４】
下記からOvTDMのみを例とし簡単に説明する。
【０００５】
まず、時間分割（以下は時分割と略称する）多重（TDM：Time Division Multiplexing）とは、デジタル通信で狭い時間持続期間を占める幾つかの信号符号に一つの広い時間持続期間を共有させる技術である。図1は、通常の時分割多重技術を示す図である。
【０００６】
図1で各被多重信号符号の時間持続期間（エンジニアリング上はタイムスロット広さ）はそれぞれT1、T2、T3、T4、...、とする。エンジニアリング上は、通常それらにタイムスロット広さを占めさせる。ΔTは最小保護タイムスロットで、実際の保護タイムスロット広さはもっと余裕がある。ΔTは使っているデマルチプレックスゲート回路の過渡時間広さプラスシステムの最大タイムジッター量より大きい。これは最もよくある時分割多重技術である。現有のほとんどのマルチチャンネルデジタル放送システム、マルチチャンネルデジタル通信等のシステムはこのような技術を採用している。
【０００７】
このような技術がデジタル通信に応用される時の最大特徴は、被多重信号符号の間は時間の上で完全にお互いに隔離され、絶対に相互干渉がなく、被多重の信号符号はいかなる制限がなく、各信号の符号持続期間（タイムスロット広さ)は異なる広さがあり、異なる通信体制に適用されることである。それらのタイムスロットはお互いに重複や交差がなければよいため、使用用途は最も広い。しかし、このような多重自身は、システムのスペクトル効率改善には完全に役立たない。
【０００８】
そのため、伝統的な観点は隣接チャンネルの間は時間領域で重複せず、隣接チャンネルの間での干渉発生を避けるが、このような技術はスペクトル効率の向上を制約している。現有技術の時分割多重技術の観点は各チャンネルの間でお互いに隔離しないだけでなく、強い相互重複があってもよいことである。図2のように、現有技術はチャンネル間の重複を新しいエンコード制約関係と見なし、当該制約関係により相応な変調と復調技術を提示しているので、重複時分割多重（OvTDM：Overlapped Time Division Multiplexing）と呼ぶ。このような技術はスペクトル効率が重複回数Kに比例し増える。
【０００９】
図3を参照されたい。重複時分割多重システムは信号発信機A01と受信機A02を含む。
【００１０】
発信機A01は重複時分割多重変調装置101と発射装置102を含む。重複時分割多重変調装置101は、入力信号シーケンスを持つ多重変調包絡波形を生成することに使われる。発射装置102は当該多重変調包絡波形を受信機A02に発射することに使われる。
【００１１】
受信機A02は受信装置201とシーケンステスト装置202を含む。受信装置201は発射装置102からの多重変調包絡波形を受信することに使われる；シーケンステスト装置202は受信された多重変調包絡波形に対し時間領域でのデータシーケンステストを行い、判決出力をすることに使われる。
【００１２】
通常は、受信機A02はさらに受信装置201とシーケンステスト装置202の間に設置される前処理装置203を含む。それは各フレーム内の同期受信デジタル信号シーケンスを形成することにアシストすることに使われる。
【００１３】
発信機A01で、入力されるデジタル信号シーケンスは重複時分割多重変調装置101により複数符号の時間領域での相互重複の発射信号を生成し、それから発射装置102により当該発射信号を受信機A02に発射する。受信機A02の受信装置201は発射装置102からの信号を受信し、前処理装置203によりシーケンステスト装置202のテスト受信に向けるデジタル信号を生成し、シーケンステスト装置202は受信信号に対し時間領域でのデータシーケンステストを行い、判決を出力する。
【００１４】
図4を参照されたい。重複時分割多重変調装置101（OvTDM変調装置）は波形生成モジュール301、変位モジュール302、乗法モジュール303と重複モジュール304を含む。
【００１５】
波形生成モジュール301は、設計パラメータにより時間領域で波形の滑らかな初期包絡波形を生成することに使われる。
【００１６】
変位モジュール302は、重複多重回数により初期包絡波形を時間領域で予定変位間隔にで変位し、各固定間隔の変位包絡波形を得ることに使われる。
【００１７】
変調モジュール305は、入力されるデジタル信号シーケンスを正負符号表現の信号符号シーケンスに変換することに使われる。
【００１８】
乗法モジュール303は、上記信号符号シーケンスを偏移後の各固定間隔の変位包絡波形と相乘し、各変調包絡波形を得ることに使われる。
【００１９】
重複モジュール304は、各変調包絡波形を時間領域で重複し、入力信号シーケンスを持つ多重変調包絡波形を得ることに使われる。
【００２０】
図5は、受信機A02の前処理装置203のブロック図である。
【００２１】
前処理装置203は同期器501、チャンネル推定器502とデジタル化処理器503を含む。その中、同期器501は受信信号に対し受信機で符号時間同期を行う；それからチャンネル推定器502はチャンネルパラメータに対し推定する；デジタル化処理器503は書くフレーム内の受信信号に対しデジタル化処理を行うことで、シーケンステスト装置のシーケンステスト受信に適合するデジタル信号シーケンスを生成する。
【００２２】
図6は、受信機A02のシーケンステスト装置202のブロック図である。
【００２３】
シーケンステスト装置20は分析ユニットストレージ601、比較器602及び複数の保留ルートストレージ603とユークリッド距離ストレージ604又は重み付けユークリッド距離ストレージ（図に示されていない）を含む。テストプロセスでは、分析ユニットストレージ601は重複時分割多重システムの複数畳み込みエンコードモデル及びトレリスグラフを作り、重複時分割多重システムの全ての状態を並べ、保存する；比較器602は分析ユニットストレージ601のトレリス（Trellis）グラフにより、受信デジタル信号と最小ユークリッド距離又は重み付け最小ユークリッド距離のルートを捜索する；保留ルートストレージ603とユークリッド距離ストレージ604又は重み付けユークリッド距離ストレージは、それぞれ比較器602から出力される保留ルートとユークリッド距離又は重み付けユークリッド距離を保存することに使われる。その中、保留ルートストレージ603とユークリッド距離ストレージ604又は重み付けユークリッド距離ストレージは、各安定状態のためにそれぞれ一つ準備する必要がある。保留ルートストレージ603の長さは優先的に4K〜5Kを選ぶ。ユークリッド距離ストレージ604又は重み付けユークリッド距離ストレージは優先的に相対距離のみを保存する。
【００２４】
OvXDMシステムで、信号発信機は信号を変調してから信号受信機に発射し、信号受信機は変調信号を受信する後に復調する。復調プロセスはデコードステップ（即ち上記シーケンステスト装置の実行するシーケンステストステップ）を含む。伝統通信システムでは、だいたいChaseアルゴリズムでデコードする。そのアルゴリズムは大量のソート運算に関わり、計算量が非常に多い。
【発明の概要】
【００２５】
本発明は重複多重に基づく一種類の変調復調方法と装置を提供し、伝統デコードでだいたいChaseアルゴリズムでデコードし、アルゴリズムプロセスは大量のソート運算に関わり、計算複雑度が高い問題を解決している。
【００２６】
本発明の第一の特徴によれば、本発明は重複多重に基づく一種類の変調方法を提供し、下記を含む：
入力情報を取得する；
入力情報に対しパリティチェックプロダクトコードエンコードを行い、因子グラフを生成する；
重複多重変調エンコードを行う；
エンコード後の信号を発射する。
【００２７】
本発明の第二の特徴によれば、本発明は重複多重に基づく一種類の復調方法をも提供し、下記を含む：
入力信号を取得する；
入力信号に対し重複多重復調デコードを行う；
因子グラフ確率伝搬法デコードを行う；
デコード結果を出力する。
【００２８】
本発明の第三の特徴によれば、本発明は重複多重に基づく一種類の変調装置をも提供し、下記を含む：
入力情報取得モジュールは、入力情報を取得することに使われる；
パリティチェックプロダクトコードエンコードモジュールは、入力情報に対しパリティチェックプロダクトコードエンコードを行うことに使われる；
重複多重変調エンコードモジュールは、重複多重変調エンコードを行うことに使われる；
信号発射モジュールは、エンコード後の信号を発射することに使われる。
【００２９】
一つの実施例では、パリティチェックプロダクトコードエンコードモジュールは下記を含む：
情報ビット記入ユニットは、入力情報をエンコード構成の情報ビットに記入することに使われる；
具体的には、情報ビット記入ユニットは入力情報シーケンスを
2019535196.tif 000003
の二次元構成で対応する情報ビットに記入することに使われる。その中、入力情報の長さ
2019535196.tif 000004
；
2019535196.tif 000005
は行数で、
2019535196.tif 000006
は列数である；
ラインコーディングユニットは、情報ビットの中の情報に対しラインコーディングを行うことに使われる；具体的には、ラインコーディングユニットは各行の第
2019535196.tif 000007
ビットの情報を現在行の前
2019535196.tif 000008
列モジューロツーサム結果としラインコーディングすることに使われる；
コラムコーディングユニットは、情報ビットの中の情報に対しコラムコーディングを行うことに使われる；具体的には、コラムコーディングユニットは各列の第
2019535196.tif 000009
ビットの情報を現在列の前
2019535196.tif 000010
行モジューロツーサム結果としコラムコーディングすることに使われる；
因子グラフ生成ユニットは、エンコード後の結果をエンコード規則により因子グラフを生成することに使われる。
【００３０】
一つの実施例では、上記エンコード構成は対角エンコード構成、二次元エンコード構成、三次元エンコード構成又は四次元エンコード構成である。
【００３１】
本発明の第四の特徴によれば、本発明は重複多重に基づく一種類の復調装置をも提供し、下記を含む：
入力信号取得モジュールは、入力信号を取得することに使われる；
重複多重復調デコードモジュールは、入力信号に対し重複多重復調デコードを行うことに使われる；
因子グラフ確率伝搬法デコードモジュールは、因子グラフ確率伝搬法デコードを行うことに使われる；
デコード結果出力モジュールは、デコード結果を出力することに使われる。
【００３２】
一つの実施例では、因子グラフ確率伝搬法デコードモジュールは下記を含む：
初期対数尤度比計算ユニットは、初期対数尤度比を計算することに使われる；
最大反復回数設定ユニットは、反復の最大回数を設定することに使われる；
検証情報更新ユニットは、検証ノードを計算し、検証情報を更新することに使われる；
情報消息更新ユニットは、変数ノードを計算し、情報メッセージを更新することに使われる；
対数尤度比更新ユニットは、すべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算することに使われる；
判決ユニットは、情報シーケンス判決をすることに使われる；
デコード結果出力ユニットは、一定のプリセット条件を満足してから、デコード結果を出力することに使われる。
【００３３】
一つの実施例では、検証情報更新ユニットは公式
2019535196.tif 000011
又は
2019535196.tif 000012
で検証ノードを計算し、検証情報を更新することに使われる；その中、
2019535196.tif 000013
は検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000014
は情報メッセージであり、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000015
は検証ノードの制約する局部コードエレメント情報セットである；
2019535196.tif 000016
は
2019535196.tif 000017
が第j個変数ノードのサブセットを含まないことを示す；
2019535196.tif 000018
は連乗運算である；
及び/又は、上記プリセット条件は最大反復回数となることである。
【００３４】
一つの実施例では、情報消息更新ユニットは公式
2019535196.tif 000019
で変数ノードを計算し、情報メッセージを更新することに使われる；
対数尤度比更新ユニットは公式
2019535196.tif 000020
ですべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算することに使われる；
その中、
2019535196.tif 000021
は発信機の発送する信号の中の発送コードワードである；
2019535196.tif 000022
は受信機の受信した入力信号の中の受信コードワードである；
2019535196.tif 000023
は変数ノードの参加する検証セットであり、
2019535196.tif 000024
は
2019535196.tif 000025
が第i個検証ノードのサブセットを含まないことを示す；
2019535196.tif 000026
は検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000027
は受信機の初期受信のチャンネル情報の対数尤度比表示形式である；
2019535196.tif 000028
は情報メッセージであり、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す；
2019535196.tif 000029
はすべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を示す。
【００３５】
本発明の提供する重複多重に基づく変調復調方法と装置の中で、変調方法では、プリコーディング構成を採用し、発送側は入力情報シーケンスに対しまずパリティチェックプロダクトコードエンコードを行い、エンコード後の結果でエンコード規則により因子グラフを生成し、重複多重変調エンコードを行い、エンコード後の信号をアンテナにより発射する。復調方法では、信号はチャンネルで伝送され、受信側はアンテナにより信号を受信する後に、まず同期、均一化などを含むデジタル信号処理を行い、それからそれに対し重複多重復調デコードを行い、最後にデコード後の結果に対し因子グラフ確率伝搬法デコードを行い、最終的にデコード後のシーケンスを得る。本発明では、プロダクトコードデコード方法を採用し、パリティチェックコードをサブコードとし、因子グラフ確率伝搬法の考えをデコード側に使う。それにより、パリティチェックプロダクトコードを採用し、その構成が簡単で、且つデコードプロセスで因子グラフを採用する方法を実現したので、運算複雑度を下げている。
【図面の簡単な説明】
【００３６】
【図１】図1は、通常の時分割多重技術を示す図である。
【図２】図2は、重複時分割多重原理を示す図である。
【図３】図3は、一種類の重複時分割多重システムの構成を示す図である。
【図４】図4は、一種類の重複時分割多重変調装置の構成を示す図である。
【図５】図5は、一種類の受信機前処理装置の構成を示す図である。
【図６】図6は、一種類の受信機シーケンステスト装置の構成を示す図である。
【図７】図7は、本発明の一種類の実施例でのパリティチェックプロダクトコードの構成図である。
【図８】図8は、本発明の一種類の実施例での双方向伝達因子グラフである。
【図９】図9は、本発明の一種類の実施例でのパリティチェックプロダクトコードマトリックスと因子グラフの対応関係を示す図である。
【図１０】図10は、本発明の一種類の実施例でのプリコーディングOvXDMシステム発送側のブロック図である。
【図１１】図11は、本発明の一種類の実施例での重複多重に基づく変調方法のフローを示す図である。
【図１２】図12は、本発明の一種類の実施例の重複多重に基づく変調方法でのパリティチェックプロダクトエンコードステップのフローを示す図である。
【図１３】図13は、本発明の一種類の実施例の重複多重に基づく変調方法での重複多重変調エンコードステップのフローを示す図である。
【図１４】図14は、Kチャンネル波形多重の原理を示す図である。
【図１５】図15は、Kチャンネル波形の符号重複プロセスの原理を示す図である。
【図１６】図16は、本発明の一種類の実施例での基于重複多重の復調方法のフローを示す図である。
【図１７】図17は、本発明の一種類の実施例の重複時分割多重に基づく復調方法でのデコードステップのフローを示す図である。
【図１８】図18は、K＝3時の重複時分割多重システムの入力- 出力関係のトリーグラフである。
【図１９】図19は、ノード状態遷移図である。
【図２０】図20は、K=3時のOvTDMシステムTrellisグラフである。
【図２１】図21は、本発明の一種類の実施例の重複多重に基づく復調方法での因子グラフ確率伝搬法デコードステップのフローを示す図である。
【図２２】図22は、本発明の一種類の実施例での重複多重に基づく変調装置のモジュールを示す図である。
【図２３】図23は、本発明の一種類の実施例の重複多重に基づく変調装置でのパリティチェックプロダクトコードエンコードモジュールのユニットを示す図である。
【図２４】図24は、本発明の一種類の実施例での重複多重に基づく復調装置のモジュールを示す図である。
【図２５】図25は、本発明の一種類の実施例の重複多重に基づく復調装置での因子グラフ確率伝搬法デコードモジュールのユニットを示す図である。
【発明を実施するための形態】
【００３７】
まず、説明する必要のあるのは、本発明の提供する重複多重（OvXDM）に基づく変調復調方法と装置では、Xは任意領域を示すことである。例えば時間T、スペースS、周波数F、コード分割C、混合Hなど。説明のために、本発明の実施例は主としては重複時分割多重（OvTDM）と重複周波数分割多重（OvFDM）を例とし説明する。本分野の技術者にとっては、本発明が保護を要求する変調復調方法と装置は同様にほかの領域の重複多重技術に応用できると分かるのである。
【００３８】
OvXDMシステムでは、信号発信機は信号を変調する後に信号受信機に発射し、信号受信機は変調信号を受信する後に復調する。復調プロセスはデコードステップを含む。重複多重変調復調技術の研究では、本分野の技術者はみんな伝統的なデコード方法を採用している。伝統のデコード方法では、だいたいChaseアルゴリズムを使う。そのアルゴリズムは大量のソート運算に関わり、計算量が非常に多い。
【００３９】
しかし、現在のデコード方法は上記問題が存在しているが、これらの方法は広く使われているので、本分野の技術者は全面的にこのような方法を受け取り、コストを掛けより良いデコード方法を探していない現状である。
【００４０】
一般の状況では、簡単なOvXDM(Xは時間T、周波数F、スペースS、コード領域C又は混合H等を示す)システムは受信側で波形デコードを完成する後に、プロセス全体が終了する。より複雑な状況では、OvXDMシステムと普通の伝統通信技術と結び合わせることで、システム全体の性能を向上する。例えばカスケードOvXDMシステム、プリコーディングOvXDMシステムなど。而誤り訂正コードは優れる誤り訂正能力があり、システム全体の性能を上げ、エラーレートを下げるので、大部分の技術者は誤り訂正コードをOvXDMシステムに応用する。
【００４１】
プロダクトコードは是ブロック構成の誤り訂正コードであり、デコード側で反復の考えを導入することで、現在流行っているTurboプロダクトコードとなり、つまりTPCエンコードである。このようなエンコードは現在の通信システムで非常に幅広く応用されている。伝統の通信システムでは、だいたいChaseアルゴリズムでデコードしている。そのアルゴリズムは大量のソート運算に関わり、計算量が非常に多い。
【００４２】
本発明の提出するプロダクトコードは、とても簡単なパリティチェックコードをサブコードとし、非常に柔軟で便利にコード長さ制御とマッチングを行える。デコード側では、采因子グラフ確率伝搬法の考えに基づく反復デコード方法を採用し、その操作が柔軟で便利である。本発明は因子グラフ確率伝搬法に基づくデコードの考えを、パリティチェックプロダクトコードのデコード方法に導入している。
【００４３】
よく使われる誤り訂正コードはプロダクトコード（Turbo Product Code、TPCコード）、低密度パリティチェックコード（Low Density Parity Check Code、LDPCコード）を含む。本発明はパリティチェックプロダクトコードを例とする。そのエンコード構成は図7のようである。そのエンコード構成はとても簡潔で、行列サブコードは同じコード長さを選べるし、異なるコード長さを選べる。構成は柔軟で便利であり、且つコードレートを容易に調整できる。
【００４４】
このエンコード構成により、このエンコードによる因子グラフを容易に得られる。図8に示されているようである。その中、図の下側のノードは変数ノードで、数量はエンコードマトリックスブロックの中のコード長さである。図の上側のノードは検証ノードで、数量は検証ビットの長さである。OvXDMシステムに応用すると、対応するエンコードマトリックスと双方向因子グラフの関係は図9のようである。パリティチェックプロダクトコードからなる関連因子グラフにより、この因子グラフは大きな周囲長さがあり、確率伝搬法の方法によるデコードにも適用されることが分かる。確率伝搬法デコードする時に、対数領域での操作によりメッセージ放送を行える。
【００４５】
以下に、具体的な実施方法と図面とを結び合わせ、本発明に対しさらに詳しい説明を行う。
【００４６】
実施例一
本実施例では、仮に入力情報シーケンスはxとし、情報長さ
2019535196.tif 000030
=100、
2019535196.tif 000031
=10、
2019535196.tif 000032
=10、重複回数K=5、チェビシェフウィンドウを多重波形とし、変調方法はBPSKで、プリコーディングはTPC誤り訂正コードを例とする。図10のプリコーディング構成を採用し、システム処理プロセスは：発送側は入力情報シーケンスに対しまずTPCエンコードを行い、それからOvXDMエンコードを行い、エンコード後の信号をアンテナにより発射する。信号はチャンネルにより伝送され、受信側はアンテナにより信号を受信する後に、まず同期、均一化などのプロセスを含むデジタル信号処理を行う。それからそれに対しOvXDMデコードを行い、最後にデコード後の結果をTPCデコードし、即ち本発明の因子グラフ確率伝搬法デコード方法を使い、最終的にデコード後のシーケンスがx'であると得る。
【００４７】
よく使われるOvXDM(Xは時間T、周波数F、スペースS、コード領域C又は混合Hなどを示す)デコードは最大尤度シーケンスデコードアルゴリズム（例えばビタビデコード）、最大事後確率アルゴリズム（例えばBCJRアルゴリズム）、MAPアルゴリズム、Log_MAPアルゴリズムなどを含む。下記から主としてはパリティチェックプロダクトコードのエンコードプロセス及び対応する因子グラフ確率伝搬法デコード方法を説明する。
【００４８】
まず、本実施例のこれから使う変数に対し定義する。
2019535196.tif 000033
を発信コードワードとし、
2019535196.tif 000034
を受信コードワードとする。
2019535196.tif 000035
：検証ノードの制約する局部コードエレメント情報セット。
2019535196.tif 000036
は
2019535196.tif 000037
が第j個変数ノードのサブセットを含まないことを示す。
2019535196.tif 000038
：変数ノードの参加する検証セット。
2019535196.tif 000039
は
2019535196.tif 000040
が第i個検証ノードのサブセットを含まないことを示す。
2019535196.tif 000041
：初期受信のチャンネル情報の対数尤度比の表現形式。
2019535196.tif 000042
：情報メッセージで、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す。
2019535196.tif 000043
：情報メッセージで、すべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を示す。
2019535196.tif 000044
：検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す。
【００４９】
図11を参照されたい。本実施例は重複多重に基づく一種類の変調方法を提供している。下記ステップを含む：
ステップ1.1：入力情報を取得する。情報、つまり図4で入力される信号デジタルシーケンスを入力する。
【００５０】
ステップ1.2：入力情報に対しパリティチェックプロダクトコードエンコードを行う。
【００５１】
ステップ1.3：重複多重変調エンコードを行う。
【００５２】
ステップ1.4：エンコード後の信号を発射する。図4で生成した多重変調包絡波形を発射信号とし受信側に発射する。
【００５３】
図12を参照されたい。ステップ1.2は下記サブステップを含む：
サブステップ1-1：入力情報をエンコード構成の情報ビットに記入する。その中、エンコード構成は対角エンコード構成、二次元エンコード構成、三次元エンコード構成、四次元エンコード構成又はより高い次元のエンコード構成であってもよい。本実施例では二次元エンコード構成を例とし説明する。
【００５４】
具体的には、入力情報シーケンスを
2019535196.tif 000045
の二次元構成で対応する情報ビットに記入する。
2019535196.tif 000046
は行数で、
2019535196.tif 000047
は列数である。
【００５５】
サブステップ1-2：情報ビットの中の情報に対しラインコーディングを行う。
【００５６】
具体的には、各行の第
2019535196.tif 000048
ビット情報を現在行の前
2019535196.tif 000049
列モジューロツーサム結果としてラインコーディングする。
【００５７】
サブステップ1-3：情報ビットの中の情報に対しコラムコーディングを行う。
【００５８】
具体的には、各列の第
2019535196.tif 000050
ビット情報を現在列の前
2019535196.tif 000051
行モジューロツーサム結果としてコラムコーディングする。
【００５９】
エンコード後のマトリックスの大きさは
2019535196.tif 000052
で、そのTPCのコードレートは
2019535196.tif 000053
であり、本実施例では0.8264である。
【００６０】
サブステップ1-4：エンコード後の結果に対しエンコード規則により因子グラフを生成する。
【００６１】
図13のように、ステップ1.3は下記サブステップを含む：
サブステップ2.1：設計パラメータにより時間領域での初期包絡波形h(t)を生成する。
【００６２】
初期包絡波形を生成時に、ユーザの設計パラメータ入力により、実際のシステムでシステム性能指標により柔軟に配置することを実現する。
【００６３】
ある実施例では、初期包絡波形のサイドローブ減衰が確定された場合、設計パラメータは初期包絡波形の長さLを含む。例えば初期包絡波形がバートレット包絡波形である場合である。
【００６４】
ある実施例では、設計パラメータは初期包絡波形の長さLとサイドローブ減衰rを含む。例えば初期包絡波形がチェビシェフ包絡波形である場合である。
【００６５】
もちろん、初期包絡波形がほかの形式である場合、相応な初期包絡波形の特徴により設計パラメータを確定してもよい。
【００６６】
サブステップ2.2：重複多重回数Kにより初期包絡波形を相応な領域（本実施例では時間領域とする）で予定変位間隔で変位し、各固定間隔の変位包絡波形h(t-i*△T)を得る。
【００６７】
その中、変位間隔は時間間隔△Tで、時間間隔△Tは：△T=L/K。この時、信号の符号広さ△Tである。
【００６８】
また、△Tがシステムサンプリングレートの逆数より小さくないことを保証する必要がある。
【００６９】
iのバリューは入力符号長さNと関連し、且つiは0〜N-1の整数とする。例えば、N=8の時に、iは0〜7の整数とする。
【００７０】
サブステップ2.3：ステップ1.2のエンコード後に得られた信号デジタルシーケンスを正負符号表現の信号符号シーケンスに変換する。
【００７１】
具体的には、デジタル信号シーケンスの中の0を+Aに変換し、1を-Aに変換し、Aを0でない任意数とすることで、正負符号シーケンスを得る。例えば、Aを1とする時に、入力される{0,1}ビットシーケンスをBPSK（Binary Phase Shift Keying、位相偏移変調）により変調し{+1、-1}符号シーケンスに変換する。
【００７２】
サブステップ2.4：変換後の信号符号シーケンスxi（本実施例でxi={+1 +1 -1 -1 -1 +1 -1 +1}）と各固定間隔の変位包絡波形h(t-i*△T)を相乘し、各変調包絡波形xi h(t-i*△T)を得る。
【００７３】
サブステップ2.5：各変調包絡波形xi h(t-i*△T)を相応な領域（本実施例では時間領域とする）で重複し、入力信号シーケンスを持つ多重変調包絡波形、即ち発送する信号を得る。
【００７４】
発送する信号は下記のように示される：
2019535196.tif 000054
ステップ1.4でエンコード後の信号を発射するのは、得られる多重変調包絡波形を発送信号とし発射することである。
【００７５】
その故、本実施例では、Aのバリューが1の場合、重複後の出力符号（出力する信号符号シーケンス）は：s(t) = {+1 +2 +1 -1 -3 -1 -1 +1}とする。
【００７６】
図14は、Kチャンネル波形多重の原理を示す図であり、平行四辺形の形状となっている。その中、各行は一つの発送しようとする符号xiと早々時刻の包絡波形h(t-i*△T)の相乘後に得られる発送待ち信号波形xih(t-i*△T)を示す。a0〜ak-1は各ウィンドウ関数波形（包絡波形）に対しK回セグメントを分けてから得られる各部分の係数値を示す。具体的には振幅値に関する係数である。
【００７７】
入力されるデジタル信号シーケンスを正負符号シーケンスに変換する時に、入力されるデジタル信号シーケンスの中の0,1を±Aに変換し、Aのバリューは0でない任意数であることで、正負符号シーケンスを得る。例えば、Aのバリューが1の場合、入力される{0,1}ビットシーケンスをBPSK変調により{+1、-1}符号シーケンスに変換することで、正負符号シーケンスを得る。そのため、図15はKチャンネル波形の符号重複プロセスの原理を示す図である。図15の重複プロセスでは、第1行左側の三つの数字は第1個入力符号+1を示し、第2行左側の三つの数字は第2個入力符号+1を示し、第3行左側の三つの数字は第3個入力符号-1を示し、第1行の真ん中の三つの数字は第4個入力符号-1を示し、第2行の真ん中の三つの数字は第5個入力符号-1を示し、第3行の真ん中の三つの数字は第6個入力符号+1を示し、第1行右側の三つの数字は第7個入力符号-1を示し、第2行右側の三つの数字は第8個入力符号+1を示す。そのため、三つの波形が重複する後に、得られる出力符号は｛+1 +2 +1 -1 -3 -1 -1 +1｝である。
【００７８】
もちろん、入力符号の長さはほかの数値である場合、図14と図15に示される方法により重複し、出力符号を得る。
【００７９】
説明する必要のあるのは、ステップ1.3の重複多重変調エンコードは上記方法以外に、現有技術の任意一種類実行可能な方法を採用してもよい。
【００８０】
実施例二
上記実施例一の提供する重複多重に基づく変調方法により、本実施例はそれと相応に重複多重に基づく一種類の復調方法を提供している。図16を参照されたい。重複多重に基づく復調方法は下記ステップを含む：
ステップ3.1：入力信号を取得する。入力信号は図4で発射される多重変調包絡波形信号である。
【００８１】
ステップ3.2：入力信号に対し重複多重復調デコードを行う。
【００８２】
ステップ3.3：因子グラフ確率伝搬法デコードを行う。
【００８３】
ステップ3.4：デコード結果を出力する。
【００８４】
図17を参照されたい。本実施例では、ステップ3.2は具体的に下記サブステップを含む：
(1)まず受信信号に対し搬送波同期、フレーム同期、符号時間同期などを含む同期を行う。
【００８５】
(2)サンプリング定理により、各フレーム内の受信信号に対しデジタル化処理を行う。
【００８６】
(3)受信された波形に対し波形発送時間間隔により切断する。
【００８７】
(4)一定のデコードアルゴリズムにより切割後の波形をデコードする。デコードプロセスは図18を参照されたい。それは K=3の場合、重複時分割入力-出力関係図である。図19はノード状態遷移図である。図20は K=3の場合、OvTDMシステム のTrellisグラフである。ビタビ（Viterbi）デコードは畳み込みコードで最も幅広く使われている方法である。それの基本考えはTrellisグラフのすべてのルートを遍歴し、Trellisグラフ状態遷移プロセスで各状態に到達する幾つかの支線と正確ルートの距離比較により、距離最小のルートしか保留しない。比較選別により正確ルートの推定を入手し、デコードを実現する。
【００８８】
説明する必要のあるのは、ステップ3.2での重複多重復調デコードは上記方法を採用してもよく、現有技術のほかの方法を採用してもよいことである。
【００８９】
図21を参照されたい。仮にデジタル信号処理により得られるノイズ分散は
2019535196.tif 000055
とし、OvXDMデコード後のシーケンスは
2019535196.tif 000056
とすると、ステップ3.3は下記サブステップを含む：
サブステップ3-1：初期対数尤度比を計算する。
【００９０】
具体的には、対数尤度比の結果は一般的に変調方法、OvXDMシステムのデコード方法（即ち重複多重復調デコード方法）などと関連し、最終的な目的はそのソフトバリューの取得である。
【００９１】
サブステップ3-2：反復の最大回数を設定する。
【００９２】
具体的には、本実施例は反復回数を50とする。反復回数が大きいほど、その最終のデコードシーケンスは理論シーケンスに近づけるが、大きすぎると、それと相応にアルゴリズムの複雑度を上げてしまう。
【００９３】
サブステップ3-3：検証ノードを計算し、検証情報を更新する。
【００９４】
具体的には、公式
2019535196.tif 000057
により各行と各列の検証ノード并を計算し更新する。その中、
2019535196.tif 000058
は連乗運算である。
【００９５】
サブステップ3-4：変数ノードを計算し、情報メッセージを更新する。
【００９６】
具体的には、公式
2019535196.tif 000059
により各行と各列の情報メッセージを更新する。
【００９７】
サブステップ3-5：すべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算する
2019535196.tif 000060
。
【００９８】
具体的には、公式
2019535196.tif 000061
により新しい情報メッセージを計算する。
【００９９】
サブステップ3-6：判決をする。
【０１００】
具体的には、本実施例はハード判決の方法で、BPSK変調方法を採用する。仮にここで対応する変調マッピングは1->-1、0->+1とすれば、即ち
2019535196.tif 000062
。
【０１０１】
サブステップ3-7：一定のプリセット条件を満足してから、デコード結果を出力する。本実施例では、プリセット条件は最大反復回数となることである。
【０１０２】
最大反復回数となっているかを判断する。満足すれば、今回のデコードが終了し、デコードシーケンスを出力する。そうでなければ、サブステップ3-3に戻り計算しなおす。この時の情報メッセージ
2019535196.tif 000063
は前回反復運算の結果である。
【０１０３】
実施例三
図22を参照されたい。実施例一の提供する重複多重に基づく変調方法により、本実施例はそれと相応に重複多重に基づく一種類の変調装置を提供している。それは入力情報取得モジュールA11、パリティチェックプロダクトコードエンコードモジュールA12、重複多重変調エンコードモジュールA13と信号発射モジュールA14を含む。
【０１０４】
入力情報取得モジュールA11は入力情報を取得することに使われる。
【０１０５】
パリティチェックプロダクトコードエンコードモジュールA12は入力情報に対しパリティチェックプロダクトコードエンコードを行うことに使われる。
【０１０６】
重複多重変調エンコードモジュールA13は重複多重変調エンコードを行うことに使われる。
【０１０７】
信号発射モジュールA14はエンコード後の信号を発射することに使われる。
【０１０８】
図23を参照されたい。本実施例では、パリティチェックプロダクトコードエンコードモジュールA12は、情報ビット記入ユニットA21、ラインコーディングユニットA22、コラムコーディングユニットA23と因子グラフ生成ユニットA24を含む。
【０１０９】
情報ビット記入ユニットA21は入力情報をエンコード構成の情報ビットに記入することに使われる。その中、エンコード構成は対角エンコード構成、二次元エンコード構成、三次元エンコード構成、四次元エンコード構成又はより高い次元のエンコード構成であってもよい。本実施例では二次元エンコード構成を例とし説明する。
【０１１０】
具体的には、情報ビット記入ユニットは入力情報シーケンスを
2019535196.tif 000064
の二次元構成で対応する情報ビットに記入することに使われる。その中、入力情報長さ
2019535196.tif 000065
。
2019535196.tif 000066
は行数、
2019535196.tif 000067
は列数である。
【０１１１】
ラインコーディングユニットA22は、情報ビットの中の情報に対しラインコーディングを行うことに使われる。具体的には、ラインコーディングユニットは、各行の第
2019535196.tif 000068
ビット情報を現在行の前
2019535196.tif 000069
列モジューロツーサムの結果としラインコーディングすることに使われる。
【０１１２】
コラムコーディングユニットA23は、情報ビットの中の情報に対しコラムコーディングを行うことに使われる。具体的には、コラムコーディングユニットは、各列の第
2019535196.tif 000070
ビット情報を現在列の前
2019535196.tif 000071
行モジューロツーサムの結果としコラムコーディングすることに使われる。
【０１１３】
因子グラフ生成ユニットA24は、エンコード後の結果をエンコード規則により因子グラフを生成することに使われる。
【０１１４】
本実施例の提供する重複多重に基づく変調装置は、実施例一の提供する重複多重に基づく変調方法と対応しているので、ここでその原理については再度説明しない。
【０１１５】
実施例四
図24を参照されたい。実施例二の提供する重複多重に基づく復調方法により、本実施例はそれと相応に重複多重に基づく一種類の復調装置を提供している。それは入力信号取得モジュールB11、重複多重復調デコードモジュールB12、因子グラフ確率伝搬法デコードモジュールB13とデコード結果出力モジュールB14を含む。
【０１１６】
入力信号取得モジュールB11用于入力信号を取得する。
【０１１７】
重複多重復調デコードモジュールB12は、入力信号に対し重複多重復調デコードを行うことに使われる。
【０１１８】
因子グラフ確率伝搬法デコードモジュールB13は、因子グラフ確率伝搬法デコードを行うことに使われる。
【０１１９】
デコード結果出力モジュールB14は、デコード結果を出力することに使われる。
【０１２０】
図25を参照されたい。本実施例では、因子グラフ確率伝搬法デコードモジュールB13は、初期対数尤度比計算ユニットB21、最大反復回数設定ユニットB22、検証情報更新ユニットB23、情報メッセージ更新ユニットB24、対数尤度比更新ユニットB25、判決ユニットB26とデコード結果出力ユニットB27を含む。
【０１２１】
初期対数尤度比計算ユニットB21は、初期対数尤度比を計算することに使われる。
【０１２２】
最大反復回数設定ユニットB22は、反復の最大回数を設定することに使われる。
【０１２３】
検証情報更新ユニットB23は、検証ノードを計算し、検証情報を更新することに使われる。
【０１２４】
情報メッセージ更新ユニットB24は、変数ノードを計算し、情報メッセージを更新することに使われる。
【０１２５】
対数尤度比更新ユニットB25は、すべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算することに使われる。
【０１２６】
判決ユニットB26は、判決をすることに使われる。
【０１２７】
デコード結果出力ユニットB27は、一定のプリセット条件を満足してから、デコード結果を出力することに使われる。上記プリセット条件は最大反復回数となることである。
【０１２８】
具体的には、本実施例はハード判決の方法で、BPSK変調方法を採用する。仮にここで対応する変調マッピングは1->-1、0->+1とすると、即ち
2019535196.tif 000072
。
【０１２９】
本実施例では、検証情報更新ユニットB23は公式
2019535196.tif 000073
又は
2019535196.tif 000074
で検証ノードを計算し、検証情報を更新することに使われる。その中、
2019535196.tif 000075
は検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す。
2019535196.tif 000076
は情報メッセージであり、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す。
2019535196.tif 000077
は検証ノードの制約する局部コードエレメント情報セットである。
2019535196.tif 000078
は
2019535196.tif 000079
が第j個変数ノードのサブセットを含まないことを示す。
2019535196.tif 000080
は連乗運算である。
【０１３０】
情報メッセージ更新ユニットB24は公式
2019535196.tif 000081
で変数ノードを計算し、情報メッセージを更新することに使われる。
【０１３１】
対数尤度比更新ユニットB25は、公式
2019535196.tif 000082
ですべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を計算することに使われる。
【０１３２】
その中、
2019535196.tif 000083
は発信機の発送する信号の中の発送コードワードである。
2019535196.tif 000084
は受信機の受信した入力信号の中の受信コードワードである。
2019535196.tif 000085
は変数ノードの参加する検証セットであり、
2019535196.tif 000086
は
2019535196.tif 000087
が第i個検証ノードのサブセットを含まないことを示す。
2019535196.tif 000088
は検証情報であり、第j個変数ノードを除くほかの変数ノードの情報提供状況での対数尤度比を示す。
2019535196.tif 000089
は受信機の初期受信のチャンネル情報の対数尤度比表示形式である。
2019535196.tif 000090
は情報メッセージであり、第i個検証ノードを除くほかの検証ノードの情報提供状況での対数尤度比を示す。
2019535196.tif 000091
はすべての情報ビットの関連検証ノードの情報提供状況での対数尤度比を示す。
【０１３３】
本実施例の提供する重複多重に基づく復調装置は、実施例二の提供する重複多重に基づく復調方法と対応しているので、ここでその原理に対し再度説明しない。
【０１３４】
重複多重に基づく一種類の変調復調方法と装置、変調方法では、プリコーディング構成を採用し、発送側は入力情報シーケンスに対しまずパリティチェックプロダクトコードエンコードを行い、エンコード後の結果でエンコード規則により因子グラフを生成し、それから重複多重変調エンコードを行い、エンコード後の信号をアンテナで発射する。復調方法では、信号はチャンネルにより伝送され、受信側はアンテナにより信号を受信すると、まずデジタル信号処理を行う。それは同期、均一化等のプロセスを含む。その後、それに対し重複多重復調デコードを行い、最後にデコード後の結果に対し因子グラフ確率伝搬法デコードを行い、最終的にデコード後のシーケンスを得る。本発明では、プロダクトコードデコード方法を採用し、パリティチェックコードをサブコードとし、因子グラフ確率伝搬法の考えをデコード側に使う。それにより、パリティチェックプロダクトコードを採用し、その構成が簡単で、且つデコードプロセスで因子グラフを採用する方法を実現したので、運算複雑度を下げている。
【０１３５】
説明する必要のあるのは、本発明の実施例の提供する重複多重に基づく変調復調方法と装置はモバイル通信、衛星通信、マイクロウェーブ視界通信、散乱通信、大気圏光通信、赤外線通信、水中音響通信等の無線通信システムに応用でき、大容量の無線伝送にも小容量の小型無線システムにも応用できることである。
【０１３６】
本分野の技術者としては、上記実施方法で各種類の方法の全部又は一部のステップはプログラムで関連ハードウエアを制御することにより完成でき、当該プログラムはパソコン読取可能な保存媒質に保存できることが理解できる。保存媒質は、読取専用メモリ、ランダムアクセスメモリ、磁気ディスク又はCDなどを含む。
【０１３７】
上記説明内容は具体的な実施方法と結び合わせ、本発明をさらに詳しく説明したものであり、本発明の具体的な実施はこれらの説明に限られると見なしてはならない。本発明の所属技術分野の一般技術者にとっては、本発明の構想を逸脱しない前提で、また様々な演繹や差し替えを行うことができる。
【図１】
2019535196.tif 000119
【図２】
2019535196.tif 000120
【図３】
2019535196.tif 000121
【図４】
2019535196.tif 000122
【図５】
2019535196.tif 000123
【図６】
2019535196.tif 000124
【図７】
2019535196.tif 000125
【図８】
2019535196.tif 000126
【図９】
2019535196.tif 000127
【図１０】
2019535196.tif 000128
【図１１】
2019535196.tif 000129
【図１２】
2019535196.tif 000130
【図１３】
2019535196.tif 000131
【図１４】
2019535196.tif 000132
【図１５】
2019535196.tif 000133
【図１６】
2019535196.tif 000134
【図１７】
2019535196.tif 000135
【図１８】
2019535196.tif 000136
【図１９】
2019535196.tif 000137
【図２０】
2019535196.tif 000138
【図２１】
2019535196.tif 000139
【図２２】
2019535196.tif 000140
【図２３】
2019535196.tif 000141
【図２４】
2019535196.tif 000142
【図２５】
2019535196.tif 000143
【国際調査報告】
2019535196.tif 000144
2019535196.tif 000145
2019535196.tif 000146
2019535196.tif 000147