Data creation apparatus, light control apparatus, data creation method, and data creation program

Abstract

An iterative Fourier transform unit of a modulation pattern calculation apparatus performs a Fourier transform on a waveform function including an intensity spectrum function and a phase spectrum function, performs a replacement of a temporal intensity waveform function based on a desired waveform after the Fourier transform, and then performs an inverse Fourier transform. The iterative Fourier transform unit performs the replacement using a result of multiplying a function representing the desired waveform by a coefficient. The coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

Claims

1. A data creation method for creating, by a computer, data for controlling a spatial light modulator in an optical system for generating, from input light, output light having a temporal intensity waveform different from that of the input light, the method comprising: generating a phase spectrum function to be used for creating data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; creating the data on the basis of the intensity spectrum function and the phase spectrum function; providing a control signal based on the created data to the spatial light modulator; and controlling the spatial light modulator on the basis of the control signal, wherein the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication, and the coefficient has a value with which a time integral value of a result of multiplying a function, including a difference between the temporal intensity waveform function after the Fourier transform and the function after the multiplication, by a first weight function having a larger weight value when the intensity is higher at each time of the function before the multiplication, is minimized.

2. The data creation method according to claim 1, wherein generating the phase spectrum function to be used for creating the data by repeatedly performing the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data is created on the basis of the generated phase spectrum function by the repeatedly performing of the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement.

3. The data creation method according to claim 1, wherein the first weight function includes a function obtained by multiplying the function by another coefficient.

4. The data creation method according to claim 1, wherein the first replacement is performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function has a larger weight value when the intensity is higher at each time of the function before the multiplication.

5. A data creation method for creating, by a computer, data for controlling a spatial light modulator in an optical system for generating, from input light, output light having a temporal intensity waveform different from that of the input light, the method comprising: generating a phase spectrum function to be used for creating data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; creating the data on the basis of the intensity spectrum function and the phase spectrum function; providing a control signal based on the created data to the spatial light modulator; and controlling the spatial light modulator on the basis of the control signal, wherein the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication, the first replacement is performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function has a larger weight value when the intensity is higher at each time of the function before the multiplication.

6. The data creation method according to claim 5, wherein generating the phase spectrum function to be used for creating the data by repeatedly performing the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data is created on the basis of the generated phase spectrum function by the repeatedly performing of the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement.

7. A computer including a processor and a memory that stores instructions, the instructions, when executed by the processor perform the method comprising: generating a phase spectrum function to be used for creating data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; creating the data on the basis of the intensity spectrum function and the phase spectrum function; providing a control signal based on the created data to the spatial light modulator; and controlling the spatial light modulator on the basis of the control signal, wherein the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication, and the coefficient has a value with which a time integral value of a result of multiplying a function, including a difference between the temporal intensity waveform function after the Fourier transform and the function after the multiplication, by a first weight function having a larger weight value when the intensity is higher at each time of the function before the multiplication, is minimized.

8. The computer according to claim 7, wherein generating the phase spectrum function to be used for creating the data by repeatedly performing the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data is created on the basis of the generated phase spectrum function by the repeatedly performing of the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement.

9. The computer according to claim 7, wherein the first weight function includes a function obtained by multiplying the function by another coefficient.

10. The computer according to claim 7, wherein the first replacement is performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function has a larger weight value when the intensity is higher at each time of the function before the multiplication.

11. A computer including a processor and a memory that stores instructions, the instructions, when executed by the processor perform the method comprising: generating a phase spectrum function to be used for creating data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; creating the data on the basis of the intensity spectrum function and the phase spectrum function; providing a control signal based on the created data to the spatial light modulator; and controlling the spatial light modulator on the basis of the control signal, wherein the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication, the first replacement is performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function has a larger weight value when the intensity is higher at each time of the function before the multiplication.

12. The computer according to claim 11, wherein generating the phase spectrum function to be used for creating the data by repeatedly performing the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data is created on the basis of the generated phase spectrum function by the repeatedly performing of the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement.

13. A non-transitory computer-readable medium storing a program, the program, when executed by a computer, performing a data creation method comprising: generating a phase spectrum function to be used for creating data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; creating the data on the basis of the intensity spectrum function and the phase spectrum function; providing a control signal based on the created data to the spatial light modulator; and controlling the spatial light modulator on the basis of the control signal, wherein the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication, and the coefficient has a value with which a time integral value of a result of multiplying a function, including a difference between the temporal intensity waveform function after the Fourier transform and the function after the multiplication, by a first weight function having a larger weight value when the intensity is higher at each time of the function before the multiplication, is minimized.

14. The non-transitory computer-readable medium according to claim 13, wherein generating the phase spectrum function to be used for creating the data by repeatedly performing the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data is created on the basis of the generated phase spectrum function by the repeatedly performing of the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement.

15. The non-transitory computer-readable medium according to claim 13, wherein the first weight function includes a function obtained by multiplying the function by another coefficient.

16. The non-transitory computer-readable medium according to claim 13, wherein the first replacement is performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function has a larger weight value when the intensity is higher at each time of the function before the multiplication.

17. A non-transitory computer-readable medium storing a program, the program, when executed by a computer, performing a data creation method comprising: generating a phase spectrum function to be used for creating data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; creating the data on the basis of the intensity spectrum function and the phase spectrum function; providing a control signal based on the created data to the spatial light modulator; and controlling the spatial light modulator on the basis of the control signal, wherein the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication, the first replacement is performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function has a larger weight value when the intensity is higher at each time of the function before the multiplication.

18. The non-transitory computer-readable medium according to claim 17, wherein generating the phase spectrum function to be used for creating the data by repeatedly performing the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data is created on the basis of the generated phase spectrum function by the repeatedly performing of the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a diagram schematically illustrating a configuration of a light control apparatus according to an embodiment.

(2) FIG. 2 is a diagram illustrating a configuration of an optical system included in the light control apparatus.

(3) FIG. 3 is a diagram illustrating a modulation plane of an SLM.

(4) FIG. 4A and FIG. 4B are (A) a diagram illustrating, as an example, a spectrum waveform (a spectrum phase and a spectrum intensity) of single pulse shaped input light, and (B) a diagram illustrating a temporal intensity waveform of the input light.

(5) FIG. 5A and FIG. 5B are (A) a diagram illustrating, as an example, a spectrum waveform (a spectrum phase and a spectrum intensity) of output light when rectangular wave shaped phase spectrum modulation is applied in the SLM, and (B) a diagram illustrating a temporal intensity waveform of the output light.

(6) FIG. 6 is a diagram schematically illustrating a hardware configuration example of a modulation pattern calculation apparatus.

(7) FIG. 7 is a diagram illustrating a phase spectrum calculation procedure by a normal iterative Fourier method as a first comparative example.

(8) FIG. 8 is a diagram illustrating a phase spectrum calculation procedure by an iterative Fourier method described in Non Patent Document 2 as a second comparative example.

(9) FIG. 9 is a diagram illustrating a phase spectrum calculation procedure described in Patent Document 1.

(10) FIG. 10A and FIG. 10B are diagrams for visually explaining a process of a process number (4) in FIG. 9.

(11) FIG. 11 is a diagram illustrating a phase spectrum calculation procedure by an improved iterative Fourier method in an embodiment.

(12) FIG. 12A to FIG. 12C are diagrams for visually explaining a process of a process number (4) in FIG. 11.

(13) FIG. 13 is a diagram illustrating an example of a weight function We(t) when Target.sub.0(t) is a multi pulse.

(14) FIG. 14 is a diagram illustrating a calculation procedure in an iterative Fourier transform unit of an intensity spectrum design unit.

(15) FIG. 15 is a flowchart illustrating a modulation pattern calculation method.

(16) FIG. 16 is a graph illustrating an output temporal waveform when a phase spectrum function Ψ.sub.IFTA(ω) and an intensity spectrum function A.sub.IFTA(ω) are calculated using a method described in Patent Document 1.

(17) FIG. 17 is a graph illustrating an output temporal waveform when a phase spectrum function Ψ.sub.IFTA(ω) and an intensity spectrum function A.sub.IFTA(ω) are calculated using a method according to an embodiment.

(18) FIG. 18 is a graph illustrating a result of improving the uniformity of peak values by the calculated phase spectrum function Ψ.sub.IFTA(ω) and intensity spectrum function A.sub.IFTA(ω).

(19) FIG. 19 is a graph illustrating an output temporal waveform when a phase spectrum function Ψ.sub.IFTA(ω) and an intensity spectrum function A.sub.IFTA(ω) are calculated using a method according to an embodiment.

(20) FIG. 20 is a diagram illustrating a phase spectrum calculation procedure by an iterative Fourier method according to a modification.

(21) FIG. 21 is a diagram illustrating an intensity spectrum calculation procedure by an iterative Fourier method according to a modification.

(22) FIG. 22A to FIG. 22C are diagrams for visually explaining a process of a process number (5) in an embodiment.

(23) FIG. 23A to FIG. 23C are diagrams for visually explaining a process of a process number (5) in a modification.

DETAILED DESCRIPTION

(24) Hereinafter, embodiments of a data creation apparatus, a light control apparatus, a data creation method, and a data creation program will be described in detail with reference to the accompanying drawings. In the description of the drawings, the same elements will be denoted by the same reference symbols, without redundant description.

(25) FIG. 1 is a diagram schematically illustrating a configuration of a light control apparatus 1A according to an embodiment. FIG. 2 is a diagram illustrating a configuration of an optical system 10 included in the light control apparatus 1A. The light control apparatus 1A according to the present embodiment generates, from input light La, output light Ld having an arbitrary temporal intensity waveform different from that of the input light La. As illustrated in FIG. 1, the light control apparatus 1A includes a light source 2, an optical system 10, and a modulation pattern calculation apparatus (data creation apparatus) 20.

(26) The light source 2 outputs the input light La input to the optical system 10. The light source 2 is a laser light source such as a solid-state laser light source and a fiber laser light source, for example, and the input light La is coherent pulse light, for example. The optical system 10 has an SLM 14, and receives a control signal SC for controlling each pixel of the SLM 14 from the modulation pattern calculation apparatus 20 by the SLM 14. The optical system 10 converts the input light La from the light source 2 into the output light Ld having the arbitrary temporal intensity waveform. A modulation pattern is data for controlling the SLM 14, and is data in which intensities of a complex amplitude distribution or intensities of a phase distribution are output in a file. The modulation pattern is, for example, a computer-generated hologram (CGH).

(27) As illustrated in FIG. 2, the optical system 10 includes a diffraction grating 12, a lens 13, an SLM 14, a lens 15, and a diffraction grating 16. The diffraction grating 12 is a dispersive element in the present embodiment, and is optically coupled to the light source 2. The SLM 14 is optically coupled to the diffraction grating 12 via the lens 13. The diffraction grating 12 spectrally disperses the input light La for each wavelength component. In addition, as the dispersive element, other optical components such as a prism may be used instead of the diffraction grating 12. Further, the dispersive element may be of a reflection type or a transmission type.

(28) The input light La is obliquely incident on the diffraction grating 12, and is spectrally dispersed into a plurality of wavelength components. Light Lb including the plurality of wavelength components is focused for each wavelength component by the lens 13, and forms an image on a modulation plane of the SLM 14. The lens 13 may be a convex lens made of a light transmitting member or a concave mirror having a concave light reflection surface.

(29) The SLM 14 simultaneously performs phase modulation and intensity modulation of the light Lb to generate the output light Ld having an arbitrary temporal intensity waveform different from that of the input light La. The SLM 14 may perform only the intensity modulation. The SLM 14 is of a phase modulation type, for example. In one example, the SLM 14 is of a liquid crystal on silicon (LCOS) type. Further, the SLM 14 may be an intensity modulation type SLM such as a digital micromirror device (DMD). Further, the SLM 14 may be of a reflection type or a transmission type.

(30) FIG. 3 is a diagram illustrating a modulation plane 17 of the SLM 14. As illustrated in FIG. 3, in the modulation plane 17, a plurality of modulation regions 17a are arranged along a certain direction A, and each modulation region 17a extends in a direction B intersecting with the direction A. The direction A is a dispersing direction by the diffraction grating 12. The modulation plane 17 functions as a Fourier transform plane, and each corresponding wavelength component after the dispersion is incident on each of the plurality of modulation regions 17a. The SLM 14 modulates a phase and an intensity of each incident wavelength component, independently from the other wavelength components, in each modulation region 17a. In addition, because the SLM 14 in the present embodiment is of the phase modulation type, the intensity modulation is realized by a phase pattern (phase image) presented on the modulation plane 17.

(31) Each wavelength component of modulated light Lc modulated by the SLM 14 is focused at a point on the diffraction grating 16 by the lens 15. At this time, the lens 15 functions as a focusing optical system for focusing the modulated light Lc. The lens 15 may be a convex lens made of a light transmitting member or a concave mirror having a concave light reflection surface. Further, the diffraction grating 16 functions as a combining optical system and combines the respective wavelength components after the modulation. That is, by the lens 15 and the diffraction grating 16, the plurality of wavelength components of the modulated light Lc are focused and combined with each other and become the output light Ld.

(32) A region (spectrum domain) before the lens 15 and a region (time domain) behind the diffraction grating 16 are in a Fourier transform relation with each other, and the phase modulation and the intensity modulation in the spectrum domain affect the temporal intensity waveform in the time domain. Therefore, the output light Ld has a desired temporal intensity waveform different from that of the input light La in accordance with the modulation pattern of the SLM 14.

(33) Here, as an example, FIG. 4A illustrates a spectrum waveform (a spectrum phase G11 and a spectrum intensity G12) of the single pulse shaped input light La, and FIG. 4B illustrates the temporal intensity waveform of the input light La. Further, as an example, FIG. 5A illustrates a spectrum waveform (a spectrum phase G21 and a spectrum intensity G22) of the output light Ld when rectangular wave shaped phase spectrum modulation is applied in the SLM 14, and FIG. 5B illustrates a temporal intensity waveform of the output light Ld. In FIG. 4A and FIG. 5A, a horizontal axis shows a wavelength (nm), a left vertical axis shows an intensity value (arbitrary unit) of the intensity spectrum, and a right vertical axis shows a phase value (rad) of the phase spectrum. Further, in FIG. 4B and FIG. 5B, a horizontal axis shows a time (femtosecond) and a vertical axis shows a light intensity (arbitrary unit).

(34) In this example, a single pulse of the input light La is converted into a double pulse with high-order light as the output light Ld by giving a rectangular wave shaped phase spectrum waveform to the output light Ld. In addition, the spectra and the waveforms illustrated in FIG. 4A, FIG. 4B, FIG. 5A, and FIG. 5B are examples, and the temporal intensity waveform of the output light Ld can be controlled in various shapes by combinations of various phase spectra and intensity spectra.

(35) FIG. 1 is referred to again. The modulation pattern calculation apparatus 20 is a computer having a processor including, for example, a personal computer, a smart device such as a smart phone and a tablet terminal, and a cloud server. The modulation pattern calculation apparatus 20 is electrically coupled to the SLM 14, calculates a phase modulation pattern to bring the temporal intensity waveform of the output light Ld closer to a desired waveform, and provides a control signal SC including the phase modulation pattern to the SLM 14.

(36) The modulation pattern calculation apparatus 20 of the present embodiment presents, on the SLM 14, a phase pattern including a phase pattern for phase modulation that gives a phase spectrum for obtaining the desired waveform to the output light Ld and a phase pattern for intensity modulation that gives an intensity spectrum for obtaining the desired waveform to the output light Ld. For this purpose, the modulation pattern calculation apparatus 20 includes an arbitrary waveform input unit 21, a phase spectrum design unit 22, an intensity spectrum design unit 23, and a modulation pattern generation unit (data generation unit) 24.

(37) That is, the processor of the computer provided in the modulation pattern calculation apparatus 20 realizes a function of the arbitrary waveform input unit 21, a function of the phase spectrum design unit 22, a function of the intensity spectrum design unit 23, and a function of the modulation pattern generation unit 24. The respective functions may be realized by the same processor or may be realized by different processors.

(38) FIG. 6 is a diagram schematically illustrating a hardware configuration example of the modulation pattern calculation apparatus 20. As illustrated in FIG. 6, the modulation pattern calculation apparatus 20 can be physically configured as a normal computer including a processor (CPU) 201, a main storage device such as a ROM 202 and a RAM 203, an input device 204 such as a keyboard, a mouse, and a touch screen, an output device 205 such as a display (including a touch screen), a communication module 206 such as a network card for transmitting/receiving data to/from other devices, an auxiliary storage device 207 such as a hard disk, and the like.

(39) The processor 201 of the computer can realize the above respective functions (the arbitrary waveform input unit 21, the phase spectrum design unit 22, the intensity spectrum design unit 23, and the modulation pattern generation unit 24) by a modulation pattern calculation program (data creation program). Therefore, the modulation pattern calculation program operates the processor 201 of the computer as the arbitrary waveform input unit 21, the phase spectrum design unit 22, the intensity spectrum design unit 23, and the modulation pattern generation unit 24 in the modulation pattern calculation apparatus 20.

(40) The modulation pattern calculation program is stored in a storage device (storage medium) inside or outside the computer such as the auxiliary storage device 207, for example. The storage device may be a non-transitory recording medium. Examples of the recording medium include a recording medium such as a flexible disk, a CD, and a DVD, a recording medium such as a ROM, a semiconductor memory, a cloud server, and the like.

(41) The arbitrary waveform input unit 21 receives the desired temporal intensity waveform input from an operator. The operator inputs information (for example, a pulse width, a pulse number, or the like) on the desired temporal intensity waveform to the arbitrary waveform input unit 21. The information on the desired temporal intensity waveform is given to the phase spectrum design unit 22 and the intensity spectrum design unit 23. The phase spectrum design unit 22 calculates the phase spectrum of the output light Ld suitable for realizing the given desired temporal intensity waveform. The intensity spectrum design unit 23 calculates the intensity spectrum of the output light Ld suitable for realizing the given desired temporal intensity waveform.

(42) The modulation pattern generation unit 24 calculates a phase modulation pattern (for example, a computer-generated hologram) to give the phase spectrum obtained in the phase spectrum design unit 22 and the intensity spectrum obtained in the intensity spectrum design unit 23 to the output light Ld. In addition, the control signal SC including the calculated phase modulation pattern is provided to the SLM 14, and the SLM 14 is controlled on the basis of the control signal SC.

(43) Here, a method of calculating a phase spectrum and an intensity spectrum corresponding to a desired temporal intensity waveform will be described in detail. The desired temporal intensity waveform is represented as a function in the time domain, and the phase spectrum and the intensity spectrum are represented as functions in the frequency domain. Therefore, the phase spectrum and the intensity spectrum corresponding to the desired temporal intensity waveform are obtained by an iterative Fourier transform based on the desired temporal intensity waveform. In the method described below, the phase spectrum and the intensity spectrum are calculated using an iterative Fourier transform method. Therefore, as illustrated in FIG. 1, the phase spectrum design unit 22 has an iterative Fourier transform unit 22a. Further, the intensity spectrum design unit 23 has an iterative Fourier transform unit 23a.

(44) However, in a normal iterative Fourier method, a solution may be led to a local solution, and the solution is not necessarily an optimal solution. Further, there are problems in that it takes long time to converge the calculation (that is, the number of iterations increases), and the intensity loss increases. FIG. 7 illustrates a phase spectrum calculation procedure by the normal iterative Fourier method, as a first comparative example for the present embodiment.

(45) First, an initial intensity spectrum function A.sub.0(ω) and a phase spectrum function Ψ.sub.n=0(ω) to be functions of a frequency ω are prepared (process number (1) in the drawing). In one example, the intensity spectrum function A.sub.0(ω) and the phase spectrum function Ψ.sub.n=0(ω) represent the intensity spectrum and the phase spectrum of the input light La, respectively. Next, a waveform function (a) in the frequency domain including the intensity spectrum function A.sub.0(ω) and the phase spectrum function Ψ.sub.n(ω) is prepared (process number (2) in the drawing).
√{square root over (A.sub.0(ω))}exp{iΨ.sub.n(ω)}  (a)
A subscript n represents after an n-th Fourier transform process. Before a first Fourier transform process, the initial phase spectrum function Ψ.sub.n=0(ω) described above is used as the phase spectrum function Ψ.sub.n(ω). i is an imaginary number.

(46) Next, a Fourier transform from the frequency domain to the time domain is performed on the function (a) (arrow A1 in the drawing). As a result, a waveform function (b) in the frequency domain including a temporal intensity waveform function b.sub.n(t) is obtained (process number (3) in the drawing).
√{square root over (b.sub.n(t))}exp{iΘ.sub.n(t)}  (b)

(47) Next, the temporal intensity waveform function b.sub.n(t) included in the function (b) is replaced by Target.sub.0(t) based on the desired waveform (process numbers (4) and (5) in the drawing).
b.sub.n(t)Target.sub.0(t)  (c)
√{square root over (Target.sub.0(t))}exp{iΘ.sub.n(t)}  (d)

(48) Next, an inverse Fourier transform from the time domain to the frequency domain is performed on the function (d) (arrow A2 in the drawing). As a result, a waveform function (e) in the frequency domain including an intensity spectrum function B.sub.n(ω) and the phase spectrum function Ψ.sub.n(ω) is obtained (process number (6) in the drawing).
√{square root over (B.sub.n(ω))}exp{iΨ.sub.n(ω)}  (e)

(49) Next, to constrain the intensity spectrum function B.sub.n(ω) included in the function (e), this is replaced by the initial intensity spectrum function A.sub.0(ω) (process number (7) in the drawing).
B.sub.n(ω)A.sub.0(ω)  (f)

(50) Subsequently, the above processes (1) to (7) are repeatedly performed a plurality of times, so that the phase spectrum shape represented by the phase spectrum function Ψ.sub.n(ω) in the waveform function can be brought closer to a phase spectrum shape corresponding to the desired temporal intensity waveform. A phase spectrum function Ψ.sub.IFTA(ω) to be finally obtained is used for calculating the modulation pattern.

(51) For the normal iterative Fourier method described above, a process for preventing a solution from being led to a local solution is included in a method described in Non Patent Document 2. FIG. 8 illustrates a phase spectrum calculation procedure by an iterative Fourier method (hereinafter, referred to as IFTA-Fienup) described in Non Patent Document 2, as a second comparative example for the present embodiment. In addition, in FIG. 8, because the processes (1) to (3) and (6) and (7) are the same as those in the method described above, the description thereof is omitted.

(52) In IFTA-Fienup, when processes (4) and (5) are executed, that is, a replacement based on a desired waveform is performed on the temporal intensity waveform function b.sub.n(t) included in the function (b) after the Fourier transform, instead of Target.sub.0(t), Target.sub.n(t) calculated by the following formula (g) is used (process numbers (4) and (5) in the drawing).

(53) $\begin{array}{cc}Targe{t}_n\left(t\right)=\{\begin{array}{c}Targe{t}_0\left(t\right)+\beta \left(Targe{t}_0\left(t\right)-b_n\left(t\right)\right)\\ 0,\mathrm{if}{\mathrm{Target}}_0\left(t\right)+\beta \left(Targe{t}_0\left(t\right)-b_n\left(t\right)\right)<0\end{array}& \left(g\right)\end{array}$$\begin{array}{cc}\sqrt{Targe{t}_n\left(t\right)}\exp \left\{i{\Theta }_n\left(t\right)\right\}& \left(h\right)\end{array}$

(54) In addition, in the above formula (g), a difference (Target.sub.0(t)−b.sub.n(t)) between the function Target.sub.0(t) representing the desired waveform and the waveform function b.sub.n(t) after the Fourier transform is multiplied by a predetermined coefficient β, and a result thereof is added to the desired waveform Target.sub.0(t), so that Target.sub.n(t) is calculated. Further, when this value is smaller than 0, Target.sub.n(t)=0 is set.

(55) However, even in IFTA-Fienup, for example, when the function Target.sub.0(t) representing the desired waveform and the waveform function b.sub.n(t) after the Fourier transform are greatly different from each other, a solution may still be led to a local solution.

(56) Therefore, in the method described in Patent Document 1, the iterative Fourier method is further improved as described below. FIG. 9 illustrates the phase spectrum calculation procedure described in Patent Document 1. First, an initial intensity spectrum function A.sub.0(ω) and a phase spectrum function Ψ.sub.n=0(ω) to be functions of a frequency ω are prepared (process number (1) in the drawing). In one example, the intensity spectrum function A.sub.0(ω) and the phase spectrum function Ψ.sub.n=0(ω) represent the intensity spectrum and the phase spectrum of the input light, respectively.

(57) Next, a waveform function (i) in the frequency domain including the intensity spectrum function A.sub.0(ω) and the phase spectrum function Ψ.sub.n(ω) is prepared (process number (2) in the drawing).
√{square root over (A.sub.0(ω))}exp{iΨ.sub.n(ω)}  (i)

(58) A subscript n represents after an n-th Fourier transform process. Before a first Fourier transform process, the initial phase spectrum function Ψ.sub.n=0(ω) described above is used as the phase spectrum function Ψ.sub.n(ω). i is an imaginary number.

(59) Next, a Fourier transform from the frequency domain to the time domain is performed on the function (1). As a result, a waveform function (j) in the frequency domain including a temporal intensity waveform function b.sub.n(t) is obtained (process number (3) in the drawing).
√{square root over (b.sub.n(t))}exp{iΘ.sub.n(t)}  (j)

(60) Next, a coefficient α with which a difference between the waveform function b.sub.n(t) after the Fourier transform and the function Target.sub.0(t) multiplied by the coefficient α (α×Target.sub.0(t)) is smaller than a difference between the waveform function b.sub.n(t) and the function Target.sub.0(t) is obtained (process number (4) in the drawing). In one example, the coefficient α with which a standard deviation σ of α×Target.sub.0(t) for the waveform function b.sub.n(t) after the Fourier transform is minimized (σ.sub.min) is derived by searching using an evaluation function shown by the following formula (k). In addition, in the formula (k), D represents the number of data points, and t.sub.s and t.sub.e represent a start point and an end point of a time axis, respectively.

(61) $\begin{array}{cc}{\sigma }_{\min }=\sqrt{\frac{1}{D}\underset{t=t_s}{\overset{t=t_e}{.Math.}}{\left(\alpha \times \mathrm{Targe}{t}_{n=0}\left(t\right)-b_n\left(t\right)\right)}^2}& \left(k\right)\end{array}$

(62) FIG. 10A and FIG. 10B are diagrams for visually explaining the process of the process number (4). FIG. 10A illustrates an example of each of the waveform function b.sub.n(t) after the Fourier transform, the function Target.sub.0(t), and α×Target.sub.0(t) after the multiplication, when the desired waveform is a protruded pulse waveform. Further, FIG. 10B illustrates a graph illustrating a calculation result of (α×Target.sub.0(t)−b.sub.n(t)).sup.2 included in the formula (k), corresponding to the example illustrated in FIG. 10A.

(63) Next, a replacement based on the desired waveform is performed on the temporal intensity waveform function b.sub.n(t) included in the function (j) after the Fourier transform (first replacement). At this time, the replacement is performed using the result (α×Target.sub.0(t)) of multiplying the function Target.sub.0(t) representing the desired waveform by the coefficient α. In one example, the replacement is performed by Target.sub.n(t) calculated by the formula (m) obtained by replacing Target.sub.0(t) in the formula (g) in IFTA-Fienup described above by α×Target.sub.0(t) (process numbers (5) and (6) in the drawing). In addition, in the formula, β is an arbitrary coefficient, and by appropriately selecting the coefficient β, it can be expected to find a better solution with a small number of iterations n and to prevent a solution from being led to a local solution.

(64) $\begin{array}{cc}Targe{t}_n\left(t\right)=\{\begin{array}{c}\alpha \times Targe{t}_0\left(t\right)+\beta \left(\alpha \times Targe{t}_0\left(t\right)-b_n\left(t\right)\right)\\ 0,\mathrm{if}\alpha \times \mathrm{Tar}ge{t}_0\left(t\right)+\beta \left(\alpha \times Targe{t}_0\left(t\right)-b_n\left(t\right)\right)<0\end{array}& \left(m\right)\end{array}$$\begin{array}{cc}\sqrt{Targe{t}_n\left(t\right)}\exp \left\{i{\Theta }_n\left(t\right)\right\}& \left(n\right)\end{array}$

(65) Next, an inverse Fourier transform from the time domain to the frequency domain is performed on the function (n) (arrow A2 in the drawing). As a result, a waveform function (o) in the frequency domain including the intensity spectrum function B.sub.n(ω) and the phase spectrum function Ψ.sub.n(ω) is obtained (process number (7) in the drawing).
√{square root over (B.sub.n(ω))}exp{iΨ.sub.n(ω)}  (o)

(66) Next, to constrain the intensity spectrum function B.sub.n(ω) included in the function (o), this is replaced by the initial intensity spectrum function A.sub.0(ω) (second replacement, process number (8) in the drawing).
B.sub.n(ω)A.sub.0(ω)  (p)

(67) Subsequently, the above processes (1) to (8) are repeatedly performed a plurality of times, so that the phase spectrum shape represented by the phase spectrum function Ψ.sub.n(ω) in the waveform function can be brought closer to a phase spectrum shape corresponding to the desired temporal intensity waveform. A phase spectrum function Ψ.sub.IFTA(ω) to be finally obtained is provided for calculating the modulation pattern.

(68) However, in the method described in Patent Document 1 described above, the evaluation function shown in the formula (k) performs uniform evaluation within a certain time region. That is, a part where the light intensity is relatively low in the temporal waveform of the output light Ld, or a part where only noise exists is evaluated with the same importance as that of a part where the light intensity is high. For this reason, the above method is effective when waveform control in the part where the light intensity is relatively low is important, however, for example, when a utility value of the temporal waveform in the part with the high light intensity is high and the temporal waveform in the part with the low light intensity can be ignored, as in laser processing or a nonlinear optical microscope, the above method is not necessarily an optimal design method.

(69) Therefore, in the iterative Fourier transform unit 22a of the phase spectrum design unit 22 in the present embodiment, an iterative Fourier method is further improved as described below. FIG. 11 illustrates a phase spectrum calculation procedure by an improved iterative Fourier method in the present embodiment. In addition, since the calculation procedure is the same as the calculation procedure illustrated in FIG. 9 in many places, the description thereof will be appropriately omitted.

(70) First, the iterative Fourier transform unit 22a performs the processes of the same process numbers (1) to (3) as those in the calculation procedure illustrated in FIG. 9. Next, the iterative Fourier transform unit 22a obtains a coefficient α having the following features (A) and (B) (process number (4) in the drawing).

(71) (A) A difference (α×Target.sub.0(t)−b.sub.n(t)) between the waveform function b.sub.n(t) after the Fourier transform and the function Target.sub.0(t) multiplied by the coefficient α is smaller than a difference (Target.sub.0(t)−b.sub.n(t)) between the waveform function b.sub.n(t) and the function Target.sub.0(t). Specifically, a time integral value of the difference (α×Target.sub.0(t)−b.sub.n(t)) is smaller than a time integral value of the difference (Target.sub.0(t)−b.sub.n(t)).
(B) At each time of the function Target.sub.0(t), a ratio of the difference (α×Target.sub.0(t)−b.sub.n(t)), that is, a ratio of the difference (α×Target.sub.0(t)−b.sub.n(t)) based on the intensity value of the function Target.sub.0(t) as a reference decreases when the intensity increases.

(72) In one example, the coefficient α with which a pseudo standard deviation σ of α×Target.sub.0(t) for the waveform function b.sub.n(t) after the Fourier transform is minimized (σ.sub.min) is derived by searching using an evaluation function shown by the following formula (q). In addition, in the formula (q), D represents the number of data points, and t.sub.s and t.sub.e represent a start point and an end point of a time axis, respectively. Further, We(t) is a first weight function.

(73) $\begin{array}{cc}{\sigma }_{\min }=\sqrt{\frac{1}{D}\underset{t=t_s}{\overset{t=t_e}{.Math.}}\mathrm{We}\left(t\right)\times {\left(\alpha \times Targe{t}_{n=0}\left(t\right)-b_n\left(t\right)\right)}^2}& \left(q\right)\end{array}$

(74) As shown in the formula (q), the evaluation function is a function including the difference (α×Target.sub.0(t)−b.sub.n(t)) between the waveform function b.sub.n(t) after the Fourier transform and the function α×Target.sub.0(t) after the multiplication, and specifically, the evaluation function includes (α×Target.sub.0(t)−b.sub.n(t)).sup.2. Further, the evaluation function includes a weight function We(t) by which the function is multiplied, and includes a time integral of a result of multiplying the function by the weight function We(t). Further, the coefficient α with which the evaluation function, that is, the time integral is minimized (σ.sub.min) is derived by searching.

(75) The weight function We(t) is a function having a larger weight value when the intensity is higher, at each time of the function Target.sub.0(t) before the multiplication. In one example, the weight function We(t) includes a function obtained by multiplying the function Target.sub.0(t) by another coefficient C.sub.1, and is represented by, for example, the following formula (r). In other words, the weight function We(t) may be determined on the basis of the function Target.sub.0(t).
We(t)=C.sub.1.Math.Target.sub.0(t)  (r)

(76) FIG. 12A to FIG. 12C are diagrams for visually explaining the process of the process number (4). FIG. 12A illustrates a weight function We(t) having the same shape as that of Target.sub.0(t) illustrated in FIG. 10A as an example of the weight function We(t). FIG. 12B is a graph illustrating a calculation result of (α×Target.sub.0(t)−b.sub.n(t)).sup.2 illustrated in FIG. 10B. FIG. 12C illustrates a graph of a result of multiplying the weight function We(t) illustrated in FIG. 12A by (α×Target.sub.0(t)−b.sub.n(t)).sup.2 illustrated in FIG. 10B.

(77) As illustrated in FIG. 12C, by multiplication of the weight function We(t), the magnitude of a difference is emphasized in a part D1 where the intensity is high in Target.sub.0(t), as compared with the other parts. Therefore, at the time of finding the coefficient α satisfying the formula (q), a value with which the difference in the part D1 is particularly reduced is found.

(78) As described above, the evaluation function shown in the formula (q) includes the weight function We(t), so that the feature (B) described above can be applied to the coefficient α. In addition, a method and a formula for applying the feature (B) to the coefficient α are not limited to the method and the formula described above. As another example, FIG. 13 illustrates an example of the weight function We(t) when Target.sub.0(t) is a multi pulse. A graph G31 in FIG. 13 shows the case where a coefficient C of the formula (r) is 1, and a graph G32 shows the case where the coefficient C of the formula (r) is 2.

(79) Then, the iterative Fourier transform unit 22a performs the processes of the same process numbers (5) to (8) as those in the calculation procedure illustrated in FIG. 9. Subsequently, the above processes (1) to (8) are repeatedly performed a plurality of times, so that a phase spectrum shape represented by the phase spectrum function Ψ.sub.n(ω) in the waveform function can be brought closer to a phase spectrum shape corresponding to the desired temporal intensity waveform. A phase spectrum function Ψ.sub.IFTA(ω) to be finally obtained is provided to the modulation pattern generation unit 24.

(80) FIG. 14 illustrates a calculation procedure in the iterative Fourier transform unit 23a of the intensity spectrum design unit 23. The iterative Fourier transform unit 23a calculates an intensity spectrum by the similar method as the calculation method by the iterative Fourier transform unit 22a described above.

(81) First, the iterative Fourier transform unit 23a prepares an initial intensity spectrum function A.sub.k=0(ω) and a phase spectrum function Ψ.sub.0(ω) in the same manner as when the phase spectrum is calculated (process number (1) in the drawing). Next, the iterative Fourier transform unit 23a prepares a waveform function (s) of a frequency domain including the intensity spectrum function A.sub.k(ω) and the phase spectrum function Ψ.sub.0(ω) (process number (2) in the drawing).
√{square root over (A.sub.k(ω))}exp{iΨ.sub.0(ω}  (s)

(82) A subscript k represents after a k-th Fourier transform process. Before a first Fourier transform process, the initial intensity spectrum function A.sub.k=0(ω) described above is used as the intensity spectrum function A.sub.k(ω). i is an imaginary number.

(83) Next, the iterative Fourier transform unit 23a performs the Fourier transform from the frequency domain to the time domain on the function (s). As a result, a waveform function (t) in the frequency domain including the temporal intensity waveform function b.sub.k(t) is obtained (process number (3) in the drawing).
√{square root over (b.sub.k(t))}exp{iΘ.sub.k(t)}  (t)

(84) Next, the iterative Fourier transform unit 23a obtains a coefficient α having the following features (C) and (D) (process number (4) in the drawing).

(85) (C) A difference (α×Target.sub.0(t)−b.sub.k(t)) between the waveform function b.sub.k(t) after the Fourier transform and the function Target.sub.0(t) multiplied by the coefficient α is smaller than a difference (Target.sub.0(t)−b.sub.k(t)) between the waveform function b.sub.k(t) and the function Target.sub.0(t). Specifically, a time integral value of the difference (α×Target.sub.0(t)−b.sub.k(t)) is smaller than a time integral value of the difference (Target.sub.0(t)−b.sub.k(t)).
(D) At each time of the function Target.sub.0(t), a ratio of the difference (α×Target.sub.0(t)−b.sub.k(t)), that is, a ratio of the difference (α×Target.sub.0(t)−b.sub.k(t)) based on the intensity value of the function Target.sub.0(t) as a reference decreases when the intensity increases.

(86) In one example, the coefficient α with which a pseudo standard deviation σ of α×Target.sub.0(t) for the waveform function b.sub.k(t) after the Fourier transform is minimized (σ.sub.min) is derived by searching using an evaluation function shown by the following formula (u). In addition, in the formula (u), D represents the number of data points, and t.sub.s and t.sub.e represent a start point and an end point of a time axis, respectively. Further, We(t) is a first weight function.

(87) $\begin{array}{cc}{\sigma }_{\min }=\sqrt{\frac{1}{D}\underset{t=t_s}{\overset{t=t_e}{.Math.}}\mathrm{We}\left(t\right)\times {\left(\alpha \times Targe{t}_{k=0}\left(t\right)-b_k\left(t\right)\right)}^2}& \left(u\right)\end{array}$

(88) As shown in the formula (u), the evaluation function is a function including a difference (α×Target.sub.0(t)−b.sub.k(t)) between the waveform function b.sub.k(t) after the Fourier transform and the function α×Target.sub.0(t) after the multiplication, and specifically, the evaluation function includes (α×Target.sub.0(t)−b.sub.k(t)).sup.2. Further, the evaluation function includes a weight function We(t) by which the function is multiplied, and includes a time integral of a result of multiplying the function by the weight function We(t). Further, the coefficient α with which the evaluation function, that is, the time integral is minimized (σ.sub.min) is derived by searching.

(89) In addition, the features and the specific examples of the weight function We(t) are the same as those in the case of calculating the phase spectrum function described above (refer to the formulas (q) and (r)). In addition, the following formula (v) is used instead of the formula (q) described above.

(90) $\begin{array}{cc}{\sigma }_{\min }=\sqrt{\frac{1}{D}\underset{t=t_s}{\overset{t=t_e}{.Math.}}\mathrm{We}\left(t\right)\times {\left(\alpha \times Targe{t}_{k=0}\left(t\right)-b_k\left(t\right)\right)}^2}& \left(v\right)\end{array}$

(91) Next, the iterative Fourier transform unit 23a performs a replacement based on the desired waveform on the temporal intensity waveform function b.sub.k(t) included in the function (t) after the Fourier transform (first replacement). At this time, the iterative Fourier transform unit 23a performs the replacement using a result (α×Target.sub.0(t)) of multiplying the function Target.sub.0(t) representing the desired waveform by the coefficient α. In one example, the function is replaced by Target.sub.k(t) calculated by the formula (w) (process numbers (5) and (6) in the drawing).

(92) $\begin{array}{cc}Targe{t}_k\left(t\right)=\{\begin{array}{c}\alpha \times Targe{t}_0\left(t\right)+\beta \left(\alpha \times Targe{t}_0\left(t\right)-b_k\left(t\right)\right)\\ 0,\mathrm{if}\alpha \times \mathrm{Tar}ge{t}_0\left(t\right)+\beta \left(\alpha \times Targe{t}_0\left(t\right)-b_k\left(t\right)\right)<0\end{array}& \left(w\right)\end{array}$$\begin{array}{cc}\sqrt{Targe{t}_k\left(t\right)}\exp \left\{i{\Theta }_k\left(t\right)\right\}& \left(x\right)\end{array}$

(93) Next, the iterative Fourier transform unit 23a performs an inverse Fourier transform from the time domain to the frequency domain on the function (x). As a result, a waveform function (y) in the frequency domain including the intensity spectrum function C.sub.k(ω) and the phase spectrum function Ψ.sub.k(ω) is obtained (process number (7) in the drawing).
√{square root over (C.sub.k(ω)}exp{iΨ.sub.k(ω)}  (y)

(94) Next, to constrain the phase spectrum function Ψ.sub.k(ω) included in the function (y), the iterative Fourier transform unit 23a replaces the function by the initial phase spectrum function Ψ.sub.0(ω) (second replacement, process number (8) in the drawing).
Ψ.sub.k(ω)Ψ.sub.0(ω)  (z)

(95) Further, the iterative Fourier transform unit 23a performs a filter process based on the intensity spectrum of the input light La on the intensity spectrum function C.sub.k(ω) in the frequency domain after the inverse Fourier transform. Specifically, a portion exceeding a cutoff intensity for each wavelength, which is determined on the basis of the intensity spectrum of the input light La, is cut from the intensity spectrum represented by the intensity spectrum function C.sub.k(ω).

(96) In one example, the cutoff intensity for each wavelength is set to be matched with the intensity spectrum of the input light La (initial intensity spectrum function A.sub.k=0(ω) in the present embodiment). In this case, as shown in the following formula (z1), at a frequency where the intensity spectrum function C.sub.k(ω) is larger than the initial intensity spectrum function A.sub.k=0(ω), a value of the initial intensity spectrum function A.sub.k=0(ω) is taken as a value of the intensity spectrum function A.sub.k(ω). Further, at a frequency where the intensity spectrum function C.sub.k(ω) is equal to or smaller than the initial intensity spectrum function A.sub.k=0(ω), a value of the intensity spectrum function C.sub.k(ω) is taken as the value of the intensity spectrum function A.sub.k(ω).

(97) $\begin{array}{cc}{A}_k\left(\omega \right)=\{\begin{array}{c}{A}_{k=0}\left(\omega \right),A_{k=0}\left(\omega \right)<C_k\left(\omega \right)\\ C_k\left(\omega \right),A_{k=0}\left(\omega \right)\ge C_k\left(\omega \right)\end{array}& \left(\mathrm{z1}\right)\end{array}$

(98) The iterative Fourier transform unit 23a replaces the intensity spectrum function C.sub.k(ω) included in the function (y) by the intensity spectrum function A.sub.k(ω) after the filter process according to the formula (z1). Further, a method for relatively changing the cutoff intensity by defining a function C′.sub.k(ω) obtained by multiplying C.sub.k(ω) by an arbitrary coefficient may be used (process number (9) in the drawing).

(99) Subsequently, the iterative Fourier transform unit 23a repeatedly performs the above processes (1) to (8) (or (1) to (9)) a plurality of times, so that the intensity spectrum shape represented by the intensity spectrum function A.sub.k(ω) in the waveform function can be brought closer to an intensity spectrum shape corresponding to the desired temporal intensity waveform. An intensity spectrum function A.sub.IFTA(ω) to be finally obtained is provided to the modulation pattern generation unit 24.

(100) FIG. 15 is a flowchart illustrating a modulation pattern calculation method realized by the modulation pattern calculation apparatus 20 described above. In addition, the modulation pattern calculation program described above causes the processor 201 (refer to FIG. 6) of the computer to execute each step included in the flowchart.

(101) As illustrated in FIG. 15, first, information on a desired temporal intensity waveform is input to the arbitrary waveform input unit 21 by the operator (input step S1). Next, a phase spectrum and an intensity spectrum to bring the temporal intensity waveform close to a desired waveform are calculated respectively in the phase spectrum design unit 22 and the intensity spectrum design unit 23 (phase spectrum calculation step S11 and intensity spectrum calculation step S21).

(102) The phase spectrum calculation step S11 includes an iterative Fourier transform step S12 by the iterative Fourier transform unit 22a. The details of the iterative Fourier transform step S12 are the same as the operation of the iterative Fourier transform unit 22a described above. The phase spectrum function Ψ.sub.IFTA(ω) to be finally obtained is provided to a subsequent modulation pattern calculation step S3. Further, the intensity spectrum calculation step S21 includes an iterative Fourier transform step S22 by the iterative Fourier transform unit 23a. The details of the iterative Fourier transform step S22 are the same as the operation of the iterative Fourier transform unit 23a described above. The intensity spectrum function A.sub.IFTA(ω) to be finally obtained is provided to the subsequent modulation pattern calculation step S3.

(103) In the modulation pattern calculation step S3 (data generation step), a modulation pattern is calculated on the basis of the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω). The modulation pattern is presented on the SLM 14.

(104) The effects obtained by the modulation pattern calculation apparatus 20, the modulation pattern calculation method, and the modulation pattern calculation program according to the present embodiment described above will be described.

(105) In the present embodiment, the phase spectrum and the intensity spectrum for bringing the temporal intensity waveform of light closer to the desired waveform are calculated in the iterative Fourier transform units 22a and 23a (or the iterative Fourier transform steps S12 and S22), respectively. At this time, when a replacement based on the desired waveform is performed on the temporal intensity waveform function b.sub.n(t) (or b.sub.k(t)) in the time domain obtained by the Fourier transform, a result of multiplying the function Target.sub.0(t) representing the desired waveform by the coefficient α is used. The coefficient α has a value with which a difference between the function α×Target.sub.0(t) after the multiplication of the coefficient α and the temporal intensity waveform function b.sub.n(t) (or b.sub.k(t)) after the Fourier transform decreases, as compared with that before the multiplication of the coefficient α.

(106) As a result, a difference (Target.sub.n(t)−b.sub.n(t) or Target.sub.k(t)−b.sub.k(t)) between the function before the replacement (that is, the temporal intensity waveform function b.sub.n(t) or b.sub.k(t) after the Fourier transform) and the function Target.sub.n(t) (or Target.sub.k(t)) after the replacement based on the desired waveform decreases, so that a solution is prevented from being led to a local solution. Therefore, according to the present embodiment, it is possible to accurately calculate the intensity spectrum and the phase spectrum for bringing the temporal waveform of light closer to the desired waveform.

(107) In particular, when the intensity spectrum function A.sub.IFTA(ω) is derived by the iterative Fourier method, the intensity spectrum after the inverse Fourier transform may exceed the intensity spectrum of the input light La depending on a wavelength region. In such a case, since realization is impossible with respect to a portion exceeding the intensity spectrum of the input light La, an intensity loss occurs with the iterative calculation, and the intensity of the temporal intensity waveform b.sub.n(t) gradually decreases. Therefore, in the conventional iterative Fourier method, a difference between the waveform Target.sub.0(t) after the replacement and the temporal intensity waveform b.sub.n(t) (or b.sub.k(t)) tends to be large, and the risk of being led to a local solution increases.

(108) Meanwhile, according to the present embodiment, since the difference between the waveform function Target.sub.n(t) (or Target.sub.k(t)) after the replacement and the temporal intensity waveform b.sub.n(t) (or b.sub.k(t)) can be reduced, it is possible to reduce the risk of being led to the local solution while suppressing the intensity loss.

(109) Further, in the present embodiment, the coefficient α by which the function Target.sub.0(t) is multiplied has a value with which a ratio of the difference (Target.sub.n(t)−b.sub.n(t) or Target.sub.k(t)−b.sub.k(t)) is smaller when the intensity is higher, at each time of the function Target.sub.0(t) before the multiplication. That is, by multiplication of the coefficient α, in a part with the high intensity in the function Target.sub.0(t), the difference (Target.sub.0(t)−b.sub.n(t) or Target.sub.k(t)−b.sub.k(t)) between the function b.sub.n(t) (or b.sub.k(t)) before the replacement and the function Target.sub.n(t) (or Target.sub.k(t)) after the replacement becomes smaller. Therefore, the phase spectrum and the intensity spectrum can be calculated with higher accuracy. For this reason, the temporal waveform in the part where the light intensity is high particularly in the temporal waveform of the output light Ld can be brought closer to the desired waveform with higher accuracy.

(110) Further, as in the present embodiment, the iterative Fourier transform units 22a and 23a (in the iterative Fourier transform steps S12 and S22) may repeatedly perform the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the modulation pattern generation unit 24 (in the modulation pattern calculation step S3) may calculate the modulation pattern on the basis of the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω) obtained after the repetition. As described above, by repeatedly performing the iterative Fourier transform, the intensity spectrum A.sub.IFTA(ω) and the phase spectrum Ψ.sub.IFTA(ω) for bringing the temporal intensity waveform of light closer to the desired waveform can be calculated with higher accuracy.

(111) Further, as shown in the above formulas (q) and (v), the coefficient α may have a value with which a time integral value of a result of multiplying a function including the difference between the temporal intensity waveform function b.sub.n(t) (or b.sub.k(t)) after the Fourier transform and the function α×Target.sub.0(t) by the weight function We(t) having a larger weight value when the intensity is higher at each time of the function Target.sub.0(t) is minimized. As a result, the difference between the function α×Target.sub.0(t) and the temporal intensity waveform function b.sub.n(t) (or b.sub.k(t)) after the Fourier transform in the part where the intensity is high in the desired temporal intensity waveform can be selectively reduced, and the temporal waveform in the part with the high light intensity can be brought closer to the desired waveform with higher accuracy.

(112) Further, as shown in the formula (r) described above, the weight function We(t) may include a function obtained by multiplying the function Target.sub.0(t) by another coefficient C. For example, with such a configuration, it is possible to appropriately determine the weight function We(t) having a larger weight value when the intensity is higher at each time of the function Target.sub.0(t).

(113) Further, the light control apparatus 1A of the present embodiment includes the light source 2 that outputs the input light La, the diffraction grating 12 that functions as the dispersive element spectrally dispersing the input light La, the SLM 14 that modulates the light Lb after the dispersion and outputs modulated light, and the lens 15 that functions as the optical system focusing the modulated light Lc. Further, the SLM 14 modulates the input light La on the basis of the data created by the modulation pattern calculation apparatus 20. According to the light control apparatus 1A, it is possible to calculate an intensity spectrum and a phase spectrum to bring a temporal waveform in a part with a high light intensity particularly in a temporal waveform of the output light Ld closer to a desired waveform with higher accuracy while preventing a solution from being led to a local solution at the time of the iterative Fourier operation.

(114) In addition, in the present embodiment, to bring the temporal intensity waveform closer to the desired waveform, both the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω) are calculated, and the modulation pattern based on these functions is presented on the SLM 14. However, it is not limited to such a configuration, and for example, only one of the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω) for bringing the temporal intensity waveform closer to the desired waveform may be calculated. In this case, a spectrum prepared (or selected) in advance may be used as the other spectrum, or the other spectrum may not be modulated with the input light La. Even in any configuration, the effects of the present embodiment can be suitably obtained.

(115) FIG. 16 and FIG. 17 are graphs illustrating output temporal waveforms when the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω) are calculated using a pulse waveform with a time width of about 135 femtoseconds as the input light La and a rectangular pulse waveform as the desired waveform. FIG. 16 illustrates the case of calculating the functions using the method (refer to FIG. 9) described in Patent Document 1, and FIG. 17 illustrates the case of calculating the functions using the method (refer to FIG. 11 and FIG. 14) of the present embodiment. Further, in these drawings, graphs G41 to G43 show cases where the time width of the rectangular pulse is set as 6 picoseconds, 7 picoseconds, and 9 picoseconds, respectively.

(116) When FIG. 16 and FIG. 17 are compared, it can be seen that a waveform of the present embodiment is smoother in a flat part (top part of the rectangular pulse) where the intensity of the rectangular pulse is high, even in any time width. Further, a standard deviation of the temporal waveform illustrated in FIG. 17 is 1.2 times to a maximum of about 3.4 times smaller than that of the temporal waveform illustrated in FIG. 16. As described above, according to the present embodiment, the temporal waveform in the part with the high light intensity particularly in the temporal waveform of the output light Ld can be brought closer to the desired waveform with higher accuracy.

(117) Further, a pulse waveform having the time width of about 135 femtoseconds is set as the input light La and a multi pulse waveform having a uniform peak value of each pulse is set as a desired waveform, and the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω) are calculated using the method of the present embodiment (refer to FIG. 11 and FIG. 14) and the method described in Patent Document 1. Here, the standard deviation of the peak value of each pulse is adopted as an evaluation value indicating the uniformity.

(118) FIG. 18 is a graph illustrating a result of improving the uniformity of peak values by the calculated phase spectrum function Ψ.sub.IFTA(ω) and intensity spectrum function A.sub.IFTA(ω). A horizontal axis represents the number of pulses. A vertical axis represents a ratio (σ2/σ1) between a standard deviation σ1 according to the method of the present embodiment and a standard deviation σ2 according to the method described in Patent Document 1 in order to show a degree of improvement of the uniformity. A value of the ratio exceeds 1 even at any number of pulses, and it is confirmed that the uniformity of the peak values of the respective pulses is improved by the method of the present embodiment.

(119) Further, FIG. 19 is a graph illustrating an output temporal waveform when a pulse waveform with a time width of about 135 femtoseconds is set as the input light La, a time shape having two parts Da and Dc of a certain constant intensity and a part Db of a constant intensity located between the two parts Da and Dc and having an intensity higher than that of the parts Da and Dc is set as the desired waveform, and the phase spectrum function Ψ.sub.IFTA(ω) and the intensity spectrum function A.sub.IFTA(ω) are calculated using the method of the present embodiment (refer to FIG. 11 and FIG. 14). In the graph, a standard deviation of an intensity variation in the part Db is improved by about 2.2 times as compared with a standard deviation of an intensity variation in the parts Da and Dc.

(120) FIG. 20 and FIG. 21 are diagrams illustrating a calculation procedure of the phase spectrum and the intensity spectrum by the iterative Fourier method in accordance with a modification of the above embodiment. A difference between the calculation procedure and the calculation procedure of the above embodiment (refer to FIG. 11 and FIG. 14) is that the coefficient β in the process number (5) is replaced by the weight function Wr(t). In the present modification, the formulas (m) and (w) in the above embodiment are replaced by the following formulas (z2) and (z3), respectively.

(121) $\begin{array}{cc}\mathrm{Targ}e{t}_n\left(t\right)=\{\begin{array}{c}\alpha \times \mathrm{Ta}rge{t}_0\left(t\right)+Wr\left(t\right).Math.\left(\alpha \times {\mathrm{Target}}_0\left(t\right)-b_n\left(t\right)\right)\\ 0,\mathrm{if}\alpha \times {\mathrm{Target}}_0\left(t\right)+Wr\left(t\right).Math.\left(\alpha \times {\mathrm{Target}}_0\left(t\right)-b_n\left(t\right)\right)<0\end{array}& \left(\mathrm{z2}\right)\end{array}$$\begin{array}{cc}\mathrm{Targe}{t}_k\left(t\right)=\{\begin{array}{c}\alpha \times Targe{t}_0\left(t\right)+Wr\left(t\right).Math.\left(\alpha \times {\mathrm{Target}}_0\left(t\right)-b_k\left(t\right)\right)\\ 0,\mathrm{if}\alpha \times \mathrm{Tar}ge{t}_0\left(t\right)+Wr\left(t\right).Math.\left(\alpha \times {\mathrm{Target}}_0\left(t\right)-b_k\left(t\right)\right)<0\end{array}& \left(\mathrm{z3}\right)\end{array}$

(122) That is, in the present modification, the first replacement is performed using the sum of a result of multiplying a function obtained by subtracting the temporal intensity waveform function b.sub.n(t) (or b.sub.k(t)) after the Fourier transform from the function {α×Target.sub.0(t)} by the weight function Wr(t) and the function {α×Target.sub.0(t)}. In addition, the weight function Wr(t) is an example of the second weight function in the present modification.

(123) The weight function Wr(t) is a function having a larger weight value when the intensity is higher, at each time of the function Target.sub.0(t). In one example, the weight function Wr(t) includes a function obtained by multiplying the function Target.sub.0(t) by another coefficient C.sub.2, and is represented by, for example, the following formula (z4). In other words, the weight function Wr(t) may be determined on the basis of the function Target.sub.0(t).
Wr(t)=C.sub.2.Math.Target.sub.0(t)  (z4)

(124) FIG. 22A to FIG. 22C and FIG. 23A to FIG. 23C are diagrams for visually explaining the process of the process number (5) in the above embodiment and the present modification, respectively.

(125) FIG. 22A and FIG. 23A illustrate an example of a calculation result of {α×Target.sub.0(t)}. FIG. 22B illustrates an example of a result of multiplying {α×Target.sub.0(t)−b.sub.n(t)} by the coefficient β (refer to the formulas (m) and (w)) of the embodiment, and FIG. 23B illustrates an example of a result of multiplying {α×Target.sub.0(t)−b.sub.n(t)} by the weight function Wr(t). Further, FIG. 22C illustrates a graph of the sum (that is, Target.sub.n(t) of the formula (m)) of FIG. 22A and FIG. 22B, and FIG. 23C illustrates a graph of the sum (that is, Target.sub.n(t) of the formula (z2)) of FIG. 23A and FIG. 23B.

(126) As is clear from comparison of FIG. 22C and FIG. 23C, by replacing the coefficient β by the weight function Wr(t), the magnitude of a difference is emphasized in a part D1 where the intensity is high in Target.sub.0(t), as compared with the other parts. Therefore, at the time of the iterative Fourier operation, a result that particularly reduces the difference in the part D1 is calculated. Thus, the temporal waveform in the part where the light intensity is high particularly in the temporal waveform of the output light Ld can be brought closer to the desired waveform with higher accuracy.

(127) The data creation apparatus, the light control apparatus, the data creation method, and the data creation program are not limited to the embodiments and the modifications described above, and a variety of other modifications can be made.

(128) For example, in the above embodiment, the case where the coefficient α and the weight function We(t) are applied to IFTA-Fienup (second comparative example) is exemplified, however, the present disclosure is not limited thereto. For example, the coefficient α and the weight function We(t) may be applied to the normal iterative Fourier method according to the first comparative example.

(129) Further, in the above embodiment, the coefficient α with which a difference between the waveform function after the Fourier transform and a result of multiplying the function representing the desired waveform by the coefficient α decreases is obtained by minimizing the standard deviation σ. However, a method for calculating the coefficient α is not limited thereto. For example, the total sum (that is, the time integral value) of the difference values between the waveform function after the Fourier transform and the result of multiplying the function representing the desired waveform by the coefficient α may be minimized.

(130) The data creation apparatus of the above embodiment is an apparatus for creating data for controlling a spatial light modulator, and is configured to include a phase spectrum design unit for generating a phase spectrum function to be used for creating the data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; and a data generation unit for creating the data on the basis of the intensity spectrum function and the phase spectrum function generated by the phase spectrum design unit, and the phase spectrum design unit performs the first replacement using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, and the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

(131) The data creation method of the above embodiment is a method for creating data for controlling a spatial light modulator, and is configured to include a phase spectrum calculation step of generating a phase spectrum function to be used for creating the data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; and a data generation step of creating the data on the basis of the intensity spectrum function and the phase spectrum function generated in the phase spectrum calculation step, and in the phase spectrum calculation step, the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, and the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

(132) The data creation program of the above embodiment is a program for creating data for controlling a spatial light modulator, and is configured to cause a computer to execute a phase spectrum calculation step of generating a phase spectrum function to be used for creating the data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the intensity spectrum function in the frequency domain after the inverse Fourier transform; and a data generation step of creating the data on the basis of the intensity spectrum function and the phase spectrum function generated in the phase spectrum calculation step, and in the phase spectrum calculation step, the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, and the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

(133) The data creation apparatus of the above embodiment is an apparatus for creating data for controlling a spatial light modulator, and is configured to include an intensity spectrum design unit for generating an intensity spectrum function to be used for creating the data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the phase spectrum function in the frequency domain after the inverse Fourier transform; and a data generation unit for creating the data on the basis of the phase spectrum function and the intensity spectrum function generated by the intensity spectrum design unit, and the intensity spectrum design unit performs the first replacement using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, and the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

(134) The data creation method of the above embodiment is a method for creating data for controlling a spatial light modulator, and is configured to include an intensity spectrum calculation step of generating an intensity spectrum function to be used for creating the data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the phase spectrum function in the frequency domain after the inverse Fourier transform; and a data generation step of creating the data on the basis of the phase spectrum function and the intensity spectrum function generated in the intensity spectrum calculation step, and in the intensity spectrum calculation step, the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, and the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

(135) The data creation program of the above embodiment is a program for creating data for controlling a spatial light modulator, and is configured to cause a computer to execute an intensity spectrum calculation step of generating an intensity spectrum function to be used for creating the data by performing a Fourier transform on a waveform function in a frequency domain including an intensity spectrum function and a phase spectrum function, performing a first replacement of a temporal intensity waveform function based on a desired temporal intensity waveform in a time domain after the Fourier transform and then performing an inverse Fourier transform, and performing a second replacement to constrain the phase spectrum function in the frequency domain after the inverse Fourier transform; and a data generation step of creating the data on the basis of the phase spectrum function and the intensity spectrum function generated in the intensity spectrum calculation step, and in the intensity spectrum calculation step, the first replacement is performed using a result of multiplying a function representing the desired temporal intensity waveform by a coefficient, and the coefficient has a value with which a difference between the function after the multiplication of the coefficient and the temporal intensity waveform function after the Fourier transform is smaller than a difference before the multiplication, and a ratio of the difference is smaller when an intensity is higher at each time of the function before the multiplication.

(136) In the above configuration, the phase spectrum design unit may repeatedly perform the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data generation unit may create the data on the basis of the phase spectrum function obtained after the repetition. Further, in the above configuration, the intensity spectrum design unit may repeatedly perform the Fourier transform, the first replacement, the inverse Fourier transform, and the second replacement, and the data generation unit may create the data on the basis of the intensity spectrum function obtained after the repetition.

(137) As described above, by repeatedly performing the iterative Fourier transform, the intensity spectrum or the phase spectrum for bringing the temporal intensity waveform of light closer to the desired waveform can be calculated with higher accuracy.

(138) In the above configuration, the coefficient may have a value with which a time integral value of a result of multiplying a function, including a difference between the temporal intensity waveform function after the Fourier transform and the function after the multiplication, by a first weight function having a larger weight value when the intensity is higher at each time of the function before the multiplication, is minimized. As a result, the difference between the function after the multiplication and the temporal intensity waveform function after the Fourier transform in the part where the intensity is high in the desired temporal intensity waveform can be selectively reduced, and the temporal waveform in the part with the high light intensity can be brought closer to the desired waveform with higher accuracy.

(139) In this case, the first weight function may include a function obtained by multiplying the function before the multiplication by another coefficient. For example, with the above configuration, it is possible to appropriately determine the first weight function having a larger weight value when the intensity is higher at each time of the function before the multiplication.

(140) In the above configuration, the first replacement may be performed using a sum of a result of multiplying a function, obtained by subtracting the temporal intensity waveform function after the Fourier transform from the function after the multiplication, by a second weight function, and the function after the multiplication, and the second weight function may have a larger weight value when the intensity is higher at each time of the function before the multiplication. As a result, the temporal waveform in the part where the light intensity is high particularly in the temporal waveform of the output light can be brought closer to the desired waveform with higher accuracy.

(141) The light control apparatus of the above embodiment is configured to include a light source for outputting input light; a dispersive element for spectrally dispersing the input light; a spatial light modulator for modulating the input light after the dispersion and outputting modulated light; and an optical system for focusing the modulated light, and the spatial light modulator modulates the input light on the basis of the data created by the data creation apparatus having the above configuration.

(142) According to the above light control apparatus, the input light is modulated on the basis of the data created by the data creation apparatus of the above configuration, so that it is possible to calculate an intensity spectrum or a phase spectrum to bring a temporal waveform in a part with a high light intensity particularly in a temporal waveform of light closer to a desired waveform with higher accuracy while preventing a solution from being led to a local solution at the time of an iterative Fourier operation.

(143) An embodiment can be used as a data creation apparatus, a light control apparatus, a data creation method, and a data creation program, which are capable of calculating an intensity spectrum or a phase spectrum to bring a temporal waveform in a part with a high light intensity particularly in a temporal waveform of light closer to a desired waveform with higher accuracy while preventing a solution from being led to a local solution at the time of an iterative Fourier operation.

(144) From the invention thus described, it will be obvious that the invention may be varied in many ways. Such variations are not to be regarded as a departure from the spirit and scope of the invention, and all such modifications as would be obvious to one skilled in the art are intended for inclusion within the scope of the following claims.