CN112668213A - Rapid simulation analysis method for multi-scale antenna array - Google Patents

Rapid simulation analysis method for multi-scale antenna array Download PDF

Info

Publication number
CN112668213A
CN112668213A CN202011241448.4A CN202011241448A CN112668213A CN 112668213 A CN112668213 A CN 112668213A CN 202011241448 A CN202011241448 A CN 202011241448A CN 112668213 A CN112668213 A CN 112668213A
Authority
CN
China
Prior art keywords
equivalent
operator
plane
current
array
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011241448.4A
Other languages
Chinese (zh)
Inventor
陶诗飞
王意钊
王昊
李男
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Xinnuo Communication Technology Co ltd
Original Assignee
Nanjing University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202011241448.4A priority Critical patent/CN112668213A/en
Publication of CN112668213A publication Critical patent/CN112668213A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention discloses a rapid simulation analysis method of a multi-scale antenna array. The method comprises the following steps: firstly, establishing an equivalent surface outside each array element structure, solving the radiation problem of the array elements by using a region decomposition method based on an equivalent principle to obtain the relation between equivalent incident electromagnetic flow and equivalent scattering electromagnetic flow, and calculating an equivalent operator S; then calculating a transfer operator T of the interaction of the equivalent surfaces of the two array elements, and expanding to the calculation of the interaction problem among the multiple array elements to obtain a specific equation for solving the equivalent radiation electromagnetic flow on the multiple equivalent surfaces; and finally, accelerating filling of the equivalent operator S of the electrically large equivalent surface by using a self-adaptive cross algorithm. The method is simple to implement, good in robustness, capable of reducing data calculation amount, saving iterative solution time and array analysis time, and capable of improving the efficiency and accuracy of electromagnetic characteristic analysis of the multi-scale array antenna.

Description

Rapid simulation analysis method for multi-scale antenna array
Technical Field
The invention belongs to the technical field of antenna array rapid simulation analysis, and particularly relates to a rapid simulation analysis method of a multi-scale antenna array.
Background
Wireless communication technology is continuously developed, antennas are an essential part of communication systems, and the performance of antennas required by communication technology is continuously improved. For the current electromagnetic problem, the integral equation method has been combined with a fast algorithm to achieve great success, and the fast and accurate solution of the electrically large complex target can be realized. However, problems in practical engineering tend to have multi-scale characteristics. For electrically large and smooth structures, a coarse uniform mesh division can be used, while for electrically small fine structures, a very fine mesh division must be used. For coarser, uniform mesh parts, the mesh size is comparable to the order of the electromagnetic wave wavelength, so that the wave physics plays a major role in this respect. For the fine mesh portion, the mesh size is much smaller than the wavelength of the electromagnetic wave, so the circuit physics plays a major role in this respect. Due to the fact that the eigenvalue distribution and the eigenvalue vector of the two physical processes are greatly different, the finally obtained impedance matrix is ill-conditioned, and the traditional iterative method is difficult to converge or even does not converge when being used for solving.
Disclosure of Invention
The invention discloses an antenna array analysis method which is suitable for rapid analysis of various multi-scale large-scale antenna arrays, simple in method, high in calculation efficiency and high in accuracy.
The technical solution for realizing the purpose of the invention is as follows: a rapid simulation analysis method of a multi-scale antenna array comprises the following steps:
step 1, establishing an equivalent surface outside each array element structure, solving the radiation problem of the array elements by using a region decomposition method based on an equivalent principle to obtain the relation between equivalent incident electromagnetic flow and equivalent scattering electromagnetic flow, and calculating an equivalent operator S;
step 2, calculating a transfer operator T of the interaction of the equivalent surfaces of the two array elements, and expanding to the calculation of the interaction problem among the multiple array elements to obtain a specific equation for solving the equivalent radiation electromagnetic flow on the multiple equivalent surfaces;
and 2, accelerating filling of the equivalent operator S of the electrically large equivalent surface by using a self-adaptive cross algorithm.
Compared with the prior art, the invention has the following remarkable advantages: (1) the self-adaptive cross algorithm is adopted, the method is not limited by a Green function expression form, the realization is simple, and the robustness is good; (2) by using the repeatability of the array unit and the translation invariance of the Green function, the same unit type only needs to be calculated once, so that the data calculation amount is reduced; (3) the fine structure is uniformly subdivided by using the equivalent surface, so that the matrix behavior and the iterative convergence effect are improved, and the iterative solution time and the array analysis time are saved; (4) the matrix compression is accelerated for the electrically large-sized structure, the array analysis efficiency is improved, and the efficient and accurate analysis of the electromagnetic characteristics of the multi-scale array antenna is realized.
Drawings
Fig. 1 is a schematic flow chart of a multi-scale antenna array rapid simulation analysis method according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and specific embodiments.
The equivalent principle is also called as the Huygens principle, and means that a field surrounding each point on a closed surface of a wave source can be regarded as a secondary wave source to radiate again to the outside of the closed surface, and an electromagnetic field of any point inside or outside the closed surface is generated by all the electromagnetic fields on the closed surface together. All wave sources are arranged in or out of the closed surface S, and an infinite closed surface S is formed at infinityAnd a passive region is formed between the two closed surfaces, the electric field and the magnetic field of any point in the passive region can be obtained by the following formula:
Figure BDA0002768487480000021
wherein equivalent current on the closed side
Figure BDA0002768487480000022
Electromagnetic current
Figure BDA0002768487480000023
G (r, r ') is a free space Green function, and r' represents a field vector of any point; define integration operator L, K:
Figure BDA0002768487480000024
then, a matrix relationship between the electromagnetic field and the equivalent surface electromagnetic current can be obtained:
Figure BDA0002768487480000025
the traditional integral equation method can directly carry out geometric modeling and grid dispersion on the target surface by utilizing an equivalent principle, and establish an equation on the target surface. For practical complex targets, the number of grids is often very large and uneven, which directly results in huge unknowns of matrix equations to be solved and difficult iterative convergence. The equivalent principle algorithm is that unknown quantity on a target is transferred to a virtual closed equivalent surface, the equivalent surface is regular and smooth, so that subdivision grids are uniform and regular, the finally obtained matrix equation is small in unknown quantity, good in performance and fast in convergence, iterative solution time and array analysis time are saved, and rapid simulation analysis of the electromagnetic characteristics of the array antenna can be achieved.
With reference to fig. 1, the present invention provides a method for rapid simulation analysis of a multi-scale antenna array, which includes the following steps:
step 1, establishing an equivalent surface outside each array element structure, solving the radiation problem of the array elements by using a region decomposition method based on an equivalent principle to obtain the relation between equivalent incident electromagnetic flow and equivalent scattering electromagnetic flow, and calculating an equivalent operator S, wherein the method specifically comprises the following steps:
step 1.1, the radiation problem of the antenna array element is specifically divided into: (1) spreading from outside to inside; (2) scattering on the internal metal; (3) spread from the inside to the outside.
(1) From outside to inside, incident field Einc,HincFirstly, the incident electromagnetic current is irradiated on the equivalent surface S, the equivalent incident electromagnetic current is generated on the equivalent surface S, a source field outside the equivalent surface is replaced, the equivalent incident electromagnetic current on the equivalent surface excites the most original incident electric field and magnetic field inside the equivalent surface, and the electric field and magnetic field outside the equivalent surface become zero, which is a zero-field equivalent principle, and the following formula can show the principle of the equivalent surfaceThe process is as follows:
Figure BDA0002768487480000031
(2) the scattering effect on the internal metal, and the electromagnetic current on the equivalent plane interior can be solved by a moment method.
(3) From inside to outside, the electromagnetic current of the inner target surface is used for solving the radiation electromagnetic current on the equivalent surface, at the moment, the radiation electromagnetic current on the equivalent surface excites a zero field inside the equivalent surface, an original radiation electromagnetic field is excited outside the equivalent surface, and the zero field equivalence principle is also adopted.
Step 1.2, the three processes can be written into a matrix equation form, and the RWG basis functions are used for dispersing the current and the magnetic current on the equivalent surface and the internal target.
Current J on the internal targetaCan be dispersed as:
Figure BDA0002768487480000032
wherein f isai(r) represents RWG basis function, jaiRepresenting the current coefficient on the internal target, Na represents the fraction.
Step 1.3, by combining the above analysis, the total matrix equation for solving the single target radiation by using the equivalence principle can be obtained as follows:
Figure BDA0002768487480000033
wherein the subscript S represents an equivalent plane,
Figure BDA0002768487480000034
and
Figure BDA0002768487480000035
for an equivalent incident current and an equivalent incident magnetic current,
Figure BDA0002768487480000036
and
Figure BDA0002768487480000037
for equivalent scattering currents and equivalent incident magnetic currents,
Figure BDA0002768487480000038
denotes the normal vector on the equivalent plane S, the subscript p denotes the object relating to the operator operation within the equivalent plane, [ L ]pp]-1The process of solving the current on the equivalent plane internal target is represented, and L and K are integral operators:
Figure BDA0002768487480000041
where G (r, r ') is a free space Green's function, k0And η are the free space wavenumber and wave impedance, respectively.
As can be seen from equation (6), the unknowns of each sub-region are transferred to the equivalent plane, and the information of each sub-region is stored in the S matrix.
The equivalent radiation electromagnetic current on the final equivalent surface can be solved by the matrix equation, and the radiation electromagnetic field outside the equivalent surface can be obtained by the radiation electromagnetic current on the equivalent surface.
Step 2, calculating a transfer operator T of the interaction of the equivalent surfaces of the two array elements, and expanding to the calculation of the interaction problem among the multiple array elements to obtain a specific equation for solving the equivalent radiation electromagnetic flow on the multiple equivalent surfaces;
step 2.1, when analyzing the electromagnetic problem of the multi-scale target or the periodically repeated target, the action between any two single targets can be replaced by the interaction between the equivalent surfaces, it is assumed that only two solution sub-areas are surrounded by the two equivalent surfaces, the electromagnetic current on the first equivalent surface will generate a radiation action on the second equivalent surface, so as to generate an additional incident current and an additional incident magnetic current on the second equivalent surface (where the additional incident current and the incident magnetic current represent that the additional incident current and the incident magnetic current generated by the radiation action of the first equivalent surface after removing the incident current and the incident magnetic current generated on the equivalent surface by the original incident electric field and the original incident magnetic field), and the current and the magnetic current on the second equivalent surface will also generate a radiation action on the first equivalent surface, therefore, extra incident current and incident magnetic current are generated on the first equivalent surface, the closer the two sub-regions are, the stronger the interaction is, but the equivalent electromagnetic current on each equivalent surface cannot generate self-action on the equivalent surface per se, namely the equivalent surface per se has no self-action process, and only has the interaction process between different equivalent surfaces, because the self-action process of the equivalent surface is actually the filling process of the S matrix in the equivalent principle algorithm. Therefore, the T-transition matrix, which represents the interaction between the equivalent surface and the equivalent surface, can be defined as follows:
Figure BDA0002768487480000042
wherein the superscript h is the equivalent plane.
Step 2.2, the S operator can be used for processing
Figure BDA0002768487480000043
Converting into equivalent radiation electromagnetic current to be solved
Figure BDA0002768487480000044
Obtaining a specific equation for solving equivalent radiation electromagnetic currents on two equivalent surfaces:
Figure BDA0002768487480000051
wherein I is an identity matrix, and I is an identity matrix,
Figure BDA0002768487480000052
for equivalently radiating electromagnetic currents
Figure BDA0002768487480000053
Figure BDA0002768487480000054
For equivalent incident electromagnetic flow
Figure BDA0002768487480000055
S11Is the S equivalent operator on the equivalent plane 1, and so on.
Step 2.3, if M equivalent surfaces are included in the problem, equation (8) can be expanded to the following form:
Figure BDA0002768487480000056
step 3, performing accelerated filling on the equivalent operator S of the electrically large-size equivalent surface by using a self-adaptive cross algorithm, which specifically comprises the following steps:
step 3.1, grouping the targets to be solved, establishing a tree structure, and dividing the original matrix into a plurality of sub-matrices with different sizes on the basis of the tree structure;
and 3.2, directly calculating the submatrices formed by the self-acting groups and the adjacent two groups by using a moment method, and quickly filling the submatrices between other non-adjacent groups by using a self-adaptive cross algorithm.
The impedance matrix obtained by the moment method through the equivalent surface subdivision solution is a full-rank matrix, but the impedance matrix comprises a plurality of low-rank sub-matrix blocks due to the property of the Green function. If the quantities to be requested on the target are divided into groups, the interactions between the groups that are further apart have a low rank characteristic and can be compressed. An acceleration algorithm based on such idea is called a low rank type compression algorithm.
For the target with larger electrical size, the region decomposition algorithm based on the equivalence principle will encounter problems, the area of the equivalent plane may be larger than the surface area of the internal metal, and finally the unknown quantity on the equivalent plane will be larger than the total unknown quantity. An Adaptive Crossover Algorithm (ACA), a widely used low rank algorithm, was therefore introduced. The ACA method is a pure algebraic method, is different from a fast multipole, is not limited by a Green function expression form, and is simple to implement and good in robustness.
It is advantageous to use the EPA algorithm to solve for electromagnetic radiation characteristics of multi-scale targets with fine structures, or periodically repeating targets. Especially when analyzing a multi-scale target with a fine structure, the fine structure can approach an original object more truly only by adopting a denser subdivision size, so that the unknown quantity is larger, the equivalent surface shape is regular, and the larger subdivision size can be adopted, so that the unknown quantity on the equivalent surface is greatly reduced compared with the unknown quantity on the surface of the internal fine structure target, the interaction between calculation regions is converted into the interaction between calculation of the equivalent surface and the equivalent surface surrounding each region, and the same unit type only needs to be calculated once by utilizing the repeatability of an array unit and the translation invariance of a Green function; for the structure with large electrical size, the matrix can be greatly compressed by using an Adaptive Cross Algorithm (ACA), and the original matrix can be approximately obtained only by sampling partial elements of the original matrix.
For a fine structure, the matrix performance is improved by utilizing the uniform subdivision of an equivalent surface, the iterative convergence effect is improved, the iterative solving time and the array analysis time are saved, the matrix compression is accelerated for an electrically large-size structure, and finally the efficient and accurate analysis of the electromagnetic characteristics of the multi-scale array antenna is realized.

Claims (4)

1. A rapid simulation analysis method of a multi-scale antenna array is characterized by comprising the following steps:
step 1, establishing an equivalent surface outside each array element structure, solving the radiation problem of the array elements by using a region decomposition method based on an equivalent principle to obtain the relation between equivalent incident electromagnetic flow and equivalent scattering electromagnetic flow, and calculating an equivalent operator S;
step 2, calculating a transfer operator T of the interaction of the equivalent surfaces of the two array elements, and expanding to the calculation of the interaction problem among the multiple array elements to obtain a specific equation for solving the equivalent radiation electromagnetic flow on the multiple equivalent surfaces;
and 3, accelerating filling of the equivalent operator S of the electrically large-size equivalent surface by using a self-adaptive cross algorithm.
2. The method for rapid simulation analysis of a multi-scale antenna array according to claim 1, wherein step 1 establishes an equivalent plane outside each array element structure, solves the radiation problem of the array elements by using a region decomposition method based on an equivalent principle, obtains the relationship between an equivalent incident electromagnetic current and an equivalent scattering electromagnetic current, and calculates an equivalent operator S as follows:
step 1.1, establishing an equivalent surface outside each array element structure;
step 1.2, solving the radiation problem of the array element by using a regional decomposition method based on an equivalent principle to obtain the relation between equivalent incident electromagnetic current and equivalent scattering electromagnetic current;
step 1.3, obtaining a total matrix equation for solving the radiation of the single target by using the equivalent principle as follows:
Figure FDA0002768487470000011
wherein the subscript S represents an equivalent plane,
Figure FDA0002768487470000012
and
Figure FDA0002768487470000013
for an equivalent incident current and an equivalent incident magnetic current,
Figure FDA0002768487470000014
and
Figure FDA0002768487470000015
are equivalent scattering currents and equivalent scattering magnetic currents,
Figure FDA0002768487470000016
denotes the normal vector on the equivalent plane S, the subscript p denotes the object relating to the operator operation within the equivalent plane, [ L ]pp]-1Representing the process of solving the current on the internal target of the equivalent plane; l and K are integral operators:
Figure FDA0002768487470000017
wherein G (r, r ') is a free space Green's function,
Figure FDA0002768487470000018
as gradient operator, k0And η are the free space wavenumber and wave impedance, respectively.
3. The method for rapid simulation analysis of a multi-scale antenna array according to claim 2, wherein the step 2 of calculating the transfer operator T of the interaction between the equivalent surfaces of the two array elements is extended to calculating the interaction problem between the plurality of array elements, so as to obtain a specific equation for solving the equivalent radiation electromagnetic flow on the plurality of equivalent surfaces, which is specifically as follows:
step 2.1, solving a transfer operator of interaction between two array elements, wherein the formula is as follows:
Figure FDA0002768487470000021
wherein the content of the first and second substances,
Figure FDA0002768487470000022
representing the equivalent radiation current and magnetic current of the array element 1;
Figure FDA0002768487470000023
representing the equivalent radiation current and magnetic current, T, of the array element 1hhRepresenting a transfer operator between equivalent surfaces of the array elements 1 and 2;
step 2.2, utilize S operator to be with
Figure FDA0002768487470000024
Converting into equivalent radiation electromagnetic current to be solved
Figure FDA0002768487470000025
Obtaining a specific equation for solving equivalent radiation electromagnetic currents on two equivalent surfaces:
Figure FDA0002768487470000026
wherein, I is an identity matrix,
Figure FDA0002768487470000027
for an equivalent radiated electromagnetic current on the equivalent plane 1,
Figure FDA0002768487470000028
for an equivalent radiated electromagnetic current on the equivalent plane 2,
Figure FDA0002768487470000029
for an equivalent incident electromagnetic flow on the equivalent plane 1,
Figure FDA00027684874700000210
for an equivalent incident electromagnetic current, S, on the equivalent plane 211Is an equivalence operator on an equivalence plane 1, S22Is the equivalence operator on the equivalence plane 2; t is12、T21Respectively, a transfer operator of the interaction between the equivalent surface 1 and the equivalent surface 2;
step 2.3, if the problem contains M equivalent surfaces, the formula (8) is expanded to the following form:
Figure FDA00027684874700000211
in the formula, SMMIs an equivalent operator on the equivalent plane M; t isM2、T2MRespectively, the equivalent plane M and the equivalent plane 2 are mutually acted transfer operators; t isM1、T1MRespectively, the transfer operators of the interaction between the equivalent plane M and the equivalent plane 1.
4. The method for rapid simulation analysis of a multi-scale antenna array according to claim 3, wherein the equivalent operator S for the electrically large equivalent surface in step 3 is filled in with an adaptive cross algorithm at an accelerated speed, specifically as follows:
step 3.1, grouping the targets to be solved, establishing a tree structure, and dividing the original matrix into a plurality of sub-matrices with different sizes on the basis of the tree structure;
and 3.2, directly calculating the submatrices of the self-acting groups and the submatrices formed by the interaction between the two adjacent groups by using a moment method, and quickly filling the submatrices between the other non-adjacent groups by using a self-adaptive cross algorithm.
CN202011241448.4A 2020-11-09 2020-11-09 Rapid simulation analysis method for multi-scale antenna array Pending CN112668213A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011241448.4A CN112668213A (en) 2020-11-09 2020-11-09 Rapid simulation analysis method for multi-scale antenna array

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011241448.4A CN112668213A (en) 2020-11-09 2020-11-09 Rapid simulation analysis method for multi-scale antenna array

Publications (1)

Publication Number Publication Date
CN112668213A true CN112668213A (en) 2021-04-16

Family

ID=75404056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011241448.4A Pending CN112668213A (en) 2020-11-09 2020-11-09 Rapid simulation analysis method for multi-scale antenna array

Country Status (1)

Country Link
CN (1) CN112668213A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115935671A (en) * 2022-12-20 2023-04-07 安徽大学 Regional decomposition electromagnetic simulation method and system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邵汉儒: "等效原理算法在三维电磁问题中的分析及其应用", 《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115935671A (en) * 2022-12-20 2023-04-07 安徽大学 Regional decomposition electromagnetic simulation method and system
CN115935671B (en) * 2022-12-20 2023-08-18 安徽大学 Regional decomposition electromagnetic simulation method and system

Similar Documents

Publication Publication Date Title
CN102129523B (en) Method for analyzing electromagnetic scattering of complex target through MDA and MLSSM
Gies et al. Particle swarm optimization for reconfigurable phase‐differentiated array design
Guo et al. Pattern synthesis of conformal array antenna in the presence of platform using differential evolution algorithm
CN112329303B (en) Array antenna electromagnetic characteristic solving method based on finite element region decomposition
CN102081690B (en) MDA (Matrix Decomposition Algorithm)-combined novel SVD (Singular Value Decomposition) method for complex circuit
CN106294894B (en) Finite element boundary integration method for rapidly analyzing electromagnetic scattering characteristics of non-uniform target
CN110442919A (en) A kind of microwave component micro discharge method for numerical simulation based on GPU architecture
Burkholder et al. Fast dual-MGS block-factorization algorithm for dense MoM matrices
CN112668213A (en) Rapid simulation analysis method for multi-scale antenna array
Li et al. Increasing radiation efficiency using antenna shape optimization approach
Raji et al. Fast optimization of sparse antenna array using numerical Green's function and genetic algorithm
CN116562143A (en) Antenna topology and parameter mixing optimization method based on normalized Gaussian network
Chu et al. A multilevel fast multipole algorithm for electrically small composite structures
Dang et al. A fast modeling algorithm for quasi-periodic array
CN106650094A (en) Method for calculating tolerance of microwave device design parameters
Wang et al. An adaptive nested complex source beam method for electromagnetic scattering of composite conducting-dielectric objects
CN104699876B (en) The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna
CN103150415A (en) Full wave analysis method for microstrip circuit
Ding et al. Boundary‐based optimization strategies for design of fragment‐type isolation structure in compact MIMO systems
Ergul et al. Fast and accurate solutions of large-scale scattering problems with parallel multilevel fast multipole algorithm
Prabhakar et al. A hybrid approach on metamaterial-loaded fractal antenna design
Zhang et al. Adaptive Multi-Group Differential Evolution Optimization for MIMO Antenna Design
Ayari On the Efficiency of the Advanced TWA Approach to the 60-GHz Microstrip Antenna Analysis for 5G Wireless Communication Systems
Hu et al. Parallel computing graphene frequency selective surface (GFSS) with large finite array using HIE-FDTD method on high performance computer
CN118152706A (en) CSIE electromagnetic scattering analysis method based on characteristic basis function

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20230605

Address after: E701, No. 388 Ruoshui Road, Suzhou Industrial Park, Suzhou Area, China (Jiangsu) Pilot Free Trade Zone, 215125, Jiangsu Province

Applicant after: Suzhou Xinnuo Communication Technology Co.,Ltd.

Address before: 210094 Xuanwu District, Jiangsu, Xiaolingwei 200, Nanjing

Applicant before: NANJING University OF SCIENCE AND TECHNOLOGY

TA01 Transfer of patent application right
RJ01 Rejection of invention patent application after publication

Application publication date: 20210416

RJ01 Rejection of invention patent application after publication