ホーム > 特許ランキング > 浜松ホトニクス株式会社
特開2024-159073データ作成方法、データ作成プログラム、及びデータ作成装置
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024159073
(43)【公開日】2024-11-08
(54)【発明の名称】データ作成方法、データ作成プログラム、及びデータ作成装置
(51)【国際特許分類】
G02F 1/01 20060101AFI20241031BHJP
【ＦＩ】
G02F1/01 B
【審査請求】未請求
【請求項の数】15
【出願形態】ＯＬ
(21)【出願番号】P 2023074832
(22)【出願日】2023-04-28
【新規性喪失の例外の表示】特許法第３０条第２項適用申請有り ウェブサイトの掲載日 令和５年１月１８日 ウェブサイトのアドレス（ＵＲＬ） ｈｔｔｐｓ：／／ｃｏｎｆｉｔ．ａｔｌａｓ．ｊｐ／ｇｕｉｄｅ／ｅｖｅｎｔ／ｌｓｊ４３／ｐｒｏｃｅｅｄｉｎｇｓ／ｌｉｓｔ 開催日 令和５年１月２０日 集会名、開催場所 一般社団法人レーザー学会学術講演会第４３回年次大会 （ウインクあいち 愛知県名古屋市中村区名駅４丁目４－３８）
(71)【出願人】
【識別番号】000236436
【氏名又は名称】浜松ホトニクス株式会社
(74)【代理人】
【識別番号】100088155
【弁理士】
【氏名又は名称】長谷川 芳樹
(74)【代理人】
【識別番号】100113435
【弁理士】
【氏名又は名称】黒木 義樹
(74)【代理人】
【識別番号】100140442
【弁理士】
【氏名又は名称】柴山 健一
(74)【代理人】
【識別番号】100174399
【弁理士】
【氏名又は名称】寺澤 正太郎
(72)【発明者】
【氏名】渡辺 向陽
(72)【発明者】
【氏名】井上 卓
【テーマコード（参考）】
2K102
【Ｆターム（参考）】
2K102AA21
2K102BA05
2K102BA21
2K102BB01
2K102BB04
2K102BC04
2K102BD09
2K102BD10
2K102DB08
2K102DD02
2K102EA21
2K102EB08
2K102EB10
2K102EB14
2K102EB20
(57)【要約】
【課題】所望の時間強度波形及び波長成分を有する光を精度良く実現するためのＳＬＭの変調パターンを得る。
【解決手段】このデータ作成方法は、空間光変調器を制御するデータを作成する方法である。このデータ作成方法は、第１ステップＳＴ１と、第２ステップＳＴ２と、第３ステップＳＴ３と、を含む。第１ステップＳＴ１では、複数の初期位相スペクトル関数を準備する。第２ステップＳＴ２では、複数の初期位相スペクトル関数それぞれを用いて、空間光変調器を制御し得る複数の予備データそれぞれを生成する。第３ステップＳＴ３では、複数の予備データの中から少なくとも一つの予備データを選択して、ＳＬＭを制御するデータとする。
【選択図】図１８
【解決手段】このデータ作成方法は、空間光変調器を制御するデータを作成する方法である。このデータ作成方法は、第１ステップＳＴ１と、第２ステップＳＴ２と、第３ステップＳＴ３と、を含む。第１ステップＳＴ１では、複数の初期位相スペクトル関数を準備する。第２ステップＳＴ２では、複数の初期位相スペクトル関数それぞれを用いて、空間光変調器を制御し得る複数の予備データそれぞれを生成する。第３ステップＳＴ３では、複数の予備データの中から少なくとも一つの予備データを選択して、ＳＬＭを制御するデータとする。
【選択図】図１８
【特許請求の範囲】
【請求項1】
空間光変調器を制御するデータを作成する方法であって、
複数の初期位相スペクトル関数を準備する第１ステップと、
前記複数の初期位相スペクトル関数それぞれを用いて、前記空間光変調器を制御し得る複数の予備データそれぞれを生成する第２ステップと、
前記複数の予備データの中から少なくとも一つの予備データを選択して前記空間光変調器を制御する前記データとする第３ステップと、
を含み、
前記第２ステップは、
強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する第１変換ステップと、
時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、前記第２波形関数から求める第２変換ステップと、
前記第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する第３変換ステップと、を含み、
前記第２ステップでは、前記第１波形関数を前記第４波形関数に置き換えながら前記複数の予備データのそれぞれに関して前記第１変換ステップ、前記第２変換ステップ、および前記第３変換ステップを繰り返し行い、繰り返し動作の最初における前記第１変換ステップにおいて前記複数の初期位相スペクトル関数それぞれを前記第１波形関数の前記位相スペクトル関数とし、繰り返し動作後に得られる前記第４波形関数の前記位相スペクトル関数に基づいて前記複数の予備データそれぞれを生成する、データ作成方法。
複数の初期位相スペクトル関数を準備する第１ステップと、
前記複数の初期位相スペクトル関数それぞれを用いて、前記空間光変調器を制御し得る複数の予備データそれぞれを生成する第２ステップと、
前記複数の予備データの中から少なくとも一つの予備データを選択して前記空間光変調器を制御する前記データとする第３ステップと、
を含み、
前記第２ステップは、
強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する第１変換ステップと、
時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、前記第２波形関数から求める第２変換ステップと、
前記第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する第３変換ステップと、を含み、
前記第２ステップでは、前記第１波形関数を前記第４波形関数に置き換えながら前記複数の予備データのそれぞれに関して前記第１変換ステップ、前記第２変換ステップ、および前記第３変換ステップを繰り返し行い、繰り返し動作の最初における前記第１変換ステップにおいて前記複数の初期位相スペクトル関数それぞれを前記第１波形関数の前記位相スペクトル関数とし、繰り返し動作後に得られる前記第４波形関数の前記位相スペクトル関数に基づいて前記複数の予備データそれぞれを生成する、データ作成方法。
【請求項2】
前記ターゲット強度スペクトログラムは、中心波長が互いに異なる複数の光パルスを含む光パルス列に関する強度スペクトログラムである、請求項１に記載のデータ作成方法。
【請求項3】
前記ターゲット強度スペクトログラムにおける前記複数の光パルス間の中心波長差を、目標とする前記複数の光パルス間の中心波長差よりも大きく設定する、請求項２に記載のデータ作成方法。
【請求項4】
前記ターゲット強度スペクトログラムにおける前記複数の光パルス間の中心波長差を、目標とする前記複数の光パルス間の中心波長差の１．１倍よりも大きく設定する、請求項２に記載のデータ作成方法。
【請求項5】
前記ターゲット強度スペクトログラムにおける前記複数の光パルス間の中心波長差を、前記空間光変調器への入力光の波長帯域を（前記光パルス列のパルス本数－１）で除算した値よりも小さく設定する、請求項２～４のいずれか１項に記載のデータ作成方法。
【請求項6】
前記第２変換ステップは、
前記第２波形関数を、強度スペクトログラム及び位相スペクトログラムに変換するステップと、
前記強度スペクトログラムを前記ターゲット強度スペクトログラムに置き換えるとともに、前記位相スペクトログラムを拘束するステップと、
置き換え後の前記強度スペクトログラム及び拘束された前記位相スペクトログラムを、前記第３波形関数に変換するステップと、
を含む、請求項１～４のいずれか１項に記載のデータ作成方法。
前記第２波形関数を、強度スペクトログラム及び位相スペクトログラムに変換するステップと、
前記強度スペクトログラムを前記ターゲット強度スペクトログラムに置き換えるとともに、前記位相スペクトログラムを拘束するステップと、
置き換え後の前記強度スペクトログラム及び拘束された前記位相スペクトログラムを、前記第３波形関数に変換するステップと、
を含む、請求項１～４のいずれか１項に記載のデータ作成方法。
【請求項7】
前記第２変換ステップは、
前記第２波形関数に対し、前記ターゲット強度スペクトログラムに対応する目標波形に基づく前記時間強度波形関数の置き換えを行うステップと、
前記ターゲット強度スペクトログラムに前記第２波形関数のスペクトログラムが近づくように前記第２波形関数を修正するステップと、
修正後の前記第２波形関数から前記第３波形関数を生成するステップと、
を含む、請求項１～４のいずれか１項に記載のデータ作成方法。
前記第２波形関数に対し、前記ターゲット強度スペクトログラムに対応する目標波形に基づく前記時間強度波形関数の置き換えを行うステップと、
前記ターゲット強度スペクトログラムに前記第２波形関数のスペクトログラムが近づくように前記第２波形関数を修正するステップと、
修正後の前記第２波形関数から前記第３波形関数を生成するステップと、
を含む、請求項１～４のいずれか１項に記載のデータ作成方法。
【請求項8】
空間光変調器を制御するデータを作成するためのプログラムであって、
複数の初期位相スペクトル関数を準備する第１ステップと、
前記複数の初期位相スペクトル関数それぞれを用いて、前記空間光変調器を制御し得る複数の予備データそれぞれを生成する第２ステップと、
前記複数の予備データの中から少なくとも一つの予備データを選択して前記空間光変調器を制御する前記データとする第３ステップと、
をコンピュータに実行させ、
前記第２ステップは、
強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する第１変換ステップと、
時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、前記第２波形関数から求める第２変換ステップと、
前記第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する第３変換ステップと、を含み、
前記第２ステップでは、前記第１波形関数を前記第４波形関数に置き換えながら前記複数の予備データのそれぞれに関して前記第１変換ステップ、前記第２変換ステップ、および前記第３変換ステップを繰り返し行い、繰り返し動作の最初における前記第１変換ステップにおいて前記複数の初期位相スペクトル関数それぞれを前記第１波形関数の前記位相スペクトル関数とし、繰り返し動作後に得られる前記第４波形関数の前記位相スペクトル関数に基づいて前記複数の予備データそれぞれを生成する、データ作成プログラム。
複数の初期位相スペクトル関数を準備する第１ステップと、
前記複数の初期位相スペクトル関数それぞれを用いて、前記空間光変調器を制御し得る複数の予備データそれぞれを生成する第２ステップと、
前記複数の予備データの中から少なくとも一つの予備データを選択して前記空間光変調器を制御する前記データとする第３ステップと、
をコンピュータに実行させ、
前記第２ステップは、
強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する第１変換ステップと、
時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、前記第２波形関数から求める第２変換ステップと、
前記第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する第３変換ステップと、を含み、
前記第２ステップでは、前記第１波形関数を前記第４波形関数に置き換えながら前記複数の予備データのそれぞれに関して前記第１変換ステップ、前記第２変換ステップ、および前記第３変換ステップを繰り返し行い、繰り返し動作の最初における前記第１変換ステップにおいて前記複数の初期位相スペクトル関数それぞれを前記第１波形関数の前記位相スペクトル関数とし、繰り返し動作後に得られる前記第４波形関数の前記位相スペクトル関数に基づいて前記複数の予備データそれぞれを生成する、データ作成プログラム。
【請求項9】
空間光変調器を制御するデータを作成する装置であって、
複数の初期位相スペクトル関数を記憶する記憶部と、
前記複数の初期位相スペクトル関数それぞれを用いて、前記空間光変調器を制御し得る複数の予備データそれぞれを生成する予備データ生成部と、
前記複数の予備データの中から少なくとも一つの予備データを選択して前記空間光変調器を制御する前記データとするデータ選択部と、
を備え、
前記予備データ生成部は、
強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する第１変換部と、
時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、前記第２波形関数から求める第２変換部と、
前記第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する第３変換部と、を有し、
前記予備データ生成部は、前記第１波形関数を前記第４波形関数に置き換えながら前記複数の予備データのそれぞれに関して前記第１変換部、前記第２変換部、および前記第３変換部の動作を繰り返し行い、繰り返し動作の最初における前記第１変換部において前記複数の初期位相スペクトル関数それぞれを前記第１波形関数の前記位相スペクトル関数とし、繰り返し動作後に得られる前記第４波形関数の前記位相スペクトル関数に基づいて前記複数の予備データそれぞれを生成する、データ作成装置。
複数の初期位相スペクトル関数を記憶する記憶部と、
前記複数の初期位相スペクトル関数それぞれを用いて、前記空間光変調器を制御し得る複数の予備データそれぞれを生成する予備データ生成部と、
前記複数の予備データの中から少なくとも一つの予備データを選択して前記空間光変調器を制御する前記データとするデータ選択部と、
を備え、
前記予備データ生成部は、
強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する第１変換部と、
時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、前記第２波形関数から求める第２変換部と、
前記第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する第３変換部と、を有し、
前記予備データ生成部は、前記第１波形関数を前記第４波形関数に置き換えながら前記複数の予備データのそれぞれに関して前記第１変換部、前記第２変換部、および前記第３変換部の動作を繰り返し行い、繰り返し動作の最初における前記第１変換部において前記複数の初期位相スペクトル関数それぞれを前記第１波形関数の前記位相スペクトル関数とし、繰り返し動作後に得られる前記第４波形関数の前記位相スペクトル関数に基づいて前記複数の予備データそれぞれを生成する、データ作成装置。
【請求項10】
前記ターゲット強度スペクトログラムは、中心波長が互いに異なる複数の光パルスを含む光パルス列に関する強度スペクトログラムである、請求項９に記載のデータ作成装置。
【請求項11】
前記ターゲット強度スペクトログラムにおける前記複数の光パルス間の中心波長差が、目標とする前記複数の光パルス間の中心波長差よりも大きく設定される、請求項１０に記載のデータ作成装置。
【請求項12】
前記ターゲット強度スペクトログラムにおける前記複数の光パルス間の中心波長差が、目標とする前記複数の光パルス間の中心波長差の１．１倍よりも大きく設定される、請求項１０に記載のデータ作成装置。
【請求項13】
前記ターゲット強度スペクトログラムにおける前記複数の光パルス間の中心波長差が、前記空間光変調器への入力光の波長帯域を（前記光パルス列のパルス本数－１）で除算した値よりも小さく設定される、請求項１０～１２のいずれか１項に記載のデータ作成装置。
【請求項14】
前記第２変換部は、
前記第２波形関数を、強度スペクトログラム及び位相スペクトログラムに変換する部分と、
前記強度スペクトログラムを前記ターゲット強度スペクトログラムに置き換えるとともに、前記位相スペクトログラムを拘束する部分と、
置き換え後の前記強度スペクトログラム及び拘束された前記位相スペクトログラムを、前記第３波形関数に変換する部分と、
を含む、請求項９～１２のいずれか１項に記載のデータ作成装置。
前記第２波形関数を、強度スペクトログラム及び位相スペクトログラムに変換する部分と、
前記強度スペクトログラムを前記ターゲット強度スペクトログラムに置き換えるとともに、前記位相スペクトログラムを拘束する部分と、
置き換え後の前記強度スペクトログラム及び拘束された前記位相スペクトログラムを、前記第３波形関数に変換する部分と、
を含む、請求項９～１２のいずれか１項に記載のデータ作成装置。
【請求項15】
前記第２変換部は、
前記第２波形関数に対し、前記ターゲット強度スペクトログラムに対応する目標波形に基づく前記時間強度波形関数の置き換えを行う部分と、
前記ターゲット強度スペクトログラムに前記第２波形関数のスペクトログラムが近づくように前記第２波形関数を修正する部分と、
修正後の前記第２波形関数から前記第３波形関数を生成する部分と、
を含む、請求項９～１２のいずれか１項に記載のデータ作成装置。
前記第２波形関数に対し、前記ターゲット強度スペクトログラムに対応する目標波形に基づく前記時間強度波形関数の置き換えを行う部分と、
前記ターゲット強度スペクトログラムに前記第２波形関数のスペクトログラムが近づくように前記第２波形関数を修正する部分と、
修正後の前記第２波形関数から前記第３波形関数を生成する部分と、
を含む、請求項９～１２のいずれか１項に記載のデータ作成装置。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、データ作成方法、データ作成プログラム、及びデータ作成装置に関する。
【背景技術】
【0002】
特許文献１及び非特許文献１には、空間光変調器（Spatial Light Modulator：ＳＬＭ）を用いてスペクトル位相及び／又はスペクトル強度を変調することにより、光パルスを成形する技術が開示されている。非特許文献１では、所望の光パルス波形を得るためのスペクトル位相及びスペクトル強度を、反復フーリエ法を用いて算出している。特許文献１では、時間強度波形を構成する光の波長成分（周波数成分）を制御している。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２０１８－０３６４８６号公報
【非特許文献】
【0004】
【非特許文献1】M. Hacker, G. Stobrawa, T. Feurer, “Iterative Fourier transformalgorithm for phase-only pulse shaping”, Optics Express, Vol. 9, No. 4, pp.191-199,13 August 2001
【非特許文献2】Olivier Ripoll, Ville Kettunen, Hans Peter Herzig, “Review ofiterative Fouriertransform algorithms for beam shaping applications”, OpticalEngineering, Vol. 43, No. 11, pp.2549-2556, November 2004
【発明の概要】
【発明が解決しようとする課題】
【0005】
例えば超短パルス光といった種々の光の時間波形を制御するための技術として、光パルスのスペクトル位相及びスペクトル強度をＳＬＭによって変調するものがある。このような技術では、光の時間強度波形を所望の波形に近づけるためのスペクトル位相及びスペクトル強度を算出し、そのスペクトル位相及びスペクトル強度を光に与えるための変調パターンをＳＬＭに呈示させる。また、時間強度波形の形状の制御に加えて、時間強度波形を構成する光の波長成分（周波数成分）を制御することもできる（例えば特許文献１を参照）。複数の光パルスを含む出力光を生成する場合、複数の光パルス毎に波長を異ならせることによって、分散計測装置、レーザ加工装置、超高速撮像カメラ、テラヘルツ波発生装置など様々な装置への応用が可能となる。このような技術において、所望の時間強度波形及び波長成分を有する光を精度良く実現するための変調パターンを得ることは重要である。
【0006】
本開示は、所望の時間強度波形及び波長成分を有する光を精度良く実現するためのＳＬＭの変調パターンを得ることができるデータ作成方法、データ作成プログラム、及びデータ作成装置を提供することを目的とする。
【課題を解決するための手段】
【0007】
［１］一形態のデータ作成方法は、ＳＬＭを制御するデータを作成する方法である。このデータ作成方法は、第１ステップと、第２ステップと、第３ステップと、を含む。第１ステップでは、複数の初期位相スペクトル関数を準備する。第２ステップでは、複数の初期位相スペクトル関数それぞれを用いて、空間光変調器を制御し得る複数の予備データそれぞれを生成する。第３ステップでは、複数の予備データの中から少なくとも一つの予備データを選択して、ＳＬＭを制御する上記データとする。第２ステップは、第１変換ステップと、第２変換ステップと、第３変換ステップと、を含む。第１変換ステップでは、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する。第２変換ステップでは、時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、第２波形関数から求める。第３変換ステップでは、第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する。第２ステップでは、第１波形関数を第４波形関数に置き換えながら、複数の予備データのそれぞれに関して第１変換ステップ、第２変換ステップ、および第３変換ステップを繰り返し行う。その際、繰り返し動作の最初における第１変換ステップにおいて、複数の初期位相スペクトル関数それぞれを、第１波形関数の位相スペクトル関数とする。また、繰り返し動作後に得られる第４波形関数の位相スペクトル関数に基づいて、複数の予備データそれぞれを生成する。
【0008】
［２］一形態のデータ作成プログラムは、ＳＬＭを制御するデータを作成するためのプログラムである。このデータ作成プログラムは、第１ステップと、第２ステップと、第３ステップと、をコンピュータに実行させる。第１ステップでは、複数の初期位相スペクトル関数を準備する。第２ステップでは、複数の初期位相スペクトル関数それぞれを用いて、空間光変調器を制御し得る複数の予備データそれぞれを生成する。第３ステップでは、複数の予備データの中から少なくとも一つの予備データを選択して、ＳＬＭを制御する上記データとする。第２ステップは、第１変換ステップと、第２変換ステップと、第３変換ステップと、を含む。第１変換ステップでは、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する。第２変換ステップでは、時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、第２波形関数から求める。第３変換ステップでは、第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する。第２ステップでは、第１波形関数を第４波形関数に置き換えながら、複数の予備データのそれぞれに関して第１変換ステップ、第２変換ステップ、および第３変換ステップを繰り返し行う。その際、繰り返し動作の最初における第１変換ステップにおいて、複数の初期位相スペクトル関数それぞれを、第１波形関数の位相スペクトル関数とする。また、繰り返し動作後に得られる第４波形関数の位相スペクトル関数に基づいて、複数の予備データそれぞれを生成する。
【0009】
［３］一形態のデータ作成装置は、ＳＬＭを制御するデータを作成する装置である。データ作成装置は、記憶部と、予備データ生成部と、データ選択部と、を備える。記憶部は、複数の初期位相スペクトル関数を記憶する。予備データ生成部は、複数の初期位相スペクトル関数それぞれを用いて、空間光変調器を制御し得る複数の予備データそれぞれを生成する。データ選択部は、複数の予備データの中から少なくとも一つの予備データを選択して空間光変調器を制御するデータとする。予備データ生成部は、第１変換部と、第２変換部と、第３変換部と、を有する。第１変換部は、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第１波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の第２波形関数に変換する。第２変換部は、時間強度波形関数及び時間位相波形関数を含み、予め作成されたターゲット強度スペクトログラムに対応する時間領域の第３波形関数を、第２波形関数から求める。第３変換部は、第３波形関数を、強度スペクトル関数及び位相スペクトル関数を含む周波数領域の第４波形関数に変換する。予備データ生成部は、第１波形関数を第４波形関数に置き換えながら複数の予備データのそれぞれに関して第１変換部、第２変換部、および第３変換部の動作を繰り返し行う。その際、予備データ生成部は、繰り返し動作の最初における第１変換部において、複数の初期位相スペクトル関数それぞれを、第１波形関数の位相スペクトル関数とする。予備データ生成部は、繰り返し動作後に得られる第４波形関数の位相スペクトル関数に基づいて、複数の予備データそれぞれを生成する。
【0010】
一般的に、ＳＬＭに呈示する変調パターンを算出するためには、まず一つの初期位相スペクトル関数を用意し、繰り返し計算によって、位相スペクトル関数を、目的とする時間強度波形及び波長成分を実現可能な値に徐々に近づける。しかしながら、最終的に得られるスペクトル位相及びスペクトル強度は、初期位相スペクトル関数によって異なる。その結果、その位相スペクトル関数に基づく変調パターンをＳＬＭに呈示させたとき、得られる光の時間強度波形及び波長成分も、初期位相スペクトル関数によって異なる。そこで、上記［１］～［３］のデータ作成方法、データ作成プログラム、及びデータ作成装置では、予め複数の初期値（初期位相スペクトル関数）を準備する。そして、第１変換ステップないし第３変換ステップ、または第１変換部ないし第３変換部において、初期位相スペクトル関数ごとに繰り返し計算を行い、得られた複数の予備データの中から、ＳＬＭを制御するデータとして少なくとも一つの予備データを選択する。これにより、所望の時間強度波形及び波長成分を有する光を精度良く実現するための変調パターンを含むデータを得ることができる。
【0011】
［４］上記［１］～［３］のデータ作成方法、データ作成プログラム、又はデータ作成装置において、ターゲット強度スペクトログラムは、中心波長が互いに異なる複数の光パルスを含む光パルス列に関する強度スペクトログラムであってもよい。その場合、中心波長が互いに異なる複数の光パルスを含む光パルス列の時間強度波形及び波長成分の精度を高めることができる。
【0012】
［５］上記［４］のデータ作成方法、データ作成プログラム、又はデータ作成装置において、ターゲット強度スペクトログラムにおける複数の光パルス間の中心波長差を、目標とする複数の光パルス間の中心波長差よりも大きく設定してもよい。或いは、［６］上記［４］のデータ作成方法、データ作成プログラム、又はデータ作成装置において、ターゲット強度スペクトログラムにおける複数の光パルス間の中心波長差を、目標とする複数の光パルス間の中心波長差の１．１倍よりも大きく設定してもよい。本発明者の知見によれば、繰り返し計算によって得られるデータにより得られる複数の光パルス間の中心波長差は、ターゲット強度スペクトログラムの基となった複数の光パルス間の中心波長差よりも小さくなる傾向がある。従って、ターゲット強度スペクトログラムにおける複数の光パルス間の中心波長差を上記のように設定することにより、複数の光パルス間の中心波長差をより精度良く実現することができる。
【0013】
［７］上記［４］～［６］のうちいずれかのデータ作成方法、データ作成プログラム、又はデータ作成装置において、ターゲット強度スペクトログラムにおける複数の光パルス間の中心波長差を、ＳＬＭへの入力光の波長帯域を（光パルス列のパルス本数－１）で除算した値よりも小さく設定してもよい。ターゲット強度スペクトログラムにおける複数の光パルス間の中心波長差をこのように設定することにより、複数の光パルス間の中心波長差をより精度良く実現することができる。
【0014】
［８］上記［１］、［２］、及び［４］～［７］のうちいずれかのデータ作成方法又はデータ作成プログラムにおいて、第２変換ステップは、第２波形関数を、強度スペクトログラム及び位相スペクトログラムに変換するステップと、強度スペクトログラムをターゲット強度スペクトログラムに置き換えるとともに、位相スペクトログラムを拘束するステップと、置き換え後の強度スペクトログラム及び拘束された位相スペクトログラムを、第３波形関数に変換するステップと、を含んでもよい。同様に、上記［３］～［７］のうちいずれかのデータ作成装置において、第２変換部は、第２波形関数を、強度スペクトログラム及び位相スペクトログラムに変換する部分と、強度スペクトログラムをターゲット強度スペクトログラムに置き換えるとともに、位相スペクトログラムを拘束する部分と、置き換え後の強度スペクトログラム及び拘束された位相スペクトログラムを、第３波形関数に変換する部分と、を含んでもよい。例えばこれらのデータ作成方法、データ作成プログラム、又はデータ作成装置によって、所望の時間強度波形及び波長成分を有する光を精度良く実現するための複数の予備データを得ることができる。
【0015】
［９］上記［１］、［２］、及び［４］～［７］のうちいずれかのデータ作成方法又はデータ作成プログラムにおいて、第２変換ステップは、第２波形関数に対し、ターゲット強度スペクトログラムに対応する目標波形に基づく時間強度波形関数の置き換えを行うステップと、ターゲット強度スペクトログラムに第２波形関数のスペクトログラムが近づくように第２波形関数を修正するステップと、修正後の第２波形関数から第３波形関数を生成するステップと、を含んでもよい。同様に、上記［３］～［７］のうちいずれかのデータ作成装置において、第２変換部は、第２波形関数に対し、ターゲット強度スペクトログラムに対応する目標波形に基づく時間強度波形関数の置き換えを行う部分と、ターゲット強度スペクトログラムに第２波形関数のスペクトログラムが近づくように第２波形関数を修正する部分と、修正後の第２波形関数から第３波形関数を生成する部分と、を含んでもよい。例えばこれらのデータ作成方法、データ作成プログラム、又はデータ作成装置によって、所望の時間強度波形及び波長成分を有する光を精度良く実現するための複数の予備データを得ることができる。
【発明の効果】
【0016】
本開示によるデータ作成方法、データ作成プログラム、及びデータ作成装置によれば、所望の時間強度波形及び波長成分を有する光を精度良く実現するためのＳＬＭの変調パターンを得ることができる。
【図面の簡単な説明】
【0017】
【図12】図１２（ａ）は、初期スペクトル関数の一例として、初期強度スペクトル関数及び初期位相スペクトル関数を模式的に示すグラフである。図１２（ｂ）は、第１波形関数からフーリエ変換された、１周目における第２波形関数の時間強度波形関数及び時間位相波形関数を模式的に示すグラフである。
【図13】図１３（ａ）は、図１２（ｂ）に示される第２波形関数から変換された強度スペクトログラムを示す図である。図１３（ｂ）は、図１２（ｂ）に示される第２波形関数から変換された位相スペクトログラムを示す図である。図１３（ｃ）は、ターゲット強度スペクトログラムの一例を示す図である。図１３（ｄ）は、拘束された位相スペクトログラムを示す図である。
【図14】図１４（ａ）は、図１３（ｃ）及び図１３（ｄ）に示される強度スペクトログラム及び位相スペクトログラムから逆ＳＴＦＴされた、１周目における第３波形関数の時間強度波形関数及び時間位相波形関数を模式的に示すグラフである。図１４（ｂ）は、図１４（ａ）に示される第３波形関数から逆フーリエ変換された、１周目における第４波形関数の強度スペクトル関数及び位相スペクトル関数を模式的に示すグラフである。
【図15】図１５（ａ）は、図１４（ｂ）ののちの３周目（ｎ＝３）における、第１波形関数の強度スペクトル関数及び位相スペクトル関数を模式的に示すグラフである。図１５（ｂ）は、図１５（ａ）に示される第１波形関数からフーリエ変換された、３周目における第２波形関数の時間強度波形関数及び時間位相波形関数を模式的に示すグラフである。
【図16】図１６（ａ）及び（ｂ）それぞれは、図１５（ｂ）に示される第２波形関数から変換された、３周目における強度スペクトログラム及び位相スペクトログラムそれぞれを示す図である。図１６（ｃ）は、３周目における強度スペクトログラムを示す図である。図１６（ｄ）は、３周目における拘束された位相スペクトログラムを示す図である。
【図17】図１７（ａ）は、図１６（ｃ）及び図１６（ｄ）に示される強度スペクトログラム及び位相スペクトログラムから逆ＳＴＦＴされた、３周目における第３波形関数の時間強度波形関数及び時間位相波形関数を模式的に示すグラフである。図１７（ｂ）は、図１７（ａ）に示される第３波形関数から逆フーリエ変換された、３周目における第４波形関数の強度スペクトル関数及び位相スペクトル関数を模式的に示すグラフである。
【図20】図２０（ａ）は、ターゲット波形関数の一例として、図６に示された一つの光パルスについてのターゲット波形関数の強度スペクトル関数及び位相スペクトル関数を示すグラフである。図２０（ｂ）は、時間領域の波形関数の一例として、図２０（ａ）に示されるターゲット波形関数から生成された時間強度波形関数及び時間位相波形関数を示すグラフである。
【図21】図２１（ａ）は、強度スペクトログラムの一例として、図２０（ｂ）に示される波形関数から生成された強度スペクトログラムを示す図である。図２１（ｂ）は、一例として、３つの光パルスについての強度スペクトログラムを互いに重畳して得られる強度スペクトログラムを示す図である。
【図23】図２３は、ターゲット強度スペクトログラムを用いて得られたデータを空間光変調器に呈示することにより実際に得られる光の、複数のパルス間の中心波長差と、複数のパルスのピーク強度のばらつきとの関係を、複数の初期位相スペクトル関数それぞれについてプロットしたグラフである。
【図24】図２４は、ターゲット強度スペクトログラムを用いて得られたデータを空間光変調器に呈示することにより実際に得られる光の、複数のパルス間の中心波長差と、複数のパルスのピーク強度のばらつきとの関係を、複数の初期位相スペクトル関数それぞれについてプロットしたグラフである。
【発明を実施するための形態】
【0018】
以下、添付図面を参照しながら実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。
【0019】
図１は、本開示の一実施形態に係る光制御装置１の構成を概略的に示す図である。本実施形態の光制御装置１は、入力光Ｌａから、該入力光Ｌａとは異なる任意の時間強度波形を有する出力光Ｌｂを生成する。図１に示されるように、光制御装置１は、光源２、光波形制御部３、及びデータ作成装置２０を備える。
【0020】
光源２は、光波形制御部３に入力される入力光Ｌａを出力する。入力光Ｌａは、例えば連続光である。或いは、光源２は例えば固体レーザ光源等のレーザ光源であり、入力光Ｌａ（初期の光パルス）は、例えばコヒーレントな単一の光パルスである。光源２は、例えばフェムト秒レーザであり、一実施例ではＬＤ直接励起型Ｙｂ：ＹＡＧパルスレーザである。光パルスの時間波形は例えばガウス関数状である。光パルスの半値全幅（ＦＷＨＭ）は、例えば１０ｆｓ～１００００ｆｓの範囲内であり、一例では１００ｆｓである。入力光Ｌａは、或る程度の帯域幅を有し、連続する複数の波長成分を含む。一実施例では、入力光Ｌａの帯域幅は１０ｎｍであり、入力光Ｌａの中心波長は８００ｎｍである。
【0021】
光波形制御部３は、空間光変調器（ＳＬＭ）１４を有しており、データ作成装置２０からの制御信号ＳＣをＳＬＭ１４に受ける。光波形制御部３は、光源２からの入力光Ｌａを、任意の時間強度波形を有する出力光Ｌｂに変換する。制御信号ＳＣは、ＳＬＭ１４を制御するためのデータ、すなわち複素振幅分布の強度あるいは位相分布の強度を含むデータに基づいて生成される。データは、例えば、計算機合成ホログラム（Computer-Generated Holograms(CGH)）である。
【0022】
図２は、光波形制御部３の構成例を示す図である。光波形制御部３は、回折格子１２、レンズ１３、ＳＬＭ１４、レンズ１５、及び回折格子１６を有する。回折格子１２は本実施形態における分光素子であり、光源２と光学的に結合されている。ＳＬＭ１４はレンズ１３を介して回折格子１２と光学的に結合されている。回折格子１２は、入力光Ｌａに含まれる複数の波長成分を、波長毎に空間的に分離する。なお、分光素子として、回折格子１２に代えてプリズム等の他の光学部品を用いてもよい。
【0023】
入力光Ｌａは、回折格子１２に対して斜めに入射し、複数の波長成分に分光される。この複数の波長成分を含む光Ｌ１は、レンズ１３によって波長成分毎に集光され、ＳＬＭ１４の変調面に結像される。レンズ１３は、光透過部材からなる凸レンズであってもよく、凹状の光反射面を有する凹面鏡であってもよい。
【0024】
ＳＬＭ１４は、入力光Ｌａを出力光Ｌｂに変換するために、回折格子１２から出力された複数の波長成分の位相を相互にずらす。そのために、ＳＬＭ１４は、データ作成装置２０から制御信号ＳＣを受けて、光Ｌ１の位相変調と強度変調とを同時に行う。なお、ＳＬＭ１４は、位相変調のみ、または強度変調のみを行ってもよい。ＳＬＭ１４は、例えば位相変調型である。一実施例では、ＳＬＭ１４はＬＣＯＳ（Liquid crystal on silicon）型である。なお、図には透過型のＳＬＭ１４が示されているが、ＳＬＭ１４は反射型であってもよい。また、空間光変調器１４は、位相変調型の空間光変調器に限らず、ＤＭＤ（Digital Micro Mirror Device）などの強度変調型の空間光変調器、或いは位相－強度変調型の空間光変調器であってもよい。
【0025】
図３は、ＳＬＭ１４の変調面１７を示す図である。図３に示すように、変調面１７には、複数の変調領域１７ａが或る方向ＡＡに沿って並んでおり、各変調領域１７ａは方向ＡＡと交差する方向ＡＢに延びている。方向ＡＡは、回折格子１２による分光方向である。この変調面１７はフーリエ変換面として働き、複数の変調領域１７ａのそれぞれには、分光後の対応する各波長成分が入射する。ＳＬＭ１４は、各変調領域１７ａにおいて、入射した各波長成分の位相及び強度を他の波長成分から独立して変調する。なお、本実施形態のＳＬＭ１４は位相変調型であるため、強度変調は、変調面１７に呈示される位相パターン（位相画像）によって実現される。
【0026】
ＳＬＭ１４によって変調された変調光Ｌ２の各波長成分は、レンズ１５によって回折格子１６上の一点に集められる。レンズ１５及び回折格子１６は、変調光Ｌ２を集光する光学系として機能する。レンズ１５は、光透過部材からなる凸レンズであってもよく、凹状の光反射面を有する凹面鏡であってもよい。また、回折格子１６は合波光学系として機能し、変調後の各波長成分を合波する。すなわち、これらのレンズ１５及び回折格子１６により、変調光Ｌ２の複数の波長成分は互いに集光・合波されて、出力光Ｌｂとなる。ＳＬＭ１４が反射型である場合、レンズ１３及び１５は共通のレンズによって構成されてもよく、回折格子１２及び１６は共通の回折格子によって構成されてもよい。
【0027】
レンズ１５よりも前の領域（スペクトル領域）と、回折格子１６よりも後ろの領域（時間領域）とは、互いにフーリエ変換の関係にあり、スペクトル領域における位相変調は、時間領域における時間強度波形に影響する。従って、出力光Ｌｂは、ＳＬＭ１４の変調パターンに応じた、入力光Ｌａとは異なる所望の時間強度波形を有することとなる。ここで、図４（ａ）は、一例として、単パルス状の入力光Ｌａのスペクトル波形（スペクトル位相Ｇ１１及びスペクトル強度Ｇ１２）を示し、図４（ｂ）は、該入力光Ｌａの時間強度波形を示す。また、図５（ａ）は、一例として、ＳＬＭ１４において矩形波状のスペクトル位相変調を与えたときの出力光Ｌｂのスペクトル波形（スペクトル位相Ｇ２１及びスペクトル強度Ｇ２２）を示し、図５（ｂ）は、該出力光Ｌｂの時間強度波形を示す。図４（ａ）及び図５（ａ）において、横軸は波長（ｎｍ）を示し、左の縦軸は強度スペクトルの強度値（任意単位）を示し、右の縦軸はスペクトル位相の位相値（ｒａｄ）を示す。また、図４（ｂ）及び図５（ｂ）において、横軸は時間（フェムト秒）を表し、縦軸は光強度（任意単位）を表す。この例では、矩形波状の位相スペクトル波形を出力光Ｌｂに与えることにより、入力光Ｌａのシングルパルスが、出力光Ｌｂとして高次光を伴うダブルパルスに変換されている。なお、図５に示されるスペクトル及び波形は一つの例であって、様々なスペクトル位相及びスペクトル強度の組み合わせにより、出力光Ｌｂの時間強度波形を様々な形状に整形することができる。
【0028】
出力光Ｌｂが複数の光パルスを含む光パルス列である場合、光パルス列は、入力光Ｌａを構成するスペクトルを複数の波長帯域に分け、それぞれの波長帯域を用いて生成したシングルパルス群であってもよい。この場合、複数の波長帯域の境界では、互いに重なり合う部分があってもよい。このような光パルス列は「帯域制御バーストパルス」と称される。
【0029】
図６は、帯域制御バーストパルスの例を示す図である。この例では、３つの光パルスＰｂ1～Ｐｂ3からなる光パルス列Ｐｂが示されている。図６（ａ）は、スペクトログラムであって、横軸に時間、縦軸に波長を示しており、光強度を色の濃淡で表している。図６（ｂ）は、光パルス列Ｐｂの時間波形を表している。各光パルスＰｂ1～Ｐｂ3の時間波形は例えばガウス関数状である。
【0030】
図６（ａ）及び図６（ｂ）に示すように、３つの光パルスＰｂ1～Ｐｂ3のピーク同士は時間的に互いに離れており、３つの光パルスＰｂ1～Ｐｂ3の伝搬タイミングは互いにずれている。言い換えると、一の光パルスＰｂ1に対して別の光パルスＰｂ2が時間遅れを有しており、該別の光パルスＰｂ2に対して更に別の光パルスＰｂ3が時間遅れを有している。但し、隣り合う光パルスＰｂ1，Ｐｂ2（又はＰｂ2，Ｐｂ3）の裾部分同士が互いに重なっていてもよい。隣り合う光パルスＰｂ1，Ｐｂ2（又はＰｂ2，Ｐｂ3）の時間間隔（ピーク間隔）は、例えば１０ｆｓ～１００００ｆｓの範囲内であり、一例では２０００ｆｓである。また、各光パルスＰｂ１～Ｐｂ３のＦＷＨＭは、例えば１０ｆｓ～５０００ｆｓの範囲内であり、一例では３００ｆｓである。
【0031】
図６（ｃ）は、３つの光パルスＰｂ1～Ｐｂ3を合成したスペクトルを表している。図６（ｃ）に示すように３つの光パルスＰｂ1～Ｐｂ3を合成したスペクトルは単一のピークを有するが、図６（ａ）を参照すると３つの光パルスＰｂ1～Ｐｂ3の中心波長は互いにずれている。図６（ｃ）に示す単一のピークを有するスペクトルは、ほぼ入力光Ｌａのスペクトルと同じである。
【0032】
隣り合う光パルスＰｂ1，Ｐｂ2（又はＰｂ2，Ｐｂ3）のピーク波長間隔は、入力光Ｌａのスペクトル帯域幅によって定まり、概ね半値全幅の２倍の範囲内である。一例では、入力光Ｌａのスペクトル帯域幅が１０ｎｍの場合、ピーク波長間隔は５ｎｍである。具体例として、入力光Ｌａの中心波長が８００ｎｍである場合、３つの光パルスＰｂ1～Ｐｂ3のピーク波長はそれぞれ７９５ｎｍ、８００ｎｍ、及び８０５ｎｍであることができる。
【0033】
図７は、帯域制御されていないバーストパルスの例を示す図である。この例では、３つの光パルスＰｄ1～Ｐｄ3からなる光パルス列Ｐｄが示されている。図７（ａ）は、図６（ａ）と同様に、スペクトログラムであって、横軸に時間、縦軸に波長を示しており、光強度を色の濃淡で表している。図７（ｂ）は、光パルス列Ｐｄの時間波形を表している。図７（ｃ）は、３つの光パルスＰｄ1～Ｐｄ3を合成したスペクトルを表している。
【0034】
図７（ａ）～（ｃ）に示すように、３つの光パルスＰｄ1～Ｐｄ3のピーク同士は時間的に互いに離れているが、３つの光パルスＰｄ1～Ｐｄ3の中心波長は互いに一致している。本実施形態の光波形制御部３は、このような光パルス列Ｐｄに限らず、図６に示されたような、各光パルスＰｂ1～Ｐｂ3の中心波長が互いに異なる光パルス列Ｐｂを生成することができる。
【0035】
再び図１を参照する。データ作成装置２０は、ＳＬＭ１４と通信可能に接続されており、出力光Ｌｂの時間強度波形を所望の波形に近づけるための位相変調パターンに関するデータを作成し、該データを含む制御信号ＳＣをＳＬＭ１４に提供する。本実施形態のデータ作成装置２０は、所望の波形を得る為のスペクトル位相及びスペクトル強度を出力光Ｌｂに与える位相変調用の位相パターンをＳＬＭ１４に呈示させる。そのために、データ作成装置２０は、記憶部２１と、スペクトログラム設定部２２と、予備データ生成部２３と、データ選択部２４とを有する。
【0036】
なお、データ作成装置２０は、制御信号ＳＣをＳＬＭ１４に提供する装置とは別に設けられてもよい。その場合、データ作成装置２０は、制御信号ＳＣをＳＬＭ１４に提供する装置に、位相変調パターンに関するデータを提供する。データ作成装置２０は、例えば、パーソナルコンピュータ；スマートフォン、タブレット端末などのスマートデバイス；あるいはクラウドサーバなどのプロセッサを有するコンピュータであってよい。
【0037】
図８は、データ作成装置２０のハードウェアの構成例を概略的に示す図である。図８に示すように、データ作成装置２０は、物理的には、プロセッサ（ＣＰＵ）６１、ＲＯＭ６２及びＲＡＭ６３等の主記憶装置、キーボード、マウス及びタッチスクリーン等の入力デバイス６４、ディスプレイ（タッチスクリーン含む）等の出力デバイス６５、他の装置との間でデータの送受信を行うためのネットワークカード等の通信モジュール６６、ハードディスク等の補助記憶装置６７などを含む、通常のコンピュータとして構成され得る。
【0038】
コンピュータのプロセッサ６１は、データ作成プログラムによって、上記のデータ作成装置２０の各機能を実現することができる。故に、データ作成プログラムは、コンピュータのプロセッサ６１を、データ作成装置２０におけるスペクトログラム設定部２２、予備データ生成部２３、及びデータ選択部２４として動作させる。データ作成プログラムは、コンピュータの内部または外部の記憶装置（記憶媒体）に記憶される。記憶装置は、非一時的記録媒体であってもよい。記録媒体としては、フレキシブルディスク、ＣＤ、ＤＶＤ等の記録媒体、ＲＯＭ等の記録媒体、半導体メモリ、クラウドサーバ等が例示される。ＲＯＭ６２または補助記憶装置６７は、記憶部２１を構成する。
【0039】
記憶部２１は、出力光Ｌｂの時間強度波形を所望の波形に近づけるための位相変調パターンに関するデータを作成する際に用いられる、複数の初期位相スペクトル関数を記憶する。初期位相スペクトル関数については、後に詳述する。
【0040】
スペクトログラム設定部２２は、ターゲット強度スペクトログラムを作成する。ターゲット強度スペクトログラムに関するデータは、予備データ生成部２３に与えられる。予備データ生成部２３は、記憶部２１に記憶された複数の初期位相スペクトル関数それぞれを用いて、与えられたターゲット強度スペクトログラムの実現に適した複数の位相スペクトル関数を算出する。そして、予備データ生成部２３は、複数の位相スペクトル関数それぞれに基づいて、光波形制御部３において出力光Ｌｂに与えるための複数の位相変調パターン（例えば、計算機合成ホログラム）それぞれを算出し、複数の位相変調パターンそれぞれに関する、ＳＬＭ１４を制御し得る複数の予備データそれぞれを生成する。
【0041】
データ選択部２４は、予備データ生成部２３により生成された複数の予備データの中から少なくとも一つの予備データを選択して、その予備データを、ＳＬＭ１４を制御するデータに決定する。そして、決定されたデータを含む制御信号ＳＣが、ＳＬＭ１４に提供される。ＳＬＭ１４は、制御信号ＳＣに基づいて制御される。データ選択部２４では、操作者の操作によってデータが選択されてもよいし、所定のルールに従って自動的にデータが選択されてもよい。
【0042】
図９は、予備データ生成部２３の内部構成を示すブロック図である。図９に示されるように、予備データ生成部２３は、フーリエ変換部２５（第１変換部）、関数変換部２６（第２変換部）、逆フーリエ変換部２７（第３変換部）、及びデータ生成部２８を有する。関数変換部２６は、スペクトログラム変換部２６１、スペクトログラム置換部２６２、スペクトログラム逆変換部２６３を含む。予備データ生成部２３は、以下に説明する算出方法を用いて、変調パターンの基になる位相スペクトル関数を算出する。図１０は、予備データ生成部２３における位相スペクトル関数の計算手順を示すブロック図である。図１１は、予備データ生成部２３における位相スペクトル関数の計算手順を数式により示す図である。
【0043】
まず、初期スペクトル関数Ａ１、すなわち周波数ωの関数である初期強度スペクトル関数Ａ0（ω）及び初期位相スペクトル関数Φ0（ω）を用意する。一例では、初期強度スペクトル関数Ａ0（ω）は入力光Ｌａの強度スペクトルを表すが、これに限られない。また、初期位相スペクトル関数Φ0（ω）は、記憶部２１に記憶された複数の初期位相スペクトル関数Φ0（ω）の中から順に選択される。複数の初期位相スペクトル関数Φ0（ω）は、例えば乱数を用いて自動的に生成され得る。初期位相スペクトル関数Φ0（ω）の個数は、例えば１００個以上である。データ作成装置２０は、そのような複数の初期位相スペクトル関数Φ0（ω）を生成する部分、すなわち初期位相スペクトル関数生成部を更に備えてもよい。また、予め多数の初期位相スペクトル関数Φ0（ω）を記憶部２１に記憶させておき、操作者がその中から一部の初期位相スペクトル関数Φ0（ω）を選択してもよい。その場合、記憶部２１に記憶される初期位相スペクトル関数Φ0（ω）は、入力光Ｌａのパラメータ（中心波長、帯域幅、スペクトル形状、スペクトル位相など）、目標とする光パルス列Ｐｂのパラメータ（パルス本数、パルスの時間間隔、各パルスの中心波長、各パルスの帯域幅、各パルスのピーク強度比など）といった種々の条件ごとにグループ分けして記録されてもよい。その場合、操作者が望む光パルス列Ｐｂに適した初期位相スペクトル関数Φ0（ω）のグループを選択し得る。
【0044】
図１２（ａ）は、初期スペクトル関数Ａ１の一例として、初期強度スペクトル関数Ａ0（ω）及び初期位相スペクトル関数Φ0（ω）を模式的に示すグラフである。図１２（ａ）において、グラフＧ３１は初期強度スペクトル関数Ａ0（ω）を表し、グラフＧ３２は初期位相スペクトル関数Φ0（ω）を表す。横軸は波長を示し、縦軸は強度スペクトル関数の強度値または位相スペクトル関数の位相値を示す。
【0045】
第１波形関数Ａ２は、下記の数式（１）によって示される。但し、ｎは繰り返し番号（ｎ＝１，２，・・・，Ｎ）であり、ｉは虚数である。第１周目（ｎ＝１）においては、初期スペクトル関数Ａ１を、第１波形関数Ａ２とする。すなわち、繰り返し第１周目の第１波形関数Ａ２のΦ1（ω）＝Φ0（ω）とする。
予備データ生成部２３のフーリエ変換部２５は、第１波形関数Ａ２に対して周波数領域から時間領域へのフーリエ変換を行う（図中の矢印Ｂ１、第１変換ステップ）。これにより、下記の数式（２）によって示される、時間強度波形関数ａn（ｔ）及び時間位相波形関数φn（ｔ）を含む、時間領域の第２波形関数Ａ３が得られる。
図１２（ｂ）は、図１２（ａ）に示される第１波形関数Ａ２からフーリエ変換された、１周目における第２波形関数Ａ３の時間強度波形関数ａ1（ｔ）及び時間位相波形関数φ1（ｔ）を模式的に示すグラフである。図１２（ｂ）において、グラフＧ４１は時間強度波形関数ａ1（ｔ）を表し、グラフＧ４２は時間位相波形関数φ1（ｔ）を表す。横軸は時間を示し、縦軸は時間強度波形関数の強度値または時間位相波形関数の位相値を示す。
【数1】
予備データ生成部２３のフーリエ変換部２５は、第１波形関数Ａ２に対して周波数領域から時間領域へのフーリエ変換を行う（図中の矢印Ｂ１、第１変換ステップ）。これにより、下記の数式（２）によって示される、時間強度波形関数ａn（ｔ）及び時間位相波形関数φn（ｔ）を含む、時間領域の第２波形関数Ａ３が得られる。
【数2】
図１２（ｂ）は、図１２（ａ）に示される第１波形関数Ａ２からフーリエ変換された、１周目における第２波形関数Ａ３の時間強度波形関数ａ1（ｔ）及び時間位相波形関数φ1（ｔ）を模式的に示すグラフである。図１２（ｂ）において、グラフＧ４１は時間強度波形関数ａ1（ｔ）を表し、グラフＧ４２は時間位相波形関数φ1（ｔ）を表す。横軸は時間を示し、縦軸は時間強度波形関数の強度値または時間位相波形関数の位相値を示す。
【0046】
続いて、関数変換部２６は、予め作成された強度スペクトログラムＡ４３（ターゲット強度スペクトログラム）に対応する、時間強度波形関数ａ'n（ｔ）及び時間位相波形関数φ'n（ｔ）を含む時間領域の第３波形関数Ａ５を、第２波形関数Ａ３から求める（第２変換ステップ）。以下、具体的な関数変換部２６の動作を説明する。
【0047】
関数変換部２６のスペクトログラム変換部２６１は、第２波形関数Ａ３に対して短時間フーリエ変換（Short-Time Fourier Transform；ＳＴＦＴ）を行う（図中の矢印Ｂ２）。これにより、下記の数式（３）によって示される強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２が得られる。
図１３（ａ）及び図１３（ｂ）それぞれは、図１２（ｂ）に示される第２波形関数Ａ３から変換された強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２それぞれを示す図である。なお、図１３（ａ）及び図１３（ｂ）において横軸は時間を示し、縦軸は波長を示す。また、スペクトログラムの値は、図の明暗によって示されており、明るいほどスペクトログラムの値が大きい。この例では、強度スペクトログラムＡ４１において単一の光パルスＰｂ0が現れている。
【数3】
図１３（ａ）及び図１３（ｂ）それぞれは、図１２（ｂ）に示される第２波形関数Ａ３から変換された強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２それぞれを示す図である。なお、図１３（ａ）及び図１３（ｂ）において横軸は時間を示し、縦軸は波長を示す。また、スペクトログラムの値は、図の明暗によって示されており、明るいほどスペクトログラムの値が大きい。この例では、強度スペクトログラムＡ４１において単一の光パルスＰｂ0が現れている。
【0048】
なお、第２波形関数Ａ３を強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２に変換するためのスペクトログラム変換部２６１における処理は、ＳＴＦＴに限られず、他の処理であってもよい。時間波形をスペクトログラムに変換する処理は、ＳＴＦＴを含めて時間－周波数変換と呼ばれる。時間－周波数変換では、時間波形のような複合信号に対して、周波数フィルタ処理または数値演算処理（窓関数をずらしながら乗算して、各々の時間に対してスペクトルを導出する処理）を施し、時間、周波数、信号成分の強さ（強度スペクトル）からなる３次元情報を生成する。本実施形態では、その変換結果（時間、周波数、強度スペクトル）を「スペクトログラム」と定義する。時間－周波数変換としては、ＳＴＦＴのほか、ウェーブレット変換（ハールウェーブレット変換、ガボールウェーブレット変換、メキシカンハットウェーブレット変換、モルレーウェーブレット変換）などがある。
【0049】
続いて、関数変換部２６のスペクトログラム置換部２６２は、強度スペクトログラムＡ４１を予め作成された強度スペクトログラムＡ４３（ターゲット強度スペクトログラム）に置き換えるとともに、位相スペクトログラムＡ４２を拘束する（図中の矢印Ｂ３）。強度スペクトログラムＡ４３は、スペクトログラム設定部２２（図１を参照）から提供される。位相スペクトログラムＡ４２を拘束するとは、位相スペクトログラムＡ４２を変化させない（そのまま残す）ことを意味する。これにより、上記の数式（３）は、下記の数式（４）に置き換えられる。ＴＳＧ（ω，ｔ）はターゲット強度スペクトログラム関数である。
図１３（ｃ）は、強度スペクトログラムＡ４３の一例を示す図である。図１３（ｄ）は、拘束された位相スペクトログラムＡ４２を示す図である。図１３（ｃ）及び図１３（ｄ）において、横軸は時間を示し、縦軸は波長を示す。スペクトログラムの値は、図の明暗によって示されており、明るいほどスペクトログラムの値が大きい。この例では、強度スペクトログラムＡ４３に、互いに時間差を有し中心波長が互いに異なる３つの光パルスＰｂ1，Ｐｂ2，Ｐｂ3を含む光パルス列Ｐｂが含まれている。拘束された位相スペクトログラムＡ４２は、図１３（ｂ）と全く同一である。
【数4】
図１３（ｃ）は、強度スペクトログラムＡ４３の一例を示す図である。図１３（ｄ）は、拘束された位相スペクトログラムＡ４２を示す図である。図１３（ｃ）及び図１３（ｄ）において、横軸は時間を示し、縦軸は波長を示す。スペクトログラムの値は、図の明暗によって示されており、明るいほどスペクトログラムの値が大きい。この例では、強度スペクトログラムＡ４３に、互いに時間差を有し中心波長が互いに異なる３つの光パルスＰｂ1，Ｐｂ2，Ｐｂ3を含む光パルス列Ｐｂが含まれている。拘束された位相スペクトログラムＡ４２は、図１３（ｂ）と全く同一である。
【0050】
続いて、関数変換部２６のスペクトログラム逆変換部２６３は、強度スペクトログラムＡ４３及び位相スペクトログラムＡ４２に対して逆ＳＴＦＴを行う（図中の矢印Ｂ４）。なお、ここでもスペクトログラム変換部２６１と同様に、ＳＴＦＴ以外の時間－周波数変換を用いてもよい。これにより、下記の数式（５）によって示される、時間強度波形関数ａ'n（ｔ）及び時間位相波形関数φ'n（ｔ）を含む時間領域の第３波形関数Ａ５が得られる。但し、１周目ではｎ＝１である。
図１４（ａ）は、図１３（ｃ）及び図１３（ｄ）に示される強度スペクトログラムＡ４３及び位相スペクトログラムＡ４２から逆ＳＴＦＴされた、１周目における第３波形関数Ａ５の時間強度波形関数ａ'1（ｔ）及び時間位相波形関数φ'1（ｔ）を模式的に示すグラフである。図１４（ａ）において、グラフＧ５１は時間強度波形関数ａ'1（ｔ）を表し、グラフＧ５２は時間位相波形関数φ'1（ｔ）を表す。横軸は時間を示し、縦軸は時間強度波形関数の強度値または時間位相波形関数の位相値を示す。
【数5】
図１４（ａ）は、図１３（ｃ）及び図１３（ｄ）に示される強度スペクトログラムＡ４３及び位相スペクトログラムＡ４２から逆ＳＴＦＴされた、１周目における第３波形関数Ａ５の時間強度波形関数ａ'1（ｔ）及び時間位相波形関数φ'1（ｔ）を模式的に示すグラフである。図１４（ａ）において、グラフＧ５１は時間強度波形関数ａ'1（ｔ）を表し、グラフＧ５２は時間位相波形関数φ'1（ｔ）を表す。横軸は時間を示し、縦軸は時間強度波形関数の強度値または時間位相波形関数の位相値を示す。
【0051】
続いて、逆フーリエ変換部２７は、第３波形関数Ａ５に対して時間領域から周波数領域への逆フーリエ変換を行う（図中の矢印Ｂ５、第３変換ステップ）。これにより、下記の数式（６）によって示される、強度スペクトル関数Ａ'n（ω）及び位相スペクトル関数Φ'n（ω）を含む周波数領域の第４波形関数Ａ６が得られる。但し、１周目ではｎ＝１である。
図１４（ｂ）は、図１４（ａ）に示される第３波形関数Ａ５から逆フーリエ変換された、１周目における第４波形関数Ａ６の強度スペクトル関数Ａ'1（ω）及び位相スペクトル関数Φ'1（ω）を模式的に示すグラフである。図１４（ｂ）において、グラフＧ６１は強度スペクトル関数Ａ'1（ω）を表し、グラフＧ６２は位相スペクトル関数Φ'1（ω）を表す。横軸は波長を示し、縦軸は強度スペクトル関数の強度値または位相スペクトル関数の位相値を示す。
【数6】
図１４（ｂ）は、図１４（ａ）に示される第３波形関数Ａ５から逆フーリエ変換された、１周目における第４波形関数Ａ６の強度スペクトル関数Ａ'1（ω）及び位相スペクトル関数Φ'1（ω）を模式的に示すグラフである。図１４（ｂ）において、グラフＧ６１は強度スペクトル関数Ａ'1（ω）を表し、グラフＧ６２は位相スペクトル関数Φ'1（ω）を表す。横軸は波長を示し、縦軸は強度スペクトル関数の強度値または位相スペクトル関数の位相値を示す。
【0052】
その後、予備データ生成部２３は、第１波形関数Ａ２の強度スペクトル関数Ａ0（ω）を拘束しつつ、第１波形関数Ａ２の位相スペクトル関数Φ1（ω）を第４波形関数の位相スペクトル関数Φ'1（ω）に置き換える（すなわち第１波形関数Ａ２のΦ2（ω）＝Φ'1（ω）とする。図中の矢印Ｂ６を参照）。そして、上述したフーリエ変換部２５、スペクトログラム変換部２６１、スペクトログラム置換部２６２、スペクトログラム逆変換部２６３、及び逆フーリエ変換部２７の動作を、強度スペクトログラムＡ４１とターゲット強度スペクトログラムＡ４３との一致度合いを示す評価値が収束するまでＮ回にわたって繰り返す。このように、予備データ生成部２３では、第１波形関数Ａ２の強度スペクトル関数Ａ0（ω）を拘束しつつ、（ｎ＋１）周目の第１波形関数Ａ２の位相スペクトル関数Φn+1（ω）をｎ周目の第４波形関数の位相スペクトル関数Φ'n（ω）に置き換えながら、フーリエ変換部２５、スペクトログラム変換部２６１、スペクトログラム置換部２６２、スペクトログラム逆変換部２６３、及び逆フーリエ変換部２７がこの順で繰り返し動作する。
【0053】
図１５（ａ）は、図１４（ｂ）ののちの３周目（ｎ＝３）における、第１波形関数Ａ２の強度スペクトル関数Ａ0（ω）（グラフＧ７１）及び位相スペクトル関数Φ3（ω）（グラフＧ７２）を模式的に示すグラフである。位相スペクトル関数Φ3（ω）の波形が１周目から変化していることがわかる。図１５（ｂ）は、図１５（ａ）に示される第１波形関数Ａ２からフーリエ変換された、３周目における第２波形関数Ａ３の時間強度波形関数ａ3（ｔ）（グラフＧ８１）及び時間位相波形関数φ3（ｔ）（グラフＧ８２）を模式的に示すグラフである。図１６（ａ）及び（ｂ）それぞれは、図１５（ｂ）に示される第２波形関数Ａ３から変換された、３周目における強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２それぞれを示す図である。互いに時間差を有し中心波長が互いに異なる３つの光パルスＰｂ1～Ｐｂ3が生成され始めていることがわかる。図１６（ｃ）は、３周目における強度スペクトログラムＡ４３を示す図であって、これは図１３（ｃ）に示される１週目の強度スペクトログラムＡ４３と同一である。図１６（ｄ）は、３周目における拘束された位相スペクトログラムＡ４２を示す図であって、これは図１６（ｂ）の位相スペクトログラムＡ４２と同一である。図１７（ａ）は、図１６（ｃ）及び図１６（ｄ）に示される強度スペクトログラムＡ４３及び位相スペクトログラムＡ４２から逆ＳＴＦＴされた、３周目における第３波形関数Ａ５の時間強度波形関数ａ'3（ｔ）（グラフＧ９１）及び時間位相波形関数φ'3（ｔ）（グラフＧ９２）を模式的に示すグラフである。図１７（ｂ）は、図１７（ａ）に示される第３波形関数Ａ５から逆フーリエ変換された、３周目における第４波形関数Ａ６の強度スペクトル関数Ａ'3（ω）（グラフＧ１０１）及び位相スペクトル関数Φ'3（ω）（グラフＧ１０２）を模式的に示すグラフである。
【0054】
上記のような繰り返し動作によって、強度スペクトログラムＡ４１が強度スペクトログラムＡ４３に次第に近づくように、位相スペクトル関数Φ'n（ω）が修正される。最終的に、第４波形関数Ａ６に含まれる位相スペクトル関数Φ'n（ω）が、所望のスペクトル位相解Φresult（ω）となる。このスペクトル位相解Φresult（ω）が、データ生成部２８に提供される。
【0055】
データ生成部２８は、予備データ生成部２３において算出されたスペクトル位相解Φresult（ω）に基づくスペクトル位相及び／又はスペクトル強度を出力光Ｌｂに与えるための位相変調パターン（例えば、計算機合成ホログラム）に関する予備データを算出する。
【0056】
予備データ生成部２３は、上記のような繰り返し計算および予備データの算出を、記憶部２１に記憶されている複数の初期位相スペクトル関数Φ0（ω）のそれぞれについて行う。予備データ生成部２３は、繰り返し計算および予備データの算出を、一つの初期位相スペクトル関数Φ0（ω）ずつ順次行ってもよいし、複数の初期位相スペクトル関数Φ0（ω）について並行して行ってもよい。
【0057】
算出された複数の予備データは、それらによって生成される光パルス列Ｐｂのパラメータ（中心波長差、ピーク強度のばらつきなど）ごとにグループ分けして記録されてもよい。その場合、データ選択部２４において、操作者が望む光パルス列Ｐｂに適したデータを選択し易くなる。
【0058】
図１８は、本実施形態に係るデータ作成方法を示すフローチャートである。このデータ作成方法は、ＳＬＭ１４を制御するデータを作成する方法であって、上述したデータ作成装置２０を用いて好適に実現される。また、本実施形態に係るデータ作成プログラムは、このデータ作成方法の各ステップをコンピュータに実行させるためのプログラムである。
【0059】
図１８に示されるように、まず、複数の初期位相スペクトル関数Φ0（ω）を準備する（第１ステップＳＴ１）。前述したように、複数の初期位相スペクトル関数Φ0（ω）は、乱数を用いて生成されてもよい。準備された複数の初期位相スペクトル関数Φ0（ω）は、記憶部２１に記憶される。
【0060】
次に、予備データ生成部２３が、複数の初期位相スペクトル関数Φ0（ω）のそれぞれを用いて、ＳＬＭ１４を制御し得る複数の予備データそれぞれを、予備データ生成部２３により生成する（第２ステップＳＴ２）。この第２ステップＳＴ２は、第１変換ステップＳＴ２１と、第２変換ステップＳＴ２２と、第３変換ステップＳＴ２３と、データ生成ステップＳＴ２４とを含む。第１変換ステップＳＴ２１では、フーリエ変換部２５が、強度スペクトル関数Ａ0（ω）及び位相スペクトル関数Φn（ω）を含む周波数領域の第１波形関数Ａ２を、時間強度波形関数ａn（ｔ）及び時間位相波形関数φn（ｔ）を含む時間領域の第２波形関数Ａ３に変換する。第２変換ステップＳＴ２２では、関数変換部２６が、時間強度波形関数ａ'n（ｔ）及び時間位相波形関数φ'n（ｔ）を含み、予め作成されたターゲット強度スペクトログラムＡ４３に対応する時間領域の第３波形関数Ａ５を、第２波形関数Ａ３から求める。第２変換ステップＳＴ２２の詳細は、前述した関数変換部２６の動作のとおりである。第３変換ステップＳＴ２３では、逆フーリエ変換部２７が、第３波形関数Ａ５を、強度スペクトル関数Ａ'n（ω）及び位相スペクトル関数Φ'n（ω）を含む周波数領域の第４波形関数Ａ６に変換する。第２ステップＳＴ２では、第１波形関数Ａ２を第４波形関数Ａ６に置き換えながら、第１変換ステップＳＴ２１、第２変換ステップＳＴ２２、および第３変換ステップＳＴ２３を繰り返し行う。また、この繰り返し動作の最初における第１変換ステップＳＴ２１においては、複数の初期位相スペクトル関数Ａ0（ω）のうち一つを、第１波形関数Ａ２の位相スペクトル関数とする。そして、データ生成ステップＳＴ２４では、データ生成部２８が、繰り返し動作後に得られる第４波形関数Ａ６の位相スペクトル関数Φ'n（ω）に基づいて、予備データを生成する。第２ステップＳＴ２では、複数の初期位相スペクトル関数Ａ0（ω）の全てについて、上記の繰り返し計算を行い、データ生成ステップＳＴ２４において予備データを生成する（ステップＳＴ２５）。こうして、複数の初期位相スペクトル関数Ａ0（ω）にそれぞれ対応する複数の予備データが得られる。
【0061】
その後、データ選択部２４が、複数の予備データの中から少なくとも一つの予備データを選択して、ＳＬＭ１４を制御するデータとする（第３ステップＳＴ３）。
【0062】
以上に説明した本実施形態のデータ作成装置２０、データ作成方法、及びデータ作成プログラムによって得られる効果について説明する。一般的に、ＳＬＭ１４に呈示する変調パターンを算出するためには、まず一つの初期位相スペクトル関数Φ0（ω）を用意し、繰り返し計算によって、位相スペクトル関数Φn（ω）を、目的とする時間強度波形及び波長成分を実現可能な値に徐々に近づける。しかしながら、最終的に得られる位相スペクトル関数Φresult（ω）は、初期位相スペクトル関数Φ0（ω）によって異なる。その結果、その位相スペクトル関数Φresult（ω）に基づく変調パターンをＳＬＭ１４に呈示させたとき、得られる光の時間強度波形及び波長成分も、初期位相スペクトル関数Φ0（ω）によって異なる。そこで、本実施形態では、予め複数の初期位相スペクトル関数Φ0（ω）を準備する。そして、第１変換ステップＳＴ２１ないし第３変換ステップＳＴ２３において、初期位相スペクトル関数Φ0（ω）ごとに繰り返し計算を行い、得られた複数の予備データの中から、ＳＬＭ１４を制御するデータとして少なくとも一つの予備データを選択する。これにより、所望の時間強度波形及び波長成分を有する光を精度良く実現するための変調パターンを含むデータを得ることができる。
【0063】
本実施形態のように、ターゲット強度スペクトログラムＡ４３は、中心波長が互いに異なる複数の光パルスＰｂ1～Ｐｂ3を含む光パルス列Ｐｂに関する強度スペクトログラムであってもよい。その場合、中心波長が互いに異なる複数の光パルスＰｂ1～Ｐｂ3を含む光パルス列Ｐｂの時間強度波形及び波長成分の精度を高めることができる。
【0064】
本実施形態のように、第２変換ステップＳＴ２２は、第２波形関数Ａ３を、強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２に変換するステップと、強度スペクトログラムＡ４１をターゲット強度スペクトログラムＡ４３に置き換えるとともに、位相スペクトログラムＡ４２を拘束するステップと、置き換え後の強度スペクトログラムＡ４３及び拘束された位相スペクトログラムＡ４２を、第３波形関数Ａ５に変換するステップと、を含んでもよい。同様に、関数変換部２６は、第２波形関数Ａ３を、強度スペクトログラムＡ４１及び位相スペクトログラムＡ４２に変換する部分、すなわちスペクトログラム変換部２６１と、強度スペクトログラムＡ４１をターゲット強度スペクトログラムＡ４３に置き換えるとともに、位相スペクトログラムＡ４２を拘束する部分、すなわちスペクトログラム置換部２６２と、置き換え後の強度スペクトログラムＡ４３及び拘束された位相スペクトログラムＡ４２を、第３波形関数Ａ５に変換する部分、すなわちスペクトログラム逆変換部２６３と、を含んでもよい。例えばこのような構成により、所望の時間強度波形及び波長成分を有する光を精度良く実現するための複数の予備データを得ることができる。
【0065】
［ターゲット強度スペクトログラムの作成］
続いて、スペクトログラム設定部２２におけるターゲット強度スペクトログラムＡ４３の作成手順について説明する。図１９は、スペクトログラム設定部２２の機能的構成を示すブロック図である。スペクトログラム設定部２２は、互いに時間差を有し中心波長が互いに異なる複数の光パルス（例えば上述した光パルスＰｂ1～Ｐｂ3）を含む光パルス列Ｐｂに関するターゲット強度スペクトログラムＡ４３を作成する。図１９に示されるように、スペクトログラム設定部２２は、波形関数設定部１０１と、フーリエ変換部１０２と、スペクトログラム変換部１０３と、ターゲット生成部１０４と、を備える。
続いて、スペクトログラム設定部２２におけるターゲット強度スペクトログラムＡ４３の作成手順について説明する。図１９は、スペクトログラム設定部２２の機能的構成を示すブロック図である。スペクトログラム設定部２２は、互いに時間差を有し中心波長が互いに異なる複数の光パルス（例えば上述した光パルスＰｂ1～Ｐｂ3）を含む光パルス列Ｐｂに関するターゲット強度スペクトログラムＡ４３を作成する。図１９に示されるように、スペクトログラム設定部２２は、波形関数設定部１０１と、フーリエ変換部１０２と、スペクトログラム変換部１０３と、ターゲット生成部１０４と、を備える。
【0066】
波形関数設定部１０１は、強度スペクトル関数及び位相スペクトル関数を含む周波数領域のターゲット波形関数を、光パルス毎に設定する。例えば、波形関数設定部１０１は、光パルスＰｂ1についてターゲット波形関数を設定し、光パルスＰｂ2について別のターゲット波形関数を設定し、光パルスＰｂ3について更に別のターゲット波形関数を設定する。ターゲット波形関数のパラメータとしては、例えば下記のものが挙げられる。
ａ）各光パルスの強度スペクトル関数の形状（例えばガウシアン）
ｂ）各光パルスの強度スペクトル関数のスペクトルエネルギー量
ｃ）各光パルスの強度スペクトル関数の帯域幅（半値全幅）λｓ
ｄ）各光パルスの強度スペクトル関数の中心波長
ｅ）各光パルスの位相スペクトル関数
なお、各光パルスのスペクトル位相が周波数に対して線形である場合、その線形関数の傾きは、各光パルスの時間シフト量に相当する。図２０（ａ）は、ターゲット波形関数の一例として、図６に示された光パルスＰｂ2についてのターゲット波形関数の強度スペクトル関数Ｇ１３１及び位相スペクトル関数Ｇ１３２を示すグラフである。図２０（ａ）において、横軸は波長（ｎｍ）を示し、縦軸は強度スペクトル関数Ｇ１３１の強度値（任意単位）及び位相スペクトル関数Ｇ１３２の位相値（ｒａｄ）を示す。この例は、時間的に等間隔である３つの光パルスＰｂ1～Ｐｂ3のうち中央に位置する光パルスＰｂ2についてのターゲット波形関数であるので、位相スペクトル関数Ｇ１３２の傾き（すなわち時間シフト量）はゼロとされている。
ａ）各光パルスの強度スペクトル関数の形状（例えばガウシアン）
ｂ）各光パルスの強度スペクトル関数のスペクトルエネルギー量
ｃ）各光パルスの強度スペクトル関数の帯域幅（半値全幅）λｓ
ｄ）各光パルスの強度スペクトル関数の中心波長
ｅ）各光パルスの位相スペクトル関数
なお、各光パルスのスペクトル位相が周波数に対して線形である場合、その線形関数の傾きは、各光パルスの時間シフト量に相当する。図２０（ａ）は、ターゲット波形関数の一例として、図６に示された光パルスＰｂ2についてのターゲット波形関数の強度スペクトル関数Ｇ１３１及び位相スペクトル関数Ｇ１３２を示すグラフである。図２０（ａ）において、横軸は波長（ｎｍ）を示し、縦軸は強度スペクトル関数Ｇ１３１の強度値（任意単位）及び位相スペクトル関数Ｇ１３２の位相値（ｒａｄ）を示す。この例は、時間的に等間隔である３つの光パルスＰｂ1～Ｐｂ3のうち中央に位置する光パルスＰｂ2についてのターゲット波形関数であるので、位相スペクトル関数Ｇ１３２の傾き（すなわち時間シフト量）はゼロとされている。
【0067】
フーリエ変換部１０２は、複数の光パルスそれぞれのターゲット波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の波形関数に変換する。図２０（ｂ）は、時間領域の波形関数の一例として、図２０（ａ）に示されるターゲット波形関数から生成された時間強度波形関数Ｇ１４１及び時間位相波形関数Ｇ１４２を示すグラフである。図２０（ｂ）において、横軸は時間（ｆｓ）を示し、縦軸は時間強度波形関数Ｇ１４１の強度値（任意単位）及び時間位相波形関数Ｇ１４２の位相値（ｒａｄ）を示す。
【0068】
スペクトログラム変換部１０３は、複数の光パルスそれぞれの時間領域の波形関数から強度スペクトログラムを生成する。図２１（ａ）は、強度スペクトログラムの一例として、図２０（ｂ）に示される波形関数から生成された強度スペクトログラムを示す図である。なお、強度スペクトログラムの作成方法及びスペクトログラムの定義は、前述した予備データ生成部２３のスペクトログラム変換部２６１における説明と同様である。また、スペクトログラム変換部１０３においても、ＳＴＦＴに限らず、他の時間－周波数変換（例えばウェーブレット変換）を用いてよい。
【0069】
ターゲット生成部１０４は、複数の光パルスそれぞれの強度スペクトログラムを互いに重畳して強度スペクトログラムＡ４３を生成する。図２１（ｂ）は、一例として、３つの光パルスＰｂ1～Ｐｂ3についての強度スペクトログラムを互いに重畳して得られる強度スペクトログラムＡ４３を示す図である。なお、ターゲット生成部１０４は、複数の光パルスそれぞれの強度スペクトログラムを重畳したものに対して、補正のための係数を乗算してもよい。補正のための係数は、例えば、作成したターゲット強度スペクトログラムＡ４３のスペクトル強度分布を、入力光Ｌａの光パルスのスペクトル強度分布に近づけるための係数である。
【0070】
図２２は、強度スペクトログラムＡ４３を作成するための方法を示すフローチャートである。本実施形態のターゲット強度スペクトログラム作成プログラムは、下記の各ステップをコンピュータに実行させる。
【0071】
まず、波形関数設定ステップＳ１では、複数の光パルスのうち一の光パルスについて、強度スペクトル関数及び位相スペクトル関数を含む周波数領域のターゲット波形関数を設定する。フーリエ変換ステップＳ２では、上記一の光パルスのターゲット波形関数を、時間強度波形関数及び時間位相波形関数を含む時間領域の波形関数に変換する。スペクトログラム変換ステップＳ３では、フーリエ変換ステップＳ２において生成された時間領域の波形関数から強度スペクトログラムを生成する。ターゲット生成ステップＳ４では、スペクトログラム変換ステップＳ３において生成された強度スペクトログラムをターゲット強度スペクトログラムに重畳する。以上の波形関数設定ステップＳ１、フーリエ変換ステップＳ２、スペクトログラム変換ステップＳ３、及びターゲット生成ステップＳ４を、光パルスの本数と同じ繰り返し数にわたって繰り返す（ステップＳ５）。これにより、強度スペクトログラムＡ４３が作成される。なお、この例に限らず、例えば波形関数設定ステップＳ１、フーリエ変換ステップＳ２、及びスペクトログラム変換ステップＳ３のそれぞれにおいて複数の光パルスについての処理をまとめて行ったのち、生成された複数の強度スペクトログラムをターゲット生成ステップＳ４において重畳してもよい。
【0072】
以上に説明したスペクトログラム設定部２２によれば、光パルス列Ｐｂに含まれる複数のパルス毎に中心波長を異ならせるためのターゲット強度スペクトログラムＡ４３を好適に作成することができる。
【0073】
図２３～図２５は、ターゲット強度スペクトログラムＡ４３を用いて得られたデータをＳＬＭ１４に呈示することにより実際に得られる光の、複数のパルス間の中心波長差と、複数のパルスのピーク強度のばらつきとの関係を、複数の初期位相スペクトル関数Φ0（ω）それぞれについてプロットしたグラフである。図２３～図２５において、縦軸は複数のパルス間の中心波長差の平均値（ｎｍ）を表し、横軸は複数のパルスのピーク強度の標準偏差（任意単位）を表す。図２３は、ターゲット強度スペクトログラムＡ４３を作成する際に、複数のパルス間の中心波長差を２ｎｍに設定した場合を示す。図２４及び図２５は、ターゲット強度スペクトログラムＡ４３を作成する際に、複数のパルス間の中心波長差を、それぞれ２．５ｎｍ及び３ｎｍに設定した場合を示す。なお、複数のパルス間の中心波長差の設定値は、図中において星印により示されている。
【0074】
図２３～図２５を参照すると、初期位相スペクトル関数Φ0（ω）に依存して、複数のパルス間の中心波長差およびピーク強度のばらつきが変化することがわかる。そして、中心波長差が或る値になるときにピーク強度のばらつきが最小となり、その値から中心波長差が離れるほど、ピーク強度のばらつきが大きくなることがわかる。
【0075】
加えて、図２３～図２５を参照すると、ピーク強度のばらつきが最小となるときの複数のパルス間の中心波長差の値が、ターゲット強度スペクトログラムＡ４３を作成する際に設定された中心波長差の値よりも有意に小さいことがわかる。例えば、図２３に示されるように、ターゲット強度スペクトログラムＡ４３を作成する際に設定された中心波長差が２ｎｍである場合、ピーク強度のばらつきが最小となるときの複数のパルス間の中心波長差の値は１．７ｎｍとなっており、設定された中心波長差（２ｎｍ）よりも１５％小さい。また、図２４に示されるように、ターゲット強度スペクトログラムＡ４３を作成する際に設定された中心波長差が２．５ｎｍである場合、ピーク強度のばらつきが最小となるときの複数のパルス間の中心波長差の値は２．１ｎｍとなっており、設定された中心波長差（２．５ｎｍ）よりも１６％小さい。また、図２５に示されるように、ターゲット強度スペクトログラムＡ４３を作成する際に設定された中心波長差が３ｎｍである場合、ピーク強度のばらつきが最小となるときの複数のパルス間の中心波長差の値は２．３ｎｍとなっており、設定された中心波長差（３ｎｍ）よりも２３％小さい。
【0076】
このことから、スペクトログラム設定部２２がターゲット強度スペクトログラムＡ４３を作成する際には、ターゲット強度スペクトログラムＡ４３における複数の光パルス間の中心波長差を、目標とする複数の光パルス間の中心波長差よりも大きく（例えば、目標とする複数の光パルス間の中心波長差の１．１倍よりも大きく）設定すればよいことがわかる。目標とする複数の光パルス間の中心波長差は、例えば、図８に示された入力デバイス６４を通じて操作者により入力される。ターゲット強度スペクトログラムＡ４３における複数の光パルス間の中心波長差をこのように設定することにより、光パルス列Ｐｂにおける複数の光パルス間の中心波長差をより精度良く実現することができる。
【0077】
なお、このようにターゲット強度スペクトログラムＡ４３における複数の光パルス間の中心波長差を目標とする中心波長差よりも大きくする場合、その中心波長差は、ＳＬＭ１４への入力光Ｌａの波長帯域を（光パルス列Ｐｂのパルス本数－１）で除算した値よりも小さいことが好ましい。例えば、光パルス列Ｐｂのパルス本数が３である場合、ターゲット強度スペクトログラムＡ４３における複数の光パルス間の中心波長差は、入力光Ｌａの波長帯域の１／２よりも小さいことが好ましい。ターゲット強度スペクトログラムＡ４３における複数の光パルス間の中心波長差をこのように設定することにより、光パルス列Ｐｂにおける複数の光パルス間の中心波長差をより精度良く実現することができる。一例では、入力光Ｌａのスペクトルの波長帯域幅が１０ｎｍであり、光パルス列Ｐｂのパルス本数が３である場合、中心波長差の設定値は５ｎｍが上限となる。
【0078】
［変形例］
続いて、スペクトログラム設定部２２及び予備データ生成部２３の変形例について説明する。図２６は、一変形例に係るターゲット設定部３５及び予備データ生成部３０の内部構成を示すブロック図である。図２６に示されるように、予備データ生成部３０は、フーリエ変換部３１（第１変換部）、関数変換部３２（第２変換部）、逆フーリエ変換部３３（第３変換部）、及びデータ生成部３４を有する。関数変換部３２は、関数置換部３２１及び波形関数修正部３２２を含む。また、ターゲット設定部３５は、フーリエ変換部３５１及びスペクトログラム修正部３５２を有する。これらの各構成要素の機能については、後に詳述する。
続いて、スペクトログラム設定部２２及び予備データ生成部２３の変形例について説明する。図２６は、一変形例に係るターゲット設定部３５及び予備データ生成部３０の内部構成を示すブロック図である。図２６に示されるように、予備データ生成部３０は、フーリエ変換部３１（第１変換部）、関数変換部３２（第２変換部）、逆フーリエ変換部３３（第３変換部）、及びデータ生成部３４を有する。関数変換部３２は、関数置換部３２１及び波形関数修正部３２２を含む。また、ターゲット設定部３５は、フーリエ変換部３５１及びスペクトログラム修正部３５２を有する。これらの各構成要素の機能については、後に詳述する。
【0079】
図２７は、予備データ生成部３０における位相スペクトル関数の計算手順を示す図である。まず、周波数ωの関数である初期強度スペクトル関数Ａ0（ω）及び初期位相スペクトル関数Φ0（ω）を用意する（図中の処理番号（１））。一例では、初期強度スペクトル関数Ａ0（ω）は入力光Ｌａの強度スペクトルを表すが、これに限られない。また、初期位相スペクトル関数Φ0（ω）は、上記実施形態と同様に、記憶部２１に記憶された複数の初期位相スペクトル関数Φ0（ω）の中から順に選択される。
【0080】
次に、初期強度スペクトル関数Ａ0（ω）及び初期位相スペクトル関数Φ0（ω）を含む周波数領域の第１波形関数（７）を用意する（処理番号（２－ａ））。但し、ｉは虚数である。
【数7】
【0081】
続いて、フーリエ変換部３１は、上記関数（７）に対して周波数領域から時間領域へのフーリエ変換を行う（図中の矢印Ｃ１）。これにより、時間強度波形関数ａ0（ｔ）及び時間位相波形関数φ0（ｔ）を含む時間領域の第２波形関数（８）が得られる（フーリエ変換ステップ、図中の処理番号（３））。
【数8】
【0082】
続いて、関数変換部３２の関数置換部３２１は、次の数式（９）に示されるように、時間強度波形関数ｂ0（ｔ）に、目標とする時間波形を表す時間強度波形関数Ｔａｒｇｅｔ0（ｔ）を代入する（処理番号（４－ａ））。
【数9】
【0083】
続いて、関数変換部３２の関数置換部３２１は、次の数式（１０）に示されるように、時間強度波形関数ａ0（ｔ）を時間強度波形関数ｂ0（ｔ）で置き換える。すなわち、上記関数（８）に含まれる時間強度波形関数ａ0（ｔ）を、目標とする時間波形に基づく時間強度波形関数Ｔａｒｇｅｔ0（ｔ）に置き換える（関数置換ステップ、図中の処理番号（５））。
【数10】
【0084】
続いて、関数変換部３２の波形関数修正部３２２は、置き換え後の第２波形関数（１０）のスペクトログラムが、所望の波長帯域に従って予め生成されたターゲットスペクトログラムに近づくように、第２波形関数を修正する。まず、置き換え後の第２波形関数（１０）に対して時間－周波数変換を施すことにより、第２波形関数（１０）をスペクトログラムＳＧ0,k（ω，ｔ）に変換する（図中の処理番号（５－ａ））。添え字ｋは、第ｋ回目の変換処理を表す。
【0085】
また、所望の波長帯域に従って予め生成されたターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）をターゲット設定部３５から読み出す。このターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）は、目標とする時間波形（時間強度波形とそれを構成する周波数成分）と概ね同値であり、処理番号（５－ｂ）のターゲットスペクトログラム関数において生成される。
【0086】
次に、関数変換部３２の波形関数修正部３２２は、スペクトログラムＳＧ0,k（ω，ｔ）とターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）とのパターンマッチングを行い、類似度（どの程度一致しているか）を調べる。本実施形態では、類似度を表す指標として、評価値を算出する。そして、続く処理番号（５－ｃ）では、得られた評価値が、所定の終了条件を満たすか否かの判定を行う。条件を満たせば処理番号（６）へ進み、満たさなければ処理番号（５－ｄ）へ進む。処理番号（５－ｄ）では、第２波形関数に含まれる時間位相波形関数φ0（ｔ）を任意の時間位相波形関数φ0,k（ｔ）に変更する。時間位相波形関数を変更した後の第２波形関数は、ＳＴＦＴなどの時間－周波数変換により再びスペクトログラムに変換される。以降、上述した処理番号（５－ａ）～（５－ｃ）が繰り返し行われる。こうして、スペクトログラムＳＧ0,k（ω，ｔ）がターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）に次第に近づくように、第２波形関数が修正される（波形関数修正ステップ）。
【0087】
その後、逆フーリエ変換部３３は、修正後の第２波形関数に対して逆フーリエ変換を行い（図中の矢印Ｃ２）、周波数領域の第３波形関数（１１）を生成する（逆フーリエ変換ステップ、図中の処理番号（６））。
この第３波形関数（１１）に含まれる位相スペクトル関数Φ0,k（ω）が、最終的に得られる所望の位相スペクトル関数ΦTWC-TFD（ω）となる。この位相スペクトル関数ΦTWC-TFD（ω）が、データ生成部３４に提供される。
【数11】
この第３波形関数（１１）に含まれる位相スペクトル関数Φ0,k（ω）が、最終的に得られる所望の位相スペクトル関数ΦTWC-TFD（ω）となる。この位相スペクトル関数ΦTWC-TFD（ω）が、データ生成部３４に提供される。
【0088】
データ生成部３４は、位相スペクトル関数ΦTWC-TFD（ω）により示されるスペクトル位相を入力光Ｌａに与えるための位相変調パターン（例えば、計算機合成ホログラム）を算出する（データ生成ステップ）。
【0089】
ここで、図２８は、ターゲット設定部３５におけるターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）の生成手順の一例を示す図である。ターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）は、目標とする時間波形（時間強度波形とそれを構成する周波数成分（波長帯域成分））を示すので、ターゲットスペクトログラムの作成は、周波数成分（波長帯域成分）を制御するために極めて重要な工程である。図２８に示されるように、ターゲット設定部３５は、まずスペクトル波形（初期の強度スペクトル関数Ａ0（ω）及び初期の位相スペクトル関数Φ0（ω））、並びに所望の時間強度波形関数Ｔａｒｇｅｔ0（ｔ）を入力する。また、所望の周波数（波長）帯域情報を含む時間関数ｐ0（ｔ）を入力する（処理番号（１））。
【0090】
次に、ターゲット設定部３５は、例えば一般的な反復フーリエ変換法、或いは非特許文献１または２に記載された方法を用いて、時間強度波形関数Ｔａｒｇｅｔ0（ｔ）を実現するための位相スペクトル関数ΦIFTA（ω）を算出する（処理番号（２））。
【0091】
続いて、ターゲット設定部３５は、先に得られた位相スペクトル関数ΦIFTA（ω）を利用した反復フーリエ変換法により、時間強度波形関数Ｔａｒｇｅｔ0（ｔ）を実現するための強度スペクトル関数ＡIFTA（ω）を算出する（処理番号（３））。ここで、図２９は、強度スペクトル関数ＡIFTA（ω）を算出する手順の一例を示す図である。
【0092】
まず、初期の強度スペクトル関数Ａk=0（ω）及び位相スペクトル関数Ψ0（ω）を用意する（図中の処理番号（１））。次に、強度スペクトル関数Ａk（ω）及び位相スペクトル関数Ψ0（ω）を含む周波数領域の波形関数（１２）を用意する（図中の処理番号（２））。
添え字ｋは、第ｋ回目のフーリエ変換処理後を表す。最初（第１回目）のフーリエ変換処理の前においては、強度スペクトル関数Ａk（ω）として上記の初期強度スペクトル関数Ａk=0（ω）が用いられる。ｉは虚数である。
【数12】
添え字ｋは、第ｋ回目のフーリエ変換処理後を表す。最初（第１回目）のフーリエ変換処理の前においては、強度スペクトル関数Ａk（ω）として上記の初期強度スペクトル関数Ａk=0（ω）が用いられる。ｉは虚数である。
【0093】
続いて、上記関数（１２）に対して周波数領域から時間領域へのフーリエ変換を行う（図中の矢印Ｃ３）。これにより、時間強度波形関数ｂk（ｔ）を含む周波数領域の波形関数（１３）が得られる（図中の処理番号（３））。
【数13】
【0094】
続いて、上記関数（１３）に含まれる時間強度波形関数ｂk（ｔ）を、所望の波形に基づく時間強度波形関数Ｔａｒｇｅｔ0（ｔ）に置き換える（図中の処理番号（４）、（５））。
【数14】
【数15】
【0095】
続いて、上記関数（１５）に対して時間領域から周波数領域への逆フーリエ変換を行う（図中の矢印Ｃ４）。これにより、強度スペクトル関数Ｃk（ω）及び位相スペクトル関数Ψk（ω）を含む周波数領域の波形関数（１６）が得られる（図中の処理番号（６））。
続いて、上記関数（１６）に含まれる位相スペクトル関数Ψk（ω）を拘束するため、初期の位相スペクトル関数Ψ0（ω）に置き換える（図中の処理番号（７－ａ））。
【数16】
続いて、上記関数（１６）に含まれる位相スペクトル関数Ψk（ω）を拘束するため、初期の位相スペクトル関数Ψ0（ω）に置き換える（図中の処理番号（７－ａ））。
【数17】
【0096】
また、逆フーリエ変換後の周波数領域における強度スペクトル関数Ｃk（ω）に対し、入力光Ｌａの強度スペクトルに基づくフィルタ処理を行う。具体的には、強度スペクトル関数Ｃk（ω）により表される強度スペクトルのうち、入力光Ｌａの強度スペクトルに基づいて定められる波長毎のカットオフ強度を超える部分をカットする。一例では、波長毎のカットオフ強度は、入力光Ｌａの強度スペクトル（例えば初期の強度スペクトル関数Ａk=0（ω））と一致するように設定される。その場合、次の数式（１８）に示されるように、強度スペクトル関数Ｃk（ω）が強度スペクトル関数Ａk=0（ω）よりも大きい周波数では、強度スペクトル関数Ａk（ω）の値として強度スペクトル関数Ａk=0（ω）の値が取り入れられる。また、強度スペクトル関数Ｃk（ω）が強度スペクトル関数Ａk=0（ω）以下である周波数では、強度スペクトル関数Ａk（ω）の値として強度スペクトル関数Ｃk（ω）の値が取り入れられる（図中の処理番号（７－ｂ））。
上記関数（１６）に含まれる強度スペクトル関数Ｃk（ω）を、上記数式（１８）によるフィルタ処理後の強度スペクトル関数Ａk（ω）に置き換える。
【数18】
上記関数（１６）に含まれる強度スペクトル関数Ｃk（ω）を、上記数式（１８）によるフィルタ処理後の強度スペクトル関数Ａk（ω）に置き換える。
【0097】
以降、上記の処理（１）～（７－ｂ）を繰り返し行うことにより、波形関数中の強度スペクトル関数Ａk（ω）が表す強度スペクトル形状を、所望の時間強度波形に対応する強度スペクトル形状に近づけることができる。最終的に、強度スペクトル関数ＡIFTA（ω）が得られる。
【0098】
再び図２８を参照する。以上に説明した処理番号（２）、（３）における位相スペクトル関数ΦIFTA（ω）及び強度スペクトル関数ＡIFTA（ω）の算出によって、これらの関数を含む周波数領域の第３波形関数（１９）が得られる（図中の処理番号（４））。
ターゲット設定部３５のフーリエ変換部３５１は、上の波形関数（１９）をフーリエ変換する。これにより、時間領域の第４波形関数（２０）が得られる（図中の処理番号（５））。
【数19】
ターゲット設定部３５のフーリエ変換部３５１は、上の波形関数（１９）をフーリエ変換する。これにより、時間領域の第４波形関数（２０）が得られる（図中の処理番号（５））。
【数20】
【0099】
ターゲット設定部３５のスペクトログラム修正部３５２は、時間－周波数変換により第４波形関数（２０）をスペクトログラムＳＧIFTA（ω，ｔ）に変換する（処理番号（６））。そして、処理番号（７）では、所望の周波数（波長）帯域情報を含む時間関数ｐ0（ｔ）を基にスペクトログラムＳＧIFTA（ω，ｔ）を修正することにより、ターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）を生成する。例えば、２次元データにより構成されるスペクトログラムＳＧIFTA（ω，ｔ）に現れる特徴的パターンを部分的に切り出し、時間関数ｐ0（ｔ）を基に当該部分の周波数成分の操作を行う。以下、その具体例について詳細に説明する。
【0100】
例えば、所望の時間強度波形関数Ｔａｒｇｅｔ0（ｔ）として時間間隔が２ピコ秒であるトリプルパルスを設定した場合について考える。このとき、スペクトログラムＳＧIFTA（ω，ｔ）は、図７（ａ）に示されるような結果となる。仮に光パルス列Ｐｂの時間強度波形のみを制御したい（単にトリプルパルスを得たい）場合には、スペクトログラムＳＧIFTA（ω，ｔ）に含まれるパルスＰｂ1，Ｐｂ2，Ｐｂ3の波長帯域を操作する必要はない。しかし、各パルスの波長帯域を制御したい場合には、これらのパルスＰｂ1，Ｐｂ2，Ｐｂ3の操作が必要となる。すなわち、図６（ａ）に示されるように、波長軸（縦軸）に沿った方向に各パルスＰｂ1，Ｐｂ2，Ｐｂ3を互いに独立して移動させる。このような各パルスの波長帯域の変更は、時間関数ｐ0（ｔ）を基に行われる。
【0101】
例えば、パルスＰｂ2のピーク波長を８００ｎｍで据え置き、パルスＰｂ1及びＰｂ3のピーク波長がそれぞれ－２ｎｍ、＋２ｎｍだけ平行移動するように時間関数ｐ0（ｔ）を記述するとき、スペクトログラムＳＧIFTA（ω，ｔ）は、図６（ａ）に示されるターゲットスペクトログラムＴａｒｇｅｔＳＧ0（ω，ｔ）に変化する。例えばスペクトログラムにこのような処理を施すことによって、時間強度波形の形状を変えずに、各パルスの波長帯域が任意に制御されたターゲットスペクトログラムを作成することができる。
【0102】
本開示によるデータ作成装置、データ作成方法、及びデータ作成プログラムは、上述した実施形態に限られるものではなく、他に様々な変形が可能である。例えば、上記実施形態では、互いに時間差を有し中心波長が互いに異なる複数の光パルスを含む光パルス列を生成する場合について主に説明したが、互いに時間差を有し中心波長が互いに等しい複数の光パルスを含む光パルス列を生成する場合であっても本開示によるデータ作成装置、データ作成方法、及びデータ作成プログラムは有効である。また、単一の光パルスの時間波形を成形する場合など、光パルス列以外の光の波形制御においても本開示によるデータ作成装置、データ作成方法、及びデータ作成プログラムは有効である。
【符号の説明】
【0103】
１…光制御装置，２…光源，３…光波形制御部，１２…回折格子，１３…レンズ，１４…空間光変調器（ＳＬＭ），１５…レンズ，１６…回折格子，１７…変調面，１７ａ…変調領域，２０…データ作成装置，２１…記憶部，２２…スペクトログラム設定部，２３，３０…予備データ生成部，２４…データ選択部，２５，３１…フーリエ変換部（第１変換部），２６，３２…関数変換部（第２変換部），２７，３３…逆フーリエ変換部（第３変換部），２８，３４…データ生成部，３５…ターゲット設定部，６１…プロセッサ，６２…ＲＯＭ，６３…ＲＡＭ，６４…入力デバイス，６５…出力デバイス，６６…通信モジュール，６７…補助記憶装置，１０１…波形関数設定部，１０２…フーリエ変換部，１０３，２６１…スペクトログラム変換部，１０４…ターゲット生成部，２６２…スペクトログラム置換部，２６３…スペクトログラム逆変換部，３２１…関数置換部，３２２…波形関数修正部，３５１…フーリエ変換部，３５２…スペクトログラム修正部，Ａ１…初期スペクトル関数，Ａ２…第１波形関数，Ａ３…第２波形関数，Ａ５…第３波形関数，Ａ６…第４波形関数，Ａ４１…強度スペクトログラム，Ａ４２…位相スペクトログラム，Ａ４３…強度スペクトログラム（ターゲット強度スペクトログラム），ＡＡ，ＡＢ…方向，Ｌ１…光，Ｌ２…変調光，Ｌａ…入力光，Ｌｂ…出力光，Ｐｂ，Ｐｄ…光パルス列，Ｐｂ０…光パルス，Ｐｂ１～Ｐｂ３，Ｐｄ１～Ｐｄ３…光パルス。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】