• 検索結果がありません。

時空間ターボ符号化変調方式の符号間干渉通信路への適用

N/A
N/A
Protected

Academic year: 2021

シェア "時空間ターボ符号化変調方式の符号間干渉通信路への適用"

Copied!
5
0
0

読み込み中.... (全文を見る)

全文

(1)

時空間ターボ符号化変調方式の符号間干渉通信路への適用

兼田一幸

**

正本利行

***

荻原春生

****

A Space-Time Turbo Trellis Coded scheme with Decision Aided Equalizer for Intersymbol Interference

Kazuyuki KANEDA, Toshiyuki SHOHON and Haruo OGIWARA

1.はじめに

ディジタルデータの高速伝送の要求に伴って,無 線通信路でも多量のディジタルデータを伝送する必 要性が大きくなってきている。この無線通信路では,

マルチパス伝播に起因する符号間干渉が生じるため,

高速でディジタル伝送を行うと長いデータ伝送区間 に渡って特性劣化が現れる。

このような符号間干渉通信路に対して,

Lindskog

Paulraj [1]

は送信ダイバシチ効果

[2]

を得る特性 改善手法を提案している。しかしながら,この方式 では,ダイバシチ受信機において最尤系列推定を仮 定しているため,符号間干渉を受ける信号が長くな ると,その信号を最尤推定するための演算量が多く なり実現困難となる問題がある。そこで,この点を 解決するために,本論文では,この符号間干渉通信 路に判定補助等化器を用いた繰り返し等化器を用い ることを提案する。繰り返し等化器としては,ター ボ符号化変調方式(

Turbo-TCM

)を用い,送信デー タ系列の推定と干渉通信路の等化を同時に実現する と共に,伝送帯域を拡大させずにデータ伝送速度の 向上を図る。

この提案方式では一度誤りが生じるとその誤り がフィードバックし,さらなる誤りが生じる誤り伝 搬現象が発生し,誤り特性が改善しない問題が生じ た。そこで,この問題を解決するために,等化信号 のスムージングを行うことを提案し,誤り率特性を 改善する。

本論文ではまず,

2

節で

Lindskog

の送信ダイバ

シチを説明し,

3

節で送信ダイバシチを実現するタ ーボ符号化変調を用いた伝送システムを提案する。

4

節では,等化システムのスムージングを説明し,

5

節でそれらの特性評価を示す。

6

節は結論である。

2.

Lindskog

の送信ダイバシチシステム[1]

図1に

Lindskog

の送信ダイバシチシステムを示

す。まず,伝送信号

S(t)

は長さ

N

の2つの系列

s

1

(t)

s

2

(t)

に分けられる。分けられた情報はダイバシチ 効果を実現するため,表

1

の伝送手順に従い,

2

つ のアンテナからビット系列を伝送する。第一フレー ムでは,アンテナ1から

s

1

(t)

の信号を送信し,アン テナ2から

s

2

(t)

の信号を送信する。同様に,第

2

フ レームでは,アンテナ1から

-s

2

*(N-t)

の信号を送信 し,アンテナ

2

から

s

1

*(N-t)

の信号を送信する。

*

の 記号は複素共役を示す。

伝送した信号は符号間干渉通信路に入る。図

2

に その通信路モデルを示す。ここでは通信路のインパ ルス応答を

h

1

(t)

h

2

(t)

として記述しており,通信 路のインパルス応答の長さを両アンテナとも

q

個と している。また,雑音

n(t)

は白色ガウス雑音を仮定 している。

受信側でもフレーム構成を用い,

r

1 フレームでは 次式に従って受信信号

R

1

(t)

を作成する。

R

1

(t) = s

1

(t) h

1

(t) + s

2

(t) h

2

(t) +n

1

(t), (1)

ただし,

はたたみ込みを示す。また,

r

2 フレーム では

,

時間反転回路と複素共役回路を用いて,受信信 号

R

2

(t)

を以下のように作成する。

R

2

(t) = s

1

(t)  h

2*

(-t) - s

2

(t)  h

1*

(-t) + n

2

(t). (2)

但し,

h

1*

(-t)

h

2*

(-t)

は時間順序を逆順に入れ替え 複素共役を取った通信路のインパルス応答を示す。

*

原稿受付 平成

24

10

23

** 佐世保工業高等専門学校 電子制御工学科

*** 香川工業高等専門学校 通信ネットワーク工学科

****元長岡技術科学大学 電気系

(2)

これらの式をベクトル表記用いて以下のように表す ことができる。

R(t) =H s(t) +n(t), (3)

ここで,

R(t)

s(t)

及び

n(t)

は以下である。

R(t) =

 

 ) (

) (

2 1

t R

t

R , s(t)= 

 

 ) (

) (

2 1

t s

t

s , n(t)=

 

 ) (

) (

2 1

t n

t n

また,通信路行列

H

は次式を表現している。

H= 

 

 ) * ( )

* (

) ( )

(

1 2

2 1

t h t h

t h t

h .

マッチドフィルタ

H

Hはこの受信信号

R(t)

から次式 を用いて信号

~ s ( t )

を取り出す。

)

~ ( t

s = H

H

R(t) = H

H

H s(t) + H

H

n(t)

=[h

1

(t) h

1

*(-t)+h

2

(t) h

2

*(-t)] I s(t)+ v(t),(4)

但し,

~ s ( t )

H

H

v(t)

は次式となる。

)

~ ( t

s =  

 

) 2 (

~ ~ 1 ( ) t s

t

s , H

H

=

 

 

) 1 ( )

* ( 2

) 2 ( )

* ( 1

t h t h

t h t

h ,

v(t)= 

 

 ) (

) (

2 1

t v

t

v .

得られた信号

~ ( )

1

t

s

~ ( )

2

t

s

は最尤系列推定器に入 り,ここで最終推定信号系列

s ˆ

1

( t )

s ˆ

2

( t )

を推定す る。最後に,コンバイナで,送信信号と対応した最 終推定系列

S ˆ ( t )

を作成する。

H

Sep arat or

Time

H

Reversal S(t)

s

1

(t)

s

2

(t)

MLSE MLSE s

1

'(t)

s

2

'(t) r

1

or r

2

frame

R

1

(t) R

2

(t) r

1

or r

2

frame

Complex Conjugation

Complex Conjugation

Reversal Time

Complex Conjugation

Reversal Time

-1

s

2

(t) s

2

(t) s

1

(t) s

1

(t)

~

~

Antenna 1

Antenna 2

S(t)

C om bi ne r

S(t)

図 1 Lindskog の送信ダイバシチシステム

表1 伝送フレーム構成

Transmission order

0 1 2 3

Frame number

r

1

r

2

r

1

r

2

Antenna 1

s

1

( t ) -s

2*

( N - t ) s

1

( t ) -s

2*

( N - t )

Antenna 2

s

2

( t ) s

1*

( N - t ) s

2

( t ) s

1*

( N - t )

h

1

(t) n(t)

h

2

(t) s

1

'(t) s

2

'(t)

R(t)

図 2 干渉通信路モデル

3.Turbo-TCM を用いたダイバシチシステム

Lindskog

の送信ダイバシチシステムの最尤系列

推定器の代わりに,判定補助等化器と

Turbo-TCM

を用いた等化システムを適用する。図3にその判定 補助等化器を用いた

Turbo-TCM

等化システムを示 す。

ま ず 、 入 力 シ ン ボ ル

d(k)  {0, 1,…, 3}

Turbo-TCM

の畳み込み符号器に入り,ここで符号化

される。この畳み込み符号器は終端処理を行なって いる。符号化されたパリティビット

p(k)

は入力デー タ

d(k)

と一緒にシンボル

x(k)  {0, 1,

・・・

, 7}

に対 応付けられ,その後

8PSK

信号

S(k)

に写像される。

8PSK

信号は

Lindskog

の送信ダイバシチに入り,

2 節と同じ信号処理が行われる。この信号処理によ り信号

s

1

 ( t )

s

2

 ( t )

を作成し,

2

つのアンテナから 送信する。干渉通信路通過後の受信信号は

1

つのア ンテナで受信される。通信路のインパルス応答の値 は,独立なガウス変数を想定し,同じ分散の値を持 つと仮定する。さらに一つの伝送ブロック中ではそ のインパルス応答は変化しないと仮定する。受信信 号

R(t)

は送信ダイバシチの受信回路に入り,そこで,

信号

~ s

1

( t )

~ s

2

( t )

に変換され,コンバイナに入る。コ ンバイナは送信信号

S(t)

と対応した信号

S ~ ( t )

を出力 する。

次に,この信号

S ~ ( t )

Turbo

等化器に入る。初回 は,符号間干渉を削減するために,線形等化器(

LE

) を用いて等化を行う。

2

回目以後は,

Turbo-TCM

と 判定補助等化器(

DAE

)を組み合わせて,等化を行

(3)

いながら系列の推定を行う。それぞれの復号器は信 頼度として対数事後確率比を用いる。復号器1は受 信系列

y(k)

から次式の対数事後確率比

L

1

(x(n)|Y)

を計 算する。

 

 

 

 

Y p

Y Y p

| )) ( ), ( (

| ) ( )) ( ), ( log ( )

| ) (

1 ( P d k k 0

k x k k d k P

x

L , (5)

但し,

Y

は 受信信号

y(k)

の系列を示している。ま

た,同様に,復号器2はインタリーブされた受信信 号系列から対数事後確率比

L

2

(x(n)|Y)

を計算する。

判定補助等化器は,まず,この復号器の出力にお ける信頼度を用いて,符号間干渉のレプリカを作成 する。次に,このレプリカを受信信号から抜き去り,

新しい受信信号を作成する。そして,この受信信号 に対して,

Turbo-TCM

を用いて送信シンボル

d(k)

の 信頼度を再計算する。これを繰り返して行い,ター ボ等化器は最終的に送信信号

d(n)

を決定する。ただ し,初回の繰り返しでは,送信信号

d(k)

の信頼度が 得られないので,線形等化器を用いて受信信号

y(t)

の推定を行う。

8-PSK mapper

Enc 1

Enc 2

Decision

Lindskog's transmit diversity scheme Deinter

leaver

p(2)

(k)

p(1)

(k)

p(k) d(k)

S(k) Antenna 1

s1'(t)

Antenna 2 s2'(t)

H H

R1(t) R2(t) r1 or r2 frame

Time

*

Reversal ~

s~1(t)

s2(t)

C om bi ne r

DAE Dec

1

first

LE

else

Dec

2

Deinter leaver

inter leaver

inter leaver inter leaver Feedback

signal S(k)~

d(k) y(k)

Deinter leaver

Deinter leaver y(n) L1(x(k)|Y)

L2(x(k)|Y) ES(k)

図3 提案システム (Enc は畳み込み符号器,Dec は復号器,

LE は線形等化器,DAE は判定補助等化器を示す。)

D D D D ... D

(k)

...

ES

+

+

h

'

h

'

h

'

h

'

h

'

(k-1) (k-q)

(k+1) (k+q-1)

(k+q)

S(k) ~ y(k)

(k-q) (k-1)

(k+1) (k+q-1)

(k+q)

replica

ES ES

ES ES

ES

図4 判定補助等化器の構成

3.1 判定補助等化器

図4に判定補助等化器を示す。この判定補助等化 器は推定した信号

S ˆ ( k )

から符号間干渉を作成する。

シンボル

D

は遅延素子を示す。

h(k )

は式(4)の インパルス応答を示しており、次式で計算する。

(k )

h  = [h

1

(k) h

1

*(-k)+h

2

(k) h

2

*(-k)]. (6)

本論文では,これらのインパルス応答

h

1

(k)

h

2

(k)

は パイロット信号等を用いて既に分かっており、受信 機側で既知と仮定する。また,入力される伝送信号 の推定値

E

S

(k)

は以下の式で計算する。

 

 

 7

0 exp( ( ( ) | ))

7 0 exp( ( ( ) | )) ( )

) (

x L x k

x L x k Φ q

S k

E Y

Y , (7)

但し,

(q )

8PSK

信号とし,

L(x(k)|Y)

は対数事後 確率比を示す。符号化システムがパリティビットを 交互に用いているので,それに併せて式

(7)

L(x(k)|Y)

は以下の式を用いる。

L(x(k)|Y)=



 

number.

even is : )

| ) ( 2 (

number, odd

is

: )

| ) (

1

(

k k

x L

k k

x L

Y

Y (8)

3.2 スムージング

提案システムを用いて誤り率特性を求めた結果,

誤り率特性が繰り返しと共に悪くなることがある ことが分かった。これは,一度誤りが復号器で生じ ると,判定補助等化器で符号間干渉を抜き去る際に さらなる誤りが生じ,それを用いて復号を繰り返し 行うため、誤り伝搬が生じて特性が劣化するためと 考えられる。図5にその伝送ブロックに対する等化 信号の変動の様子を示す。この図では白丸で等化出 力における位相信号の変動を示している。この図か ら,位相の変動が繰り返しと共に次第に増加してい ることが分かる。

そこで,この変動に対して等化信号のスムージン グ処理を行うことを提案する。図6にその構成を示 す。このシステムでは,等化信号

y(t)

を既に得られ ている等化信号

y(t-1)

に加え,それを等化された信 号

y (t )

として用いる。この提案方式を用いたシミュ レーション結果を図5の黒丸で示す。この図より,

スムージングが位相変動を押さえ,推定された位相 が正しい位相値になることが分かる。

(4)

1 2 3 4 5 6 7 8 9 10

−3

−2

−1 0 1 2 3

Iteration number

Estimated phase [rad]

−

: Without smoothing : Smoothing : Correct phase

図5 等化システムの位相変動(送信位相が 0°の場合)

Decision

DAE LE

else

Trubo TCMdecoder

Feedback signal

S(k)~ d(k)

First Hard

Transmit diversity receiver

D

y(k) y(k)

図6 スムージング処理の構成

4.シミュレーション結果

この提案システムを用いて特性評価を行った結 果を図7に示す。シミュレーション条件としては,

干渉波の数を9とし直接波の数を1とした。また,

同じ長さの符号間干渉が各アンテナに生じるとし,

干渉波の電力も同じとした。符号器は

(11, 02, 04)

8

[5]

Ungerbock

符号とし,符号器内の疑似ランダム インターリーバの長さを

1023

とした。線形等化器 は,タップ数

22

個の

RLS

アルゴリズムを用いるも のとし,トレーニング長は

200

とした。この

RLS

アルゴリズムにおいては,等化誤差を通信路の雑音 の増加と捉えることにし,計算機シミュレーション では等化の平均

2

乗誤差を計算しそれを雑音の電 力に加えた。

図7に提案システムの誤り率特性を実線で示す。

また,図7の点線はスムージングを用いない場合の 特性を示す。この図から,スムージングを用いるこ とで特性が改善できていることが分かる。

10

回の 繰り返しで比較すると,

10

-5において,スムージン グを用いない場合に比べて平均

CNR

の改善量は6 d

B

であった。

0 10 20

10

−8

10

−6

10

−4

10

−2

10

0

ISI 10(direct+9delay waves) Training length = 200 symbols

Iteration = 1 time 2 times 3 times 4 times 10 times

Average CNR[dB]

BER

図7 提案 Turbo-TCM システムのビット誤り率特性

次にダイバシチ効果を検討する。図8にその結果 を示す。図8の実線は受信アンテナを

2

本用いたダ イバチシ構成の誤り率特性を示し,破線は1つのア ンテナを用いた場合の誤り率特性を示す。シミュレ ーション条件は図7の条件と同じとしている。この 図より、得られたダイバシチ効果は大きく,

10

回の 繰り返しで比較すると、誤り率

10

-4において,平均

CNR

が6d

B

ほど改善できていることが分かる。

次に,初回の線形等化器で用いる

RLS

アルゴリズ ムのトレーニング長の効果を検討する。図

9

にその トレーニングの長さの変化に対する誤り率特性をシ ミュレーションで求めた結果を示す。この図では平 均

CNR

12d

B

で評価を行った。シミュレーション 条件は図7と同じとした。このシミュレーションで は,伝送ブロック毎に毎回トレーニングを行い,そ の都度タップの更新を行っている。この図から,ト レーニング長さ

100

以上で特性が大きく改善してお り,この程度の長さのトレーニングが必要であるこ とが分かる。伝送ビット長が

1023

なので,トレーニ ング長さ

100

は全体の約

10

%となる。長い伝送ブロ ックを用いれば,トレーニング長に対する全体の割 合はさらに小さくなりスループットの向上が期待で きる。

(5)

0 10 20 10

−8

10

−6

10

−4

10

−2

10

0

ISI 10(direct+9delay waves) Training length = 200 symbols

1 time 2 times 3 times 4 times 10 times

Average CNR[dB]

BER

Iteration = diversity

No diversity With

8 ダイバシチの有無に対する特性比較

0 100 200

10

−8

10

−6

10

−4

10

−2

10

0

BER

Training length of RLS algorithm[symbol]

Avarege CNR=12[dB]

iteration = 1 time

iteration = 2 times iteration = 3 times iteration = 10 times

図9 トレーニング長さに対する誤り率特性

5.まとめ

本論文では,長い符号間干渉通信路に対する送信 ダイバシチを実現するために,

Lindskog

の送信ダ イ バ シ チ シ ス テ ム に 判 定 補 助 等 化 器 を 用 い た

Turbo-TCM

等化器を適用することを提案した。

Lindskog

の送信ダイバシチは最尤系列推定器を仮

定しており,干渉が長くなればなるほど,提案した 本方式の方が装置規模を小さく実現できる。また、

提案システムでは誤り伝搬が生じたため,等化信号 のスムージング処理を用いることを提案した。その 結果,繰り返し毎に特性を大きく改善することがで

きた。さらに,ダイバシチ効果についても検討し,

提案システムを用いることで,誤り率

10

-4で,平均

CNR

6dB

ほど特性を改善できることを確認した。

最後に,線形等化器のトレーニング長の効果を検討 した。その結果,長さ

100

程度で特性が大きく向上 することを確認した。この長さは伝送長さが長くな ればさらに受け入れられるものになる。

参考文献

1) E. Lindskog and A. Paulraj, “A transmit diver- sity scheme for channels with intersymbol in- terference,” in IEEE Int. Conf. commun.

(ICC2000), no. 1, 307-311, June 2000.

2) S. B. Alamounti, “A simple transmit diversity technique for wireless communications,” IEEE Journal on Select Areas in Communication, vol. 16, no. 8, pp. 1451-1458, Oct. 1998.

3) Z.-N. Wu and J. M. Cioffi, “Low-complexity iterative decoding with decision-aided equali- zation for magnetic recording channel,” IEEE J.

Selected Areas Commun., vol. 19, no. 4, pp.

699-708, Apr. 2001.

4) P. Robertson and T. W¨orz, “Coded modulation scheme employing turbo code,” Electronics Letters, vol. 31, no. 18, pp.1546-1547, 1995.

5) G. Ungerbock, “Channel coding with multi-

level/phase signals,” IEEE Transactions on In-

formation Theory, vol. IT-28, pp. 56-67, Jan

1982.

参照

関連したドキュメント

[r]

あらまし MPEG は Moving Picture Experts Group の略称であり, ISO/IEC JTC1 におけるオーディオビジュアル符号化標準の

In this paper, we consider the discrete deformation of the discrete space curves with constant torsion described by the discrete mKdV or the discrete sine‐Gordon equations, and

古物営業法第5条第1項第6号に規定する文字・番号・記号 その他の符号(ホームページのURL)

アクセス道路の多重化・道路の補強 工事中 通信設備の増強(衛星電話の設置等) 完了 環境モニタリング設備等の増強・モニタリングカーの増設 完了 高台への緊急時用資機材倉庫の設置※

現行の HDTV デジタル放送では 4:2:0 が採用されていること、また、 Main 10 プロファイルおよ び Main プロファイルは Y′C′ B C′ R 4:2:0 のみをサポートしていることから、 Y′C′ B

システムの許容範囲を超えた気海象 許容範囲内外の判定システム システムの不具合による自動運航の継続不可 システムの予備の搭載 船陸間通信の信頼性低下

16 単列 GIS配管との干渉回避 17 単列 DG連絡ダクトとの干渉回避 18~20 単列 電気・通信ケーブル,K排水路,.