JP2000227450A - Electro-magnetic field analyzer, electro-magnetic field analyzing method, and computer-readable recording medium of recorded program for executing it by computer - Google Patents

Abstract

PROBLEM TO BE SOLVED: To reduce a genegated FDTD(Finite Difference Time Domain) error to conduct highly precise electro-magnetic field analysis, by using an FDTD method. SOLUTION: This electro-magnetic field analyzer for analyzing a propagating condition of an electro-magnetic field using an FDTD method is provided with an incident wave information inputting part 101 for inputting information as to an incident wave, an incident wave calculating part 102 for calculating the incident wave based on the information input by the information inputting part 101, an electromagnetic field calculating part 105 for calculating the electro- magnetic field based on the incident wave calculated by the incident wave calculating part 102, and an electro-magnetic field calculation-outputting part 106 for outputting a result calculated by the electro-magnetic field calculating part 105.

Description

DETAILED DESCRIPTION OF THE INVENTION

[0001]

[0001] The present invention relates to an FDTD (Fi
Nite Difference Time Doma
The present invention relates to an electromagnetic field analysis device that analyzes the propagation of an electromagnetic field using an in (finite difference time domain) method, an electromagnetic field analysis method, and a computer-readable recording medium that records a program for causing a computer to execute the method.

[0002]

2. Description of the Related Art If it is possible to logically predict the state in which an electromagnetic field such as light or radio waves propagates in the air, and is reflected or refracted and transmitted through a material surface, an antenna, an optical component, a laser device, Guidance can be given to the design of microscopes, integrated circuits, etc. Since it is known that the electromagnetic field follows Maxwell's equation, necessary information can be obtained by solving Maxwell's equation of a given system.

In the FDTD method, the way a wave travels from time to time is calculated from the arrangement of matter in a three-dimensional space and the wavelength and direction of an incident wave by numerically solving Maxwell's equations using a computer. How to here,
One of the difficult problems when using the FDTD method is a problem of inputting an incident wave.

The simplest method of inputting an incident wave is a method called a hard source. This is because the electromagnetic field in the entire analysis area is set to “0 (zero)” as an initial state, and the electric field or magnetic field on the wave incident surface is forcibly changed so that the electromagnetic field propagates in the analysis area. Is what you do.

The term "forcibly changing" means that a variable representing the electromagnetic field on the incident surface is overwritten with the value of the incident wave during the execution of the program. This is the same as the operation of generating a wave by waving the end of the rope of a skipping rope that is stretched horizontally, as shown in FIG. 10, and only the end point is forcibly controlled. After that, the waves naturally occur and propagate.

A problem that can be considered in such a hard source method is that if a substance that reflects a wave (light) is present in the analysis area, the reflected wave arrives at the incident surface. Then, the sum of the incident wave and the reflected wave should be the electromagnetic field on the incident surface, but this is forcibly overwritten by the incident wave. This state is shown in FIG.

[0007] As a result, the incident surface and the surface inside one mesh are not aligned, and a large spatial inclination occurs. This causes the electromagnetic field at the next time to be disturbed. This is repeatedly performed, and propagates one after another in the analysis area, and eventually the whole is filled with completely different (random) values.

In order to solve the above-mentioned problems in the hard source method, a “global / scattering field” method is used. FIG. 12 shows an outline of the “whole field / scattering field” method.
In FIG. 12, another area is provided in the analysis area, and this area is defined as an entire field 1201. In addition, the whole venue 12
The outside of 01 is a scattered field 1202.

Although both the incident wave and the scattered wave exist in the region of the whole field 1201, the interface 1 between the whole field 1201 and the scattered field exists.
At 203, the interface 1203 becomes the entrance surface by subtracting the contribution from the incident wave component. Although an electric field E is defined at the interface 1203 in FIG. 12, this may be a magnetic field.

The magnetic field H ₁ should be updated using the electric field on both sides thereof, that is, the difference between E ₁ and E _2. Here, in order to make H ₁ only a scattered wave component, the incident wave from E ₂ must be updated. Update H ₁ with the remaining E ₂ ′ minus the component instead. H ₁
Is updated using the electric field (E ₁ and E ₂ ′) of only the scattered wave component, and therefore has only the scattered wave component.

[0011] Conversely, to update the E ₂ using H ₁ 'plus the incident wave components instead of H _1. E ₂ is the whole field (incident wave +
Since it is updated using the magnetic field (H ₁ and H ₁ ′) of the scattered wave, the electric field becomes the whole field. By doing so, as a result, the incident wave can be made to enter the whole field, and the scattered wave can be transmitted through the incident surface and exit to the outside.

As described above, in the “whole field / scattering field” method, since the whole field is included in the scattering field and the incident wave component is added or subtracted at the interface, the value of the incident wave at the interface (the electric field and the Magnetic field) must be known. Usually, a simple wave such as a sine wave or a pulse wave is used as the incident wave.

FIG. 13 is an explanatory diagram showing an ideal calculation result when a state in which a plane wave propagates in a system in which no scatterer exists, that is, in a vacuum, is calculated using the “whole field / scattering field” method. It is. FIG. 13A is an explanatory diagram showing a model of the whole field and the scattering field, and FIG.
FIG. 13A shows the electromagnetic field distribution between “A” and “B” in FIG. 13A, and FIG. 13C shows the propagation direction, that is, between “C” and “D” in FIG. 3 shows an electromagnetic field distribution.

By definition, a plane wave is uniform in the plane of incidence and is a sine wave in the direction of propagation. Since there is no scatterer, there is no scattered wave, and the wave intensity in the scattered field is “0 (zero)”.

FIG. 14 is an explanatory diagram showing a calculation result obtained by actually calculating plane wave propagation by the FDTD method. In FIG. 14, it can be seen that light propagates in the z direction and reaches z = 500. FIG. 15 shows this in contour display. This makes it possible to see errors that were not clear in FIG. In the scattered field, in this case, the wave is "0
(Zero) ", but the error is causing a wave.

FIG. 16 is an explanatory diagram showing a cross-sectional view of FIG. 15 cut in the x direction at the point z = 350. As can be seen from FIG. 16, not only the scattering field is not “0 (zero)”, but also the waves of the entire field are not uniform. Thus, in the “global / scattered field” method,
It can be seen that the incident surface does not interfere with the scattered waves, but the error once generated is not swept out.

In the FDTD method, the electric field and the magnetic field are updated every time step. The electric field will be updated from the difference between the magnetic fields on both sides, and the magnetic field will be updated from the difference between the electric fields on both sides. Therefore, as a result, the electric field and the magnetic field satisfy the Maxwell equation, and the light is propagated as a wave at the speed of light.

In order for the calculation not to diverge, the size of the time step (step time) must be smaller than a certain value (Courant number). In this time step, light travels about 0.5 mesh in one step.

[0019]

However, in the above conventional method, as shown in FIG. 17, after the incident wave is incident from the incident surface, the electric field and the magnetic field are updated every time step. Changes in one time step
The mesh will propagate.

For this reason, a non-zero value is generated in the electromagnetic field before the arrival of the incident wave, and propagates forward at a speed higher than the speed of light. In principle, as shown in FIG. 18, the speed of the wave itself (for example, the moving speed of the wave peak)
Even if is equal to the speed of light, the tip of the wave is distorted and this error propagates.

In FIG. 18, the FDTD calculation value has reached a region where the sine wave as a logical value has not yet reached, and the whole field / scattering field interface at “E”-“F” in FIG. There is a problem that the operation shown in FIG. 12 cannot be performed even though a wave already exists in the field.

This is an unavoidable error that occurs from the first step of the calculation, does not attenuate during the calculation, and is included as an error component in the final calculation result. This error is about 10 compared to the incident wave amplitude.
^Since it is about ^-3, it is a fatal problem when dealing with weak light whose incident light is only about 10 ^{-3 of the} incident light as in near-field light analysis.

The present invention solves the above-mentioned problems of the conventional method, reduces the FDTD error generated by using the FDTD method, and performs an electromagnetic field analysis with high accuracy and an electromagnetic field analysis method. It is another object of the present invention to provide a computer-readable recording medium on which a program for causing a computer to execute the method is recorded.

[0024]

Means for Solving the Problems The above-mentioned problems are solved,
In order to achieve the above object, an electromagnetic field analysis apparatus according to the first aspect of the present invention provides an FDTD (Finite Difference).
In an electromagnetic field analyzer that analyzes a propagation state of an electromagnetic field using an e-time domain method, an input unit that inputs information regarding an incident wave, and an incident wave calculation that calculates the incident wave based on the information input by the input unit Means, an electromagnetic field calculating means for calculating an electromagnetic field based on the incident wave calculated by the incident wave calculating means, and an output means for outputting a result calculated by the electromagnetic field calculating means. .

According to the first aspect of the present invention, the incident wave calculation means accurately reproduces the error of the incident wave, and the value including the error is used as the incident wave, so that the electromagnetic field calculation means cancels the FDTD error. be able to.

Further, the electromagnetic field analysis method according to the second aspect of the present invention provides an FDTD (Finite Difference).
In an electromagnetic field analysis method for analyzing a propagation state of an electromagnetic field using a Time Domain method, a first step of inputting information on an incident wave and a second step of calculating an incident wave based on the information input in the first step A step of calculating an electromagnetic field based on the incident wave calculated in the first step, and a fourth step of outputting a result calculated in the third step. .

According to the second aspect of the present invention, the error of the incident wave is accurately reproduced in the second step, and the value including the error is used as the incident wave, thereby canceling the FDTD error in the third step. Can be.

Further, the storage medium according to the invention of claim 3 is
By recording a program for causing a computer to execute the method described in claim 2, the program can be machine-readable, whereby the operation of claim 2 can be realized by a computer.

[0029]

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS Referring to the accompanying drawings, preferred embodiments of an electromagnetic field analyzing apparatus, an electromagnetic field analyzing method and a computer-readable recording medium for recording a program for causing a computer to execute the method will be described below. The form will be described in detail.

First, the functional configuration of the electromagnetic field analyzer according to the present embodiment will be described. FIG. 1 is a block diagram functionally showing the configuration of the electromagnetic field analyzer according to the present embodiment. In FIG. 1, the electromagnetic field analysis apparatus includes an incident wave information input unit 101, an incident wave calculation unit 102, a structure information input unit 103, a calculation parameter input unit 104, an electromagnetic field calculation unit 105, and an electromagnetic field calculation output unit 106. It is a configuration including.

The incident wave information input unit 101 inputs information on an incident wave. Specifically, the input device 20 described later
4 is used to input data. The information on the incident wave includes information on the wavelength and direction of the incident wave.

The incident wave calculator 102 calculates the incident wave based on the information on the incident wave input from the incident wave information input unit 101. For example, a sine wave is input as an incident wave, and the input sine wave is
Calculate using the method.

The structure information input unit 103 inputs structure information. Specifically, data is input by the input device 204 described later. As the structural information, there is information on the arrangement of substances in a three-dimensional space.

The calculation parameter input unit 104 inputs calculation parameters which are parameters set when calculating an electromagnetic field. Specifically, data is input by the input device 204 described later. The calculation parameters are arbitrarily set in accordance with the incident wave information and the structure information to derive an optimum calculation result.

The electromagnetic field calculation unit 105 is calculated by the incident wave calculation unit 102 based on the structure information input by the structure information input unit 103 and the calculation parameters input by the calculation parameter input unit 104 using the FDTD method. Calculate the electromagnetic field for the sinusoidal wave.

The electromagnetic field calculation output section 106 outputs the result calculated by the electromagnetic field calculation section 105. In particular,
The calculated result may be displayed as a numerical value or a graph on the display device 205 described later, or may be printed on a printer.

The incident wave calculation unit 102 and the electromagnetic field calculation unit 105 are provided in a memory 201 or a recording medium 203 described later.
The function of each unit is realized by the arithmetic device 202 executing the instruction processing according to the instructions described in the program recorded in the program.

Next, the hardware configuration of the electromagnetic field analyzer according to the embodiment of the present invention will be described. FIG. 2 is an explanatory diagram illustrating a hardware configuration of the electromagnetic field analysis device according to the present embodiment.

In FIG. 2, the electromagnetic field analyzing apparatus according to the present embodiment has a memory 201 storing a program or data, an arithmetic unit 202 for controlling the whole or performing processing such as calculation, and a program or data. A removable recording medium (read / write device) 203,
An input device 204 for inputting data and a display device 205 for displaying calculation results and the like are shown. Also, 200
Indicates a bus for connecting the above components.

Next, the principle of the FDTD method will be described. As described above, the FDTD method numerically solves Maxwell's equations using a computer to determine the arrangement of matter in a three-dimensional space and the wavelength and direction of incident waves.
This is a method of calculating the way a wave travels at an hourly time.

The Maxwell equation has the following form.

(Equation 1) Here, μ indicates magnetic permeability, ε indicates dielectric constant, and σ indicates electrical conductivity. E indicates an electric field, H indicates a magnetic field,
∇ × indicates the rotation of the vector.

It is important to note that the Maxwell equation is:

(Equation 2) The point is that. That is, if the spatial distribution of electric and magnetic fields at a certain time is known, its spatial change (slope)
From this we can see how the electric and magnetic fields change for the next time. Values such as magnetic permeability can be obtained if the substance existing at that location is known. That is, if it is determined whether or not the substance is arranged in the three-dimensional space, the coefficients A, B, and C can be calculated.

FIG. 3 is an explanatory diagram showing a spatial arrangement of an electric field and a magnetic field. In the FDTD method, as shown in FIG. 3, the analysis target space is divided into fine meshes, and an electric field and a magnetic field are arranged alternately spatially. Calculates the temporal change of the magnetic field (or electric field) sandwiched between.

FIG. 4 is a flowchart of a program for obtaining the above temporal change by using the FDTD method. (FIG. 3
4 and FIG. 4 are quoted from "Electromagnetic Field and Antenna Analysis by FDTD Method" (Toru Uno, Corona, ISBN 4-339-00689-0). In the flowchart of FIG. 4, first, parameters, scatterers, wave sources, and the like are initialized (step S401). Next, an electric field at T = 0 is calculated (step S402), and an absorption boundary condition is set (step S403).

Subsequently, the calculation of the magnetic field at T = T + Δt / 2 is performed (step S403). Next, in step S405, a steady state, that is, T ≧ T
it is determined whether the _max, if not T ≧ T _max (step S405 negative), the process proceeds to step S402, further subjected to an electric field calculation of the T = T + Δt / 2 (step S402), and the same Repeat the process.

If it is determined in step S405 that T ≧ _Tmax (Yes in step S405), the calculation result is output (step S406), thereby ending a series of processes.

In other words, first, the coefficients A, B, and C are obtained, the initial states of the electric field and the magnetic field are input, and thereafter, the electric field is updated using the gradient of the magnetic field, and the new gradient of the electric field is calculated. The process of updating the magnetic field by using the same is repeated, and when the electromagnetic field in the analysis region is in a steady state, all the processes are completed.

FIG. 5 is an explanatory view schematically showing a one-dimensional FDTD. In FIG. 5, the electric field is indicated by E and the magnetic field is H
Indicated by The superscript indicates the time, and the subscript indicates the position. It is assumed that the electric field E and the magnetic field H are known at all times at a certain time (time = 0).

Based on the gradient of the electric field E, the amount of change in the magnetic field H at the next time can be calculated. For example, the magnetic field H ₁ is the amount of change from the difference to the next time the electric field E ₁ and E ₂ of its spatial neighboring be calculated. The same applies to H ₂ .

Next, the amount of change in the electric field E to the next time can be calculated from the difference between the magnetic fields H on both sides thereof. By repeating this operation, the time change of the electromagnetic field can be obtained every moment.

Next, the principle of the electromagnetic field analyzer according to the present embodiment will be described. In the electromagnetic field analysis apparatus according to the present embodiment, when the incident wave electromagnetic field is used in the “whole field / scattered field” method, not the sine wave itself, but light propagating by the FDTD method when a sine wave is inserted is used. .

That is, as shown in FIG. 6, a virtual one-dimensional FDTD region 603 is provided outside the analysis region (entire field) 601 (for example, in the scattering field 602), and no scatterer exists there. , So that the incident light simply propagates.

As described above, even if a sine wave is incident, the FD
Although the tip of the wave is distorted due to the characteristics of the TD method, the virtual one-dimensional FDTD region 603 accurately reproduces the error. By using the value including this error as the incident wave, the FDTD error in the analysis region of the whole field 601 is canceled.

FIG. 7 is an explanatory diagram showing the results of plane wave propagation calculation by the electromagnetic field analyzer according to the present embodiment. Although the wave distribution at the same time as FIG. 7 shown in the description of the prior art is shown,
It can be seen that the error has been reduced to such an extent that the error cannot be seen. Also, when FIG. 8 is similarly compared with FIG. 16, it can be seen that it is clearly improved. The error is the accuracy of the calculator 1
It decreased to about 0-13.

Next, the procedure of the electromagnetic field analysis processing will be described. FIG. 9 is a flowchart showing a processing procedure of an analysis process of the electromagnetic field analysis device according to the present embodiment. In the flowchart of FIG. 9, first, structure information, incident wave information, and calculation parameters are input (step S90).
1).

Next, the magnetic field (whole field and scattered field) is updated (step S902), the incident wave magnetic field (virtual area) is updated (step S903), and the magnetic field (whole field / scattered field interface) is updated (step S903). S904) is performed respectively.

Next, the time is updated (step S90).
5) Then, update of the electric field (whole field and scattered field) (step S906), update of the electric field (virtual area) of the incident wave (step S907), and update of the magnetic field (whole field / scattered field interface) (step S906) S908) is performed.

Further, the time is updated (step S90).
9) After that, it is determined whether or not it is in a steady state (step S910). Here, when it is not in the steady state (No at Step S910), Step S902 is performed.
The process proceeds to steps S902 to S910.
Is repeated until a steady state is reached.

If a steady state is reached in step S910 (Yes in step S910), the calculated electromagnetic field is output (step S911), and all the processing ends.

As described above, according to the present embodiment, instead of the theoretical value of the incident wave necessary to give the incident wave by the FDTD method, the FDTD in the virtual FDTD region is used.
By using the value calculated by the D method, the error can be reduced by about 10
Digit reduction can be achieved. This is extremely important in near-field light analysis or the like that requires high-precision calculations.

As an application field of the electromagnetic field analysis according to the present embodiment, for example, the design of an optical waveguide or an optical fiber can be considered. More specifically, regarding the input,
These are the structure and material inside the waveguide and the fiber, or the wavelength and intensity of the light used. As for the output, it is a clue to know the state of light propagation in the waveguide, the fiber, and the intensity at the output section.

As another application field, an antenna,
The design of a personal wireless communication device such as a mobile phone, a car navigation system, and the like can be considered. More specifically, the input refers to the structure and material of the antenna or the wavelength and intensity of the electromagnetic wave used. As for the output, it is a clue to know the intensity and direction of the electromagnetic wave transmitted or received from the antenna.

As another field of application, medical treatment, in particular, the design of a heat treatment system for cancer and the like can be considered. More specifically, the input is the tissue distribution in the patient or the wavelength, intensity, and direction of incident light. As for the output, it is a clue to know the light intensity distribution in the patient.

Further, other fields of application include biological research,
In particular, research on the mechanism of the human retina is considered. More specifically, the input is the tissue structure of the retina or the wavelength, intensity, and direction of incident light. As for the output, it is a clue to know the state of light propagation inside the retina.

As another application field, the design of an integrated circuit and the like can be considered. More specifically, the input is the structure of the multilayer integrated circuit board or the voltage and frequency used. As for the output, it is a clue to know how the electromagnetic field propagates between the boards.

As another application field, the design of a fighter or the like can be considered. More specifically, regarding the input,
The shape and material of the fighter, or the radar intensity and frequency. As for the output, it is a clue to know the reflection propagation intensity and direction.

The electromagnetic field analysis method described in the present embodiment is realized by executing a prepared program on a computer such as a personal computer or a workstation. This program includes hard disk, floppy disk, CD-ROM, M
The program is recorded on a computer-readable recording medium such as O and DVD, and is executed by being read from the recording medium by the computer. This program can be distributed via the recording medium and a network such as the Internet.

[0068]

As described above, according to the first aspect of the present invention, when analyzing the propagation state of an electromagnetic field by using the FDTD method, input means for inputting information relating to an incident wave, and input information by the input means. Incident wave calculating means for calculating an incident wave based on the information obtained, an electromagnetic field calculating means for calculating an electromagnetic field based on the incident wave calculated by the incident wave calculating means, and a result calculated by the electromagnetic field calculating means. Output means for outputting, the incident wave calculation means accurately reproduces the error of the incident wave, and the value including the error is used as the incident wave, so that the electromagnetic field calculation means cancels the FDTD error. Since the FDTD method can be used, an FDTD error generated by using the FDTD method can be reduced, and an electromagnetic field analyzer capable of performing a highly accurate electromagnetic field analysis can be obtained.

According to the second aspect of the present invention, FDTD
When analyzing the propagation state of the electromagnetic field using the method, a first step of inputting information regarding the incident wave, and a second step of calculating the incident wave based on the information input in the first step,
A third step of calculating an electromagnetic field based on the incident wave calculated in the first step; and a fourth step of outputting a result calculated in the third step. Is accurately reproduced and the value including the above error is used as the incident wave, so that the FDTD error can be canceled by the third step. Therefore, the FDTD error generated by using the FDTD method is reduced, There is an effect that an electromagnetic field analysis method capable of performing an accurate electromagnetic field analysis is obtained.

According to the invention of claim 3, according to claim 2,
By recording a program for causing a computer to execute the method described in (1), the program becomes machine-readable, thereby obtaining a recording medium capable of realizing the operation of claim 2 by a computer. Play.

[Brief description of the drawings]

FIG. 1 is a block diagram functionally showing a configuration of an electromagnetic field analyzing apparatus according to an embodiment of the present invention.

FIG. 2 is an explanatory diagram showing a hardware configuration of the electromagnetic field analysis device according to the present embodiment.

FIG. 3 is an explanatory diagram showing a spatial arrangement of an electric field and a magnetic field.

FIG. 4 is a flowchart of a program for obtaining the above temporal change using the FDTD method.

FIG. 5 is an explanatory diagram showing an outline of a one-dimensional FDTD.

FIG. 6 is an explanatory diagram modeling and showing a whole field and a scattered field of the electromagnetic field analyzer according to the present embodiment.

FIG. 7 is an explanatory diagram showing a plane wave propagation calculation result of the electromagnetic field analyzer according to the present embodiment.

FIG. 8 is an explanatory diagram showing a cross-sectional view of FIG. 7 cut in the x direction at a point of z = 350.

FIG. 9 is a flowchart showing a processing procedure of an analysis process of the electromagnetic field analysis device according to the present embodiment.

FIG. 10 is an explanatory diagram modeling and showing a principle of forcibly changing an electromagnetic field on an incident surface.

FIG. 11 is another explanatory diagram modeling and showing the principle of forcibly changing the electromagnetic field on the incident surface.

FIG. 12 is an explanatory diagram showing an outline of a conventional “whole field / scattering field” method.

FIG. 13 is an explanatory diagram showing an ideal calculation result when a state in which a plane wave propagates in a vacuum is calculated using the “whole field / scattering field” method.

FIG. 14 is an explanatory diagram showing a calculation result obtained by actually calculating plane wave propagation by the FDTD method.

FIG. 15 is an explanatory diagram showing FIG. 14 as contour lines.

FIG. 16 is an explanatory diagram showing a cross-sectional view of FIG. 15 cut in the x direction at z = 350.

FIG. 17 is an explanatory diagram showing a propagation state of an incident wave for each lapse of time.

FIG. 18 is an explanatory view showing a model of a whole field and a scattered field in a conventional “whole field / scattered field” method.

[Explanation of symbols]

 Reference Signs List 101 incident wave information input unit 102 incident wave calculation unit 103 structural information input unit 104 calculation parameter input unit 105 electromagnetic field calculation unit 106 electromagnetic field calculation output unit 200 bus 201 memory 202 arithmetic unit 203 recording medium 204 input device 205 display device

 ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Tokuo Chiba 1-8-1, Nakase, Mihama-ku, Chiba-shi Inside Seiko Instruments Inc. (72) Inventor Nobuyuki Kasama 1-8-8, Nakase, Mihama-ku, Chiba-shi, Chiba Inside Iko Instruments Co., Ltd. (72) Inventor Kenji Kato 1-8-1, Nakase, Mihama-ku, Chiba City, Chiba Prefecture Inside Inko Instruments Co., Ltd. (72) Takashi Shinawa 1-8-1, Nakase, Nakase, Mihama-ku, Chiba Prefecture, Sico Instruments Co., Ltd. F-term (reference) 5B056 AA00 BB04 BB95 HH01 9A001 GG11 KK15

Claims (3)

[Claims] 1. An FDTD (Finite Differ)
An electromagnetic field analysis apparatus for analyzing a propagation state of an electromagnetic field by using an input time domain method, an input means for inputting information about an incident wave, and an incident wave calculating means for calculating an incident wave based on the information input by the input means An electromagnetic field, comprising: an electromagnetic field calculation unit that calculates an electromagnetic field based on the incident wave calculated by the incident wave calculation unit; and an output unit that outputs a result calculated by the electromagnetic field calculation unit. Analysis device. 2. An FDTD (Finite Differ)
a first step of inputting information about an incident wave, and a second step of calculating an incident wave based on the information input in the first step. And a third step of calculating an electromagnetic field based on the incident wave calculated in the second step; and a fourth step of outputting a result calculated in the third step. Electromagnetic field analysis method. 3. A computer-readable recording medium on which a program for causing a computer to execute the method according to claim 2 is recorded.