Mutual reciprocity-based system-level electromagnetic coupling fast analysis method, device and equipment
By modeling the cable coupling model of the electronic system as a multi-port antenna system, and utilizing the fast Fourier transform and reciprocity theorem, the problems of low efficiency and insufficient accuracy in electromagnetic coupling response calculation in the prior art are solved, and efficient and accurate solutions for omnidirectional and fully polarized coupling responses are achieved.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- NAT UNIV OF DEFENSE TECH
- Filing Date
- 2026-05-08
- Publication Date
- 2026-06-02
AI Technical Summary
Existing technologies for solving the electromagnetic coupling response of electronic systems under electromagnetic radiation, especially at high frequencies and with complex structures, suffer from low computational efficiency and insufficient accuracy, and cannot accurately determine the coupling terminal response under the worst operating conditions.
A system-level electromagnetic coupling fast analysis method based on reciprocity is adopted. By modeling the cable coupling model of the electronic system as a multi-port antenna system, the induced voltage is calculated by using fast Fourier transform and frequency pointer, combined with antenna reciprocity theorem and impedance matrix, and then converted to the time domain by inverse Fourier transform, so as to achieve efficient solution of omnidirectional and fully polarized coupling response.
It achieves efficient and accurate calculations in the frequency and time domains, significantly reduces computational overhead, and can objectively determine the worst-case coupling terminal response of an electronic system, making it suitable for complex electronic systems.
Smart Images
Figure CN122133360A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of electromagnetic coupling analysis technology, and in particular to a system-level rapid electromagnetic coupling analysis method, apparatus, and equipment based on reciprocity. Background Technology
[0002] Under electromagnetic radiation, energy can couple into the various sub-modules of an electronic system through antennas, cables, and other means. The induced voltage or current generated can easily interfere with the normal operation of the electronic system and even cause damage to the system hardware. Therefore, accurately solving the electromagnetic coupling response of an electronic system under electromagnetic radiation, especially determining the system-level electromagnetic coupling terminal response under the worst operating conditions, has become a key core link in the electromagnetic compatibility design process of electronic systems.
[0003] Currently, the industry mainly uses transmission line theory and numerical methods to solve the coupling response caused by external electromagnetic radiation. Among them, the transmission line theory method establishes a coupling response model under distributed source excitation by extracting characteristic parameters of the transmission line. However, this method is an approximate model and does not consider the influence of many undesirable factors of the transmission line under high-frequency conditions, resulting in insufficient accuracy of the coupling response calculation results at high frequencies. Moreover, this type of transmission line model itself is not applicable to the electromagnetic coupling analysis of complex electronic systems.
[0004] Numerical methods are mainly divided into two categories: frequency domain solutions and time domain solutions. Currently, the industry largely relies on commercial electromagnetic calculation software for modeling and simulation to solve for coupled responses. While this method is applicable to modeling complex electronic systems, electromagnetic waves can enter from various directions in space and couple with the electronic system. To calculate the maximum induced voltage or current generated by electromagnetic waves with different incident directions and polarizations in the system, it is necessary to repeat the electromagnetic simulation operation for each parameter combination. For complex electronic systems, the computation time of a single electromagnetic simulation often reaches tens of minutes. A large number of repeated simulations will result in a huge computational load, leading to extremely low overall computational efficiency. Summary of the Invention
[0005] Therefore, it is necessary to provide a system-level electromagnetic coupling rapid analysis method, apparatus, and equipment based on reciprocity that can efficiently and accurately solve the worst electromagnetic coupling response under different electromagnetic wave incident and polarization conditions to address the above-mentioned technical problems.
[0006] A fast system-level electromagnetic coupling analysis method based on reciprocity, the method comprising: Step 1: Model the cable coupling model of the electronic system as a multi-port antenna system; Step 2: Perform a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave to obtain the complex frequency sequence of the incident electromagnetic wave, and establish a frequency pointer according to the target frequency band. Step 3: Set the simulation frequency range, perform a transmission mode simulation on the multi-port antenna system, sequentially excite a single port while keeping the other ports open, and obtain the far-field radiation direction complex vector, port complex current, and impedance matrix at different frequencies for each port at the frequency point determined by the frequency pointer. Step 4: Based on the antenna reciprocity theorem, according to the far-field radiation direction complex vector and the port complex current, calculate the equivalent open-circuit voltage complex vector generated by the electromagnetic wave at each port at each frequency point under any incident direction and polarization state according to the frequency pointer. Step 5: Calculate the induced voltage complex spectrum sequence of each port load based on the port load impedance matrix, the impedance matrix, and the equivalent open-circuit voltage complex vector; Step 6: Substitute the induced voltage complex spectrum sequence into the incident electromagnetic wave complex frequency sequence and perform inverse Fourier transform to obtain the induced voltage time domain sequence of each port load. Step 7: Perform steps 4 to 6 for each of the preset groups of different incident directions and polarization states to obtain the time-domain sequence of the induced voltage for each port load; calculate the voltage peak value of each induced voltage time-domain sequence, and select the maximum value from the voltage peak values as the worst coupling terminal response of the electronic system.
[0007] On the other hand, a system-level electromagnetic coupling rapid analysis device based on reciprocity is also provided, comprising: The system modeling module is used to model the cable coupling model of an electronic system as a multi-port antenna system. The frequency pointer generation module is used to perform a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave to obtain the complex frequency sequence of the incident electromagnetic wave, and to establish a frequency pointer according to the target frequency band. The parameter acquisition module is used to set the simulation frequency range, perform transmission mode simulation on the multi-port antenna system, sequentially excite a single port while keeping the other ports open, and acquire the far-field radiation direction complex vector, port complex current, and impedance matrix at different frequency points for each port at the frequency point determined by the frequency pointer. The equivalent open-circuit voltage complex vector calculation module is used to calculate the equivalent open-circuit voltage complex vector generated by electromagnetic waves at each port at each frequency point under any incident direction and polarization state based on the antenna reciprocity theorem, according to the far-field radiation direction complex vector and the port complex current, and according to the frequency pointer. The induced voltage complex frequency sequence calculation module is used to calculate the induced voltage complex spectrum sequence of each port load based on the port load impedance matrix, the impedance matrix and the equivalent open-circuit voltage complex vector; The time-domain sequence generation module is used to substitute the induced voltage complex spectrum sequence into the incident electromagnetic wave complex frequency sequence and perform inverse Fourier transform to obtain the induced voltage time-domain sequence of each port load. The worst-case response determination module is used to execute the process from the equivalent open-circuit voltage complex vector calculation module to the time-domain sequence generation module for each preset group of different incident directions and polarization states to obtain the induced voltage time-domain sequence of each port load; calculate the voltage peak value of each induced voltage time-domain sequence, and select the maximum value from the voltage peak values as the worst-case coupling terminal response of the electronic system.
[0008] In another aspect, a computer device is also provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps of the above-described system-level electromagnetic coupling rapid analysis method based on reciprocity.
[0009] Compared with existing technologies, the system-level electromagnetic coupling rapid analysis method, apparatus, and device based on reciprocity provided by this invention have the following beneficial effects: 1. By performing a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave and establishing a frequency pointer based on the target frequency band, the precise definition and orderly management of the frequency domain calculation range are achieved, avoiding the waste of resources caused by invalid calculations across the entire frequency band. At the same time, it ensures that the frequency domain data is processed uniformly in a fixed order, thereby improving the standardization and accuracy of subsequent induced voltage sequence calculations.
[0010] 2. By simulating the emission mode, the complex vector of the far-field radiation direction, the complex current of the port, and the full-frequency impedance matrix of each port at the frequency point determined by the frequency pointer are obtained. There is no need to repeatedly perform electromagnetic simulation for each set of incident directions and polarization states. Moreover, when the time domain sequence of the incident electromagnetic wave changes, the calculation can be called repeatedly, which significantly reduces the computational overhead and shortens the solution time. It is especially suitable for complex electronic systems.
[0011] 3. Based on the antenna reciprocity theorem, the equivalent open-circuit voltage complex vector of each port is directly calculated from the complex vector of the far-field radiation direction and the complex current of the port. Then, the impedance matrix and the port load impedance matrix are combined to complete the accurate solution of the complex frequency sequence of the induced voltage, which fully preserves the correlation between the incident angle of the electromagnetic wave, the polarization state and the coupling characteristics of the system ports.
[0012] 4. After substituting the complex frequency sequence of the induced voltage into the complex frequency sequence of the incident electromagnetic wave and performing an inverse Fourier transform, the time-domain sequence of the induced voltage of each port load is directly obtained. This truly extends the field-line coupling analysis based on the reciprocity theorem from the frequency domain to the time domain, and can directly reflect the change law of the coupling voltage of the cable port with time. It can be used for coupling energy calculation and transient interference assessment.
[0013] 5. By performing equivalent open-circuit voltage calculation, induced voltage complex frequency sequence solution, and time-domain sequence restoration on multiple incident directions and polarization states, and then extracting peak values from each time-domain sequence and selecting the global maximum value, the worst-case coupling terminal response of the electronic system can be determined objectively and comprehensively, avoiding the result deviation caused by selecting only local operating conditions. Attached Figure Description
[0014] To more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention, and those skilled in the art can obtain other related drawings based on these drawings without creative effort.
[0015] Figure 1 This is a flowchart illustrating the rapid analysis method for system-level electromagnetic coupling based on reciprocity in Example 1; Figure 2 This is a schematic diagram of the electromagnetic coupling model in Example 1; Figure 3 This is a structural block diagram of the system-level electromagnetic coupling rapid analysis device based on reciprocity in Example 2; Figure 4 This is a diagram of the internal structure of the computer device in Example 3.
[0016] The objectives, features, and advantages of this invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation
[0017] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative effort are within the scope of protection of the present invention.
[0018] It is understood that the technical solutions of the various embodiments of the present invention can be combined with each other, but only if they are based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such combination of technical solutions does not exist and is not within the scope of protection claimed by the present invention.
[0019] The embodiments of the present invention will now be described in detail with reference to the accompanying drawings.
[0020] Example 1 like Figure 1 As shown, this embodiment provides a fast system-level electromagnetic coupling analysis method based on reciprocity, including the following steps: Step 1: Model the cable coupling model of the electronic system as a multi-port antenna system.
[0021] It can be understood that the cable coupling model of the electronic system is an electromagnetic coupling scenario model consisting of the electronic system and the cable under electromagnetic wave irradiation. Since the cable system has linear electromagnetic characteristics, it is equivalent to a multi-port antenna system consisting of several ports, so as to use the antenna transmit and receive reciprocity characteristics to solve the coupling response and provide a model basis for subsequent reciprocity calculation.
[0022] Step 2: Perform a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave to obtain the complex spectrum sequence of the incident electromagnetic wave, and establish a frequency pointer according to the target frequency band.
[0023] It is understandable that a fast Fourier transform is performed on the time-domain signal sequence of the incident electromagnetic wave to convert the time-domain signal to the frequency domain to obtain the complex spectrum sequence of the incident electromagnetic wave; the frequency points corresponding to the target frequency band are extracted and arranged in ascending order to form a frequency pointer used for calibrating the calculated frequency points, so as to achieve precise control of the target frequency band.
[0024] Step 3: Set the simulation frequency range, perform a transmission mode simulation on the multi-port antenna system, sequentially excite a single port while keeping the other ports open, and obtain the complex vector of the far-field radiation direction, the complex current of the port, and the impedance matrix at different frequencies for each port at the frequency point determined by the frequency pointer.
[0025] It is understandable that the simulation frequency range can be greater than or equal to the frequency pointer range, so that it can be reused when the target frequency band of the analysis changes for different incident electromagnetic waves. A full-wave simulation of the transmission mode is performed on the multi-port antenna system, with each port excited sequentially and the remaining ports open. Using the frequency pointer as the target frequency, the radiation pattern, port complex current, and port complex voltage are obtained, and the impedance matrix at different frequencies is constructed accordingly.
[0026] Step 4: Based on the antenna reciprocity theorem, according to the complex vector of the far-field radiation direction and the complex current of the port, calculate the equivalent open-circuit voltage complex vector generated by the electromagnetic wave at each port at each frequency point under any incident direction and polarization state according to the frequency pointer.
[0027] It is understandable that the equivalent open-circuit voltage is an equivalent voltage source obtained by field excitation coupling based on the reciprocity theorem. This step, based on the antenna reciprocity theorem, uses the far-field radiation direction complex vector and port complex current obtained from the transmission mode as inputs to calculate the equivalent open-circuit voltage generated by electromagnetic waves coupled at each port under any incident direction and polarization state, and forms the equivalent open-circuit voltage complex vector.
[0028] Step 5: Calculate the complex spectrum sequence of induced voltage for each port load based on the port load impedance matrix, the impedance matrix, and the equivalent open-circuit voltage complex vector.
[0029] It can be understood that the induced voltage complex spectrum sequence is the frequency domain voltage sequence obtained by coupling of each port. In this step, based on the system impedance matrix, the port load impedance matrix and the equivalent open-circuit voltage complex vector, the induced current and induced voltage of each port load are solved by matrix operation, and then the induced voltage complex spectrum sequence is obtained by sorting according to the frequency pointer.
[0030] Step 6: Substitute the induced voltage complex spectrum sequence into the incident electromagnetic wave complex spectrum sequence and perform inverse Fourier transform to obtain the induced voltage time-domain sequence of each port load.
[0031] It can be understood that the induced voltage complex spectrum sequence is filled into the corresponding frequency band position of the incident electromagnetic wave complex spectrum sequence, and the electromagnetic wave complex spectrum sequence at non-target frequency band positions is assigned a value of 0 to construct a complete induced voltage complex spectrum sequence; then the frequency domain response is converted into a time domain response waveform through inverse Fourier transform to obtain the induced voltage time domain sequence of each port load.
[0032] Step 7: Perform steps 4 to 6 for each of the preset groups of different incident directions and polarization states to obtain the time-domain sequence of the induced voltage for each port load; calculate the voltage peak value of each induced voltage time-domain sequence, and select the maximum value from the voltage peak values as the worst coupling terminal response of the electronic system.
[0033] It is understandable that the worst coupling response is the maximum coupling voltage under all incident conditions. This step iterates through all preset incident directions and polarization state combinations, repeats steps 4 to 6 to obtain each group of time-domain sequences, extracts the peak values of each group and takes the global maximum value as the worst coupling terminal response of the system.
[0034] The aforementioned system-level electromagnetic coupling rapid analysis method based on reciprocity extends the reciprocity theorem from the frequency domain to the time domain by using multi-port antenna modeling, time-frequency transformation, frequency pointer control, single-transmission simulation, reciprocity theorem calculation, inverse Fourier time-domain reconstruction, and full-condition scanning optimization. It can solve the omnidirectional and full-polarization coupling response with only one simulation, solving the technical problems of long simulation time, frequency-domain-only calculation, and inaccurate worst-case results in traditional methods.
[0035] In the specific implementation process of step 1, the electromagnetic coupling model composed of various cables and structures in the electronic system is abstracted into a model with... A multi-port antenna system, where each cable port corresponds to one port of the antenna system, thus transforms the field-line coupling problem into a multi-port antenna transmit-receive coupling problem. For example... Figure 2 As shown, the electronic system is located at the origin of the coordinate system and is subjected to along the line Electromagnetic wave radiation incident in a certain direction, among which The azimuth angle of the incident electromagnetic wave. The polar angle of the incident electromagnetic wave. It is the polarization angle of the electromagnetic wave.
[0036] This step transforms a complex cable system into an equivalent multi-port antenna system, overcoming the limitation that analytical methods are only applicable to simple cables and making them suitable for modeling complex electronic systems.
[0037] In the specific implementation of step 2, firstly, for any known incident electromagnetic wave time-domain signal sequence... The corresponding complex spectrum sequence of the incident electromagnetic wave is obtained by using Fast Fourier Transform (FFT). .
[0038] Suppose the target frequency band contains a total of sequence indices of 1. arrive of If there are 10 spectral values, then the sequence of the complex frequency sequence from which the required frequency band is calculated is extracted as the target frequency band sequence. The expression for the target frequency band sequence is: ; Arrange the frequency points corresponding to the target frequency band sequence from smallest to largest to form a frequency pointer, represented as: ; In the formula, This represents the subsequence corresponding to the target frequency band in the complex spectrum vector sequence of the incident electromagnetic wave's electric field; Indicates frequency point; This represents the unextracted complex spectrum sequence; This represents the extracted complex spectrum sequence; Indicates the first frequency The corresponding complex spectrum vector of the electric field of the incident electromagnetic wave; Indicates the first frequency band within the target frequency band One frequency; Indicates the first The frequencies corresponding to a complex spectrum sequence; Indicates the total frequency band sequence length, and is the index of the total sequence, containing a total of [number of bands]. One value; Indicates the length of the target frequency band sequence, and is the index of the target frequency band sequence, containing a total of Each spectrum value, One frequency value; Indices representing frequency domain sequences, containing a total of One value, ; This indicates the frequency pointer.
[0039] This step achieves precise control of the target frequency band through time-frequency transformation and frequency pointer, reducing invalid calculations and improving the efficiency and standardization of frequency domain processing.
[0040] In the specific implementation of step 3, it is assumed that the cable coupling model of the electronic system has Each port and the model material is a linear material, which will include The cable model of each port is considered as The first antenna, then the antenna system's first antenna... Each port is powered by a voltage source. Excitation, the internal impedance of the voltage source is , No. The port complex voltage and port complex current of each port are respectively , , No. The self-impedance of the antenna at each port is The mutual coupling between the individual antennas is controlled by a current-controlled voltage source. To express.
[0041] Set up the simulation monitor according to the frequency pointer frequency, set up current and voltage monitoring ports on each port, and perform transmission mode simulation. Obtain the first transmission mode by sequentially exciting a single port and opening the rest of the ports. Complex vector of far-field radiation direction under port excitation, complex voltage at the port With port complex current and the When one port is used, the port complex voltage of the other ports is... After the simulation is completed, the above files will be saved in txt file format.
[0042] The complex vector expression for the far-field radiation direction is: ; In the formula, Indicates the first Complex vector of far-field radiation direction when excited at each port; Indicates frequency point; Represents the imaginary unit; Wave number is the amount of phase change of an electromagnetic wave per unit length in space. This represents the radial distance from the observation point to the radiation source; A unit vector representing the polar angle direction in spherical coordinates; A unit vector representing the azimuth direction in a spherical coordinate system; Representing radiation pattern Quantity; Representing radiation pattern Quantity.
[0043] Based on the port complex voltage and port complex current obtained from simulation, according to the frequency pointer The frequency values are solved sequentially column by column, at the frequency point. The impedance matrix at time The series ultimately solves for different frequency points. The impedance matrix at time , the first line, number The expression for the impedance value of the column is: ; In the formula, Representing the impedance matrix No. line, number The impedance value of the column, Indicates the launch mode number The port complex voltage corresponding to each port; Indicates the launch mode number The port complex current corresponding to the excitation of each port. It is worth noting that in the impedance value expression, With the Each port corresponds to one port. With the Each port corresponds to one. Furthermore, calculations must be performed column-by-column; for example, when calculating the first column, the denominator is always... That is, the port complex current corresponding to the first port excitation in the transmit mode.
[0044] Arrange the impedance values by rows and columns to obtain the expression for the impedance matrix: .
[0045] This step only requires performing a single transmit mode simulation for each port to obtain the basic parameters for full-condition calculation, eliminating the need for repeated simulations and significantly reducing the computational load and solution time.
[0046] In the specific implementation of step 4, theoretical calculations are first performed based on the antenna reciprocity theorem, and the antenna radiation field satisfies... ,in, Indicates pointing from the origin to space The unit vector at a point can be calculated as follows: ; Indicates pointing from the origin to space The vector at the point, Representing a spatial point The coordinates of the point in the rectangular coordinate system; Representing a rectangular coordinate system The unit vector of direction; Represents the distance from the origin to a point in space. The length at the location; This represents the phase change of the radiation field. The complex form representing the spatial radiation field. Radiation field. With the direction of radiation Related, in spherical coordinates, can be represented as In the form of.
[0047] It is understandable that the transmit and receive mode parameters of the same antenna are reciprocal. Based on this theorem, and in combination with the transmit mode parameters, the reciprocity of the first incident direction and polarization state can be calculated. The equivalent open-circuit voltage of each port is expressed as: ; In the formula, Indicates the first Each port at a frequency of The equivalent open-circuit complex voltage value at that time; Represents the imaginary unit; Indicates frequency as The corresponding wavelength; Indicates free-space wave impedance; Indicates the first Each port at a frequency of The complex current value at that time; Indicates the first The frequency when each port is excited is The complex vector of the far-field radiation direction; Indicates the first frequency points The corresponding complex spectrum vector of the electric field of the incident electromagnetic wave; Indicates the polar angle of the incident electromagnetic wave; Indicates the azimuth angle of the incident electromagnetic wave; This represents the polarization angle of the incident electromagnetic wave.
[0048] The equivalent open-circuit voltages of each port are integrated into a complex vector of equivalent open-circuit voltages: ; This step utilizes the reciprocity theorem to achieve the equivalent transformation from field excitation to port complex voltage source, fully preserving the correlation between incident angle, polarization, and coupling characteristics. The calculation results have clear physical meaning and high accuracy.
[0049] In the specific implementation of step 5, a diagonal port load impedance matrix is established: ; In the formula, Represents the port load impedance matrix; Indicates the first The impedance value of the load connected to each port.
[0050] Combined with impedance matrix With equivalent open-circuit voltage complex vector First calculate the induced current vector. The expression for calculating the induced current in the induced current vector is: .
[0051] Then calculate the induced voltage vector. The expression for calculating the induced voltage in the induced voltage vector is as follows: .
[0052] Arrange the induced voltages in frequency pointer order to obtain the first... Complex spectrum sequence of induced voltage of port load: ; In the formula, Indicates the first The complex spectrum sequence of induced voltages corresponding to each port load; Indicates the first The port at the A sequence, i.e., a frequency of The complex spectrum value at time.
[0053] This step uses matrix operations to accurately solve the frequency domain coupled response, and relies on the frequency pointer to organize the data sequence, ensuring the consistency and traceability of the calculation.
[0054] In the specific implementation of step 6, the incident electromagnetic wave complex spectrum sequence Complex values outside the calculated frequency band are set to 0, and the induced voltage complex spectrum sequence is used. Replace the complex values within the calculated frequency band to construct the first The complete complex spectrum sequence of the induced voltage of each port load: ; In the formula, Indicates the first The complete induced voltage complex spectrum sequence corresponding to each port load.
[0055] Then, the complete induced voltage complex spectrum sequence Performing an inverse Fourier transform yields the time-domain sequence of the induced voltage at each port load, expressed as: ; In the formula, Indicates the first Time-domain sequence of induced voltage corresponding to each port load; Represents the imaginary unit; Indices representing frequency domain sequences, containing a total of One value, .
[0056] This step extends the coupling analysis from the frequency domain to the time domain, and can directly output the time-domain waveform of the port induced voltage, providing data support for coupling energy calculation and transient interference assessment.
[0057] In the specific implementation of step 7, the polar angles of all preset incident electromagnetic waves are traversed. Azimuth With polarization angle For each set of incident parameters, steps 4 to 6 are performed to obtain the corresponding induced voltage time-domain sequence; the voltage peak value under a single set of incident conditions is extracted using the ratio method. ; In the formula, Indicates the first The peak voltage of each induced voltage time-domain sequence at each port under a set of incident directions and polarization states; Indicates the polar angle of the incident electromagnetic wave; Indicates the azimuth angle of the incident electromagnetic wave; This represents the polarization angle of the incident electromagnetic wave.
[0058] The maximum value among all voltage peaks is selected as the worst-case coupling termination response of the electronic system, expressed as: ; In the formula, Indicates the first The worst-case coupling terminal response for each port.
[0059] Preferably, a Python script is used to implement parametric scanning and automatic peak extraction, thereby improving computational efficiency.
[0060] This step obtains the global peak value by scanning parameters under full operating conditions, avoiding the result deviation caused by local operating condition simulation, and objectively obtaining the worst-case coupled terminal response at the system level.
[0061] It should be understood that, although this embodiment Figure 1 The steps are shown sequentially as indicated by the arrows, but they are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order in which these steps are performed; they can be executed in other orders. Figure 1 At least some of the steps in the process may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least some of the sub-steps or stages of other steps.
[0062] Example 2 Based on the reciprocity-based system-level electromagnetic coupling rapid analysis method in Embodiment 1, this embodiment discloses a reciprocity-based system-level electromagnetic coupling rapid analysis device, such as... Figure 3 As shown, the system-level electromagnetic coupling rapid analysis device based on reciprocity includes: a system modeling module 801, a frequency pointer generation module 802, a parameter acquisition module 803, an equivalent open-circuit voltage complex vector calculation module 804, an induced voltage complex frequency sequence calculation module 805, a time-domain sequence generation module 806, and a worst-case response determination module 807, wherein: System modeling module 801 is used to model the cable coupling model of the electronic system as a multi-port antenna system; The frequency pointer generation module 802 is used to perform a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave to obtain the complex spectrum sequence of the incident electromagnetic wave, and to establish a frequency pointer according to the target frequency band.
[0063] The parameter acquisition module 803 is used to set the simulation frequency range, perform transmission mode simulation on the multi-port antenna system, sequentially excite a single port while keeping the other ports open, and acquire the far-field radiation direction complex vector, port complex current, and impedance matrix at different frequency points for each port at the frequency point determined by the frequency pointer.
[0064] The equivalent open-circuit voltage complex vector calculation module 804 is used to calculate the equivalent open-circuit voltage complex vector generated at each port of the electromagnetic wave at each frequency point under any incident direction and polarization state based on the antenna reciprocity theorem, according to the complex vector of the far-field radiation direction and the port complex current, and according to the frequency pointer.
[0065] The induced voltage complex frequency sequence calculation module 805 is used to calculate the induced voltage complex frequency sequence of each port load based on the port load impedance matrix, the impedance matrix and the equivalent open-circuit voltage complex vector.
[0066] The time-domain sequence generation module 806 is used to substitute the induced voltage complex spectrum sequence into the incident electromagnetic wave complex spectrum sequence and perform inverse Fourier transform to obtain the induced voltage time-domain sequence of each port load.
[0067] The worst-case response determination module 807 is used to execute the process from the equivalent open-circuit voltage complex vector calculation module to the time-domain sequence generation module for each preset group of different incident directions and polarization states to obtain the time-domain sequence of the induced voltage of each port load; calculate the voltage peak value of each induced voltage time-domain sequence, and select the maximum value from the voltage peak value as the worst-case coupling terminal response of the electronic system.
[0068] In this embodiment, the specific working process and working principle of the system modeling module 801, frequency pointer generation module 802, parameter acquisition module 803, equivalent open-circuit voltage complex vector calculation module 804, induced voltage complex frequency sequence calculation module 805, time-domain sequence generation module 806, and worst-case response determination module 807 are the same as those in Embodiment 1, and therefore will not be described again in this embodiment. Each unit module can be implemented entirely or partially through software, hardware, or a combination thereof. Each unit module can be embedded in or independent of the processor in a computer device in hardware form, or stored in the memory of a computer device in software form, so that the processor can call and execute the operations corresponding to each of the above unit modules.
[0069] Example 3 like Figure 4 The diagram illustrates a computer device disclosed in this embodiment, including a transmitter, a receiver, a memory, and a processor. The transmitter is used to send instructions and data, the receiver is used to receive instructions and data, the memory is used to store computer execution instructions, and the processor is used to execute the computer execution instructions stored in the memory to implement the method in Embodiment 1 above.
[0070] It is important to note that the aforementioned memory can be either standalone or integrated with the processor. When the memory is set up independently, the terminal device also includes a bus for connecting the memory and the processor.
[0071] Example 4 This embodiment discloses a computer-readable storage medium storing computer-executable instructions. When a processor executes the computer-executable instructions, it implements the method in Embodiment 1 above.
[0072] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium, and when executed, it can include the processes of the embodiments of the above methods. Any references to memory, storage, databases, or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), dual data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.
[0073] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.
[0074] The embodiments described above are merely examples of several implementations of the present invention, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of the invention. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of the present invention, and these modifications and improvements all fall within the scope of protection of the present invention.
Claims
1. A system-level rapid analysis method for electromagnetic coupling based on reciprocity, characterized in that, The method includes: Step 1: Model the cable coupling model of the electronic system as a multi-port antenna system; Step 2: Perform a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave to obtain the complex spectrum sequence of the incident electromagnetic wave, and establish a frequency pointer according to the target frequency band. Step 3: Set the simulation frequency range, perform a transmission mode simulation on the multi-port antenna system, sequentially excite a single port while keeping the other ports open, and obtain the far-field radiation direction complex vector, port complex current, and impedance matrix at different frequencies for each port at the frequency point determined by the frequency pointer. Step 4: Based on the antenna reciprocity theorem, according to the far-field radiation direction complex vector and the port complex current, calculate the equivalent open-circuit voltage complex vector generated by the electromagnetic wave at each port at each frequency point under any incident direction and polarization state according to the frequency pointer. Step 5: Calculate the induced voltage complex spectrum sequence of each port load based on the port load impedance matrix, the impedance matrix, and the equivalent open-circuit voltage complex vector; Step 6: Substitute the induced voltage complex spectrum sequence into the incident electromagnetic wave complex spectrum sequence and perform inverse Fourier transform to obtain the induced voltage time domain sequence of each port load. Step 7: Perform steps 4 to 6 for each of the preset groups of different incident directions and polarization states to obtain the time-domain sequence of the induced voltage for each port load; calculate the voltage peak value of each induced voltage time-domain sequence, and select the maximum value from the voltage peak values as the worst coupling terminal response of the electronic system.
2. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to claim 1, characterized in that, In step 2, a frequency pointer is established based on the target frequency band, including: for any known incident electromagnetic wave time-domain signal sequence The corresponding complex spectrum sequence of the incident electromagnetic wave is obtained by using the Fast Fourier Transform. ; Suppose the target frequency band contains a total of sequence indices of 1. arrive of If there are 10 spectral values, then the target frequency band sequence can be extracted from the complex frequency sequence of the incident electromagnetic wave, and the expression is: ； Arrange the frequency points corresponding to the target frequency band sequence from smallest to largest to form a frequency pointer, represented as: ； In the formula, This represents the subsequence corresponding to the target frequency band in the complex spectrum vector sequence of the incident electromagnetic wave's electric field; Indicates frequency point; This represents the unextracted complex spectrum sequence; This represents the extracted complex spectrum sequence; Indicates the first Frequency The corresponding complex spectrum vector of the electric field of the incident electromagnetic wave; Indicates the first frequency band within the target frequency band One frequency; Indicates the first The frequencies corresponding to a complex spectrum sequence; Indicates the total frequency band sequence length; Indicates the length of the target frequency band sequence; Indicates the index of the frequency domain sequence; This indicates the frequency pointer.
3. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to claim 1, characterized in that, In step 4, in the equivalent open-circuit voltage complex vector, the first... The expression for calculating the equivalent open-circuit voltage of each port is: ； In the formula, Indicates the first Each port at a frequency of The equivalent open-circuit complex voltage value at that time; Represents the imaginary unit; Indicates frequency as The corresponding wavelength; Indicates free-space wave impedance; Indicates the first Each port at a frequency of The complex current value at that time; Indicates the first The frequency when each port is excited is The complex vector of the far-field radiation direction; Indicates the first frequency points The corresponding complex spectrum vector of the electric field of the incident electromagnetic wave; Indicates the polar angle of the incident electromagnetic wave; Indicates the azimuth angle of the incident electromagnetic wave; This represents the polarization angle of the incident electromagnetic wave.
4. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to any one of claims 1 to 3, characterized in that, In step 5, the time-domain sequence of the induced voltage of each port load is represented as follows: ； In the formula, Indicates the first The complex spectrum sequence of induced voltages corresponding to each port load; Indicates the first The port in the first There are sequences, i.e., frequencies of The complex spectrum value at time.
5. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to claim 4, characterized in that, In step 6, the induced voltage complex spectrum sequence is substituted into the incident electromagnetic wave complex spectrum sequence to obtain the complete induced voltage complex spectrum sequence on different port loads, as follows: ； In the formula, Indicates the first The complete complex spectrum sequence of induced voltages corresponding to each port load; Indicates the total frequency band sequence length; Indicates the length of the target frequency band sequence.
6. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to claim 5, characterized in that, In step 6, an inverse Fourier transform is performed on the complete induced voltage complex spectrum sequence to obtain the induced voltage time-domain sequence of each port load, as shown below: ； In the formula, Indicates the first Time-domain sequence of induced voltage corresponding to each port load; Represents the imaginary unit; Indices representing frequency domain sequences, containing a total of One value, .
7. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to claim 6, characterized in that, The peak voltage of each induced voltage time-domain sequence is calculated using the following expression: ； In the formula, Indicates the first The peak voltage of each induced voltage time-domain sequence at each port under a set of incident directions and polarization states; Indicates the polar angle of the incident electromagnetic wave; Indicates the azimuth angle of the incident electromagnetic wave; This represents the polarization angle of the incident electromagnetic wave.
8. The system-level electromagnetic coupling rapid analysis method based on reciprocity according to claim 7, characterized in that, The maximum value selected from the voltage peaks is taken as the worst-case coupling termination response of the electronic system, expressed as: ； In the formula, Indicates the first The worst-case coupling terminal response for each port.
9. A system-level electromagnetic coupling rapid analysis device based on reciprocity, characterized in that, The device includes: The system modeling module is used to model the cable coupling model of an electronic system as a multi-port antenna system. The frequency pointer generation module is used to perform a fast Fourier transform on the time-domain signal sequence of the incident electromagnetic wave to obtain the complex spectrum sequence of the incident electromagnetic wave, and to establish a frequency pointer according to the target frequency band. The parameter acquisition module is used to set the simulation frequency range, perform transmission mode simulation on the multi-port antenna system, sequentially excite a single port while keeping the other ports open, and acquire the far-field radiation direction complex vector, port complex current, and impedance matrix at different frequency points for each port at the frequency point determined by the frequency pointer. The equivalent open-circuit voltage complex vector calculation module is used to calculate the equivalent open-circuit voltage complex vector generated by electromagnetic waves at each port at each frequency point under any incident direction and polarization state based on the antenna reciprocity theorem, according to the far-field radiation direction complex vector and the port complex current, and according to the frequency pointer. The induced voltage complex frequency sequence calculation module is used to calculate the induced voltage complex spectrum sequence of each port load based on the port load impedance matrix, the impedance matrix and the equivalent open-circuit voltage complex vector; The time-domain sequence generation module is used to substitute the induced voltage complex spectrum sequence into the incident electromagnetic wave complex spectrum sequence and perform an inverse Fourier transform to obtain the induced voltage time-domain sequence of each port load. The worst-case response determination module is used to execute the process from the equivalent open-circuit voltage complex vector calculation module to the time-domain sequence generation module for each preset group of different incident directions and polarization states to obtain the induced voltage time-domain sequence of each port load; calculate the voltage peak value of each induced voltage time-domain sequence, and select the maximum value from the voltage peak values as the worst-case coupling terminal response of the electronic system.
10. A computer device, comprising a memory and a processor, characterized in that, The memory stores a computer program, and when the processor executes the computer program, it implements the steps of the system-level electromagnetic coupling fast analysis method based on reciprocity as described in any one of claims 1 to 8.