Rapid calculation and analysis method for electromagnetic interference between radio systems

Abstract

The embodiment of the invention provides a method for quickly calculating and analyzing electromagnetic interference between radio systems, which comprises the following steps of: circularly selecting interference equipment and interfered equipment from an equipment list in sequence; judging whether the interfering device and the interfered device belong to the same device group, if so, judging that no interference exists between the two groups of devices, jumping out of the cycle, and carrying out the next cycle; if not, entering the next step; determining a preposition parameter value; determining conditions of sideband interference, main frequency interference and harmonic interference between the two selected devices; determining a stray interference condition existing between the two selected devices; judging whether the circulation is finished or not, and if yes, displaying a calculation result; and if not, jumping out of the cycle, and carrying out the next cycle. According to the embodiment of the invention, the rationality of the whole analysis and calculation process is improved, the calculation speed and the calculation precision of the whole calculation and analysis process are improved, and the interference condition between electronic systems is represented more fully and comprehensively.

Description

Technical Field

[0001] The invention relates to the technical field of electromagnetic interference calculation and analysis, in particular to a method for rapid calculation and analysis of electromagnetic interference between radio systems. Background Art

[0002] Electromagnetic compatibility refers to the ability of electronic and electrical equipment or systems to work normally in the expected electromagnetic environment according to the design requirements. It is an important technical performance of electronic and electrical equipment or systems, which represents the degree of compatibility between equipment or systems coexisting in the same electromagnetic environment. With the development of modern science and technology and the increasingly widespread application of various electrical and electronic equipment, the electromagnetic environment in space is becoming increasingly complex. Especially in a local space area where a large number of radar, communication, navigation and other radio systems are deployed, there are many types of radio transceiver systems and dense electromagnetic signals. Serious mutual interference may occur between the radio systems in this local area, resulting in reduced performance or even failure.

[0003] Therefore, when a large number of electronic devices are deployed in regional scenarios and large platforms, electromagnetic compatibility calculation and analysis is required. In particular, before performing precise electromagnetic calculations, rapid analysis is required to improve the speed of solution iteration. However, the current calculation process is relatively cumbersome. Under the premise of deploying a large number of electronic devices, the performance requirements of the computer are very high, and the calculation speed is difficult to meet the needs of rapid analysis. At the same time, the current calculation method is not comprehensive for the relevant parameters used to characterize interference, and it is difficult for engineers to quickly grasp the specific situation of electromagnetic interference through the calculation results. Therefore, there is an urgent need for a rapid calculation and analysis method for electromagnetic interference between radio systems to assist engineers in better and faster electromagnetic compatibility optimization design. Summary of the invention

[0004] In view of the above problems existing in the prior art, an embodiment of the present invention provides a method for quickly calculating and analyzing electromagnetic interference between radio systems to solve the above technical problems existing in the prior art.

[0005] The embodiment of the present invention provides a method for rapid calculation and analysis of electromagnetic interference between radio systems, comprising the following steps:

[0006] Step S1, selecting the interfering device dev1 and the interfered device dev2 from the device list in order and reading the attribute values ​​of the selected devices;

[0007] Step S2, determine whether the interfering device dev1 and the interfered device dev2 belong to the same device group.

[0008] If yes, it is determined that there is no interference between the two groups of devices, no calculation and analysis results are generated, and the current loop is exited and returned to step S1 for the next loop;

[0009] If not, proceed to step S3;

[0010] Step S3, determining the values ​​of the pre-parameters, wherein the pre-parameters include the coordinate difference (dx, dy) of dev2 relative to dev1, the azimuth angle angRhoTR of the vector pointing from dev1 to dev2, the azimuth angle angRhoRT of the vector pointing from dev2 to dev1, the distance r between the antennas of the two devices, the elevation scanning angle theta, the azimuth first sidelobe relative to the main lobe level difference deltaGainT of the transmitting device, and the azimuth first sidelobe relative to the main lobe level difference deltaGainR of the receiving device;

[0011] Step S4, determining the sideband interference, main frequency interference and harmonic interference between the two selected devices;

[0012] Step S5, determining the spurious interference between the two selected devices;

[0013] Step S7, determine whether the cycle is finished,

[0014] If yes, the calculation result is displayed;

[0015] If not, then exit this loop and proceed to the next loop.

[0016] Compared with the prior art, the beneficial effect of a method for rapid calculation and analysis of electromagnetic interference between radio systems provided by an embodiment of the present invention is that: the embodiment of the present invention optimizes the design of the calculation and analysis process, classifies and merges the relevant interference characterization parameters, thereby improving the rationality of the entire analysis and calculation process, effectively reducing the amount of repeated calculations of intermediate calculation results, and improving the calculation speed and calculation accuracy of the entire calculation and analysis process, so that the interference situation between electronic systems is more fully and comprehensively characterized, and engineers can obtain the calculation and analysis results faster, more comprehensively and more clearly, which greatly simplifies the analysis difficulty and analysis workload of engineers, and is more friendly to the subsequent electromagnetic compatibility optimization design performed by engineers. DETAILED DESCRIPTION

[0017] In order to enable those skilled in the art to better understand the technical solution of the present invention, the present invention is described in detail below in conjunction with specific implementation methods.

[0018] Various aspects and features of the present application are described herein.

[0019] These and other characteristics of the present application will become apparent from the following description of a preferred form of embodiment given as non-limiting example.

[0020] It should also be understood that, although the present application has been described with reference to some specific examples, those skilled in the art will be able to realize many other equivalent forms of the present application that have the features described in the claims and are therefore within the scope of protection defined thereby.

[0021] The above and other aspects, features and advantages of the present application will become more apparent in view of the following detailed description.

[0022] Specific embodiments of the present application are described hereinafter; however, it should be understood that the embodiments applied for are merely examples of the present application, which may be implemented in a variety of ways. Well-known and / or repeated functions and structures are not described in detail to determine the true intent based on the user's historical operations and to avoid unnecessary or redundant details that make the present application unclear. Therefore, the specific structural and functional details applied for herein are not intended to be limiting, but only serve as the basis and representative basis for the claims to teach those skilled in the art to use the present application in a variety of ways with substantially any suitable detailed structure.

[0023] This specification may use the phrases "in one embodiment," "in another embodiment," "in a further embodiment," or "in other embodiments," which may all refer to one or more of the same or different embodiments according to the present application.

[0024] The principles and features of the present invention are described below. The examples are only used to explain the present invention and are not used to limit the scope of the present invention. The preferred embodiments of the present invention are further described in detail below:

[0025] The embodiment of the present invention provides a method for rapid calculation and analysis of electromagnetic interference between radio systems, comprising the following steps:

[0026] Step S1, selecting the interfering device dev1 and the interfered device dev2 from the device list in order and reading the attribute values ​​of the selected devices;

[0027] Step S2, determine whether the interfering device dev1 and the interfered device dev2 belong to the same device group.

[0028] If yes, it is determined that there is no interference between the two groups of devices, no calculation and analysis results are generated, and the current loop is exited and returned to step S1 for the next loop;

[0029] If not, proceed to step S3;

[0030] Step S3, determining the values ​​of the pre-parameters, wherein the pre-parameters include the coordinate difference (dx, dy) of dev2 relative to dev1, the azimuth angle angRhoTR of the vector pointing from dev1 to dev2, the azimuth angle angRhoRT of the vector pointing from dev2 to dev1, the distance r between the antennas of the two devices, the elevation scanning angle theta, the azimuth first sidelobe relative to the main lobe level difference deltaGainT of the transmitting device, and the azimuth first sidelobe relative to the main lobe level difference deltaGainR of the receiving device;

[0031] Step S4, determining the sideband interference, main frequency interference and harmonic interference between the two selected devices;

[0032] Step S5, determining the spurious interference between the two selected devices;

[0033] Step S7, determine whether the cycle is finished,

[0034] If yes, the calculation result is displayed;

[0035] If not, then exit this loop and proceed to the next loop.

[0036] Among them, the equipment list of step S1 can be a list of all electronic devices involved in the calculation area, or it can be a list of some electronic devices screened out by engineers after setting certain limiting conditions based on their own experience or the design parameters of the equipment itself (such as characteristic parameters such as frequency, power, waveform, etc. of the interference source, sensitivity parameters such as the operating frequency range and anti-interference ability of the interfered equipment, parameters such as coupling paths such as conduction coupling, radiation coupling, capacitive coupling and inductive coupling), etc. Of course, the list of the interfering equipment dev1 and the list of the interfered equipment dev2 can be consistent or inconsistent. At the same time, the loop selection method in step S1 can be traversed and selected according to the initial order of the list, or it can be traversed and selected after reordering according to certain attribute parameter setting importance conditions;

[0037] Normally, two devices belonging to the same device group (i.e., the interfering device dev1 and the interfered device dev2) have undergone electromagnetic interference compatibility calculation analysis and optimization design in the optimization design stage within the device group. Therefore, in this step, the two devices in the device group are no longer subjected to repeated calculation analysis to avoid excessive redundant calculations and improve the speed of the entire fast calculation analysis.

[0038] In step S3, the pre-parameters required in the subsequent steps are directly calculated separately to ensure that the subsequent steps can directly call the relevant parameters when they need to be referenced, without the problem of repeated calculation in different steps or different parameter calculations in the same step. At the same time, such an operation method can also optimize the calculation logic of the parameters, shorten the expression length of each parameter as much as possible, and make the parameter calculation more readable and more robust.

[0039] In one embodiment, the step S4 specifically includes:

[0040] Step S41, determining the maximum harmonic number maxXbNum;

[0041] Step S42, cycle the harmonic number xbNum from 1 to the maximum harmonic number maxXbNum;

[0042] Step S43, determine whether the xbNum harmonic level exists.

[0043] If yes, proceed to step S44;

[0044] If not, then jump out of this loop and enter the next loop;

[0045] Step S44, determining the values ​​of interference-related parameters, wherein the interference-related parameters include the lower limit of the harmful frequency band harmFS, the upper limit of the harmful frequency band harmFE, whether interference exists harmExiSt, the lower limit of the interference frequency band FS, the upper limit of the interference frequency band FE, the interference calculation frequency harmFreq, the transmission sideband attenuation tBDSJ and the reception sideband attenuation rBDSJ, wherein harmExist is a Boolean type parameter, when its value is 1, it indicates that interference exists, and when its value is 0, it indicates that interference does not exist;

[0046] Step S45, judging whether interference exists according to interference-related parameters,

[0047] If not, then jump out of this loop and enter the next loop;

[0048] If so, determine the value of the interference situation result parameter, the interference situation result parameter includes the interference band WaveBand, the interference type type, the spurious number zsID and the ratio of the interference bandwidth to the receiving bandwidth ibVsRb, wherein the value of the interference type type is "main frequency", "harmonic" or "sideband", and the specific judgment logic is that if xbNum≠1, the interference type is "harmonic", and xbNum can also be used to characterize the number of harmonic interferences, that is, xbNum subharmonic interference. If xbNum=1, subsequent judgment is required to clarify the specific interference type, that is, a. If dev1FS≤harmFreq≤dev1FE, the interference type is "main frequency", b. Otherwise, the interference type is "sideband", and when the spurious number zslD is 0, it indicates that the interference type is not spurious interference. Since the judgment of the sideband interference, main frequency interference and harmonic interference is performed in this step, zslD needs to be assigned a value of 0 in this step;

[0049] Step S46, determining whether the interfering device dev1 and the interfered device dev2 are at the same position on the xy plane.

[0050] If yes, the field strength fldStrength of the interfering device at the interfered device and the coupling value recPower are assigned to infinity∞;

[0051] If not, determine the gain dev1Gt of the interfering device, the gain dev2Gr of the interfered device, the polarization mismatch loss polLoss, the field strength fldStrength of the interfering device at the interfered device, and the coupling value recPower;

[0052] Step S47, determining the interference illumination state and intensity.

[0053] In one embodiment, step S5 includes:

[0054] Step S51, determining the number n of spurious levels;

[0055] Step S52, reading the ppth spurious level from 1 to n in a loop;

[0056] Step S53, determining the values ​​of interference-related parameters, wherein the interference-related parameters include the lower limit of the harmful frequency band harmFS, the upper limit of the harmful frequency band harmFE, whether interference exists harmExiSt, the lower limit of the interference frequency band FS, the upper limit of the interference frequency band FE, the interference calculation frequency harmFreq and the receiving sideband attenuation rBDSJ, wherein harmExist is a Boolean type parameter, when its value is 1, it indicates that interference exists, and when its value is 0, it indicates that interference does not exist;

[0057] Step S54: judging whether interference exists according to interference-related parameters.

[0058] If not, then jump out of this loop and enter the next loop;

[0059] If yes, determine the interference situation result parameter value, the interference situation result parameter includes the interference band WaveBand, the interference type type, the spurious signal zslD and the ratio of the interference bandwidth to the receiving bandwidth ibVsRb. Since this step is used to determine whether there is spurious interference, the value of the interference type type in this step is "spurious";

[0060] Step S56, determine whether the interfering device dev1 and the interfered device dev2 are in the same position on the xy plane.

[0061] If yes, the field strength fldStrength of the interfering device at the interfered device and the coupling value recPower are assigned to infinity∞;

[0062] If not, determine the gain devlGt of the interfering device, the gain dev2Gr of the interfered device, the transmission sideband attenuation tBDSJ, the polarization mismatch loss polLoss, the field strength fldStrength of the interfering device at the interfered device, and the coupling value recPower;

[0063] Step S57, determining the interference illumination state and intensity.

[0064] Although the method has the problem of displaying the calculation and analysis results in a pop-up window, in order to facilitate subsequent inquiries, in one embodiment, the method further includes between step S5 and step S7:

[0065] Step S6, export the calculation results and store them to the specified location. The calculation results are exported in the form of Excel and Word, and the calculation results all use the interfered device as the export index. When performing the export operation, it will first query whether the attachment save location has been selected. If so, the calculation result will be generated as an attachment and the attachment will be uploaded to the specified location. If not, it will be determined whether the matrix position has been selected. If so, the electromagnetic interference matrix will be generated at the specified location and the method will end. Otherwise, the method will end directly.

[0066] Taking into account, according to the requirements of GJB151A "Requirements for Electromagnetic Emissions and Sensitivity of Military Equipment and Subsystems", for harmonics above the third harmonic, the harmonic suppression of the equipment should be above 80dB. At this time, it is considered that it will not cause interference to other equipment in the far field. Therefore, in this fast calculation and analysis method, only the interference of the second and third harmonics to other equipment is analyzed, that is, the value of the maximum harmonic number maxXbNum is 3.

[0067] Example 1

[0068] 5.4.2.2.1 Parameter Description

[0069] The parameters involved are mainly divided into two parts, one is the input attribute parameters, and the other is the parameters that need to be obtained through attribute parameter calculation. Among them, the leftmost column in each table of the main process calculation logic is basically the parameters that need to be obtained through calculation.

[0070] Transmitter device parameters (need to be manually entered by the user)

[0071]

[0072]

[0073]

[0074] Receive device parameters (need to be entered manually by the user)

[0075]

[0076]

[0077] Gain parameter (needs to be manually entered by the user)

[0078]

[0079] Other parameters that need to be input include the maximum harmonic number maxXbNum, the attenuation per 100 MHz of the sideband extension BDSJp100MHz and the comprehensive judgment parameter ZongParameter.

[0080] 5.4.2.2.2 Main Process Calculation Logic

[0081] 1. Determine whether the interfering device and the interfered device belong to the same device group. If so, there is no interference and no subsequent calculation is required;

[0082] 2. Otherwise, calculate a series of pre-parameters

[0083]

[0084]

[0085]

[0086] 3. Cycle harmonic number xbNum from 1 to the maximum harmonic number maxXbNum

[0087] a. Determine whether the xbNum harmonic level exists. If not, skip this cycle calculation.

[0088] b. Otherwise, calculate a series of parameters

[0089]

[0090]

[0091] c. If interference exists, then

[0092] i. Calculate other parameters

[0093]

[0094]

[0095] i. If the two devices have the same position on the xy plane (xyDis<10 -6 )

[0096] Calculation parameters Variable Name formula Field strength of the interfering device at the interfered device fldStrength ∞ Coupling value recPower ∞

[0097] ii. If the two devices are not in the same location, then

[0098]

[0099]

[0100] iii. Determine the radiation status and intensity of interference. For the logic, see the section on determining the radiation status and intensity of interference.

[0101] iv. Generate a calculation result instance and fill the calculation result into the attributes of the instance model

[0102] v. Fill the calculation results and calculation input parameters in the export file according to the format. For details on the logic, see the "Export File" section.

[0103] 4. Loop to read the ppth spurious level. The number of spurious levels is n, then loop pp from 1 to n

[0104] a. Calculate a series of parameters

[0105]

[0106]

[0107] b. If interference exists, then

[0108] i. Calculation parameters

[0109]

[0110]

[0111] i. If the two devices have the same position on the xy plane (xyDis<10 -6 )

[0112]

[0113] ii. If the two devices are not in the same location, then

[0114]

[0115]

[0116] iii. Determine the radiation status and intensity of interference. For the logic, see the section on determining the radiation status and intensity of interference.

[0117] iv. Generate a calculation result instance and fill the calculation result into the attributes of the instance model

[0118] v. Fill the calculation results and calculation input parameters in the export file according to the format. For details on the logic, see the "Export File" section.

[0119] 5.4.2.2.3 Finding Device Gain Logic

[0120] Input parameters: device model, device theta, device frequency F, harmnum (xbNum)

[0121] Output: Device gain

[0122] 1. Based on the device model attributes, multiple device gain instances will be found within the range

[0123] 2. Determine the frequency ifreq in the gain table corresponding to F

[0124] a. If the frequency attribute value of these device gain instances is the same, ifreq is this value

[0125] b. If multiple device gain instances have multiple frequency values, sort them from small to large, such as q, e, z, to determine the device frequency, and arbitrarily select two adjacent frequency points m and n.

[0126] i. If m≤F≤n

[0127] 1. If the distance between m and F is closer - (nF)> (Fm), then take the gain corresponding to m, ifreq = m

[0128] 2. If n and F are closer, (nF) < (Fm), then take the gain corresponding to n, ifreq = n

[0129] 3. If the distances between the two are the same, (nF) = (Fm), then take the gain corresponding to n, ifreq = n

[0130] ii. Otherwise, the interrupt error gain cannot be found

[0131] 3. Filter out all device gain instances whose frequency attribute value is ifreq

[0132] c. Sort the gain instances in ascending order according to the pitch scan angle attribute value, where the two adjacent values ​​are theta1 and theta2 respectively.

[0133] i. If theta1≤theta≤theta2

[0134]

[0135] 2. Equipment gain = theta1 corresponding gain + (theta2 corresponding gain - theta1 corresponding gain) × thetaWeight

[0136] ii. Otherwise, the minimum value of the pitch scan angle is mintheta, and the maximum value is maxtheta.

[0137] 1.thetaFound=0

[0138] 2. Determine whether mintheta or maxtheta corresponds to the gain

[0139] a. If |theta-mintheta|<|theta-maxtheta|, then the device gain

[0140] gain=mintheta

[0141] b. Otherwise, the device gain gain = maxtheta

[0142] 4. If harmNum>1, then the device gain gain=gain-10

[0143] 5.4.2.2.4 Determine the interference exposure status and intensity

[0144] 1. Determine the elevation status (FYillStat)

[0145]

[0146]

[0147] 2. Determine the azimuth illumination status 1 (FWillStat1)

[0148]

[0149] 3. Master-to-master reclassification of azimuth illumination status, calculation of azimuth illumination status 2 (FWillStat2) and coupling amount 2 (recPower2)

[0150] a. If FWillStat1 is "Master", then

[0151]

[0152]

[0153] b. Otherwise

[0154] Calculation parameters Variable Name formula Azimuth illumination state 2 FWillStat2 FWillStat2 is empty Coupling value 2 recPower2 recPower2 is also empty

[0155] 4. Determination of interference degree

[0156] a. Determine whether the coupling value is less than the burnout level (recPower <dev2SH)

[0157] i. If not, calculate the average power corresponding to the coupling value

[0158] 1. If less than the burnout level, GRCD = "intensity", judged = 1

[0159] 2. Otherwise GRCD = "burned"

[0160] ii. If yes, then

[0161] 1. Determine whether the azimuth illumination status 1 is master-to-master (FWillStat1 is "master-master")

[0162] a. If not, determine whether the coupling value is less than the saturation level

[0163] (recPower <dev2BH)

[0164] If yes, then judged = 0

[0165] If not, GRCD = "strength", judged = 1

[0166] b. If yes, determine whether the coupling value is less than the saturation level

[0167] (recPower <dev2BH)

[0168] i. If yes, then judged = 0

[0169] ii. If not, determine whether the azimuth irradiation state only has master-to-master (FWFWillStat2 is empty)

[0170] 1. If so, GRCD = "Intensity", judged = 1

[0171] 2. If not (there is a secondary), determine whether the coupling value 2 is less than the saturation level (recPower2 < dev2BH)

[0172] a. If so, judged = 0

[0173] b. If not, GRCD = "Intensity", judged = 1

[0174] 2. Determine whether GRCD has not been set yet (judged == 0)

[0175] a. If so, determine the interval where the coupling amount is located

[0176] i. If recPower ≥ dev2LMD + 20, then GRCD = "Moderate"

[0177] ii. If dev2LMD ≤ recPower < dev2LMD + 20, then GRCD = "Weak"

[0178] iii. In other cases, GRCD = "None"

[0179] b. If not, no operation

[0180] 5. If GRCD ≠ "None", calculate the parameters

[0181]

[0182] 5.4.2.3 Export files

[0183] 1. Export two excel files

[0184] a. The first one

[0185] i. The name is 'Detailed Interference Situation'. If the name is repeated, automatically add a number to the name for renaming

[0186] ii. Table headers and content: (Only export interfering devices and sensitive devices with interference)

[0187]

[0188]

[0189]

[0190] b. The second

[0191] i. The name is 'Disturbance Brief'. If the name is repeated, it will be automatically renamed by adding a number after the name.

[0192] ii. Sort by the interfered devices:

[0193] 1. According to GRCD ranking: "burned" > "strong" > "moderate" > "weak"

[0194] 2. Sort by interference type: main frequency > sideband > harmonic > spurious

[0195] iii. Table header and content description

[0196]

[0197]

[0198] 2. Export a word document

[0199] a. Page Setup

[0200]

[0201] b. Output multiple tables

[0202] i. One table for each interfered device

[0203] ii. Table header and content:

[0204]

[0205]

[0206] The above embodiments are only exemplary embodiments of the present invention and are not intended to limit the present invention. The protection scope of the present invention is defined by the claims. Those skilled in the art may make various modifications or equivalent substitutions to the present invention within the essence and protection scope of the present invention, and such modifications or equivalent substitutions shall also be deemed to fall within the protection scope of the present invention.

Claims

1. A method for rapid calculation and analysis of electromagnetic interference between radio systems, characterized in that: The steps include: Step S1, selecting the interfering device dev1 and the interfered device dev2 from the device list in order and reading the attribute values ​​of the selected devices; Step S2, determine whether the interfering device dev1 and the interfered device dev2 belong to the same device group. If yes, it is determined that there is no interference between the two groups of devices, no calculation and analysis results are generated, and the current loop is exited and returned to step S1 for the next loop; If not, proceed to step S3; Step S3, determining the values ​​of the pre-parameters, wherein the pre-parameters include the coordinate difference (dx, dy) of dev2 relative to dev1, the azimuth angle angRhoTR of the vector pointing from dev1 to dev2, the azimuth angle angRhoRT of the vector pointing from dev2 to dev1, the distance r between the antennas of the two devices, the elevation scanning angle theta, the azimuth first sidelobe relative to the main lobe level difference deltaGainT of the transmitting device, and the azimuth first sidelobe relative to the main lobe level difference deltaGainR of the receiving device; Step S4, determining the sideband interference, main frequency interference and harmonic interference between the two selected devices; Step S5, determining the spurious interference between the two selected devices; Step S7, determine whether the cycle is finished, If yes, the calculation result is displayed; If not, then exit this loop and proceed to the next loop.

2. The method for rapid calculation and analysis of electromagnetic interference between radio systems according to claim 1, characterized in that: Between step S5 and step S7, the following is also included: Step S6, exporting the calculation results and storing them in a designated location.

3. The method for rapid calculation and analysis of electromagnetic interference between radio systems according to claim 2, characterized in that: The calculation results are exported in the form of Excel and Word.

4. The method for rapid calculation and analysis of electromagnetic interference between radio systems according to claim 2, characterized in that: The calculation results all use the interfered device as the derived index.

5. The method for rapid calculation and analysis of electromagnetic interference between radio systems according to claim 1, characterized in that: The step S4 comprises: Step 541, determining the maximum harmonic number maxXbNum; Step S42, cycle the harmonic number xbNum from 1 to the maximum harmonic number maxXbNum; Step S43, determine whether the xbNum harmonic level exists. If yes, proceed to step S44; If not, then jump out of this loop and enter the next loop; Step S44, determining the values ​​of interference-related parameters, wherein the interference-related parameters include the lower limit of the harmful frequency band harmFS, the upper limit of the harmful frequency band harmFE, whether interference exists harmExiSt, the lower limit of the interference frequency band FS, the upper limit of the interference frequency band FE, the interference calculation frequency harmFreq, the transmission sideband attenuation tBDSJ and the reception sideband attenuation rBDSJ, wherein harmExist is a Boolean type parameter, when its value is 1, it indicates that interference exists, and when its value is 0, it indicates that interference does not exist; Step S45, judging whether interference exists according to interference-related parameters, If not, then jump out of this loop and enter the next loop; If yes, determine the interference situation result parameter value, wherein the interference situation result parameter includes the interference band WaveBand, the interference type type, the spurious signal zslD and the ratio of the interference bandwidth to the receiving bandwidth ibVsRb; Step S46, determining whether the interfering device dev1 and the interfered device dev2 are at the same position on the xy plane. If yes, the field strength fld Strength of the interfering device at the interfered device and the coupling value rec Power are assigned to infinity∞; If not, then determine the values ​​of the gain dev1Gt of the interfering device, the gain dev2Gr of the interfered device, the polarization mismatch loss polLoss, the field strength fldStrength of the interfering device at the interfered device, and the coupling value recPower; Step S47, determining the interference illumination state and intensity.

6. The method for rapid calculation and analysis of electromagnetic interference between radio systems according to claim 1, characterized in that: The step S5 comprises: Step S51, determining the number n of spurious levels; Step S52, reading the ppth spurious level from 1 to n in a loop; Step S53, determining the values ​​of interference-related parameters, wherein the interference-related parameters include harm FS, harm FE, whether interference exists, harm Exist, FS, FE, FE, harm Freq, and rBDSJ, wherein harm Exist is a Boolean parameter, and when its value is 1, it indicates that interference exists, and when its value is 0, it indicates that interference does not exist; Step S54: judging whether interference exists according to interference-related parameters. If not, then jump out of this loop and enter the next loop; If yes, determine the interference situation result parameter value, wherein the interference situation result parameter includes the interference band WaveBand, the interference type type, the spurious signal zslD and the ratio of the interference bandwidth to the receiving bandwidth ibVsRb; Step S56, determine whether the interfering device dev1 and the interfered device dev2 are in the same position on the xy plane. If yes, the field strength fld Strength of the interfering device at the interfered device and the coupling value rec Power are assigned to infinity∞; If not, determine the gain dev1Gt of the interfering device, the gain dev2Gr of the interfered device, the transmission sideband attenuation tBDSJ, the polarization mismatch loss polLoSs, the field strength fldStrength of the interfering device at the interfered device, and the coupling value recPower; Step S57, determining the interference illumination state and intensity.

7. The method for rapid calculation and analysis of electromagnetic interference between radio systems according to claim 5, characterized in that: The value of the maximum harmonic order maxXbNum is 3.

Images