Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH052928B2 - - Google Patents
[go: Go Back, main page]

JPH052928B2 - - Google Patents

Info

Publication number
JPH052928B2
JPH052928B2 JP7354787A JP7354787A JPH052928B2 JP H052928 B2 JPH052928 B2 JP H052928B2 JP 7354787 A JP7354787 A JP 7354787A JP 7354787 A JP7354787 A JP 7354787A JP H052928 B2 JPH052928 B2 JP H052928B2
Authority
JP
Japan
Prior art keywords
pulse
light
measured
optical
intensity
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.)
Expired - Lifetime
Application number
JP7354787A
Other languages
Japanese (ja)
Other versions
JPS63238525A (en
Inventor
Kazunori Naganuma
Juichi Noda
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NTT Inc
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP7354787A priority Critical patent/JPS63238525A/en
Priority to US07/093,653 priority patent/US4792230A/en
Publication of JPS63238525A publication Critical patent/JPS63238525A/en
Publication of JPH052928B2 publication Critical patent/JPH052928B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J11/00Measuring the characteristics of individual optical pulses or of optical pulse trains

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は短光パルス評価方法、さらに詳しくは
現存の光検出器の応答時間程度乃至それ以下の時
間幅をもつて高速に変化するが故に、これら光検
出器をもつてしては、このパルス波形およびパル
ス各部における光の波長の変化、又はそれと同値
の光の周波数の変化が詳細に測定評価できないよ
うな短光パルスについて、その光パルスが一定周
期で繰り返すパルス列をなしている場合に、パル
スの強度波形および上記瞬間周波数の変化の積分
に相当するパルスの位相波形の両者を詳細に測定
評価できる方法に関するものである。
[Detailed Description of the Invention] [Field of Industrial Application] The present invention relates to a short optical pulse evaluation method, and more specifically, to a short optical pulse evaluation method, which changes rapidly with a time width comparable to or less than the response time of existing photodetectors. With these photodetectors, it is possible to detect short optical pulses in which the pulse waveform and changes in the wavelength of light in each part of the pulse, or changes in the frequency of light equivalent to the pulse waveform, cannot be measured and evaluated in detail. The present invention relates to a method that allows detailed measurement and evaluation of both the pulse intensity waveform and the pulse phase waveform corresponding to the integral of the change in instantaneous frequency when the pulse train forms a pulse train that repeats at a constant period.

〔従来の技術およびその問題点〕[Conventional technology and its problems]

第2図は短光パルス測定評価装置の第1の従来
例を示す概略図である。被測定光パルス光束は、
半透鏡24によつて反射された一部の光束がスペ
クトル測定系に振り分けられ、残部が強度相関波
形測定系に入射供給される。スペクトル測定系は
設定波長を可変させる波長送り装置(走査機構)
27が装着された高分解能分光器26、光束を分
光器26の入射スリツトに結焦させる絞り込みレ
ンズ25および光検出器28、検出器28の出力
電圧を増幅する積分増幅器29によつて構成され
る。
FIG. 2 is a schematic diagram showing a first conventional example of a short optical pulse measurement and evaluation device. The measured optical pulse flux is
A part of the light beam reflected by the semi-transparent mirror 24 is distributed to the spectrum measurement system, and the remaining part is supplied to the intensity correlation waveform measurement system. The spectrum measurement system is a wavelength sending device (scanning mechanism) that changes the set wavelength.
The spectrometer 27 is equipped with a high-resolution spectrometer 26, a diaphragm lens 25 that focuses a light beam onto an input slit of the spectrometer 26, a photodetector 28, and an integrating amplifier 29 that amplifies the output voltage of the detector 28.

強度相関波形測定系においては、入射光束は半
透鏡1により2光束に分枝され、その1光束は固
定されたプリズム2によつて反射され、他の1光
束は移動台4によつて光軸方向に移動可能なプリ
ズム3によつて反射され、以上2つの反射光束は
半透鏡1上で合波され、再び1光束となる。この
半透鏡および2つのプリズム2,3からなる光学
系はマイケルソン干渉計なる名称を有し、入射前
の光電場を時間tの関数として、E(t)と書くと
き、再合波後の光電場は〔E(t)+E(t−τ)〕/
2と表される。ここに、遅延時間τはプリズム2
を通る光路とプリズム3を通る光路との間の相対
的光路長差lを光速度Cで除した量であり、プリ
ズム3が移動台4によつて可動であることから、
τを自由に変化させ得る。この再合波後、光束は
レンズ5によつて2次高調波光発生能を有する結
晶7中に焦点を結び、〔E(t)+E(t−τ)〕2に、
その電界が比例する2次高調波光(入射光の半分
の波長の光)の発生をみる。結晶7からの出射光
はレンズ6によつて光検出器9に絞り込まれる。
この際、光学フイルタ8によつて2次高調波光の
みが通過され、光検出器9に達する。光検出器9
の応答時間は、被測定光パルスひいては結晶7よ
り発生した2次高調波パルスの持続時間に比して
極めて長いので2次高調波光パルスの波形を忠実
に検出できず、光検出器9の出力電圧パルスの波
形はなまつたものとなるが、パルス波形の時間積
分値すなわち面積は2次高調波光パルスのエネル
ギに比例したものとなる。光検出器9の出力電圧
パルスは、大きな時定数Tを持つた積分増幅器1
0により前記出力電圧パルスの面積に比例した一
定の電圧に変換される。
In the intensity correlation waveform measurement system, an incident light beam is branched into two light beams by a semi-transparent mirror 1, one of which is reflected by a fixed prism 2, and the other beam is reflected by a movable table 4 on the optical axis. The light beam is reflected by a prism 3 that is movable in the direction, and the above two reflected light beams are combined on a semi-transparent mirror 1 to become one light beam again. This optical system consisting of a semi-transparent mirror and two prisms 2 and 3 is called a Michelson interferometer, and when the optical electric field before incidence is written as E(t) as a function of time t, the optical electric field after recombination is expressed as E(t). The optical electric field is [E(t)+E(t-τ)]/
It is expressed as 2. Here, the delay time τ is the prism 2
It is the amount obtained by dividing the relative optical path length difference l between the optical path passing through and the optical path passing through the prism 3 by the speed of light C, and since the prism 3 is movable by the movable stage 4,
τ can be freely changed. After this recombination, the light beam is focused by the lens 5 into the crystal 7 which has the ability to generate second harmonic light, and becomes [E(t)+E(t-τ)] 2 .
Look at the generation of second-order harmonic light (light with half the wavelength of the incident light) whose electric field is proportional. The light emitted from the crystal 7 is focused onto a photodetector 9 by a lens 6.
At this time, only the second harmonic light is passed through the optical filter 8 and reaches the photodetector 9. Photodetector 9
The response time of is extremely long compared to the duration of the optical pulse to be measured and the second harmonic pulse generated by the crystal 7, so the waveform of the second harmonic optical pulse cannot be faithfully detected, and the output of the photodetector 9 Although the waveform of the voltage pulse becomes blunt, the time integral value, that is, the area of the pulse waveform becomes proportional to the energy of the second harmonic optical pulse. The output voltage pulse of the photodetector 9 is passed through an integrating amplifier 1 having a large time constant T.
0 is converted into a constant voltage proportional to the area of the output voltage pulse.

この従来装置による短光パルス測定手順を以下
に示す。
The short optical pulse measurement procedure using this conventional device is shown below.

まず、信号切換スイツチ31を右方に倒し、
A/D変換器32をスペクトル測定系に接続す
る。この状態で波長送り装置27によつて分光器
26の中心波長を一定速度で走査しつつ、光検出
器28の出力電圧をA/D変換器32を介してコ
ンピユータ18に時系列的に読み込み記憶させ
る。こうして光パルスのスペクトルD(λ)のデ
ータ列が得られる。次に信号切換スイツチ31を
左方に倒し、A/D変換器を強度相関波形測定系
に接続する。この状態で移動台4を一定速度で動
かしつつ、積分増幅器10の出力電圧をA/D変
換器32を介してコンピユータ18に時系列的に
読み込み記憶させる。この際、積分増幅器10の
遮断周波数(時定数Tの逆数)は被測定光パルス
の中心波長λ0および移動台移動速度Vによつて定
まる周波数 f0=2V/λ0 よりも充分低く設定される。そのため、こうして
得られたデータ列は1+2G2(τ)に比例する。
ここでG2(τ)は強度相関波形であり、パルスの
強度波形I(t)を用いて、 G2(τ)=∫ I(t)-∞ I(t−τ)dt/∫ I(t)-∞ I(t)dt と表される。
First, turn the signal changeover switch 31 to the right,
The A/D converter 32 is connected to the spectrum measurement system. In this state, while scanning the center wavelength of the spectrometer 26 at a constant speed using the wavelength sending device 27, the output voltage of the photodetector 28 is read and stored in the computer 18 in time series via the A/D converter 32. let In this way, a data string of the optical pulse spectrum D(λ) is obtained. Next, the signal changeover switch 31 is turned to the left to connect the A/D converter to the intensity correlation waveform measurement system. In this state, while moving the movable table 4 at a constant speed, the output voltage of the integrating amplifier 10 is read and stored in the computer 18 in time series via the A/D converter 32. At this time, the cutoff frequency (reciprocal of time constant T) of the integrating amplifier 10 is set sufficiently lower than the frequency f 0 =2V/λ 0 determined by the center wavelength λ 0 of the optical pulse to be measured and the moving speed V of the moving table. Ru. Therefore, the data string thus obtained is proportional to 1+2G 2 (τ).
Here, G 2 (τ) is the intensity correlation waveform, and using the pulse intensity waveform I(t), G 2 (τ)=∫ I(t) -∞ I(t-τ)dt/∫ It is expressed as I(t) -∞ I(t)dt.

以上で測定段階を終え、解析段階に入る。最初
に波長λの関数として表わされたスペクトルD
(λ)を周波数ωの関数としてスペクトルD(ω)
に換算する。被測定パルスの光電場E(t)、そのフ
ーリエ変換をE〓(ω)とすると、 |E〓(ω)|=√() に従つて、D(ω)の1/2乗からE〓(ω)の絶対値
|E〓(ω)|が得られる。後は E〓(ω)=|E〓(ω)|exp〔iφ(ω)〕 (i:虚数単位) と書くときのパルス成分の位相φ(ω)が求まれ
ば、被測定光パルスの電場E(t)を完全に再機築で
きる。そこで、φ(ω)をスペクトルの中心波長
ω0の回りに展開し、 φ(ω)=φ0+(ω−ω0)p +(ω−ω02q+(ω−ω03r この際の展開係数、q、r(φ0は物理的意味を持
たない任意定数であり、またpはスペクトル中心
周波数ω0に対する補正を与えるのみで、パルス
の強度波形に影響を持たないので、この方法で求
めることのできるのはq、rのみである)を定め
るために、 Gtr(τ)=∫-∞lEtr(t)
l2lEtr(t−τ)|2dt/∫-∞lEtr(t)l4dt の2式に従つて計算したGtr(τ)が、実測の強
度相関波形に最も良い一致をみるように最小二乗
法を用いる。これによつて、3次までの精度でφ
(ω)が近似的に求められるので、上式に従つて、
Etr(t)を計算すれば、|Etr(t)|2として被測定パル
スの強度波形が、arg(Etr(t))として被測定パル
スの位相波形がそれぞれ求められることになる。
このパルスの位相波形を時間微分すればパルス上
各点での瞬間周波数又は瞬間波長が得られる。こ
のように解析段階でφ(ω)の3次の係数までし
か求められないことから、この方法では、原理的
に任意の形のパルス強度波形、位相波形を忠実詳
細には再現しえない。これは、φ(ω)の、より
高次迄の展開係数を求めることにより改善され得
るが、計算時間の増大のため困難である。
This completes the measurement stage and begins the analysis stage. Spectrum D initially expressed as a function of wavelength λ
Spectrum D(ω) with (λ) as a function of frequency ω
Convert to . Letting the optical electric field E(t) of the pulse to be measured and its Fourier transform be E〓(ω), E〓 is obtained from D(ω) to the 1/2 power according to |E〓(ω)|=√() The absolute value of (ω) |E〓(ω)| is obtained. After finding the phase φ(ω) of the pulse component when written as E〓(ω)=|E〓(ω)|exp〔iφ(ω)〕 (i: imaginary unit), we can calculate the phase of the optical pulse to be measured. The electric field E(t) can be completely reconstructed. Therefore, φ(ω) is expanded around the center wavelength ω 0 of the spectrum, and φ(ω) = φ 0 + (ω-ω 0 ) p + (ω-ω 0 ) 2 q + (ω-ω 0 ) 3 r In this case, the expansion coefficients, q, r (φ 0 are arbitrary constants that have no physical meaning, and p only provides correction for the spectral center frequency ω 0 and has no effect on the pulse intensity waveform. , only q and r can be determined using this method), G tr (τ) = ∫ / -∞ lEtr (t)
Gtr (τ) calculated according to the following two formulas: l 2 lEtr (t - τ) | 2 dt/∫ / -∞ lEtr(t)l 4 dt is the one that best matches the measured intensity correlation waveform. The least squares method is used. As a result, φ can be obtained with accuracy up to the third order.
(ω) can be found approximately, so according to the above formula,
When Etr(t) is calculated, the intensity waveform of the pulse to be measured is obtained as |Etr(t)| 2 , and the phase waveform of the pulse to be measured is obtained as arg(Etr(t)).
By time-differentiating the phase waveform of this pulse, the instantaneous frequency or instantaneous wavelength at each point on the pulse can be obtained. As described above, since only the third-order coefficients of φ(ω) can be obtained at the analysis stage, this method cannot, in principle, faithfully reproduce arbitrary shapes of pulse intensity waveforms and phase waveforms in detail. This can be improved by finding expansion coefficients of φ(ω) up to higher orders, but this is difficult due to the increase in calculation time.

更に、スペクトルの測定のために、極めて高い
分解能が要求される分光器を使用する点に多大の
問題を内包する。一例をあげれば、0.6μmを中心
波長とする幅1psの光パルスのスペクトルの半値
全幅は概ね5Åであり、この半値全幅内で互いに
独立とみなせるような測定点を50点とろうとする
と、既に分光器に要求される分解能は0.1Å以下
となり、被測定光パルスの幅がより長くなると、
それに比例して良い分解能が必要になる。
Furthermore, there are many problems involved in using a spectrometer that requires extremely high resolution to measure spectra. For example, the full width at half maximum of the spectrum of an optical pulse with a width of 1 ps with a center wavelength of 0.6 μm is approximately 5 Å. The resolution required for the instrument becomes 0.1 Å or less, and as the width of the optical pulse to be measured becomes longer,
Proportionally good resolution is required.

上述のような仕様を満足する分光器には焦点距
離1m以上の回折格子分光器ないし掃引形フアブ
リペロエタロンがあるが、前者は大容積かつ光軸
調整が容易でなく、後者は前者に増して慎重な光
軸調整が要求され、いずれも簡便な測定が不可能
であつた。
Spectrometers that satisfy the above specifications include a diffraction grating spectrometer with a focal length of 1 m or more or a swept Fabry-Perot etalon, but the former has a large volume and is difficult to adjust the optical axis, and the latter has a larger capacity than the former. Therefore, careful optical axis adjustment was required, making easy measurement impossible.

以上のように、この従来の短光パルス測定評価
法においては、光パルスの強度相関波形測定に用
いられるマイケルソン干渉計に加えて、別個にス
ペクトル測定のために高い分解能の分光器が必要
とされ、測定評価装置の体積、価格の大半が分光
器によつて占められ、小型化、低価格化、調整の
簡便化を阻んでいた。更に、付随する測定データ
解析方法が粗い近似にとどまつており、光パルス
の詳細な評価がなし得なかつた。
As described above, in this conventional short optical pulse measurement evaluation method, in addition to the Michelson interferometer used to measure the intensity correlation waveform of the optical pulse, a separate high-resolution spectrometer is required for spectrum measurement. However, most of the volume and cost of measurement and evaluation equipment is occupied by the spectrometer, which hinders miniaturization, lower prices, and easier adjustment. Furthermore, the accompanying measurement data analysis method remains at a rough approximation, making detailed evaluation of the optical pulse impossible.

一方、分光器を用いない短光パルス測定評価装
置として特願昭61−211100記載の発明があり、そ
の概略図を第2の従来例として第3図に示す。こ
の方法は、一定周期で繰り返す短光パルスを測定
評価する方法であつて、被測定光束を半透鏡にて
2光束に分枝し、相異なる径路を経由せしめるこ
とにより、相対的光路長差を付与した後、2光束
を再び合波し、合波後光束を2次高調波光発生能
を有する結晶中に入射結焦せしめ、発生する2次
高調波光強度を光検出器にて電圧値に変換する光
学系において、前記相対的光路長差を緩慢なる一
定速度にて変化させ、出力電圧値を時系列的に記
録し、採取したデータをフーリエ解析してえられ
る基本波の干渉フリンジの周波数間隔をもつて分
布する3つのスペクトルの内、零周波数近傍のス
ペクトルを1/2乗した後、フーリエ逆変換するこ
とにより、被測定波長光パルスの強度波形の再構
築を行い、しかるのち、基本波の干渉フリンジの
周波数の2倍と周波数近傍のスペクトルを1/2乗
したものと、既求の強度波形より反復計算によつ
て被測定光パルスの位相を再構築する。いま、2
次高調波光電場を時間tの関数として、u(t)
と書くと、上記、基本波の干渉フリンジの周波数
の2倍の周波数近傍のスペクトルの1/2乗は、u
(t)のフーリエ変換u〓(ω)の絶対値|u〓(ω)|で

り、また零周波数近傍のスペクトルの1/2乗は、
強度波形I(t)=|u(t)|のフーリエ変換I(ω)
の絶対値|I〓(ω)|である。若し、u(t)がある時
刻を境にして左右対称である、いわゆる対称パル
スであるならば上記の|u〓(ω)|、|I〓(ω)|は

ルスの強度波形I(t)及び位相波形(1/2)arg(u
(t))を求めるに充分の情報を与える。これを用い
て反復計算によつて再構築されたパルスの強度波
形、位相波形は第1の従来例の如き低い次数迄の
近似とは異なり、極めて詳細な点まで精確であ
る。ところが対称パルスでない場合、これだけで
は不充分であり、パルスを正しく再構築すること
ができない。故にこの方法は対称パルスに対して
は、強度波形、位相波形を詳細忠実に再構築でき
るが、もしパルスが非対称であると誤つた結果を
与えるという欠点があつた。
On the other hand, there is an invention described in Japanese Patent Application No. 61-211100 as a short optical pulse measurement and evaluation apparatus that does not use a spectrometer, and a schematic diagram thereof is shown in FIG. 3 as a second conventional example. This method is a method for measuring and evaluating short optical pulses that repeat at a constant period.The light beam to be measured is branched into two light beams using a semi-transparent mirror and sent through different paths, thereby reducing the relative optical path length difference. After the addition, the two light beams are combined again, the combined light beam is focused into a crystal that has the ability to generate second harmonic light, and the generated second harmonic light intensity is converted into a voltage value by a photodetector. In the optical system, the relative optical path length difference is changed at a slow constant speed, the output voltage value is recorded in time series, and the frequency interval of the interference fringe of the fundamental wave is obtained by Fourier analysis of the collected data. After raising the spectrum near zero frequency to the 1/2 power among the three distributed spectra, the intensity waveform of the measured wavelength optical pulse is reconstructed by inverse Fourier transform, and then the fundamental wave is The phase of the optical pulse to be measured is reconstructed through repeated calculations using twice the frequency of the interference fringe, the spectrum near the frequency raised to the 1/2 power, and the already determined intensity waveform. Now, 2
Let us define the harmonic optical electric field as a function of time t, u(t)
Then, the 1/2 power of the spectrum near the frequency twice the frequency of the interference fringe of the fundamental wave above is u
The absolute value of the Fourier transform u〓(ω) of (t) is |u〓(ω)|, and the 1/2 power of the spectrum near zero frequency is
Fourier transform I(ω) of intensity waveform I(t)=|u(t)|
The absolute value of |I〓(ω)| is. If u(t) is a so-called symmetrical pulse that is symmetrical with respect to a certain time, the above |u〓(ω)| and |I〓(ω)| are the pulse intensity waveforms I( t) and phase waveform (1/2) arg (u
(t)). The pulse intensity waveform and phase waveform reconstructed through repeated calculations using this method are accurate down to extremely detailed points, unlike the approximation down to low orders as in the first conventional example. However, if the pulse is not symmetrical, this alone is insufficient and the pulse cannot be reconstructed correctly. Therefore, although this method can reconstruct the intensity waveform and phase waveform in detail and faithfully for symmetrical pulses, it has the disadvantage that it gives erroneous results if the pulse is asymmetrical.

以上のように、この第2の従来々の短光パルス
測定評価法においては、精度及び小型・簡便性の
点で第1の従来例の測定評価法に大きく優るもの
の、パルスが非対称であるとき生ずる原理的誤差
のため、被測定パルスが対称パルスであることが
予め確定できる場合にのみその応用が限定されて
いた。
As described above, although this second conventional short optical pulse measurement and evaluation method is significantly superior to the first conventional measurement and evaluation method in terms of accuracy, compactness, and simplicity, when the pulse is asymmetric, Due to the principle error that occurs, its application has been limited to cases where it can be determined in advance that the pulse to be measured is a symmetrical pulse.

本発明は以上の点に鑑みてなされたものであ
り、測定評価装置を巨大、高価なものとする高分
解能分光器を必要とすることなく、かつパルスが
時間的に非対称である場合を含む任意の場合に光
パルスの詳細にわたる測定評価が可能な短光パル
ス測定評価方法を提供することを目的とする。
The present invention has been made in view of the above points, and does not require a high-resolution spectrometer that makes measurement and evaluation equipment large and expensive, and can be used in any manner, including when pulses are temporally asymmetric. It is an object of the present invention to provide a short optical pulse measurement and evaluation method that allows detailed measurement and evaluation of optical pulses in the case of.

〔問題点を解決するための手段〕[Means for solving problems]

上記目的のうち、測定評価装置を簡便・安価と
する為に本発明者は、従来技術の第2の従来例の
装置をして、非対称パルスの測定を不可能ならし
めている原因を種々検討した結果、該装置で測定
される2次高調波強度の相対的光路長差依存性に
加えて、被測定光パルスのスペクトルが知れれ
ば、非対称パルスを再構築するに充分なデータの
組が得られることを見出した。しかるに、そのス
ペクトルの測定の為に高分解能分光器を導入すれ
ば、装置は一挙に巨大、高価なものとなり、発明
の目的は達成されない。この問題を解決する為に
分光器を使用しない方法として本発明者は、前記
従来例に示した装置に新たに基本波光強度を測定
するための光検出器を設備した構成を発明した。
相対的光路長差を変化させつつ、この光検出器の
出力を時系列的に記録すれば、正に被測定光パル
スのスペクトルをフーリエ逆変換した量(インタ
ーフエログラムと呼ばれる。)を測定することに
他ならず、得られたデータをフーリエ変換すれば
該スペクトルが得られるのである。
Among the above purposes, in order to make the measurement and evaluation device simple and inexpensive, the present inventor used the device of the second conventional example of the prior art and investigated various causes that made it impossible to measure asymmetric pulses. As a result, if the spectrum of the optical pulse to be measured is known in addition to the dependence of the second harmonic intensity measured by the device on the relative optical path length difference, a sufficient set of data can be obtained to reconstruct the asymmetric pulse. I found out that it can be done. However, if a high-resolution spectrometer is introduced to measure the spectrum, the apparatus will become huge and expensive, and the purpose of the invention will not be achieved. In order to solve this problem, as a method that does not use a spectrometer, the present inventor invented a configuration in which the apparatus shown in the conventional example is newly equipped with a photodetector for measuring the fundamental wave light intensity.
If the output of this photodetector is recorded in time series while changing the relative optical path length difference, it is possible to measure the inverse Fourier transform of the spectrum of the optical pulse to be measured (called an interferogram). In particular, the spectrum can be obtained by Fourier transforming the obtained data.

〔作用〕[Effect]

本発明によれば、マイケルソン干渉計によつて
相対的光路長差(遅延時間)を付与された後、重
畳された光束を2次高調波発生能を有する結晶中
に入射結焦せしめることにより2次高調波光を発
生させ、この2次高調波光および元の被測定パル
ス光束と同じ波長の基本波光の強度を各々別個の
光検出器をもつて電圧値に変換し測定する。
According to the present invention, after a relative optical path length difference (delay time) is imparted by a Michelson interferometer, the superimposed light beam is focused into a crystal having a second harmonic generation ability. A second harmonic light is generated, and the intensity of the second harmonic light and the fundamental wave light having the same wavelength as the original pulsed light flux to be measured is converted into a voltage value and measured using separate photodetectors.

基本波光の強度を観測する光検出器の出力電圧
は1+Re〔G1(τ)exp(−iω0τ)〕で表わされる。
ここでは遅延時間、ω0=2πc/λ0は被測定光パル
ス中心角周波数であり、記号Reは複素数の実数
部を採ることを意味する。パルスの光電場E0(t)
を E0(t)=Re〔E(t)exp(−iω0t)〕 ……(1) と書くとき、G1(τ)は下式により表され、電場
相関関数と呼ばれる。
The output voltage of the photodetector that observes the intensity of the fundamental wave light is expressed as 1+Re[G 1 (τ) exp (−iω 0 τ)].
Here, the delay time ω 0 =2πc/λ 0 is the center angular frequency of the optical pulse to be measured, and the symbol Re means taking the real part of a complex number. Pulse optical electric field E 0 (t)
When written as E 0 (t)=Re[E(t)exp(-iω 0 t)]...(1), G 1 (τ) is expressed by the following formula and is called an electric field correlation function.

G1(τ)=-∞ E(t)E*(t−τ)dt ……(2) このデータG1(τ)をフーリエ変換して被測定
パルスの電場E(t)のフーリエ変換E〓(ω)の絶対
値の2乗が得られる。即ち、-∞ G1(τ)exp(iωτ)dτ =|E〓(ω)|2 ……(3) 上式の右辺は光パルスのスペクトルを表わしてい
る。かくの如く分光器でなく干渉計を用いること
によつて光のスペクトルを求める手法はフーリエ
変換分光法と呼ばれる。
G 1 (τ) = -∞ E(t)E * (t-τ) dt ...(2) This data G 1 (τ) is Fourier transformed to obtain the Fourier value of the electric field E(t) of the pulse to be measured. The square of the absolute value of the transformation E〓(ω) is obtained. That is, -∞ G 1 (τ)exp(iωτ)dτ = |E〓(ω)| 2 ...(3) The right side of the above equation represents the spectrum of the optical pulse. This method of determining the spectrum of light by using an interferometer rather than a spectrometer is called Fourier transform spectroscopy.

次に、2次高調波の強度を観測する光検出器9
の出力電圧は下式に比例する。
Next, a photodetector 9 is used to observe the intensity of the second harmonic.
The output voltage of is proportional to the following formula.

1+2G2(τ)+4Re〔F1(τ)exp(−iω0τ)〕+R
e
〔F2(τ)exp(−2iω0τ)〕 ……(4) 光パルスの強度I(t)及び2次高調波光電場u(t)
を下式により定義し導入する。
1 + 2G 2 (τ) + 4Re [F 1 (τ) exp (−iω 0 τ)] + R
e
[F 2 (τ) exp (−2iω 0 τ)] ...(4) Intensity of optical pulse I(t) and second harmonic optical electric field u(t)
is defined and introduced by the following formula.

I(t)=E(t)*(t) ……(5) u(t)=E(t)E(t) ……(6) これらのI、uを用いてG2(τ)(強度相関関
数)、F2(τ)(2次高調波電場相関関数)は次の
様に表される。
I(t)=E(t) * (t) ……(5) u(t)=E(t)E(t) ……(6) Using these I and u, G 2 (τ)( intensity correlation function) and F 2 (τ) (second harmonic electric field correlation function) are expressed as follows.

G2(τ)=-∞ I(t)I(t−τ)dt ……(7) F2(τ)=-∞ u(t)u*(t−τ)dt ……(8) これらをフーリエ変換すれば、各々光パルスの強
度I(t)のフーリエ変換I〓(ω)の絶対値の2乗、
2次高調波光電場u(t)のフーリエ変換u〓(ω)
の絶対値の2乗が得られる。即ち -∞ G2(τ)exp(iωτ)dτ =|I〓(ω)|2 ……(9) -∞ F2(τ)exp(iωτ)dτ =|u〓(ω)|2 ……(10) ここで式(4)で表わされる2次高調波光強度から
平担な直流成分を引いたものをフーリエ解析する
ことを考えた。第4図左上図に式(4)の信号を、ま
た左下図にそのフーリエ解析の結果を示す。信号
上、一見して明らかな振動の周期は、被測定光の
中心波長に一致する。これは、相対的光路長差を
光速度で除し、遅延時間に置き直して考えると、
被測定光の光振動の周期である。
G 2 (τ)= -∞ I(t) I(t-τ)dt ……(7) F 2 (τ)= -∞ u(t)u * (t-τ)dt …… (8) If these are Fourier transformed, the square of the absolute value of the Fourier transform I〓(ω) of the intensity I(t) of the optical pulse,
Fourier transform u〓(ω) of second harmonic optical electric field u(t)
The square of the absolute value of is obtained. That is, -∞ G 2 (τ)exp(iωτ)dτ = |I〓(ω)| 2 ......(9) -∞ F 2 (τ)exp(iωτ)dτ = |u〓(ω) | 2 ...(10) Here, we considered performing Fourier analysis on the second harmonic light intensity expressed by equation (4) minus the flat DC component. The upper left diagram of FIG. 4 shows the signal of equation (4), and the lower left diagram shows the results of its Fourier analysis. The period of vibration that is obvious at first glance on the signal matches the center wavelength of the light to be measured. This can be calculated by dividing the relative optical path length difference by the speed of light and replacing it with the delay time.
This is the period of optical oscillation of the light to be measured.

この信号をフーリエ解析すると、左下図に見ら
れる如く、3つのピークが出現する。このうち真
中のピークは、信号に一見して現れていた振動成
分によるもので、被測定光の光振動の周波数ω0
付近にある。このω0を基本波の干渉フリンジ周
波数と称する。フーリエ解析の結果、その他に零
周波数付近、及び2ω0付近にもピークが現れる。
これら3つのピークは各々式(4)中の3つの項に由
来する。以上のように、式(4)中でG2(τ)、F1
(τ)、F2(τ)はそれらに比して充分速く変化す
る搬送周波数ω0の差をもつて現れるので、式(4)
の信号を全体としてフーリエ解析すれば、容易に
それら3者を分離することができるのである。上
記、零周波数付近のピークは式(9)の|I〓(ω)|2
を、また2ω0付近のピークを抽出し−2ω0だけ平
行移動すれば、式(10)の|u〓(ω)|2が得られる。
When this signal is subjected to Fourier analysis, three peaks appear as shown in the lower left figure. The peak in the middle is due to the vibration component that appeared at first glance in the signal, and the frequency of the optical vibration of the light to be measured is ω 0
It's nearby. This ω 0 is called the interference fringe frequency of the fundamental wave. As a result of Fourier analysis, peaks also appear near zero frequency and near 2ω 0 .
These three peaks originate from the three terms in formula (4), respectively. As mentioned above, in equation (4), G 2 (τ), F 1
(τ) and F 2 (τ) appear with a difference in carrier frequency ω 0 that changes sufficiently faster than them, so Equation (4)
By performing Fourier analysis on the signal as a whole, these three types can be easily separated. Above, the peak near zero frequency is |I〓(ω)| 2 in equation (9)
If we also extract the peak near 2ω 0 and translate it by −2ω 0 , we can obtain |u〓(ω)| 2 in equation (10).

一方、基本波光強度の信号から平担な直流成分
を引いたものをフーリエ解析すれば、第4図右に
見られる如く、ω0に中心を持つ1個のピークが
得られる。これを−ω0だけ平行移動すれば、式
(3)の|E〓(ω)|2が得られる。
On the other hand, if the fundamental wave light intensity signal minus the flat DC component is subjected to Fourier analysis, one peak centered at ω 0 will be obtained, as seen on the right side of FIG. 4. If we translate this by −ω 0 , we get the equation
In (3), |E〓(ω)| 2 is obtained.

かくして基本波光強度、2次高調波光強度の相
対的光路長差(遅延時間)依存性のデータをフー
リエ解析して|I〓|2、|u〓|2、|E〓|2が得られる

この事情を第4図に示した。
Thus, by Fourier analysis of the data on the relative optical path length difference (delay time) dependence of the fundamental wave light intensity and the second harmonic light intensity, |I〓| 2 , |u〓| 2 , and |E〓| 2 are obtained.
This situation is shown in Figure 4.

次に、本発明では、かくして得られた|I〓|、|
u〓|、|E〓|から、被測定光パルスを再構築するた
めの、反復計算法を行なう。この計算法を第5図
に示した。図中、記号F.T.はフーリエ変換を、I.
F.T.はフーリエ逆変換を、またcsqrtは複素数の
平方根ととることを表す。パルスの再構築は原理
的に電場E(t)のフーリエ変換E〓(ω)の位相を求
めることに帰着される。何とならば、そのフーリ
エ変換の絶対値|E〓(ω)|は既に測定データから
得られているから、その位相が知れればE〓が得ら
れ、フーリエ逆変換によつてEが求められる。し
かして|E|2より被測定光パルスの強度波形が、
arg(E)より被測定パルスの位相波形がそれぞれ得
られる。本発明の反復計算法はE、I、uの間に
式(5)、(6)の関係が成立つていることを利用して|
E〓|、|I〓|、|u〓|の位相を構築してゆくもので

る。即ち、先づ|E〓(ω)|に無作為な乱数の位相
を付与しE〓(ω)を作る。それをフーリエ逆変換
して得られたEをもつて、第5図左下から繰り返
し計算に入る。このEより式(5)、式(6)に従つて
I、uを計算し、更にそれらを各々フーリエ変換
する。かくして得たI〓、u〓において、それらの絶対
値|I〓|、|u〓|を測定データのフーリエ解析より
前段迄に得られている|I〓|、|u〓|に置き換える。
ここで置換前後の|I〓|又は|u〓|の差は、計算の
収束の判定に用いる。それらの差が無くなつたと
き、完全に被測定パルスが再構築されたと考えら
れる。続いてそれらのI〓、u〓をフーリエ逆変換して
I、uに戻し、かつuの絶対値を得られたIで置
き換える。こうして得られたuの複素数としての
平方根をとりEに戻す。更にそのEをフーリエ変
換したE〓の絶対値を既知の|E〓|で置き換える。
置換前後の|E〓|の差を収束の判定に用いるの
は、I〓、u〓の場合と同様である。置換後のE〓をフー
リエ逆変換すれば出発点に戻つてEが得られる。
この計算ループを一巡して得るEは、元のEに比
して被測定光パルスの電場に近いものを与える。
その近似の程度は、計算ループ中で得られる3種
のフーリエ変換の絶対値の置換に際してのそれぞ
れの差に反映されるので、それらの差が充分小さ
くなるまで、計算を繰り返す。
Next, in the present invention, the thus obtained |I〓|, |
An iterative calculation method is performed to reconstruct the optical pulse to be measured from u〓|, |E〓|. This calculation method is shown in FIG. In the figure, the symbol FT represents the Fourier transform, and the symbol I.
FT represents the inverse Fourier transform, and csqrt represents the square root of a complex number. Reconstruction of the pulse basically boils down to finding the phase of the Fourier transform E〓(ω) of the electric field E(t). The reason is that the absolute value of the Fourier transform |E〓(ω)| has already been obtained from the measured data, so if we know its phase, we can obtain E〓, and E can be found by inverse Fourier transform. . Therefore, from |E| 2 , the intensity waveform of the optical pulse to be measured is
The phase waveform of the pulse to be measured can be obtained from arg(E). The iterative calculation method of the present invention utilizes the relationship between E, I, and u in equations (5) and (6).
It constructs the topology of E〓|, |I〓|, and |u〓|. That is, first, a phase of a random number is given to |E〓(ω)| to create E〓(ω). Using E obtained by inverse Fourier transform, calculations begin repeatedly from the bottom left of Figure 5. From this E, I and u are calculated according to equations (5) and (6), and then each of them is Fourier transformed. In I〓, u〓 obtained in this way, their absolute values |I〓|, |u〓| are replaced with |I〓|, |u〓| obtained up to the previous stage through Fourier analysis of the measurement data.
Here, the difference between |I〓| or |u〓| before and after the replacement is used to determine the convergence of the calculation. When the difference between them disappears, it is considered that the pulse to be measured has been completely reconstructed. Subsequently, those I〓, u〓 are inversely Fourier transformed to return I, u, and the absolute value of u is replaced by I obtained. The square root of u obtained in this way as a complex number is taken and returned to E. Furthermore, the absolute value of E〓 obtained by Fourier transforming E is replaced with the known |E〓|.
The difference between |E〓| before and after replacement is used to determine convergence, as in the case of I〓 and u〓. If E〓 after substitution is subjected to inverse Fourier transform, E can be obtained by returning to the starting point.
The E obtained by going through this calculation loop gives an electric field closer to the electric field of the optical pulse to be measured than the original E.
The degree of approximation is reflected in the differences between the absolute values of the three Fourier transforms obtained during the calculation loop, so calculations are repeated until these differences become sufficiently small.

以上の解析法は、第1の従来例における解析法
の如き有限次元までの近似解法とは異なり、原理
的には無限次まで取り入れており、高い精度で被
測定パルスの強度、位相が求まり、かつ第2の従
来例とは異なり、対称なパルスにその有効範囲が
制限されておらず、全く任意の強度、位相波形の
パルスに対して有効である。
The above analysis method is different from the approximate solution up to finite dimensions, such as the analysis method in the first conventional example, in principle, it incorporates up to infinite orders, and the intensity and phase of the pulse to be measured can be determined with high accuracy. Moreover, unlike the second conventional example, its effective range is not limited to symmetrical pulses, and is effective for pulses of completely arbitrary intensity and phase waveform.

〔実施例〕〔Example〕

第1図に本発明による短光パルス評価方法を実
施するための装置の構成を示す。半透鏡1、固定
プリズム2及び移動台4によつて可動なプリズム
3よりなるマイケルソン干渉計を構成している。
この干渉計の相対的光路長差の較正のためにヘリ
ウムネオンレーザ11からの光が用いられてい
る。この光と平行に被測定パルス光束がマイケル
ソン干渉計に入射する。このヘリウムネオンレー
ザ光は水平面に対し45゜をなす直線偏光を持ち、
干渉計の片方の腕の途中で軸が垂直方向を向く様
に置かれた4分の1波長板19を通過する。干渉
計から出射したヘリウムネオンレーザ光は、垂直
に偏光した光と、水平に偏光した光とを分離する
偏光ビームスプリツタ20に入射し、各々の偏光
の光が光検出器21,22に導かれ、両検出器は
それぞれの光強度に比例した電気信号を出力す
る。両出力はサンプルトリガ発生器23に入力す
る。サンプルトリガ発生器はそれぞれの入力の交
流成分のみを抽出し、該交流成分が零値をよぎつ
た時刻にトリガ電圧パルスを発生する。かくして
干渉計に相対的光路長差がヘリウムネオンレーザ
の発振波長0.63μmだけ変化する間に0.25×
0.63μm毎に計4回トリガ電圧パルスの発生をみ
る。このトリガ電圧パルスは波形記憶装置17に
供給され、波形記憶装置は後述の各チヤンネルへ
の入力電圧値をトリガ電圧パルスの起つた時刻に
サンプルして、記録する。
FIG. 1 shows the configuration of an apparatus for implementing the short optical pulse evaluation method according to the present invention. A Michelson interferometer is constituted by a semi-transparent mirror 1, a fixed prism 2, and a movable prism 3 by means of a movable stage 4.
Light from a helium neon laser 11 is used to calibrate the relative optical path length difference of this interferometer. A pulsed beam to be measured enters the Michelson interferometer in parallel with this light. This helium neon laser beam has linear polarization at an angle of 45° to the horizontal plane.
The light passes through a quarter-wave plate 19 placed in the middle of one arm of the interferometer so that its axis faces vertically. The helium neon laser beam emitted from the interferometer enters a polarization beam splitter 20 that separates vertically polarized light and horizontally polarized light, and each polarized light is guided to photodetectors 21 and 22. Both detectors output electrical signals proportional to their respective light intensities. Both outputs are input to sample trigger generator 23. The sample trigger generator extracts only the AC component of each input and generates a trigger voltage pulse at the time when the AC component crosses a zero value. Thus, the relative optical path length difference in the interferometer is 0.25× while the oscillation wavelength of the helium-neon laser changes by 0.63 μm.
Observe the generation of trigger voltage pulses a total of 4 times every 0.63 μm. This trigger voltage pulse is supplied to the waveform storage device 17, and the waveform storage device samples and records the input voltage value to each channel, which will be described later, at the time when the trigger voltage pulse occurs.

被測定パルス光束は干渉計中で前述の4分の1
波長板19を通過することなく干渉計を通過し、
レンズ5によつて2次高調波光発生能を有する結
晶7(例えばヨウ素酸リチウム結晶)中で焦点を
結び、2次高調波光の発生を見る。結晶7からの
出射光はレンズ6によつて焦束光束とされ半透鏡
24によつて2分される。その一方は被測定パル
ス光束の2倍の周波数(半分の波長)を有する2
次高調波光のみが光学フイルタ8によつて抽出さ
れ、光検出器9に集光する。他方は結晶7中で2
次高調波光に変換されることなく出射した元の被
測定パルス光束と同じ波長の基本波光のみが光学
フイルタ30によつて抽出され、光検出器28に
焦光する。光検出器9及び28は各々2次高調波
光エネルギー及び、基本波光エネルギーに比例し
た面積(時間積分値)を有する電圧パルスを圧力
し、それらの電圧パルスはそれぞれ積分増幅器1
0,29によつて、前記面積に比例した一定の電
圧に変換された後、波形記憶装置17のチヤンネ
ル1、チヤンネル2に入力される。波形記憶装置
17には、測定後、記憶されたデータを読み出し
解析を行うためのコンピユータ18が接続されて
いる。以上の説明により可動プリズム3の一つの
位置、換言すると一つの相対的光路長差における
一点の測定が行われる。実際には移動台4が一定
の移動速度Vで移動しており、干渉計の光路長差
が、基準とするヘリウムネオンレーザの波長λref
の1/4ずつ変化する毎にサンプリングを行い波形
記憶装置17で記憶する。
The pulsed light flux to be measured is reduced to 1/4 of the above amount in the interferometer.
passes through the interferometer without passing through the wave plate 19,
A lens 5 focuses the light in a crystal 7 (for example, a lithium iodate crystal) capable of generating second harmonic light, and the generation of second harmonic light is observed. The light emitted from the crystal 7 is made into a focused beam by the lens 6 and divided into two by the semi-transparent mirror 24. One of them is 2 which has twice the frequency (half the wavelength) of the pulsed light flux to be measured.
Only the harmonic light is extracted by the optical filter 8 and focused on the photodetector 9. The other one is 2 in 7 crystals.
Only the fundamental wave light having the same wavelength as the original pulsed light flux to be measured, which is emitted without being converted into harmonic light, is extracted by the optical filter 30 and focused on the photodetector 28. The photodetectors 9 and 28 each apply a voltage pulse having an area (time integral value) proportional to the second harmonic optical energy and the fundamental optical energy, and these voltage pulses are applied to the integrating amplifier 1, respectively.
0 and 29, the voltage is converted into a constant voltage proportional to the area, and then input to channels 1 and 2 of the waveform storage device 17. A computer 18 is connected to the waveform storage device 17 for reading and analyzing the stored data after measurement. According to the above explanation, measurement is performed at one position of the movable prism 3, in other words, at one point at one relative optical path length difference. In reality, the moving table 4 is moving at a constant moving speed V, and the difference in optical path length of the interferometer is the wavelength λ ref of the reference helium neon laser.
Sampling is performed every time the waveform changes by 1/4 and is stored in the waveform storage device 17.

この際、積分増幅器10及び29の時定数T、
被測定パルスの繰り返し間隔P、及び移動台移動
速度Vは P<<T<<λref/(8V) ……(11) の関係を満足するように設定する。また基準とな
るレーザーの波長λrefは、被測定パルスの中心波
形λ0よりも短くなければならない。これは干渉計
の相対的光路長差がλref/4変化する毎に1回サ
ンプルが行なわれる一方、波長λ0の基本波の半分
即ちλ0/2の波長を有する2次高調波を正しく観
測するためには、相対的光路長差がλ0/2変化す
るうちに、少くとも2回サンプルを行なわねばな
らないことによる。これにより発振波長0.63μm
のヘリウムネオンレーザーを干渉計の光路長差の
基準に用いる場合、中心波長が0.63μmより長い
光パルスの測定を行なえる。波長がより短い光パ
ルスの測定には、より発振波長の短いレーザーを
基準に用いればよい。
At this time, the time constant T of the integrating amplifiers 10 and 29,
The repetition interval P of the pulse to be measured and the moving speed V of the moving table are set so as to satisfy the following relationship: P<<T<<λ ref /(8V) (11). Further, the wavelength λ ref of the reference laser must be shorter than the center waveform λ 0 of the pulse to be measured. This means that one sample is taken every time the relative optical path length difference of the interferometer changes by λ ref /4, while the second harmonic with a wavelength of half of the fundamental wave of wavelength λ 0 , i.e. λ 0 /2, is correctly sampled. This is due to the fact that in order to observe it, samples must be taken at least twice while the relative optical path length difference changes by λ 0 /2. As a result, the oscillation wavelength is 0.63μm.
When using a helium-neon laser as a reference for the optical path length difference of an interferometer, it is possible to measure optical pulses with a center wavelength longer than 0.63 μm. To measure a light pulse with a shorter wavelength, a laser with a shorter oscillation wavelength may be used as a reference.

以下の解析段階では、まず波形記憶装置17に
記憶された積分増幅器10及び29の出力電圧変
化をコンピユータに読み込む。この2組のデータ
が第4図上段に示した相関データである。
In the following analysis step, first, the output voltage changes of the integrating amplifiers 10 and 29 stored in the waveform storage device 17 are read into the computer. These two sets of data are the correlation data shown in the upper part of FIG.

このデータを各々フーリエ解析して得られたス
ペクトルを第4図下段に示す。こうして得られる
|I〓(ω)|、|u〓(ω)|、|E〓(ω)|のデー
タを用い
て、第5図に示す反復計算を行なうことにより、
パルスの強度波形、位相波形が求められる。
The spectra obtained by Fourier analysis of this data are shown in the lower part of FIG. 4. By using the data of |I〓(ω)|, |u〓(ω)|, and |E〓(ω)| obtained in this way, by performing the iterative calculation shown in Fig. 5,
The intensity waveform and phase waveform of the pulse are determined.

以上の解析方法の有効性の検証として、種々の
強度波形、位相波形を有するパルスについて計算
機を用いて想定されるデータ|I〓(ω)|、|u〓(ω

|、|E〓(ω)|を生成しそれらから上記の解析方
法によつて強度波形、位相波形を再生した結果を
第6図から第8図に示す。
As a verification of the effectiveness of the above analysis method, data |I〓(ω)|, |u〓(ω
)
|, |E〓(ω)| are generated and intensity waveforms and phase waveforms are reproduced from them using the above analysis method. The results are shown in FIGS. 6 to 8.

各図中右半分に、各々の原波形について想定さ
れるデータを示した。
The right half of each figure shows data assumed for each original waveform.

|E〓(ω)|、|u〓(ω)|については、測定デー

のスペクトル解析から得られたものを、それぞれ
−ω0、−2ω0だけ平行移動した後のものを示して
ある。
Regarding |E〓(ω)| and |u〓(ω)|, the values obtained from the spectrum analysis of the measured data are shown after parallel translation by −ω 0 and −2ω 0 , respectively.

第6図は原波形が対称の場合であり、第7図、
第8図は非対称な原強度波形、原位相波形の場合
である。これらより強度波形、位相波形ともに高
い精度で再生され、かつパルスが非対称であつて
も精度の低下が生じないことが明らかになつた。
Figure 6 shows the case where the original waveform is symmetrical, Figure 7,
FIG. 8 shows the case of asymmetric original intensity waveforms and original phase waveforms. From these results, it has become clear that both the intensity waveform and the phase waveform can be reproduced with high precision, and that even if the pulses are asymmetrical, there is no decrease in precision.

第9図に本実施例による短光パルス評価方法を
用いて半導体レーザーの発生する光パルスを測定
評価して得られた強度波形、位相波形を示す。使
用した半導体レーザーは波長1.3μmのInGaAsP
DFBレーザで、注入電流に周波数200MHzの正弦
波を重畳し、変調することによりパルス繰り返し
間隔P=5ns、パルス幅20〜30ps内外のパルスを
発生させた。移動台速度Vは3.88mm/minとし、
基本波の波長で繰り返される干渉フリンジが周波
数 2V/λ0=2×3.88/60/1.3×10-3=99.5Hz で信号に現れる条件下で測定を行つた。
FIG. 9 shows intensity waveforms and phase waveforms obtained by measuring and evaluating optical pulses generated by a semiconductor laser using the short optical pulse evaluation method according to this embodiment. The semiconductor laser used is InGaAsP with a wavelength of 1.3 μm.
Using a DFB laser, a sine wave with a frequency of 200 MHz was superimposed on the injected current and modulated to generate pulses with a pulse repetition interval P = 5 ns and a pulse width of 20 to 30 ps. The moving table speed V is 3.88 mm/min,
Measurements were carried out under conditions in which interference fringes repeated at the wavelength of the fundamental wave appeared in the signal at a frequency of 2V/λ 0 =2×3.88/60/1.3×10 -3 =99.5Hz.

ここで式(11)を満足するよう積分増幅器の時定数
Tは100μm程度に設定している。またλref/8Vは
約1msであり、被測定パルスの繰り返し間隔Pは
前述のように5nsである。
Here, the time constant T of the integrating amplifier is set to about 100 μm so as to satisfy equation (11). Further, λ ref /8V is approximately 1 ms, and the repetition interval P of the pulse to be measured is 5 ns as described above.

第9図の解析結果より、被測定パルスは急峻な
立ち上がりとなだらかな裾引きを有し、半値全幅
は23psであることがわかる。
From the analysis results in FIG. 9, it can be seen that the pulse to be measured has a steep rise and a gentle tail, and the full width at half maximum is 23 ps.

また位相波形は強度波形を積分した形に似た変
化をしていることがわかる。
It can also be seen that the phase waveform changes in a manner similar to the integrated intensity waveform.

〔発明の効果〕〔Effect of the invention〕

以上の説明より明らかなように、本発明による
短光パルス評価方法によれば、小型低廉な測定装
置を用いて、被測定光パルスの強度波形及び位相
波形を従来法に比較して高精度に求めることが可
能になる。
As is clear from the above explanation, according to the short optical pulse evaluation method according to the present invention, the intensity waveform and phase waveform of the optical pulse to be measured can be measured with higher precision than the conventional method using a small and inexpensive measuring device. It becomes possible to ask for it.

本発明により測り得る光パルスの時間幅の下限
は、該光パルスの光振動の周期の10倍程度と考え
られる。これは2次高調波光強度変化のスペクト
ル解析に際し、3つのピークが充分分離独立して
現れることが必要だからである。
The lower limit of the time width of a light pulse that can be measured by the present invention is considered to be about 10 times the period of optical oscillation of the light pulse. This is because it is necessary for the three peaks to appear sufficiently separately and independently when performing spectrum analysis of second harmonic light intensity changes.

従つて、例えば、1.3μm帯の光パルスならば測
定可能な光パルスの時間幅の下限は40fs(40×
10-15s)程度、また0.6μm帯ならば20fs程度とな
る。
Therefore, for example, for an optical pulse in the 1.3 μm band, the lower limit of the measurable time width of the optical pulse is 40 fs (40 ×
10 -15 s), or about 20fs in the 0.6μm band.

一方、本発明により測り得る光パルスの時間幅
の上限は原理的には存在せず、如何に長い光パル
スも測定可能である。実用的観点からは、本方法
は小型低廉な装置によつて実施可能な故、100ps
程度迄の光パルスの測定に利用されると考える。
On the other hand, there is no upper limit in principle to the time width of a light pulse that can be measured by the present invention, and any long light pulse can be measured. From a practical point of view, this method can be implemented using small and inexpensive equipment, so
It is thought that it can be used to measure optical pulses up to a certain degree.

また、本発明は適当な2次高調波発生能を有す
る結晶が存在する限りの波長の光パルスを測定で
きる。一例としてβ−BaB2O4(β−ホウ酸バリ
ウム)結晶を用いると、0.4〜2.5μmの波長で、更
にKTP(リン酸3水素カリウム)結晶で4.2μmの
長波長まで測定可能である。
Furthermore, the present invention can measure optical pulses of any wavelength as long as a crystal with an appropriate second harmonic generation ability exists. As an example, when a β-BaB 2 O 4 (β-barium borate) crystal is used, it is possible to measure at a wavelength of 0.4 to 2.5 μm, and with a KTP (potassium trihydrogen phosphate) crystal, it is possible to measure up to a longer wavelength of 4.2 μm.

以上、本発明は広範な時間幅、および波長の光
パルスに適用される。
As described above, the present invention is applicable to light pulses having a wide range of time widths and wavelengths.

本発明によれば、パルスが非対称であつても強
度、位相波形を詳細に求められるので、短光パル
ス光源の調整、評価に極めて有効である。
According to the present invention, the intensity and phase waveform can be obtained in detail even if the pulse is asymmetric, so it is extremely effective for adjusting and evaluating short optical pulse light sources.

また、光フアイバ伝送線路の通過前後の光パル
スを本発明により測定すれば、その強度、位相波
形の変化が詳細に観察でき、該光フアイバ線路の
特性評価を直接的かつ容易に行なうことが可能で
ある。同様のことは、一般の光学システム、光材
料についても行なえる。
Furthermore, by measuring optical pulses before and after passing through an optical fiber transmission line using the present invention, changes in their intensity and phase waveform can be observed in detail, making it possible to directly and easily evaluate the characteristics of the optical fiber line. It is. The same thing can be done for general optical systems and optical materials.

更に、本発明でのパルスの位相の詳細な観測
は、光フアイバを用いるパルス圧縮器の最適化、
調整に有力な手段を提供する。
Furthermore, detailed observation of the pulse phase in the present invention can be achieved by optimizing a pulse compressor using optical fibers,
Provide a powerful means for coordination.

この様に本発明は、短光パルスの生成、整形、
伝搬等、短光パルスに関連したあらゆる分野で広
く活用される基本的な方法である。
In this way, the present invention provides short optical pulse generation, shaping,
This is a basic method that is widely used in all fields related to short optical pulses, such as propagation.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明による短光パルス測定評価方法
の実施例の構成の概略図、第2図は短光パルス測
定評価方法の第1の従来例の構成概略図、第3図
は第2の従来例の構成概略図である。第4図は本
発明の方法により測定されたデータ及びそのフー
リエ解析より|I〓|2、|u〓|2、|E〓|2を求めた結

を表す図である。第5図は上記フーリエ解析結果
より、元の被測定光パルスを再構築するための反
復計算法を示す。第6図乃至第8図は、本発明の
実施例の解析方法の妥当性の検証に用いた種々の
原パルスの強度波形及び位相波形(実線)、解析
に使つたデータの組、及びそれらデータから実施
例の解析方法を用いて再構築されたパルスの強度
波形及び位相波形(点線)を示す図である。第9
図は本発明の実施例において実測データを解析し
て得られた被測定パルスの強度波形及び位相波形
を示す図である。 1……半透鏡、2,3……プリズム、4……移
動台、5,6……レンズ、7……2次高調波光発
生能を有する結晶、8……光学フイルタ、9……
光検出器、10……積分増幅器、11……ヘリウ
ムネオンレーザー、12,13……鏡、14……
レンズ、15……光検出器、16……増幅器、1
7……波形記憶装置、18……コンピユータ1
8、19……4分の1波長板、20……偏光ビー
ムスプリツタ、21,22……光検出器、23…
…サンプルトリガ発生器、24……半透鏡、25
……レンズ、26……分光器、27……波長送り
装置、28……光検出器、29……積分増幅器、
30……光学フイルタ、31……切換スイツチ、
32……A/D変換器。
FIG. 1 is a schematic diagram of the configuration of an embodiment of the short optical pulse measurement and evaluation method according to the present invention, FIG. 2 is a schematic diagram of the configuration of the first conventional example of the short optical pulse measurement and evaluation method, and FIG. FIG. 2 is a schematic configuration diagram of a conventional example. FIG. 4 is a diagram showing the results of |I〓| 2 , |u〓| 2 , and |E〓| 2 from the data measured by the method of the present invention and its Fourier analysis. FIG. 5 shows an iterative calculation method for reconstructing the original optical pulse to be measured from the above Fourier analysis results. Figures 6 to 8 show the intensity waveforms and phase waveforms (solid lines) of various original pulses used to verify the validity of the analysis method of the embodiment of the present invention, the data sets used for analysis, and the data. FIG. 3 is a diagram showing the intensity waveform and phase waveform (dotted line) of a pulse reconstructed from FIG. 1 using the analysis method of the example. 9th
The figure is a diagram showing an intensity waveform and a phase waveform of a pulse to be measured obtained by analyzing measured data in an example of the present invention. DESCRIPTION OF SYMBOLS 1... Half-transparent mirror, 2, 3... Prism, 4... Moving table, 5, 6... Lens, 7... Crystal having ability to generate second harmonic light, 8... Optical filter, 9...
Photodetector, 10... Integrating amplifier, 11... Helium neon laser, 12, 13... Mirror, 14...
Lens, 15...Photodetector, 16...Amplifier, 1
7...Waveform storage device, 18...Computer 1
8, 19... Quarter wavelength plate, 20... Polarizing beam splitter, 21, 22... Photodetector, 23...
...Sample trigger generator, 24...Semi-transparent mirror, 25
... Lens, 26 ... Spectrometer, 27 ... Wavelength sending device, 28 ... Photodetector, 29 ... Integrating amplifier,
30...Optical filter, 31...Selector switch,
32...A/D converter.

Claims (1)

【特許請求の範囲】 1 被測定光パルス光束を2光束に分枝し、 分枝した光路に相対的光路長差の変化を付与し
つつ前記2光束を合波し、 該合波光束を2次高調波光発生能を有する結晶
に入射せしめて2次高調波光を発生させ、 前記結晶を通過した基本波光と前記2次高調波
光との強度に比例した電気信号にそれぞれ変換
し、 前記相対的光路長差に対する基本波光および2
次高調波光の強度変化をそれぞれ記録し、 記録されたデータのフーリエ解析を行い被測定
光パルスの強度波形と位相波形とを計算により求
めることを特徴とする短光パルス測定評価方法。 2 基本波光の強度変化データのフーリエ解析よ
り、基本波光周波数ω0近傍のスペクトルとして
被測定光パルスの電場のフーリエ変換を求め、ま
た2次高調波光の強度変化データのフーリエ解析
より、2ω0近傍のスペクトルとして2次高調波光
の電場のフーリエ変換を、零周波数近傍のスペク
トルとして被測定光パルス強度のフーリエ変換を
それぞれ求め、以上求めた3者のフーリエ変換を
もとに反復計算を行う特許請求の範囲第1項記載
の短光パルス測定評価方法。
[Scope of Claims] 1. Branching the optical pulse beam to be measured into two beams, combining the two beams while imparting a change in relative optical path length difference to the branched optical path, and combining the combined beam into two beams. The second harmonic light is made incident on a crystal having the ability to generate second harmonic light, and is converted into electric signals proportional to the intensities of the fundamental wave light and the second harmonic light that have passed through the crystal, and the relative optical path is Fundamental wave light and 2 for length difference
A short optical pulse measurement and evaluation method characterized by recording each intensity change of harmonic light, performing Fourier analysis of the recorded data, and calculating the intensity waveform and phase waveform of the optical pulse to be measured. 2 From the Fourier analysis of the intensity change data of the fundamental wave light, we found the Fourier transform of the electric field of the measured optical pulse as a spectrum near the fundamental wave optical frequency ω 0 , and from the Fourier analysis of the intensity change data of the second harmonic light, we found the spectrum near 2ω 0 . A patent claim in which the Fourier transform of the electric field of the second harmonic light is determined as the spectrum of the second harmonic light, and the Fourier transform of the measured optical pulse intensity is determined as the spectrum near zero frequency, and the iterative calculation is performed based on the Fourier transforms of the three components determined above. The short optical pulse measurement and evaluation method according to item 1.
JP7354787A 1986-09-08 1987-03-27 Measuring and evaluating method for short light pulse Granted JPS63238525A (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP7354787A JPS63238525A (en) 1987-03-27 1987-03-27 Measuring and evaluating method for short light pulse
US07/093,653 US4792230A (en) 1986-09-08 1987-09-08 Method and apparatus for measuring ultrashort optical pulses

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP7354787A JPS63238525A (en) 1987-03-27 1987-03-27 Measuring and evaluating method for short light pulse

Publications (2)

Publication Number Publication Date
JPS63238525A JPS63238525A (en) 1988-10-04
JPH052928B2 true JPH052928B2 (en) 1993-01-13

Family

ID=13521368

Family Applications (1)

Application Number Title Priority Date Filing Date
JP7354787A Granted JPS63238525A (en) 1986-09-08 1987-03-27 Measuring and evaluating method for short light pulse

Country Status (1)

Country Link
JP (1) JPS63238525A (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5123646B2 (en) * 2007-11-19 2013-01-23 公立大学法人大阪府立大学 Ultrashort optical pulse intensity / phase information reconstruction method and laser processing equipment using this method
JP6081520B2 (en) * 2014-05-28 2017-02-15 インディアン インスティテュート オブ テクノロジー デリー Non-interfering phase measurement
FR3022346B1 (en) * 2014-06-16 2022-10-07 Commissariat Energie Atomique DEVICE AND METHOD FOR CHARACTERIZING A LIGHT BEAM

Also Published As

Publication number Publication date
JPS63238525A (en) 1988-10-04

Similar Documents

Publication Publication Date Title
US4792230A (en) Method and apparatus for measuring ultrashort optical pulses
CN110646805B (en) Frequency modulation continuous wave laser ranging system based on virtual sweep frequency light source
Iaconis et al. Self-referencing spectral interferometry for measuring ultrashort optical pulses
JP4951189B2 (en) Frequency conversion phase shift interferometry
Brun et al. Single-shot characterization of ultrashort light pulses
KR20110036944A (en) Interferometer and synchronization design with frequency combs
US20030025911A1 (en) Optical pulse measurement
CN105866788A (en) Distance measuring device and distance measuring device method for realizing optical sampling by adjusting resonant cavity of femtosecond laser
US5684586A (en) Apparatus for characterizing short optical pulses
Diels et al. Ultrafast diagnostics
JPH052928B2 (en)
CN120651093A (en) Doppler effect compensation method and device for optical time stretching sweep frequency interferometry
CN115655663B (en) Linewidth measurement method and system of all-fiber structure laser
CN112881797A (en) Single multipath synchronous measurement method and device based on all-fiber spectral interference
JP3607622B2 (en) Electric field cross-correlator
JP4002163B2 (en) Method and apparatus for measuring phase difference between pulses
JP2767000B2 (en) Waveguide dispersion measurement method and apparatus
Wei et al. Do we need all the frequency components of a fringe signal to obtain position information in a vertical scanning wideband interferometer?
JPH052252B2 (en)
Sarger et al. How to measure the characteristics of laser pulses
JP4296070B2 (en) Phase characteristic measuring device
JP2694201B2 (en) Dispersion measuring method and apparatus
JP3544901B2 (en) Method and apparatus for measuring time waveform of optical signal electric field
JP2660291B2 (en) Optical pulse measuring device
JPH02296166A (en) Electrical signal waveform observation device using electro-optical effect

Legal Events

Date Code Title Description
EXPY Cancellation because of completion of term
FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080113

Year of fee payment: 15