CN104272072A - 波形重构装置、波形重构系统以及波形重构方法 - Google Patents
波形重构装置、波形重构系统以及波形重构方法 Download PDFInfo
- Publication number
- CN104272072A CN104272072A CN201380010306.9A CN201380010306A CN104272072A CN 104272072 A CN104272072 A CN 104272072A CN 201380010306 A CN201380010306 A CN 201380010306A CN 104272072 A CN104272072 A CN 104272072A
- Authority
- CN
- China
- Prior art keywords
- mentioned
- optical signal
- phase frequency
- input optical
- frequency spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 50
- 238000001228 spectrum Methods 0.000 claims abstract description 345
- 230000003287 optical effect Effects 0.000 claims abstract description 226
- 230000000694 effects Effects 0.000 claims abstract description 117
- 238000004088 simulation Methods 0.000 claims abstract description 77
- 239000006185 dispersion Substances 0.000 claims description 75
- 230000008859 change Effects 0.000 claims description 48
- 239000013307 optical fiber Substances 0.000 claims description 44
- 230000001902 propagating effect Effects 0.000 claims description 20
- 238000005259 measurement Methods 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 238000009826 distribution Methods 0.000 abstract description 5
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000002123 temporal effect Effects 0.000 abstract 1
- 230000006870 function Effects 0.000 description 27
- 238000011156 evaluation Methods 0.000 description 22
- 238000002474 experimental method Methods 0.000 description 16
- 238000001069 Raman spectroscopy Methods 0.000 description 14
- 238000013459 approach Methods 0.000 description 13
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 13
- 230000008569 process Effects 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 9
- 230000004044 response Effects 0.000 description 9
- 238000004891 communication Methods 0.000 description 8
- 230000000644 propagated effect Effects 0.000 description 8
- 230000009471 action Effects 0.000 description 7
- 238000000137 annealing Methods 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- 230000000052 comparative effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 229910052691 Erbium Inorganic materials 0.000 description 1
- 244000287680 Garcinia dulcis Species 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000003292 diminished effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- UYAHIZSMUZPPFV-UHFFFAOYSA-N erbium Chemical compound [Er] UYAHIZSMUZPPFV-UHFFFAOYSA-N 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
- H04B10/075—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal
- H04B10/079—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using an in-service signal using measurements of the data signal
- H04B10/0795—Performance monitoring; Measurement of transmission parameters
- H04B10/07955—Monitoring or measuring power
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J11/00—Measuring the characteristics of individual optical pulses or of optical pulse trains
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/25—Arrangements specific to fibre transmission
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
波形重构装置(140)具备:相位频谱算出部(143),(a)按输入光信号的每个强度,设想为输入光信号具有给定的相位频谱,来进行输入光信号在光传输媒体内的传播模拟,从而算出输出光信号的计算功率谱,(b)以使计算功率谱和测量功率谱的按输入光信号的每个强度的差异变小的方式,使给定的相位频谱变化,进行传播模拟,从而探索输入光信号的相位频谱;波形重构部(144),利用探索出的相位频谱来对输入光信号的时间波形进行重构,相位频谱算出部(143)只根据非线性光学效果以及色散效果中的一方,来使给定的相位频谱变化,或者进行传播模拟。
Description
技术领域
本发明涉及重构光信号的时间波形的波形重构装置等。
背景技术
近年来,为了实现利用非线性光学效果的信息通信系统等的实用化,光信号的正确的时间波形信息的获取变得极为重要。对此,为了取得光信号的正确的时间波形信息,可利用光取样示波器、自相关器等,来测量光信号的时间波形的强度分布。其结果,根据测量出的强度分布,取得光信号的振幅信息。但是,根据测量出的强度分布,无法取得光信号的相位信息。即,只对光信号的强度分布进行测量,并不能取得光信号的时间波形的信息。
对此,已有各种各样用于取得光信号相位的方法被提出(例如,参照非专利文献1以及2)。非专利文献1以及2所记载的方法是根据时间分解分光来取得光信号的相位的方法。具体而言,非专利文献1以及2所记载的方法中,通过利用超高速时间门或者参照光源,来取得光信号的相位。然后,利用所取得的相位,重构光信号的时间波形。
现有技术文献
非专利文献
非专利文献1:D.J.Kane,R.Trebino,“Characterization of arbitraryfemtosecond pulses using frequency-resolved optical gating”,IEEEJ.Quantum Electron,Vol.29,1993,pp.571~pp.579
非专利文献2:C.Dorrer,M.Joffre,“Characterization of the spectralphase of ultrashort light pulses”,C.R.Acad.Sci.Paris,Vol.2,2001,pp.1415
发明概要
发明要解决的问题
但是,上述现有的方法需要超高速时间门或者参照光源,因此要求光信号和门的时间调整以及稳定性和SN(Signal-Noise)比的确保等非常高的技术水平。
发明内容
因此,本发明提供一种能够简单地重构光信号的时间波形的波形重构装置。
用于解决问题的手段
本发明的一形态的波形重构装置,对输入光信号的时间波形进行重构,该波形重构装置具备:输入频谱取得部,取得表示上述输入光信号的功率谱的信息;输出频谱取得部,对于多个强度的上述输入光信号,将通过对上述输入光信号在具有非线性光学效果的光传输媒体内传播之后被输出的输出光信号进行测量而获得的功率谱,分别作为测量功率谱而取得;相位频谱算出部,(a)按上述输入光信号的每个强度,设想为上述输入光信号具有给定的相位频谱,来进行上述输入光信号在上述光传输媒体内的传播模拟,从而算出上述输出光信号的计算功率谱,(b)以使上述计算功率谱和上述测量功率谱的按上述输入光信号的每个强度的差异变小的方式,使上述给定的相位频谱变化,并进行上述传播模拟,从而探索上述输入光信号的相位频谱;波形重构部,通过对由上述相位频谱算出部探索出的相位频谱和由上述输入频谱取得部取得的信息所表示的功率谱,进行频率/时间变换,从而对上述输入光信号的时间波形进行重构,上述相位频谱算出部,只根据非线性光学效果以及色散效果中的一方,来使上述给定的相位频谱变化,或者进行上述传播模拟。
根据该结构,能够利用输入光信号在具有非线性光学效果的光传输媒体内传播之后被输出的输出光信号的功率谱,来重构输入光信号的时间波形。即,如果能够测量出功率谱就能够重构输入光信号的时间波形,因此能够简单地重构光信号的时间波形。并且,能够只根据非线性光学效果以及色散效果中的一方,来使给定的相位频谱变化,或者进行传播模拟。因此,与根据或者不根据非线性光学效果以及色散效果的两者来进行传播模拟,或者使给定的相位频谱变化的情况相比,能够缩短相位频谱的探索时间。
另外,例如可以是,上述相位频谱算出部,按照在设想为上述输入光信号是在只具有非线性光学效果以及色散效果中的色散效果的媒体内传播而获得的信号的情况下应被满足的相位频谱的限制条件,使上述给定的相位频谱变化。
根据该结构,能够按照适当的限制条件来使给定的相位频谱变化。因此,在相位频谱的探索中,能够减少变量的个数,从而能够缩短探索时间。
另外,例如可以是,上述相位频谱算出部,在探索上述输入光信号的相位频谱时,在按照上述限制条件使上述给定的相位频谱变化并进行上述传播模拟之后,无论上述限制条件如何,都使上述给定的相位频谱变化并进行上述传播模拟。
根据该结构,能够在按照限制条件使给定的相位频谱变化之后,无论限制条件如何都使给定的相位频谱变化。因此,在相位频谱的探索中,能够抑制探索精度降低的同时缩短探索时间。
另外,例如可以是,上述相位频谱算出部只利用与非线性光学效果和色散效果中的非线性光学效果相关的参数,来进行上述传播模拟。
根据该结构,能够只利用与非线性光学效果和色散效果中的非线性光学效果相关的参数来进行传播模拟。因此,能够简化传播模拟,从而能够提高计算速度。
另外,例如可以是,上述相位频谱算出部,只利用与上述非线性光学效果相关的参数中的与自相位调制相关的参数,来进行上述传播模拟。
根据该结构,能够只利用与非线性光学效果相关的参数中的与自相位调制相关的参数来进行传播模拟。因此,能够进一步简化传播模拟,从而能够提高计算速度。
另外,上述光传输媒体例如可以是光纤。
根据该结构,能够简单地构成具有非线性光学效果的光传输媒体。
另外,本发明的一形态的波形重构系统,对输入光信号的时间波形进行重构,该波形重构系统具备:强度调节器,使上述输入光信号的强度变化成互不相同的多个强度;具有非线性光学效果的光传输媒体,至少该光传输媒体的与非线性光学效果相关的参数是已知的;频谱测量器,按上述输入光信号的每个强度,对由上述强度调节器进行了强度变更的上述输入光信号在上述光传输媒体内传播之后被输出的输出光信号的功率谱进行测量;上述波形重构装置。上述输出频谱取得部取得由上述频谱测量器测量出功率谱,以此作为测量功率谱。
根据该结构,如果有与非线性光学效果相关的参数是已知的光传输媒体和频谱测量器,就能够重构输入光信号的时间波形,从而能够简单地重构光信号的时间波形。并且,能够只根据非线性光学效果以及色散效果中的一方,来使给定的相位频谱变化,并进行传播模拟。因此,与根据或者不根据非线性光学效果以及色散效果的两者进行传播模拟,或者使给定的相位频谱变化的情况相比,能够缩短相位频谱的探索时间。
另外可以是,上述光传输媒体可以是光纤。
根据该结构,能够简单地构成具有非线性光学效果的光传输媒体。
在此,这些总括性或者具体性的形态可由系统、方法、集成电路、计算机程序或者计算机可读取的CD-ROM等记录介质来实现,亦可由系统、方法、集成电路、计算机程序或者记录介质的任意组合来实现。
发明效果
从以上的说明可明确看出,根据本发明的一形态,能够简单地重构光信号的时间波形。
附图说明
图1是表示实施方式1的波形重构系统的整体结构的图。
图2是表示实施方式1的波形重构装置的特征性功能结构的方框图。
图3是表示实施方式1的波形重构装置的处理动作的流程图。
图4是用于说明实施方式1的分步傅立叶法的图。
图5是表示实施方式1中由分光器测量出的输出光信号的功率谱的实验结果的图表。
图6是表示实施方式1的相位频谱探索中的评价函数值的变化的实验结果的图表。
图7是表示实施方式1中通过实验获得的功率谱和相位频谱的图。
图8是表示实施方式1中由波形重构装置重构的输入光信号的时间波形的实验结果的图。
图9是表示实施方式1中的与7种类的输入光信号的强度对应的输出光信号的功率谱的实测值和计算值的比较结果的图表。
图10是表示实施方式1中由分光器测量出的与10种类的输入光信号的强度对应的输出光信号的功率谱中的,与算出相位频谱时未利用的3种类的输入光信号的强度对应的输出光信号的功率谱的实测值和计算值的比较结果的图表。
图11是表示实施方式2的变形例的实验中使用的输入光信号的功率谱的图表。
图12是表示在实施方式2的变形例的实验中,2个强度的输入光信号在高非线性光纤中传播之后被分光器测量出的功率谱的图表。
图13是表示实施方式2的变形例的波形重构装置的实验结果的图。
图14是表示实施方式2的变形例的波形重构装置的实验结果的图。
图15是表示实施方式的变形例的波形重构装置的实验结果的图。
图16是表示计算机的硬件结构的一个例子的图。
具体实施方式
本发明的一形态的波形重构系统以及波形重构装置,利用能够引起某多频谱变化组合的时间波形是被限定的这一非线性光学效果的特征,进行输入光信号的时间波形的重构。具体是,波形重构系统以及波形重构装置,对输入光信号的强度发生变化时显现的、非线性光学效果为因的多个频谱变化的实测值和计算值进行比较,并利用该比较结果来进行输入光信号的时间波形的重构。通过使按照退火法等最佳化运演算法进行模拟而计算出的计算功率谱(计算值)接近测量功率谱(实测值)的方式,来使相位频谱变化,从而决定时间波形的重构所需要的相位频谱。
以下,关于实施方式,参照附图进行说明。在此,以下说明的实施方式均为总括性或者具体性的一个例子。以下的实施方式中给出的数值、形状、材料、结构要素、结构要素的配置以及连接形态、步骤,步骤的顺序等也都表示一个例子,并不意味权利要求的范围限定于此。另外,关于以下的实施方式中说明的结构要素中的未被记载于表示最上位概念的独立权利要求中的结构要素,视其为任意的结构要素。
(实施方式1)
<波形重构系统的结构>
图1是表示实施方式1的波形重构系统100的整体结构的图。该波形重构系统100是对输入光信号的时间波形进行重构的系统。如图1所示,波形重构系统100具备强度调节器110、高非线性光纤120、分光器130以及波形重构装置140。
强度调节器110使由光信号生成装置200生成的输入光信号的强度变化。
高非线性光纤120是具有非线性光学效果的光传输媒体的一个例子。被强度调节器110进行了强度变化之后的输入光信号,在该高非线性光纤120内传播。
高非线性光纤120是非线性光学效果的相关参数为已知的光纤。在本实施方式中,高非线性光纤120是色散效果的相关参数也为已知的光纤。具体而言,例如,除了2阶以及3阶之外,与4阶色散相关的参数也是已知的。
分光器130是频谱测量器的一个例子,通过将输出光信号分解成每个波长的光,并对分解成了每个波长的光进行O/E变换以及A/D变换,从而生成由数字值表示的功率谱。即,分光器130对输出光信号的功率谱进行测量。在此所说的输出光信号是指输入光信号在高非线性光纤120内传播之后被输出的光信号。
波形重构装置140例如由图16所示的计算机实现,进行输入光信号的时间波形重构。关于波形重构装置140,以下将参照图2进行详细说明。
光信号生成装置200生成输入光信号。具体而言,光信号生成装置200例如具备MLLD(Mode-Locked Laser Diode)、SMF(Single Mode Fiber)和EDFA(Erbium Doped Fiber Amplifier)。光信号生成装置200对于由MLLD输出的光脉冲,通过SMF进行色散补偿,并通过EDFA进行增幅。
<波形重构装置的结构>
图2是表示实施方式1的波形重构装置140的特征性功能结构的方框图,如图2所示,具备输入频谱取得部141、输出频谱取得部142、相位频谱算出部143以及波形重构部144。
输入频谱取得部141取得表示由光信号生成装置200生成的输入光信号的功率谱的信息。例如,在由光信号生成装置200生成的输入光信号的功率谱是已知的情况下,输入频谱取得部141通过读取被存放在存储单元等中的功率谱的数据,来取得表示输入光信号的功率谱的信息。相对而言,在由光信号生成装置200生成的输入光信号的功率谱是未知的情况下,输入频谱取得部141取得利用未图示的分光器等测量出的输入光信号的功率谱。功率谱是表示光信号的每个波长的光的强度的数据。
另外,表示输入光信号的功率谱的信息,并非定要表示输入光信号的功率谱本身。例如,表示输入光信号的功率谱的信息,除了输入光信号的功率谱之外,还可以表示输入光信号的自相关函数。即,表示输入光信号的功率谱的信息只要是能够获得输入光信号的功率谱的信息,可以是任何形式的信息。
输出频谱取得部142按多个强度的输入光信号,分别取得输入光信号在高非线性光纤120内传播之后被输出的输出光信号的、由分光器130测量出的功率谱,并以此作为测量功率谱。在此,测量功率谱表示输出光信号的每个波长的光的强度。
相位频谱算出部143,按输入光信号的每个强度,设想为输入光信号具有给定的相位频谱,来进行输入光信号在高非线性光纤120内的传播模拟,从而算出输出光信号的计算功率谱。并且,相位频谱算出部143,以使计算功率谱和测量功率谱的按每个强度的差异变小的方式,来使给定的相位频谱变化,并进行传播模拟,从而探索输入光信号的相位频谱。例如,相位频谱算出部143按照规定的演算法,使给定的相位频谱变化。
在此,规定的演算法是用于探索所给出的函数的最优解或者近似解的演算法。例如,规定的演算法可以是退火法、共轭方向法、共轭梯度法、遗传演算法等。另外,规定的演算法并不限定于这些演算法,可以是任何演算法。
另外,相位频谱算出部143所探索的相位频谱是表示光信号的每个波长的相位的数据。
另外,由相位频谱算出部143执行的传播模拟,例如是通过分步傅立叶法等进行的脉冲传播模拟。即,由相位频谱算出部143执行的传播模拟,是为了利用高非线性光纤120所固有的已知的参数,计算输出光信号的功率谱的光信号传播模拟。在此使用的参数包含与光学非线性效果相关的参数。在本实施方式中,参数也包含与色散效果相关的参数。
另外,给定的相位频谱是传播模拟中给出的相位频谱,是输入光信号的虚拟相位频谱。即,给定的相位频谱是传播模拟中使用的相位频谱。在获得最高评价的传播模拟中使用的给定的相位频谱,被决定为输入光信号的相位频谱。
波形重构部144通过对由相位频谱算出部143探索出的相位频谱和由输入频谱取得部141取得的信息所表示的功率谱进行频率/时间变换,来重构输入光信号的时间波形。具体是,波形重构部144例如对相位频谱和功率谱进行逆傅里叶变换,来重构时间波形。
<波形重构装置的处理动作>
以下,关于具有如上所述的结构的波形重构装置140的各种动作进行说明。
图3是表示实施方式1的波形重构装置140的动作的流程图。
首先,输入频谱取得部141取得表示由光信号生成装置200生成的输入光信号的功率谱的信息(步骤S101)。接下来,输出频谱取得部142取得由分光器130测量出的、与多个强度的输入光信号分别对应的输出光信号的功率谱,并以此作为测量功率谱(步骤S102)。
然后,相位频谱算出部143设定用于传播模拟的给定的相位频谱的初始值(步骤S103)。例如,相位频谱算出部143将任意的相位频谱设定为初始值。另外,例如,相位频谱算出部143也可以将在设想为输入光信号是预先决定的类型的脉冲时所获得的相位频谱,设定为初始值。
并且,相位频谱算出部143,按输入光信号的每个强度,设想为输入光信号具有给定的相位频谱,并进行输入光信号在高非线性光纤120内的传播模拟,从而算出输出光信号的计算功率谱(步骤S104)。关于传播模拟,详情后述。
然后,相位频谱算出部143判定是否结束相位频谱的探索(步骤S105)。例如,相位频谱算出部143根据表示所计算出的多个计算功率谱和由分光器130测量出的(在步骤S102取得的)测量功率谱的差异的值是否在规定阈值以下,来判定是否结束探索。再例如,相位频谱算出部143可以根据表示计算功率谱和测量功率谱的差异的值的变化率,来判定是否结束探索。再例如,相位频谱算出部143可以根据传播模拟的反复次数是否达到了上限次数,来判定是否结束探索。
在此,如果判定为不结束探索(步骤S105为“否”),相位频谱算出部143就以使计算功率谱和测量功率谱的差异变小的方式,使给定的相位频谱变化(步骤S106)。具体是,相位频谱算出部143例如按照退火法来设定用于传播模拟的新的给定的相位频谱。并且,再反复进行从步骤S104开始的处理。
另一方面,如果判定为结束探索(步骤S105为“是”),相位频谱算出部143就将获得了最小差值的传输模中使用的给定的相位频谱,决定为输入光信号的相位频谱(步骤S107)。
如上所述,相位频谱算出部143通过一边使给定的相位频谱变化一边反复进行传播模拟,来探索输入光信号的相位频谱,直到被判定为结束探索为止。即,相位频谱算出部143利用按每个强度对计算功率谱和测量功率谱的差异进行评价的评价函数,来探索相位频谱的最优解。
评价函数是指根据测量功率谱和计算功率谱的差异量而变化的函数。例如,作为评价函数,可以利用对于测量功率谱和计算功率谱的各频率成分值的差值多接近“0”的程度进行评价的函数。在此情况下,作为表示差异的值,例如可以采用各频率成分值的绝对差值之和或者差值平方之和。再例如,作为评价函数可以利用对于计算功率谱相对于测量功率谱的各频率成分值之比接近“1”的程度进行评价的函数。
最后,通过对输入光信号的功率谱和被决定的相位频谱进行频率/时间变换,来重构输入光信号的时间波形(步骤S108),并结束处理。
通过以上的处理,波形重构装置140能够利用输入光信号在具有非线性光学效果的光传输媒体内传播之后被输出的输出光信号的功率谱,来简单地重构输入光信号的时间波形。
<传播模拟>
以下,详细说明根据分步傅立叶法进行的传播模拟。首先,关于可对高非线性光纤120内传播的光信号造成影响的色散效果和非线性光学效果进行说明。
色散效果是指因光与物质的相互作用根据光的波长而异所造成的现象。由于色散效果,输入光信号的传播速度根据频率而变化。即,输入光信号中包含的各频率成分的相位发生偏离,而导致输入光信号的时间波形扩大。
另外,非线性光学效果是指因光(例如像超短光脉冲等强度非常强的光)和物质的相互作用而引起的非线性的多样现象。作为非线性光学效果,可举出自相位调制、自陡峭和拉曼应答等。
首先,关于自相位调制进行说明。光纤等光传输媒体的折射率,会与该传输媒体其中传播的光信号的强度成比例地稍有变化,因此会产生光信号自身的相位调制。将这样产生的相位调制称为自相位调制。
以下,关于自陡峭进行说明。自陡峭是指时间波形对称的输入光信号随着在光传输媒体中传播而其时间波形变得不对称,峰值移向后方的现象。引起自陡峭的原因在于群速度依赖于强度。
最后,关于拉曼应答进行说明。光射入物质的情况下,被散射成频率与入射光相等的强弹性散射(瑞利散射)光和,频率与入射光的频率稍有偏差的极弱的非弹性散射光。可将非弹性散射光分为,因物质中振荡的原子和离子而被散射的拉曼散射光,因物质中的音波而被散射的布里渊散射光。拉曼散射中,在入射光的强度超过阈值的情况下因受激发射而产生强的拉曼散射光的现象称之为受激拉曼散射。由于这个受激拉曼散射,在光传输媒体内,能量从光的高频率成分移向低频率成分,从而低频率成分得以强化。这种现象称之为拉曼应答。
因受到如上所述的色散效果和非线性光学效果的影响而在高非线性光纤120内传播的输入光信号的传播方程式,如式(1)所示。
在此,E表示光的电场成分,z表示高非线性光纤120内的距离。另外,D表示色散效果和损失,N表示非线性光学效果。可通过以下的式(2)以及式(3)来表示D以及N。
在此,α、β、γ表示高非线性光纤120的固有参数。具体而言,α表示与光强度损失相关的参数。另外,β表示与各阶色散效果相关的参数。另外,γ表示与非线性光学效果相关的参数。另外,T表示时间。
在式(2)的右侧,第1项表示2阶色散,第2项表示3阶色散,第3项表示4阶色散,最后的项表示在高非线性光纤120内传播的光强度的损失。另外,在式(3)的右侧,第1项表示自相位调制,第2项表示自陡峭,第3项表示拉曼应答。
在该传播方程式中,包含依赖于E本身的项。因此,难以同时计算色散项(D)和非线性项(N)。对此,在分步傅立叶法中,如图4所示,沿着光信号的传播方向,将高非线性光纤120分成虚拟的短步骤。并且,按每个步骤顺次计算色散项和非线性项,从而能够求出在高非线性光纤120中传播的光信号的近似解。
<实验结果>
以下,关于通过本实施方式的波形重构系统100获得的实验结果进行说明。
在本实验中,光信号生成装置200对于由MLLD以10GHz的周期输出的1.3ps的光脉冲,通过65m的SMF进行色散补偿之后,通过EDFA进行了增幅。光信号生成装置200将这样通过EDFA增幅之后的光脉冲,作为输入光信号而输出。
作为强度调节器110采用了可调光衰减器(VOA:Variable OpticalAttenuator)。强度调节器110通过以10.3mW为单位,使功率从17.2mW开始逐渐提高,从而将由光信号生成装置200生成的输入光信号的强度变更成10种不同的强度。并且,强度调节器110将强度变更后的输入光信号输出给高非线性光纤120。
除了2阶以及3阶色散之外,还考虑到4阶色散,作为高非线性光纤120采用了具有如表1所示参数的光纤。
(表1)
Loss:α | 1.0dB/km |
2nd order dispersion:β2 | +2.36×10-2ps2/km |
3rd order dispersion:β3 | +4.72×10-3ps3/km |
4th order dispersion:β4 | +2.90×10-4ps4/km |
Nonlinearily:γ | 15.0/W/km |
Fiber length:L | 92m |
根据由分光器130测量的、与输入光信号的10种尖强度对应的输出光信号的功率谱如图5的(a)~(j)所示。在图中,图表的横坐标表示波长(Wavelength(nm)),纵坐标强度(Intensity(a.u.))。另外,图5的(a)~(j)所示的功率谱分别表示强度调节器110的功率依次为17.2mW、27.5mW、37.8mW、48.1mW、58.4mW、68.7mW、79.0mW、89.4mW、99.7mW和110.0mW时的输出光信号的功率谱。
在本实验中,波形重构装置140所具备的相位频谱算出部143,将与7种类强度对应的(图5的(a)~(g)所示)输出光信号的功率谱的实测值和、通过进行模拟而获得的输出光信号的功率谱的计算值和绝对差值之和用作评价函数,采用退火法来使相位频谱变化,并探索评价函数的值成为最小的相位频谱。此时,作为相位频谱的初始值,采用根据由自相关求出的脉冲宽度,通过Gerchberg-Saxton法而使得与Sech型光脉冲近似时的相位频谱。通过将这种相位频谱利用为初始值,能够早期收敛评价函数的值。
图6是表示相位频谱的探索中的评价函数值的变化的实验结果的图表。在图6中,图表的横坐标表示评价函数的反复计算次数(Number ofiteration times),纵坐标表示评价函数的值(Evaluation function(a.u.))。
从图6可明显看出,评价函数的值快速收敛成了一定值。相位频谱算出部143将评价函数的值收敛为一定值时(例如,图6中的第2500次)的相位频谱决定为输入光信号的相位频谱。
图7是表示通过实验所获得的功率谱和相位频谱的图。在图7中,横坐标表示波长(Wavelength(nm))。另外,纵坐标的左侧表示强度(Intensity(a.u.)),右侧表示相位(Phase(rad))。在此,功率谱601与左侧纵坐标的强度对应,相位频谱602与右侧纵坐标的相位对应。
波形重构装置140所具备的波形重构部144,利用如图7中算出的相位频谱602和输入光信号的功率谱601,对输入光信号的时间波形进行重构。
图8是表示由波形重构装置140重构的输入光信号的时间波形的实验结果的图。在图8中,横坐标表示时间(Delay time(ps)),纵坐标表示强度(Intensity(a.u.))。
从图8可明显看出,被重构的输入光信号的时间波形,与作为一般的脉冲光被近似的Sech型或者Gaussian型的时间波形不同。即,可视为波形重构装置140重构成了比通过一般的近似而获得的时间波形更正确的时间波形。
然后,对上述实验中由相位频谱算出部143执行的模拟和算出的相位频谱的合理性进行验证。
图9是表示与7种类输入光信号的强度对应的输出光信号的功率谱的实测值和计算值的比较结果的图表。在此,功率谱的计算值是经过反复计算,评价函数的值收敛时获得的功率谱的值。另外,输入光信号的7种类强度是强度调节器110的功率为17.2mW、27.5mW、37.8mW、48.1mW、58.4mW、68.7mW以及79.0mW时的强度。
如图9的(a)~(g)所示,实测值801和计算值802大体一致。因此可知,在本实验中,正确地模拟出了对于输入光信号的时间波形反应灵敏的、通过自相位调制的功率谱的变化。
图10是表示与由分光器130测量出的10种类输入光信号的强度对应的输出光信号的功率谱中的、算出相位频谱时未被利用的3种类输入光信号的强度对应的输出光信号的功率谱的实测值和计算值的比较结果的图表。在此,功率谱的计算值是利用由相位频谱算出部143算出的输入光信号的相位频谱来进行模拟时的值。另外,算出相位频谱时未被利用的3种类输入光信号的强度是强度调节器110的功率为89.4mW、99.7mW以及110.0mW时的强度。
如10(a)~(c)所示,实测值901和计算值902大体一致。即,对于退火法中未被作为评价函数对象的输入光信号的强度,也利用由相位频谱算出部143算出的相位频谱,正确地计算了输出光信号的功率谱。由此可知,由相位频谱算出部143算出的输入光信号的相位频谱的精度高。即,利用由相位频谱算出部143算出的输入光信号的相位频谱而被重构的输入光信号的时间波形的精度也高。
如上所述,本实施方式的波形重构系统100以及波形重构装置140,能够利用输入光信号在激发自相位调制的光纤内传输之后被输出的输出光信号的功率谱,来重构输入光信号的时间波形。因此,波形重构系统100以及波形重构装置140并非一定需要时间分解频谱数据。即,波形重构系统100,为了重构光信号的时间波形,并非一定需要为了获得时间分解分光数据而必须的超高速时间门或者参照光源,因此能够简单地重构光信号的时间波形。
另外,波形重构系统100以及波形重构装置140,能够采用作为公知技术的退火法,通过在使相位频谱变化的同时反复进行传播模拟,来探索评价函数的值成为规定阈值以下的相位频谱,从而能在较短的时间内算出精度高的相位频谱。
(实施方式2)
以下,关于实施方式2进行说明。
在本实施方式中,与实施方式1的不同点在于,只根据非线性光学效果以及色散效果中的色散效果,来使给定的相位频谱变化。具体而言,波形重构装置,为了缩短相位频谱探索所需的探索时间,根据限制条件来使给定的相位频谱变化。以下,关于本实施方式的波形重构装置,以不同于实施方式1之处为重点进行说明。
另外,本实施方式的波形重构系统的结构、波形重构装置的功能结构以及处理动作,与实施方式1的图1~图3相同,在此省略图示。
在本实施方式中,相位频谱算出部143,只根据非线性光学效果以及色散效果中的色散效果,来使给定的相位频谱变化。具体是,相位频谱算出部143,根据限制条件,使给定的相位频谱变化。该限制条件是指,在设想为输入光信号是在只具有非线性光学效果以及色散效果中的色散效果的媒体内传输之后获得的信号的情况下应被满足的相位频谱的限制条件。
在此,关于限制条件进行详细说明。
设想为紧接在激光振荡之后的超短光脉冲的波形是无失真的傅立叶变换界限脉冲(TLP:Transform-Limited Pulse),并且,设想为作为重构对象的输入光信号是TLP在媒体内传播之后被输出的光信号。在此情况下,TLP按照式(1)在媒体内传播。
但是,传播TLP的媒体若是一般的媒体(即,并非具有高非线性光学性的媒体),非线性光学效果对于TLP的传播的影响就小。即,即使省略了式(1)中的N进行传播计算,计算结果也不会有大的的变换。因此,能够将输入光信号表示为TLP在只具有色散效果的媒体内传播之后被输出的光信号。在此,可通过泰勒展开,以式(4)表示媒体内的色散效果。
在式(4)中,0阶以及1阶的项只对时间波形的时间移动造成影响,而对相位并无影响。另外,在光信号在一般的光纤等媒体内的传播模拟中,考虑到4阶色散就足以满足需求。在此,能够以(5)所示的多项式函数,来表示因色散效果而失真的输入光信号的相位。
φ(ω)=α(ω-ω0)2+b(ω-ω0)3+c(ω-ω0)4 (5)
相位频谱算出部143以满足该式(5)的方式来设定给定的相位频谱。即,相位频谱算出部143按照式(5)所示的限制条件,使给定的相位频谱变化。
根据式(5),相位频谱由a、b、c的3个变量所表示。相对而言,在实施方式1中,相位频谱由每个频率成分的变量(例如512个变量)所表达。即,根据本变形例,在探索相位频谱时,比起实施方式1,能够减少变量的个数,从而能够缩小探索空间。因此,根据本变形例,能够缩短相位频谱的探索时间。
如上所述,根据本实施方式的波形重构装置140,能够根据适当的限制条件来使给定的相位频谱变化。因此,在探索相位频谱时,能够减少变量的个数,从而能够缩短探索时间。
另外,在此,作为限制条件,考虑了从2阶色散至4阶色散的情况,还可以考虑更高阶的色散。例如,在探索相位频谱时,可以使作为限制条件所考虑的色散的阶数逐步增加。从而,能够逐步扩大探索空间,从而能够缩短探索时间以及提高探索精度。
(实施方式2的变形例)
以下,关于实施方式2的变形例进行说明。
在实施方式2中,设想为输入光信号是在只具有色散效果的媒体内传播之后被输出的信号。因此,在输入光信号受到色散效果以外的影响的情况下,探索精度会降低。对此,在本变形例中,为了提高相位频谱的探索精度,将实施方式1以及实施方式2的两个方法组合起来,进行相位频谱的探索。
具体而言,相位频谱算出部143,在探索输入光信号的相位频谱时,按照限制条件来使给定的相位频谱变化并进行传播模拟之后,无论限制条件如何,都使给定的相位频谱变化并进行传播模拟。在此,限制条件与实施方式2的限制条件相同。
即,相位频谱算出部143,首先,像实施方式2所示的那样按照限制条件使给定的相位频谱变化,探索相位频谱。然后,相位频谱算出部143,将探索结果用作给定的相位频谱的初始值,像实施方式1所示的那样,无论限制条件如何都使相位频谱变化,探索相位频谱。
如上所述,根据本变形例的波形重构装置140,按照限制条件使给定的相位频谱变化之后,无论限制条件如何,都能使给定的相位频谱能变化。因此,在探索相位频谱时,能够抑制探索精度降低的同时,能够缩短探索时间。
<实验结果>
以下,关于本变形例的波形重构装置140的实验结果进行说明。
在本实验中,作为光信号生成装置200,采用了可生成中心波长为1546nm、反复频率为25MHz的超短光脉冲的光纤激光。另外,作为高非线性光纤120,采用了10m长,且具有比实施方式1大的非线性光学效果的光纤(γ=21/W/km)。另外,将生成的超短光脉冲的强度减弱到不发生非线性光学效果的程度的光信号,将其用作输入光信号。
图11是表示实施方式2的变形例的实验中采用的输入光信号的功率谱的图表。图12的(a)以及(b)是表示2个强度的输入光信号在高非线性光纤120中传播之后被分光器130测量出的功率谱的图表。如图12所示,在本实验中,对输入光信号的强度2进行等级变化。
作为评价函数,使用了对由分光器130测量出的测量功率谱和利用给定的相位频谱进行传播模拟而获得的计算功率谱的按每个频率的绝对差值之和进行正规化的值。相位频谱算出部143,以使该评价函数的值变小的方式,使给定的相位频谱变化并反复进行传播模拟,从而探索出相位频谱。此时,作为给定的相位频谱的初始值,采用了在全频率的成分值为0的相位频谱。另外,探索时的传播模拟的反复次数为1000次。
关于在这样的条件下进行的实验结果,利用图13以及图14进行说明。
图13以及图14是表示实施方式2的变形例的波形重构装置的实验结果的图。
在图13中,横坐标表示波长,纵坐标表示强度(左侧)或者相位(右侧)。另外,虚线表示功率谱,与强度的纵坐标对应。另外,实线表示相位频谱,与相位的纵坐标对应。另外,在图14中,横坐标表示反复次数,纵坐标表示评价函数。
另外,在图13以及图14中,(a)表示通过实施方式1的探索方法获得的实验结果,(b)表示通过本变形例的探索方法获得的实验结果。
另外,实施方式1的探索方法中,作为给定的相位频谱的各频率成分值,在1000次的传播模拟中全设定了互相独立的值。另外,在本变形例的探索方法中,作为给定的相位频谱的各频率成分值,在1000次的起初的52次传播模拟中,设定了可满足式(5)的各频率成分值。另外,在剩余的948次传播模拟中,作为给定的相位频谱的各频率成分值,设定了与式(5)无关的、互相独立的值。
从图13以及图14可明确看出,根据本变形例的探索方法,与实施方式1的探索方法相比,更容易收敛解,且提高了探索速度。作为理由可考虑到,由于限制条件,探索空间得以适当缩小。
另外,在本实验中,根据实施方式1的探索方法通过1000次的传播模拟所获得的评价函数的值,与根据变形例的方法通过293次的传播模拟所获得的评价函数的值相同。即,通过本变形例的探索方法,能够实现约为实施方式1的探索方法的3.4倍的高速化,从而能够缩短探索时间。
另外,如本实验所示,在给定的相位频谱的初始值与最优解相差较大的情况下,本变形例的探索方法比实施方式1的探索方法,更能够提高探索精度。作为理由可考虑是,根据限制条件能够适当地缩小探索空间,因此能够抑制收敛为局部的最优解。
如上所述,通过本实验确认到,通过按照限制条件来使给定的相位频谱变化,能够提高探索精度,并能缩短探索时间。
(实施方式3)
以下,关于实施方式3进行说明。
在本实施方式中,与实施方式1主要不同的是,只根据非线性光学效果以及色散效果中的非线性光学效果来进行传播模拟。以下,关于本实施方式的波形重构装置,以不同于实施方式1的点为中心进行说明。
在此,本实施方式的波形重构系统的结构、波形重构装置的功能结构以及处理动作,与实施方式1的图1~图3相同,因此省略图示。
在本实施方式中,相位频谱算出部143,只根据非线性光学效果以及色散效果中的非线性光学效果,来进行传播模拟。具体是,相位频谱算出部143,只利用与非线性光学效果和色散效果中的非线性光学效果相关的参数来进行传播模拟。具体而言,相位频谱算出部143,例如,省略表示式(2)中的2阶色散、3阶色散以及4阶色散的项(第1项~第3项),来进行传播模拟。即,相位频谱算出部143,根据以下的式(6)~式(8),来进行传播模拟。
另外,在非线性光学效果大的光传输媒体内的光信号传播中,非线性光学效果具有支配性,而色散效果是限定性的。即,色散效果对输入光信号在高非线性光纤120内的传播模拟造成的影响小。因此,即使只利用与非线性光学效果相关的参数来进行输入光信号在高非线性光纤120内的传播模拟,也能够抑制模拟的精度降低。
另外,通过像上述那样省略表示色散效果的项,能够在式(6)中同时计算D和N。即,相位频谱算出部143无需采用分步傅立叶法,而是通过1次的计算就能够进行传播模拟。因此,能够大幅度提高传播模拟的计算速度。
另外,高非线性光纤120具有色散效果为“0”的波长(以下,称之为“零色散波长”)。输入光信号的频率成分越是集中于零色散波长,色散效果对光信号传播的影响就越小。因此,输入光信号的中心波长和高非线性光纤120的零色散波长的差异越小,即使省略色散效果,也能够高精度地进行传播模拟。
因此,在输入光信号的中心波长是已知的情况下,优先利用具有与该中心波长近似的值的零色散波长的高非线性光纤120。由此,在只利用与非线性光学效果和色散效果中的非线性光学效果相关的参数来进行传播模拟时,能够抑制模拟结果的精度降低。
另外,实验证实,在高非线性光纤120的零色散波长为1550nm的情况下,如果输入光信号的中心波长是1520nm~1580nm,就能够有效地抑制模拟的精度降低。
如上所述,根据本实施方式的波形重构装置140,能够只利用与非线性光学效果和色散效果中的非线性光学效果相关的参数来进行传播模拟。因此,传播模拟被简化,从而能够提高计算速度。
(实施方式3的变形例)
以下,关于实施方式3的变形例进行说明。
上述实施方式3中,在考虑自相位调制、自陡峭以及拉曼应答的基础上进行了传播模拟。在表示自陡峭以及拉曼应答的项(式(8)的右侧的第2项以及第3项)中,包含时间波形的微分。时间波形的微分相当于某时间的时间波形的倾斜度。即,脉冲宽度越窄,自陡峭以及拉曼应答对传播造成的影响就越大。相对而言,在脉冲宽度宽的情况下,可认为即使省略自陡峭以及拉曼应答的项来进行传播模拟,模拟结果也不会有很大的变化。
在此,在本变形例,相位频谱算出部143,只利用与非线性光学效果相关的参数中的与自相位调制相关的参数,来进行传播模拟。具体是,相位频谱算出部143,例如利用省略了式(3)中表示自陡峭以及拉曼应答的项(第2项以及第3项)的式来进行传播模拟。即,相位频谱算出部143根据以下的式(9)~式(11),进行传播模拟。
N=iγ|E|2 (11)
如上所述,根据本变形例的波形重构装置140,比上述实施方式3更能够简化传播模拟,从而更能够提高传播模拟的计算速度。尤其是在输入光信号的脉冲宽度宽的情况下,能够抑制模拟精度降低,更显效果。另外,根据实验,如果脉冲宽度是500fs以上,能够更有效地抑制模拟的精度降低。
<实验结果>
以下,关于通过本变形例的波形重构装置140获得的实验结果进行说明。在此,采用与实施方式2的变形例的实验相同的条件进行了实验。另外,在本实验中,高非线性光纤120的零色散波长是1550nm。另外,作为探索方法,采用了实施方式2的变形例的探索方法。即,起初按照限制条件使给定的相位频谱变化。
图15是表示实施方式3的变形例的波形重构装置的实验结果的图。
在图15中,横坐标表示波长,纵坐标表示强度(左侧)以及相位(右侧)。另外,虚线表示功率谱,与强度的纵坐标对应。另外,实线表示相位频谱,与相位的纵坐标对应。另外,图15的(a)表示基于分步傅立叶法的传播模拟的实验结果。另外,图15的(b)表示本变形例的传播模拟的实验结果。
在基于分步傅立叶法的传播模拟中,根据式(1)~式(3)进行了计算。另外,在本变形例的传播模拟中,根据式(9)~式(11)进行了计算。
如图15所示,反复进行传播模拟而探索出的相位频谱,在2个传播模拟中大体相同。即,进行传播模拟时即使忽略色散效果等,探索结果也不会有不太的变化。
另外,本变形例的传播模拟比基于分步傅立叶法的传播模拟,能够更大幅度地削减计算时间。具体是,就1000次的传播模拟所需要的时间而言,基于分步傅立叶法的传播模拟需要9970ms,本变形例的传播模拟仅需要567ms。
即,通过本实验确认到,通过只考虑非线性光学效果中的自相位调制来进行传播模拟,能够抑制精度降低,并能够缩短计算时间。
以上,关于本发明的一形态的波形重构装置以及波形重构系统,基于实施方式进行了说明,但本发明并不限定于这些实施方式。只要不超出本发明的宗旨,将本领域技术人员想出的各种变形方式实施于本实施方式而成的形态,或者对不同的实施方式的结构要素进行组合而成的形态,也属于本发明的范围中。
例如,可对上述实施方式2及其变形例、上述实施方式3将其变形例进行组合。即,波形重构装置可以按照限制条件使给定的相位频谱变化,并且只利用非线性光学效果以及色散效果中的非线性光学效果的参数来反复进行传播模拟。从而,能够进一步缩短相位频谱的探索时间。
另外,在上述实施方式1~3中,波形重构系统具备高非线性光纤,此外也可以具备与光纤不同的光传输媒体。即,波形重构系统只要具备具有非线性光学效果的光传输媒体即可。
另外,能够通过用于执行这种波形重构装置的特征性结构要素所进行的处理的波形重构方法来实现本发明。另外,还能够作为用于由图16所示的计算机执行该波形重构方法的程序来实现本发明。另外,能够通过CD-ROM等记录介质、互联网等传输媒体来传送上述程序。
图16是表示计算机的硬件结构的一个例子的图。用于由计算机执行波形重构方法的程序,例如,被存储在作为计算机可读取的媒体的CD-ROM515中,通过CD-ROM装置514被读取。用于由计算机执行波形重构方法的程序,例如通过有线或者无线网络、或者广播等被传送。
计算机500具备CPU(Central Processing Unit)50、ROM(ReadOnly Memory)502、RAM(Random Access Memory)503、硬盘504和通信接口505等。
CPU501执行由CD-ROM装置514读取的程序或者通过通信接口505接收的程序。具体是,CPU501在RAM503展开通过CD-ROM装置514读取的程序或者通过通信接口505接收的程序。并且,CPU501执行在RAM503展开的程序中的被编码化的各命令。
ROM502是存储计算机500的动作所需要的程序以及数据的读取专用存储器。RAM503被用为CPU501执行程序时的工作区域。具体是,RAM503临时存储例如对执行程序时的参数等数据。硬盘504存储数据等。
通信接口505通过网络与其他计算机进行通信。总线506使CPU501、ROM502、RAM503、硬盘504、通信接口505、显示器511、键盘512、鼠标513以及CD-ROM装置514相互连接。
另外,上述各实施方式的各结构要素可由专用的硬件构成,也可以由适于执行各结构要素的软件程序实现。各结构要素可由CPU或者处理器等程序执行部读取并执行硬盘或者半导体存储器等记录介质中存储的软件程序来实现。在此,实现上述各实施方式的波形重构装置等的软件例如是如下程序。
即,该程序使计算机执行对输入光信号的时间波形进行重构的波形重构方法,该波形重构方法包括:取得表示上述输入光信号的功率谱的信息的步骤;对于多个强度的上述输入光信号,将通过对上述输入光信号在具有非线性光学效果的光传输媒体内传播之后被输出的输出光信号进行测量而获得的功率谱,分别作为测量功率谱而取得的步骤;按上述输入光信号的每个强度,设想为上述输入光信号具有给定的相位频谱,来进行上述输入光信号在上述光传输媒体内的传播模拟,从而算出上述输出光信号的计算功率谱的步骤;以使上述计算功率谱和上述测量功率谱的按上述输入光信号的每个强度的差异变小的方式,使上述给定的相位频谱变化,并进行上述传播模拟,从而探索上述输入光信号的相位频谱的步骤;对探索出的上述相位频谱和取得的上述信息所表示的功率谱,进行频率/时间变换,从而对上述输入光信号的时间波形进行重构的步骤,在探索上述输入光信号的相位频谱的步骤中,只根据非线性光学效果以及色散效果中的一方,来使上述给定的相位频谱变化,或者进行上述传播模拟。
工业实用性
本发明的一形态的波形重构装置以及波形重构系统,能够应用于光通信等研究领域中使用的超高速光脉冲时间波形测量装置。尤其是近年来,基于相位的光通信符号的利用,以及对光脉冲的状态敏感的非线性光信号处理被广泛研究或者利用,因此,能够简单且正确地测量时间波形的波形重构装置以及波形重构系统,不仅是在工业领域,在基础科学领域也极为有用。
符号说明
100 波形重构系统
110 强度调节器
120 高非线性光纤
130 分光器
140 波形重构装置
141 输入频谱取得部
142 输出频谱取得部
143 相位频谱算出部
144 波形重构部
200 光信号生成装置
Claims (10)
1.一种波形重构装置,对输入光信号的时间波形进行重构,该波形重构装置具备:
输入频谱取得部,取得表示上述输入光信号的功率谱的信息;
输出频谱取得部,对于多个强度的上述输入光信号,将通过对上述输入光信号在具有非线性光学效果的光传输媒体内传播之后被输出的输出光信号进行测量而获得的功率谱,分别作为测量功率谱而取得;
相位频谱算出部,(a)按上述输入光信号的每个强度,设想为上述输入光信号具有给定的相位频谱,来进行上述输入光信号在上述光传输媒体内的传播模拟,从而算出上述输出光信号的计算功率谱,(b)以使上述计算功率谱和上述测量功率谱的按上述输入光信号的每个强度的差异变小的方式,使上述给定的相位频谱变化,并进行上述传播模拟,从而探索上述输入光信号的相位频谱;以及
波形重构部,通过对由上述相位频谱算出部探索出的相位频谱和由上述输入频谱取得部取得的信息所表示的功率谱,进行频率/时间变换,从而对上述输入光信号的时间波形进行重构,
上述相位频谱算出部,只根据非线性光学效果以及色散效果中的一方,来使上述给定的相位频谱变化,或者进行上述传播模拟。
2.如利要求1所述的波形重构装置,
上述相位频谱算出部,按照在设想为上述输入光信号是在只具有非线性光学效果以及色散效果中的色散效果的媒体内传播而获得的信号的情况下应被满足的相位频谱的限制条件,使上述给定的相位频谱变化。
3.如权利要求2所述的波形重构装置,
上述相位频谱算出部,在探索上述输入光信号的相位频谱时,在按照上述限制条件使上述给定的相位频谱变化并进行上述传播模拟之后,无论上述限制条件如何,都使上述给定的相位频谱变化,来进行上述传播模拟。
4.如权利要求1至3的任一项所述的波形重构装置,
上述相位频谱算出部只利用与非线性光学效果和色散效果中的非线性光学效果相关的参数,来进行上述传播模拟。
5.如权利要求4所述的波形重构装置,
上述相位频谱算出部,只利用与上述非线性光学效果相关的参数中的与自相位调制相关的参数,来进行上述传播模拟。
6.如权利要求1至5的任一项所述的波形重构装置,
上述光传输媒体是光纤。
7.一种波形重构系统,对输入光信号的时间波形进行重构,该波形重构系统具备:
强度调节器,使上述输入光信号的强度变化成互不相同的多个强度;
具有非线性光学效果的光传输媒体,至少该光传输媒体的与非线性光学效果相关的参数是已知的;
频谱测量器,按上述输入光信号的每个强度,对由上述强度调节器进行了强度变更的上述输入光信号在上述光传输媒体内传播之后被输出的输出光信号的功率谱进行测量;以及
权利要求1所述的波形重构装置,
上述输出频谱取得部取得由上述频谱测量器测量出功率谱,以此作为测量功率谱。
8.如权利要求7所述的波形重构系统,
上述光传输媒体是光纤。
9.一种波形重构方法,对输入光信号的时间波形进行重构,该波形重构方法包括:
取得表示上述输入光信号的功率谱的信息的步骤;
对于多个强度的上述输入光信号,将通过对上述输入光信号在具有非线性光学效果的光传输媒体内传播之后被输出的输出光信号进行测量而获得的功率谱,分别作为测量功率谱而取得的步骤;
按上述输入光信号的每个强度,设想为上述输入光信号具有给定的相位频谱,来进行上述输入光信号在上述光传输媒体内的传播模拟,从而算出上述输出光信号的计算功率谱的步骤;
以使上述计算功率谱和上述测量功率谱的按上述输入光信号的每个强度的差异变小的方式,使上述给定的相位频谱变化,并进行上述传播模拟,从而探索上述输入光信号的相位频谱的步骤;以及
对探索出的上述相位频谱和取得的上述信息所表示的功率谱,进行频率/时间变换,从而对上述输入光信号的时间波形进行重构的步骤,
在探索上述输入光信号的相位频谱的步骤中,只根据非线性光学效果以及色散效果中的一方,来使上述给定的相位频谱变化,或者进行上述传播模拟。
10.一种程序,用于使计算机执行权利要求9所述的波形重构方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-034442 | 2012-02-20 | ||
JP2012034442A JP6025013B2 (ja) | 2012-02-20 | 2012-02-20 | 波形再構成装置、波形再構成システム及び波形再構成方法 |
PCT/JP2013/000647 WO2013125166A1 (ja) | 2012-02-20 | 2013-02-06 | 波形再構成装置、波形再構成システム及び波形再構成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104272072A true CN104272072A (zh) | 2015-01-07 |
CN104272072B CN104272072B (zh) | 2017-02-22 |
Family
ID=49005366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201380010306.9A Active CN104272072B (zh) | 2012-02-20 | 2013-02-06 | 波形重构装置、波形重构系统以及波形重构方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9337922B2 (zh) |
EP (1) | EP2818838B1 (zh) |
JP (1) | JP6025013B2 (zh) |
CN (1) | CN104272072B (zh) |
WO (1) | WO2013125166A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106768337A (zh) * | 2017-02-27 | 2017-05-31 | 中国科学院上海高等研究院 | 一种二维傅里叶变换电子光谱中的相位重构方法 |
CN107508654A (zh) * | 2017-07-04 | 2017-12-22 | 中国电子科技集团公司第五十研究所 | 一种分平面参数化数据链波形重构方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6025013B2 (ja) * | 2012-02-20 | 2016-11-16 | 国立大学法人大阪大学 | 波形再構成装置、波形再構成システム及び波形再構成方法 |
JP6290697B2 (ja) * | 2014-04-22 | 2018-03-07 | 浜松ホトニクス株式会社 | 波形計測装置およびパルス光生成装置 |
FR3106892B1 (fr) * | 2020-01-31 | 2022-02-25 | Oksenhendler Thomas | Dispositif et procédé de diagnostic de la compressibilité au foyer d’une impulsion ultra-brève à partir d’images spatio-spectrales de l’impulsion ayant subi un effet non linéaire dégénéré |
CN115208463B (zh) * | 2021-04-09 | 2024-01-26 | 西安电子科技大学 | 一种基于相位谱分析的光采样多目标测频方法 |
JP2023067514A (ja) * | 2021-11-01 | 2023-05-16 | 浜松ホトニクス株式会社 | データ作成装置、光制御装置、ターゲット強度スペクトログラム作成装置、データ作成方法、ターゲット強度スペクトログラム作成方法、データ作成プログラム、及びターゲット強度スペクトログラム作成プログラム |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070097825A1 (en) * | 2005-10-31 | 2007-05-03 | Plasmon Lms, Inc. | Waveform reconstructor for optical disk read channel |
WO2011028928A1 (en) * | 2009-09-02 | 2011-03-10 | Qualcomm Incorporated | Energy domain based peak reconstruction methods and apparatuses |
CN102141603A (zh) * | 2010-01-28 | 2011-08-03 | 西门子迈迪特(深圳)磁共振有限公司 | 一种平面回波成像方法和系统 |
US20110311223A1 (en) * | 2009-03-02 | 2011-12-22 | Tsuyoshi Konishi | Waveform reconstruction device, waveform reconstruction system, and waveform reconstruction method |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH118590A (ja) * | 1997-04-25 | 1999-01-12 | Oki Electric Ind Co Ltd | 光伝送システム及びその監視制御方法 |
US7567596B2 (en) * | 2001-01-30 | 2009-07-28 | Board Of Trustees Of Michigan State University | Control system and apparatus for use with ultra-fast laser |
JP2005315858A (ja) * | 2004-03-31 | 2005-11-10 | Sun Tec Kk | 光パルス評価装置およびインサービス光パルス評価装置 |
JP3920297B2 (ja) * | 2004-09-01 | 2007-05-30 | 富士通株式会社 | 光スイッチおよび光スイッチを利用した光波形モニタ装置 |
EP2682813B1 (en) * | 2011-03-04 | 2020-11-18 | Furukawa Electric Co., Ltd. | Optical amplifier, optical amplifying system, wavelength converter, optical amplification method, and optical communication system |
JP6025013B2 (ja) * | 2012-02-20 | 2016-11-16 | 国立大学法人大阪大学 | 波形再構成装置、波形再構成システム及び波形再構成方法 |
-
2012
- 2012-02-20 JP JP2012034442A patent/JP6025013B2/ja active Active
-
2013
- 2013-02-06 WO PCT/JP2013/000647 patent/WO2013125166A1/ja active Application Filing
- 2013-02-06 CN CN201380010306.9A patent/CN104272072B/zh active Active
- 2013-02-06 US US14/379,806 patent/US9337922B2/en active Active
- 2013-02-06 EP EP13751849.4A patent/EP2818838B1/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070097825A1 (en) * | 2005-10-31 | 2007-05-03 | Plasmon Lms, Inc. | Waveform reconstructor for optical disk read channel |
US20110311223A1 (en) * | 2009-03-02 | 2011-12-22 | Tsuyoshi Konishi | Waveform reconstruction device, waveform reconstruction system, and waveform reconstruction method |
WO2011028928A1 (en) * | 2009-09-02 | 2011-03-10 | Qualcomm Incorporated | Energy domain based peak reconstruction methods and apparatuses |
CN102141603A (zh) * | 2010-01-28 | 2011-08-03 | 西门子迈迪特(深圳)磁共振有限公司 | 一种平面回波成像方法和系统 |
Non-Patent Citations (1)
Title |
---|
A.SULLIVAN ET AL.: "Quantitative investigation of optical phase-measuring techniques for ultrashort pulse lasers", 《J.OPT.SOC.AM.B》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106768337A (zh) * | 2017-02-27 | 2017-05-31 | 中国科学院上海高等研究院 | 一种二维傅里叶变换电子光谱中的相位重构方法 |
CN106768337B (zh) * | 2017-02-27 | 2018-06-12 | 中国科学院上海高等研究院 | 一种二维傅里叶变换电子光谱中的相位重构方法 |
CN107508654A (zh) * | 2017-07-04 | 2017-12-22 | 中国电子科技集团公司第五十研究所 | 一种分平面参数化数据链波形重构方法 |
CN107508654B (zh) * | 2017-07-04 | 2020-04-28 | 中国电子科技集团公司第五十研究所 | 一种分平面参数化数据链波形重构方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2818838A4 (en) | 2015-11-11 |
WO2013125166A1 (ja) | 2013-08-29 |
EP2818838B1 (en) | 2020-04-08 |
JP2013170905A (ja) | 2013-09-02 |
US9337922B2 (en) | 2016-05-10 |
CN104272072B (zh) | 2017-02-22 |
EP2818838A1 (en) | 2014-12-31 |
US20150010303A1 (en) | 2015-01-08 |
JP6025013B2 (ja) | 2016-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104272072A (zh) | 波形重构装置、波形重构系统以及波形重构方法 | |
Churkin et al. | Numerical modelling of spectral, temporal and statistical properties of Raman fiber lasers | |
Dudley et al. | Extreme events in optics: Challenges of the MANUREVA project | |
EP3228025B1 (en) | Method for characterizing performance of a multimode fiber optical link and corresponding method for fabricating multimode optical fiber links | |
US8886037B2 (en) | Waveform reconstruction device, waveform reconstruction system, and waveform reconstruction method | |
CN105308429B (zh) | 光物理常数测量方法及光物理常数推测装置 | |
Tay et al. | Optical soliton simulation in optical fibers by OptiSystem | |
Feng et al. | Analysis of spontaneous Raman and Rayleigh scatterings in distributed Fiber Raman amplification systems based on a random distribution model | |
Liu et al. | Predicting the evolution of the mid-infrared supercontinuum in the near-zero dispersion region under weak triggering using LSTM neural networks | |
Kaliteevskiy et al. | Two-mode coupling model in a few mode fiber | |
Wang et al. | Research on fast and accurate extraction method of Brillouin frequency shift based on window-weighted centroid algorithm | |
Yao et al. | Simulation of atmospheric turbulence phase screen based on wind speed | |
Dorrer | Investigation of the spectrogram technique for the characterization of picosecond optical pulses | |
Zeytunyan et al. | Similariton-based spectral interferometry for femtosecond signal characterization | |
de Negreiros Júnior et al. | Ultrashort pulses propagation through different approaches of the Split-Step Fourier method | |
Miyake et al. | Brillouin optical correlation domain analysis based on intensity modulation | |
González-Herráez et al. | Simultaneous position-resolved mapping of chromatic dispersion and Brillouin shift along single-mode optical fibers | |
Korai et al. | Measurement of optical pulsewidth in the picosecond regime using a non-linear fiber and power meter | |
Iqbal et al. | Study & simulation of fiber nonlinearities using adaptive step size method | |
Nguyen et al. | Numerical method for simultaneous measurement of dispersion and nonlinear coefficient in optical fibers | |
Xiao et al. | Simulation of laser beam propagation in deep turbulence: Implied inner scale and scintillation | |
Kito et al. | Simplified and fast acquirable BOTDA with frequency-swept probe pulse | |
Huynh et al. | The interplay of thermal and pump fluctuations in stimulated Brillouin scattering | |
Yap et al. | Analysis of Second Order Dispersion on Free Space Optical Propagation | |
Yavtushenko et al. | Chirped pulse shaping via fiber dispersion modulation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |