京大院情報学研究科 Graduate School of Informatics, Kyoto University 第29回情報理論とその応用シンポジウム 2006年11月29日 スパースな拡散符号を持つ ＣＤＭＡマルチユーザー復調器の解析 Analysis of Sparsely-Spread CDMA via Statistical Mechanics 吉田実加　　Mika Yoshida 奈良女大人間文化　Graduate School of Humanities and Sciences, Nara Women’s University 田中利幸　　Toshiyuki Tanaka　　　　 京大院情報学研究科　Graduate School of Informatics, Kyoto University 上江洌達也　Tatsuya Uezu

発表の流れ 目的 モデル化 レプリカ解析 　　　モデル 　　　レプリカ解析 　　　レプリカ対称解 数値計算 まとめと課題

目的 １．周波数ホッピング方式のモデル化を行う →スパースな拡散符号を持つ直接拡散ＣＤＭＡモデル としてモデル化する。 　 としてモデル化する。 ２．モデルの性能評価指標をレプリカ法を用いて導出する。

モデル化 周波数ホッピング（ＦＨ）ＣＤＭＡ方式を スパースな拡散符号を持つ直接拡散（ＤＳ）ＣＤＭＡモデル としてモデル化する。 直接拡散（ＤＳ）方式 　Direct Sequence ・・・ 大システム極限でのDS-CDMA モデル解析 D. Guo and S. Verdu, Randomly spread CDMA : Asymptotic via statistical physics, IEEE Trans. Info. Theory, vol. 51, no. 6, pp. 1983-2010 2005. T. Tanaka, A statistical-mechanics approach to large-system analysis of CDMA multiuser detectors IEEE Trans. Info. Theory, vol. 48, no. 11, pp. 2888-2910 2002.

周波数ホッピングＣＤＭＡ（ＦＨ－ＣＤＭＡ）方式 Frequency Hopping Code Division Multiple Access M値シフトキーイングで変調する場合 G-1 G ユーザー数 :　K (k=1,・・・,K) 拡散符号のチップ数 : G　 (g=1,・・・,G) キャリア周波数の数: M (m=1,・・・,M) ユーザーkの拡散符号 : {skg ; g=1,・・・,G } 加法的ノイズ：{ng; g=1,・・・,G}　 受信信号 : { yg; g=1,・・・,G }

チップ数をN=MGとして１次元的に書き直す。 ２次元で表される拡散符号を チップ数をN=MGとして１次元的に書き直す。 ユーザーkの拡散符号 :{skm ;μ=1,・・・,N } N個のチップのうち、G個が0以外の値を持つ。 このとき、受信信号は次のように表される。 M が大きいとき →　　拡散符号がスパースな K ユーザー DS-CDMAモデル スパースな拡散符号を持つＤＳ－ＣＤＭＡモデルは 時間ホッピングＣＤＭＡ方式の簡単なモデルにもなっている

レプリカ解析 スパースな拡散符号を持つＤＳ－ＣＤＭＡモデルにおける 復調結果の性能を評価する指標をレプリカ法を用いて 導出する。

モデル スパースな拡散符号を持つＤＳ－ＣＤＭＡモデル モデル　 スパースな拡散符号を持つＤＳ－ＣＤＭＡモデル ユーザー数 :　K (k=1,・・・,K)　　　　　　拡散符号のチップ数 : N　 (m=1,・・・,N) ユーザーkの情報シンボル : x0k ユーザーkの拡散符号 　　 : {skm ; μ=1,・・・,N } 加法的ノイズ　　　　　 : {nμ ; μ=1,・・・,N }　 受信信号 : { yμ; μ=1,・・・,N } 拡散符号系列　｛skm ; μ=1,・・・,N ｝ の中で0以外の値を持つチップの数は ρ個とする。 キャリア周波数の数がM、ホップ回数がGの時との対応は、N=MG, G=ρ。 このとき、受信信号は次のように表される。

仮定 拡散符号の0でないチップは、N個の中からランダムにr個選ばれるとする。 拡散符号のチップの値は以下の分布関数に従って与えられるとする。 ノイズは、平均0、分散σ02 の加法的白色ガウスノイズとする。 情報シンボルは、1と-1からなるビットとし、 それぞれの値が等確率で出現する一様分布に従って発生させるとする。 以下の表記を導入

x0 の条件のもとでの y の分布は次のようになる。 ガウスノイズを仮定しているので、 x0 の条件のもとでの y の分布は次のようになる。 復調の結果得られる、x0 の推定値を x とする。 x の条件のもとでの y の分布： s2　はガウスノイズの分散の推定値 事前分布

情報シンボル x と受信情報 yとの相互情報量 情報シンボルの復調の性能を評価するため 相互情報量を計算する。 情報シンボル x と受信情報 yとの相互情報量　 ここで p0(x ,y)=p0(y|x)p0(x) , p(x,y)=p(y|x)p(x) 相互情報量 I(x;y) は次式のように分解できる H(y) ：エントロピー H(y|x) ：条件付きエントロピー

ガウスノイズのエントロピー H(n) と等しいので次のように表せる。 条件付きエントロピーは ガウスノイズのエントロピー H(n) と等しいので次のように表せる。 yのエントロピーは 解析的に計算するのは困難。 自己平均性を仮定し 大システム極限で、拡散符号S について平均をとった １ユーザー当たりのエントロピーを計算する。

大システム極限 自己平均性の仮定 エントロピーH(y) の値は拡散符号S に依存してばらつく。 ばらつきは小さくなる。 ここで自己平均性を仮定すると　　　　　は大システム極限で S について平均をとった１ユーザー当たりのエントロピーHと一致する。 Es(・) ：拡散符号S についての平均

レプリカ解析 拡散符号S について平均をとった１ユーザー当たりのエントロピーHを 大システム極限で評価する。 解析計算を続けるために、レプリカ法を用いる。 恒等式 より、Hは次のように変形する。

ここで、Xn は n個のレプリカ変数　xa=(xa1,・・・,xaK)T, a=1,・・・,n を導入する。 N個の拡散符号チップ中、r個のみ0以外の値を持つ場合のXn は、次のようになる。 ただし

ここで

極限limK→∞ とlimn→0 の順番の入れ替えが可能であると仮定して、 limK→∞ を先に評価すると Varadhanの定理より、Xnは次のように表される。 R. S. Ellis, Entropy, Large Deviations, and Statistical Mechanics, ser. A series of comprehensive studies in mathematics. New York: SpringerVerlag, 1985, vol.271..

レプリカ対称解

鞍点方程式

ここで

１ユーザー当たりのy のエントロピーの平均H は

ビット誤り率 スペクトル効率 大システム極限のスペクトル効率

数値計算 拡散符号のチップの値は次の分布関数に従う バイナリの場合を考える。 p(s) を次のようにとる。

s02=s2でのビット誤り率 ：b=1, r=2 の解析結果の数値計算結果 ：b=1 (K=N=200), r=2, サンプル数50のモンテカルロシミュレーション

まとめと課題 まとめ 周波数ホッピング方式を拡散符号がスパースな直接拡散ＣＤＭＡモデルとして表した。 スパースな拡散符号を持つ直接拡散ＣＤＭＡモデルをレプリカ法を用いて解析し、レプリカ対称解におけるスペクトル効率とビット誤り率を導出した。 ビット誤り率について、解析結果から得られた数値解とモンテカルロ法を用いたシミュレーション結果を比較し、良く一致することが分かった。

今後の課題 スペクトル効率についても数値解とシミュレーション結果の比較を行う。 バイナリ以外の拡散符号の値を与えた場合について数値計算を行う。 FH方式のシミュレーションを行い、求めた解析結果との比較を行う。

前回（ISIT2006）との違い 前回：拡散符号の0でない要素の数を平均r個として、 事前分布がガウス分布に従う場合で、 　　　　事前分布がガウス分布に従う場合で、 近似(Effective Medium Approximation)計算で 　　　　性能評価指標を導出した。 今回： 拡散符号の0でない要素の数を厳密にr個として、 　　　　バイナリの一様な事前分布の場合について解析を行った。