水中光無線通信方式の多段中継による大容量化に関する研究
代表研究者

松 元 隆 博山口大学大学院創成科学研究科 准教授

1 はじめに
四方を海に囲まれた島国の日本にとって，海洋資源の調査や開発，発掘は重要な課題となっている．現在，有人潜水調査艇等による調査，海底津波計や海底地震計等のセンサによる計測などがなされており，機器間もしくは機器と母船等への通信は，有線，もしくは超音波によって行なわれている．有線通信では敷設作業や保守管理が面倒であり，断線事故の発生が懸念される等の問題がある．一方，超音波通信は有線通信での問題を解決でき，長距離通信も可能であるが，伝搬路である水中の特性から帯域を広くとれず，音声通信程度の低速通信しか出来ないことや海中温度変化がもたらす音速分布の不均一性による影響などで通信安定性に問題がある．従って，今後の海洋開発で必要となる海底のリアルタイムモニタリング等の高速通信や多くの各種計測データを同時に伝送するための多重通信に対応できないという問題を抱えている．これらの問題解決のため，波長 380-780[nm]の可視光による無線通信が近年注目されているが，長距離通信のためにレーザーダイオードで実現する必要があることから機器の設置期間が長くなる海中調査での利用では消費電力が過大であることが課題となっている．
そこで，本研究では近距離の通信しか出来ない省電力の LED を用いた海中での光無線通信を多段中継伝送させることで，長距離伝送を可能とする新しい省電力の通信方式の実現可能性について検討を行う．ここでは更に，中継途中に並列に配置した複数の中継器を複数のアンテナに見立てた MIMO 通信技術と CDMA 方式を適用することで，受信時の SN 比を向上させ BER 特性を改善させることに加え，CDMA 方式を適用することで多重通信を実現する．

2 海中光学特性
一般に光がある媒体中を伝播するとき，その強度は散乱と吸収によって減少する．海中での散乱は，分子や粒子が光子の方向を変化させることであり，また，海中での吸収は，分子や粒子の振動による熱などの別のものに変換されることである．これらは，クロロフィルと黄色の物質を含む水分子と溶解した粒子によって引き起こされる[1],[2]．ここで，減衰係数c(Cc , λ)は，以下の式で求められる．
c(Cc , 𝜆) = 𝑎(𝐶𝑐 , 𝜆) + 𝑏(𝐶𝑐 , 𝜆)但し，Cc はクロロフィルの濃度，λは波長，𝑎(𝐶𝑐 , 𝜆)は吸収係数，𝑏(𝐶𝑐 , 𝜆)は散乱係数である.クロロフィルの濃度は水質によって異なり，特徴的な環境として，深海のような純粋な海水，透明な海水，沿岸海域の 3 つが挙げられる．例として，波長λが 550[nm]のときのそれぞれの環境におけるクロロフィルの濃度Cc と吸収係数𝑎(𝐶𝑐 , 𝜆)，散乱係数𝑏(𝐶𝑐 , 𝜆)，減衰係数c(Cc , 𝜆)を表 1 に示す[2]．
表 1．異なる水域におけるクロロフィル濃度Cc 、吸収係数𝑎(𝐶𝑐 , 𝜆)、散乱係数𝑏(𝐶𝑐 , 𝜆)、減衰係数c(Cc , 𝜆)の例
𝑎(𝐶𝑐 , 𝜆)[𝑚−1 ]𝑏(𝐶𝑐 , 𝜆)[𝑚−1 ]c(Cc , 𝜆)[𝑚−1 ]水質
Cc [mg/m3 ]
純粋な海水
0.005
0.0530.0030.056
透明な海水
0.310.069
0.080.15
沿岸の海水
0.83
0.088
0.2160.305

また，図 1 にそれぞれの環境における波長𝜆と減衰係数c(Cc , 𝜆)の関係を示す．次の節に，具体的に吸収係数𝑎(𝐶𝑐 , 𝜆)と散乱係数𝑏(𝐶𝑐 , 𝜆)について詳しく述べる．
2-1 吸収係数
吸収係数𝑎(𝐶𝑐 , 𝜆) はクロロフィルの濃度とフルボ酸、フミン酸に基づいてモデル化され，以下の式で表される．

1図 1．波長λに対する各環境での減衰係数c(Cc , 𝜆)𝑎(𝐶𝑐 , 𝜆) = 𝑎𝑤 (𝜆) + 𝑎𝑐 (𝜆)𝐶𝑐′ (𝐶𝑐 ) + 𝑎𝑓 𝐶𝑓 (𝐶𝑐 ) exp �𝑘𝑓 (𝜆)� + 𝑎ℎ 𝐶ℎ (𝐶𝑐 ) exp�𝑘ℎ (𝜆)�但し，𝑎𝑤 (𝜆)は波長𝜆における純水での吸収係数で，図 2 の関係がある[3]．また，𝑎𝑐 (𝜆)は波長𝜆におけるクロロフィルの吸収係数で文献[4]で得られ，𝑎𝑓 はフルボ酸の吸収係数で，𝐶𝑓 (𝐶𝑐 )はクロロフィルの濃度がCc のときのフルボ酸の濃度で，𝐶𝑓 (𝐶𝑐 ) = 1.74098Cc exp(0.12327𝐶𝑐 )で表される．同様に，𝑎ℎ はフミン酸の吸収係数，𝐶ℎ (𝐶𝑐 )はクロロフィルの濃度がCc のときのフミン酸の濃度で，𝐶ℎ (𝐶𝑐 ) = 0.19334𝐶𝑐 exp(0.12343𝐶𝑐 )で表される．また，𝐶𝑐′ (𝐶𝑐 )，𝑘𝑓 (𝜆)，𝑘ℎ (𝜆)はそれぞれ，𝐶𝑐′ (𝐶𝑐 ) = 0.06|𝐶𝑐 |0.0602
𝑘𝑓 (𝜆) = 0.01890(450 − 𝜆)𝑘ℎ (𝜆) = 0.01105(450 − 𝜆)である．ここで，文献[5]によればフルボ酸の吸収係数は𝑎𝑓 = 0.00728，フミン酸の吸収係数は𝑎ℎ = 0.1304である．

2-2 散乱係数
散乱係数𝑏(𝐶𝑐 , 𝜆)は純粋の散乱と小さな粒子の散乱，大きな粒子の散乱によりモデル化され，以下の式で表される．
𝑏(𝐶𝑐 , 𝜆) = 𝑏𝑤 (𝜆) + 𝑏𝑠 (𝜆)𝐶𝑠 (𝐶𝑐 ) + 𝑏𝑙 (𝜆)𝐶𝑙 (𝐶𝑐 )但し，𝑏𝑤 (𝜆)は波長𝜆における純水での散乱係数で，𝑏𝑠 (𝜆)は波長𝜆における小さな粒子による散乱係数で，𝑏𝑙 (𝜆)は波長𝜆における大きな粒子による散乱係数で，400 4.322
�𝜆400 1.7�𝑏𝑠 (𝜆) = 0.341074 �
𝜆400 0.3�
𝑏𝑙 (𝜆) = 0.341074 �𝜆で表される．また，𝐶𝑠 (𝐶𝑐 )はクロロフィルの濃度がCc のときの小さな粒子の濃度で，𝐶𝑙 (𝐶𝑐 ) はクロロフィルの濃度がCc のときの大きな粒子の濃度で，𝑏𝑤 (𝜆) = 0.005826 �2図 2．波長λに対する純水での吸収係数𝑎𝑤 (𝜆)で表される．

𝐶𝑠 (𝐶𝑐 ) = 0.01739𝐶𝑐 exp(0.11631𝐶𝑐 )𝐶𝑙 (𝐶𝑐 ) = 0.76284𝐶𝑐 exp(0.03092𝐶𝑐 )３ 海中での背景光雑音
海中での背景光雑音は，太陽光による影響が大きい．地表での太陽光の放射照度を求めるモデルとして，
1984 年に Bird 氏によって SPCTRAL2 が提案されている[6]．このモデルでは，ある時間において水平面，または傾斜平面上に澄んだ空での太陽の光の照度を計算できる．また，受信機平面が傾斜している場合は入射角，または平面の傾きと方位角を指定することで照度を計算できる．
受信機に入射する背景光雑音のパワーは，
𝐸𝑡𝑡𝑡𝑡𝑡 ⋅ 𝑅 ⋅ 𝐴 ⋅ 𝐿 ⋅ 𝑒𝑒𝑒(−𝑐(𝐶𝑐 , 𝜆)𝐷)𝑃𝑏 =𝜋で求められる．但し，𝐸𝑡𝑡𝑡𝑡𝑡 は水面に入射する全波長の太陽光の照度で，波長𝜆における照度を𝐸(𝜆)とすれば，Etotal = ∫ 𝐸(𝜆)𝑑𝑑で求められる．また，Rは太陽光の照度の水中反射率で 1.25%，Aは受光面積，Dは受信機の深さで，Lは水中の放射照度の方向依存性で，受信機の受光面が上を向いて水面に対して平行の場合は 1 で，受光面が横を向いて水面に対して垂直の場合は 0.026，受光面が下を向いて水面に対して平行の場合は 0.0037 である．
また，背景光の光子数の平均値λb は，
ηPbλb =ℎ𝑓で表される．但し，ηは量子効率，ℎはプランク定数，𝑓は光の周波数である．

４ 光 ZCZ 系列セットを用いた多段中継伝送
4-1 光 ZCZ 系列セット
Zero-correlation zone(ZCZ)系列セット[7][8]とは，周期相関特性に零相関区間(zero-correlation zone)を持つ系列のセットである．この零相関区間とは，相関特性のシフト 0 の前後のある限られた区間で自己相関特性がインパルス特性で，相互相関特性が完全に 0 となる区間である．光 ZCZ 系列セットはそのような相3

(a) 周期自己相関特性(b) 周期相互相関特性
図 3．光 ZCZ 系列の周期相関特性

関特性を持つ系列の一つで，光の on と off で情報を伝送する光 code-division multiple access(CDMA)方式に適用できるように，送信側で用いる 1 と 0 の要素からなる 2 値系列と，受信側で用いる 1 と-1 の要素からなる 2 相系列の対からなる系列セットである．
𝑗

長さ𝑁で，1 と-1 の要素からなる 2 相系列𝑎𝑁 を，
j𝑗𝑗𝑗𝑗

𝑎N = (𝑎𝑁,0 , 𝑎𝑁,1 , ⋯ , 𝑎𝑁,𝑖 , ⋯ , 𝑎𝑁,𝑁−1 ),𝑗𝑎𝑁,𝑖 ∈ {1, −1}

j,dと表記する．
但し，𝑗は系列番号で，𝑖は順序変数である．同様に，長さ𝑁で 1 と 0 の要素からなる 2 値系列𝑎�N を，
j,d𝑗,𝑑𝑗,𝑑𝑗,𝑑𝑗,𝑑𝑎�N = (𝑎�𝑁,0 , 𝑎�𝑁,1 , ⋯ , 𝑎�𝑁,𝑖 , ⋯ , 𝑎�𝑁,𝑁−1 ),𝑗,𝑑𝑎�𝑁,𝑖 ∈ {1, 0},d ∈ {1, 0}と表記する．これらの系列の対の集合𝐴を

𝑗𝑗,𝑑1,𝑑2,𝑑𝑀,𝑑2𝑀
�, �𝑎𝑁, 𝑎�𝑁�, ⋯ , �𝑎𝑁 , 𝑎�𝑁 �, ⋯ , �𝑎𝑁, 𝑎�𝑁�}
𝐴 = {�𝑎1𝑁 , 𝑎�𝑁とする．但し，𝑀は系列数である．
jj′,dこのとき，2 相系列𝑎N と 2 値系列𝑎�N の周期相関関数𝜌 𝑗𝜌 𝑗と定義したとき，この周期相関関数𝜌 𝑗𝑗′ ,𝑑𝑎𝑁 ,𝑎�𝑁 ,𝑖 ′𝑗′ ,𝑑𝑎𝑁 ,𝑎�𝑁 ,𝑖 ′𝑁−1𝑗′ ,𝑑𝑎𝑁 ,𝑎�𝑁 ,𝑖 ′𝑗を，j′ ,d= � 𝑎𝑁,𝑖 𝑎�𝑁,�𝑖+𝑖 ′ � mod Nが，𝑖=0𝑤 ; 𝑖 ′ = 0, 𝑗 = 𝑗 ′ , 𝑑 = 0,−w ; i′ = 0, j = j′ , d = 1,
𝜌 𝑗 𝑗′,𝑑 ′ = �𝑎𝑁 ,𝑎�𝑁 ,𝑖0 ; i′ = 0, j ≠ j′ ,0; 1 ≤ |i′ | ≤ 𝑍𝑍𝑍j,dとなるとき，この集合𝐴を光 ZCZ 系列セットと呼ぶ．但し，𝑍𝑍𝑍は零相関区間のサイズ，𝑤は 2 値系列𝑎�N の系
4

図 4．中継伝送の基本モデル
列要素 1 の数に等しく，
N−1j,d𝑤 = � 𝑎�N ≤ Ni=0の関係を満たす．ここでは，j = j′のときの周期相関関数を周期自己相関関数と，j ≠ j′のときの周期相関関数
を周期相互相関関数と呼ぶことにする．図 3 に周期自己相関関数と周期相互相関関数を示す．このように，
シフト 0 周辺の前後Zczの範囲が，理想的な自己相関特性，相互相関特性になっている．
4-2 光 ZCZ 系列セットを用いた多段中継伝送システムの構成
海中では，前述した光の吸収と散乱によって，光のパワーが減衰し，長距離伝送が困難である．そこで，途中に送受信機能を有する中継機を配置し，距離を延すことを考える．本稿では基本的な中継方式としてAF(Amplify and Forward)方式を適用した場合で検討を行う．図 4 に中継伝送の基本モデルを示す．送信機でj,d2 値の光 ZCZ 系列𝑎�N によって拡散変調された信号は，LED にて電気信号から光信号に変換され送信される．
j,0j,1

この際，ビットデータがd = 0の場合は𝑎�N を，d = 1の場合は𝑎�N をそれぞれ送信される．各中継局では

avalanche photodiode (APD)で変換された電気信号は，増幅され，再度 LED で光信号に変換され送信される．
j

これらの複数の中継機を何段も経て，受信機にて APD で変換された電気信号を 2 相の光 ZCZ 系列𝑎N によって逆拡散され，相関出力を正負で閾値判定し復調される．

図 5．送信機と中継機，受信機の配置例

5次に，図 5 に送信機と中継機，受信機の配置の例を示す．この例では中継数Nr が 2，並列数Np が 3 の場合である．以後，中継数Nr ，並列数Np の場合の場合を(Nr , 𝑁𝑝 )と記述する．矢印が各送信機，中継機からの光が届く経路を示している．これは LED の指向性によるもので，以後の性能評価では正面の中継機の上下 1 機の中継機，もしくは受信機に届くこととする．中継機は各段で並列に配置し，それを多段に中継することによって，経路数を増やし S/N 比を改善させる．

５ 計算機シミュレーションによる性能評価
海中光無線通信に前述の多段中継伝送を行った場合のビット誤り率（BER）特性を計算機シミュレーションで評価する．計算機シミュレーションの諸元を表 2 に示す．
表 2．計算機シミュレーションの諸元
拡散系列
系列長 N
零相関区間Zcz
試行回数
ビットレート[Mbps]
チップレート[Mcps]
光の波長λ[nm]APD 量子効率η
APD 利得GAPD イオン化係数K eff
APD 口径[mm]
バルク漏れ電流Ib [𝑛𝑛]表面漏れ電流Is [𝑛𝑛]受信機雑音温度Tr [𝐾]
受信機負荷抵抗R L [Ω]クロロフィル濃度Cc [𝑚𝑚/𝑐𝑚3 ]送信機と受信機の距離
最大水深
背景光モデル
経度[°]
緯度[°]
時差時間(hh:mm:ss)
放射輝度の方向依存性

光 ZCZ 系列
161106348
青 450，緑 530，赤 650
0.61000.021.50.11011001030
0.005, 0.31, 0.83
3030
SPCTRL2155
1910.33
12:00:00
1.0(下方向)，0.037（上方向）5-1 中継機の数Nr による性能評価
最初に，透明な海水の環境、つまりクロロフィルの濃度Cc が0.31[𝑚𝑚/𝑐𝑚3 ]の場合に，海底から海面に向けて上向けに通信した場合の中継数の違いによる BER 特性の比較を行う．ここでは，海底 30[m]の場所に送信機を，海面に受信機を設置し，その間の中継機の数Nr を変化させる．並列数Np が 1, 3, 5 の場合の BER 特性をそれぞれ図 6, 図 7, 図 8 に示す．
これらより，光の波長（色）に関わらず中継機の数Nr を増加させることで BER 特性を向上できることが確認できる．これは必要な BER を満たすための伝送距離が延ばせることを意味する．また，並列数Np に関わらず緑色（530[nm]），青色（450[nm]），赤色（650[nm]）の順で BER 特性が良いことがわかる．これは図 1 の波長λに対する減衰係数c(Cc , 𝜆)による影響と考えられる．
5-2 中継機の数Nr と並列数Np による性能評価
次に，5-1 の中継機の数Nr と並列数Np をどちらとも変化させた場合の最適な組み合わせについて検討を行う．送信機と受信機，中継機の配置，及び透明な海水の環境（クロロフィルの濃度Cc が0.31[𝑚𝑚/𝑐𝑚3 ]）は
6

図 6．並列数Np が 1 の場合の BER 特性(Cc = 0.31) 図 7．並列数Np が 3 の場合の BER 特性(Cc = 0.31)

図 8．並列数Np が 3 の場合の BER 特性(Cc = 0.31) 図 9．赤色（650[nm]）の場合の BER 特性(Cc = 0.31)

図 10．緑色（530[nm]）の場合の BER 特性(Cc = 0.31) 図 11．青色（450[nm]）の場合の BER 特性(Cc = 0.31)

5-1 と同様である．赤色（650[nm]），緑色（530[nm]），青色（450[nm]）の LED を使用した場合の BER 特性をそれぞれ図 9，図 10，図 11 に示す．
これらより，光の波長（色）に関わらず中継機の数Nr と並列数Np を増加させることで BER 特性を向上できることが確認できる．これは 5-1 と同様，必要な BER を満たすための伝送距離が延ばせることを意味する．
また，各光の波長（色）の中継機の数Nr と並列数Np の組み合わせは，(Nr , 𝑁𝑝 ) = (4,3)と(3,5)の場合は，緑色（530[nm]）と青色（450[nm]）は BER 特性がほとんど同じであるが，赤色（650[nm]）では(Nr , 𝑁𝑝 ) = (3,5)に比べて，(Nr , 𝑁𝑝 ) = (4,3)の方が BER 特性が向上できることがわかる．これは図 1 の波長λに対する減衰係数c(Cc , 𝜆)より，緑色（530[nm]）と青色（450[nm]）に比べ，赤色（650[nm]）が減衰係数c(Cc , 𝜆)が大きいため，伝送距離の増加による減衰を抑えるために中継機の数Nr を増加させる方が効果が大きいからだと考えられ
る．
7

5-3 クロロフィル濃度による性能評価
次に，クロロフィル濃度の違いによる影響について検討を行う．送信機と受信機，中継機の配置は 5-1，5-2 と同様である．また，中継機の数Nr と並列数Np の組み合わせは，(Nr , 𝑁𝑝 ) = (4,3)と(3,5)，(4,5)の 3 通りでそれぞれ比較している．赤色（650[nm]），緑色（530[nm]），青色（450[nm]）の LED を使用した場合の BER特性をそれぞれ図 12，図 13，図 14 に示す．

図 12．赤色（650[nm]）の場合の BER 特性

図 13．緑色（530[nm]）の場合の BER 特性

図 14．青色（450[nm]）の場合の BER 特性

図 15．赤色（650[nm]）で上下方向の場合の BER 特性

これらより，クロロフィル濃度が小さい，すなわち純粋な海水ほど BER 特性が向上できることが確認できる．また，中継機の数Nr と並列数Np の組み合わせ(Nr , 𝑁𝑝 ) = (4,3)と(3,5)の場合で比較した場合，図 12 より赤色（650[nm]）の場合はクロロフィルの濃度によらず(Nr , 𝑁𝑝 ) = (3,5)より，(4,3)の方が BER 特性が向上す
るが，図 13 と図 14 より緑色（530[nm]）と青色（450[nm]）の場合はクロロフィル濃度が大きいCc = 0.83の場合は赤色（650[nm]）と同様，(Nr , 𝑁𝑝 ) = (3,5)より，(4,3)の方が BER 特性が向上するが，クロロフィル濃度が小さいCc = 0.005の場合は赤色（650[nm]）とは逆に(Nr , 𝑁𝑝 ) = (4,3)より，(3,5)の方が BER 特性が向上する．これは赤色（650[nm]）の場合は，クロロフィル濃度に関係なく，減衰が大きいため，BER 特性を向上させるためには中継機の数Nr を増やしてその影響を抑えた方が良く，緑色（530[nm]）と青色（450[nm]）の場
合はクロロフィル濃度が小さい場合は減衰の影響が比較的小さく，逆に背景光の影響が大きくなるために並列数𝑁𝑝 を増やしてその影響を抑えた方がいいからだと考えられる．
5-4 通信方向の違いによる性能評価
最後に，透明な海水の環境、つまりクロロフィルの濃度Cc が0.31[𝑚𝑚/𝑐𝑚3 ]の場合の，通信方向の違いによる影響について検討を行う．送信機と受信機，中継機の配置が 5-1，5-2, 5-3 と同様に上方向の場合と，それとは逆に海面の送信機を，海底 30[m]の場所に受信機を設置し，その間の中継機の数Nr と並列数Np を変化させる．ここでは，中継機の数Nr と並列数Np の組み合わせは，(Nr , 𝑁𝑝 ) = (4,3)と(3,5)，(4,5)の 3 通りでそれぞれ比較している．赤色（650[nm]），緑色（530[nm]），青色（450[nm]）の LED を使用した場合の BER 特性をそれぞれ図 15，図 16，図 17 に示す．

8

図 16．緑色（530[nm]）で上下方向の場合の BER 特性 図 17．青色（450[nm]）で上下方向の場合の BER 特性
これらより，海面から海底への下方向への通信に比べ，海底から海面への上方向の通信の方が BER 特性が良いことが確認できる．これは下方向の場合は中継機と受信機の受光面が太陽を向き，それによる背景光の影響が大きいためである．また，中継機の数Nr と並列数Np の組み合わせ(Nr , 𝑁𝑝 ) = (4,3)と(3,5)の場合で比較した場合，図 15 より赤色（650[nm]）の場合は通信の方向によらず(Nr , 𝑁𝑝 ) = (3,5)より，(4,3)の方が BER 特性が向上するが，図 16 と図 17 より緑色（530[nm]）と青色（450[nm]）の場合は上方向であれば，赤色（650[nm]）と同様，(Nr , 𝑁𝑝 ) = (3,5)より，(4,3)の方が BER 特性が向上するが，下方向の場合，赤色（650[nm]）とは逆に(Nr , 𝑁𝑝 ) = (4,3)より，(3,5)の方が BER 特性が向上する．これは減衰の影響が比較的小さく，逆に背景光の影響が大きい場合，並列数𝑁𝑝 を増やしてその影響を抑えた方がいいからだと考えられる．

６ おわりに
本研究では，近距離の通信しか出来ない省電力の LED を用いた光 CDMA 方式に多段中継伝送を適用することで，長距離伝送を可能とする新しい省電力の通信方式を提案し，その通信性能を計算機シミュレーションによって評価した．本方式によれば，水質や使用する LED の色（波長）に関係なく，通信距離を延ばせることが明らかになった．また，中継機の配置は，水質による減衰による影響が大きいときは中継数を，背景光による影響が大きい場合は並列数を増加させることで BER 特性を改善できると考えられる．

【参考文献】[1] Vladimir I. Haltrin， “Chlorophyll-Based Model of Seawater Optical Properties,” Appl. Opt. 38，
vol. 38， no. 33， pp. 6826-6832, Nov. 1999.
[2] C. Gabriel， M.Khalighi， S. Bourennane， P. Leon， V. Rigaud， “Monte-Carlo-Based Channel
Characterization for Underwater Optical Communication System,” J. Optical Communications
and Networking， vol. 5 ， no.1， p.1-12， Jan. 2013.
[3] R. M. Pope, E. S. Fry， “Absorption spectrum(380-700nm) of pure water.II. Integrating cavity
measurements,” Applied Optics， vol.36， no.33， pp.8710-8723， Nov.1997.
[4] L. Prieur, S. Sathyendranath， “An optical classificantion of coastal and oceanic waters based on
the specific spectral absorption curves of phytoplankton pigments, dissolved organic matter, andother particulate materials,” Limnol. Oceanogr， vol. 26，no. 4， pp. 671- 689， Jan. 1981.
[5] K. L. Carde， R.G. Stewart， G. R. Harvey， and P. B. Ortner， “Marine humic and fulvic acids:
their effects on remote sensing of ocean chlorophyll，” Limnol. Oceanooogr， vol. 34，no. 1， pp.
68-81， 1989.
[6] R. E. Bird and C. Riordan,” Simple solar spectral model for direct and diffuse irradiance onhorizontal and rilted planes at the earth’s surface for cloudless atmospheres,” Journal of Climate
and Applied Meteorogy, vol. 25, no. 1, pp. 87- 97, Jan. 1986.

9

[7] S. Matsufuji, T. Matsumoto, Y. Tanada, N. Kuroyanagi, " ZCZ codes for ASK-CDMA system, "IEICE Trans. on Fundamentals of Electronics, Communications and Computer Sciences, vol. E89A, no.9, pp. 2268-2274, Sept. 2006.
[8] Y. Ohira, T. Matsumoto, H. Torii, Y. Ida, S. Matsufuji, " Code generator for an optical ZCZ
sequence with zero-correlation zone 2z, " Journal of Signal Processing, Vol. 21, No. 6, pp. 265-272,
Nov. 2017.

〈発題名Code Generator for an Optical ZCZ
Sequence with Zero-Correlation Zone 2z
BER Performance of SS System Using a
Huffman Sequence Against CW Jamming
A Compact Matched Filter Bank for an
Optical ZCZ Sequence Set withZero-Correlation Zone 2z
Modified Mutually ZCZ Set of Optical
Orthogonal Sequences
A Generation Method of a Mutually ZCZSet of Optical Orthogonal Sequences
Consisting of Binary and Ternary
Sequence Pairs
Construction of a Code Generator for anOptical ZCZ Sequence with ZeroCorrelation Zone 2z
Experimental Evaluation of Optical
ZCZ-CDMA System on the Move
Study on BER Performance of OpticalWireless Code-Shift Keying System Using
Compact Matched Filter Bank for OpticalZCZ Sequence Set
Construction of Code Generator for
Optical ZCZ Sequence Used as
Synchronization Symbol
A Generation Method of a Modified
Mutually ZCZ Set of Optical Orthogonal
Sequences
Influence of ISI on Synchronous OCDMA
System Using an Optical ZCZ Sequence
Set on a Diffuse Channel
Theoretical Analysis of BER Performanceof OCDMA System Using CompactMatched Filter Bank for an Optical ZCZ
Sequence Set
A Generation Method of a TwoDimensional Optical ZCZ Sequence with
the Smallest Zero-Correlation Zone

表資料〉

掲載誌・学会名等

発表年月

Journal of Signal Processing

Nov. 2017

IEICE Trans. Fundamentals

Jan. 2018

IEICE Trans. Fundamentals

Jan. 2018

IEICE Trans. Fundamentals

Dec. 20182017 RISP International
Workshop on Nonlinear Circuits,
Communications and Signal
Processing (NCSP'17)
2017 RISP International
Workshop on Nonlinear Circuits,
Communications and Signal
Processing (NCSP'17)
2017 International Workshop on
Smart Info-Media Systems in
Asia (SISA2017)2017 International Workshop on
Smart Info-Media Systems in
Asia (SISA2017)2017 International Workshop on
Smart Info-Media Systems in
Asia (SISA2017)2017 International Workshop on
Smart Info-Media Systems in
Asia (SISA2017)
8th International Workshop on
Signal Design and Its
Applications in Communications(IWSDA'17)

March 2017

March 2017

Sept. 2017

Sept. 2017

Sept. 2017
Sept. 2017

Sept. 2017

20th International Symposium
on Wireless Personal Multimedia Dec. 2017
Communications (WPMC2017)2019 RISP International
Workshop on Nonlinear Circuits,
Communications and Signal
Processing (NCSP2019)10March 2019