CN112668213A - Rapid simulation analysis method for multi-scale antenna array - Google Patents
Rapid simulation analysis method for multi-scale antenna array Download PDFInfo
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 24
- 238000004088 simulation Methods 0.000 title claims abstract description 13
- 238000000034 method Methods 0.000 claims abstract description 32
- 230000005855 radiation Effects 0.000 claims abstract description 29
- 230000003993 interaction Effects 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 18
- 238000012546 transfer Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims 1
- 238000004613 tight binding model Methods 0.000 description 6
- 230000005672 electromagnetic field Effects 0.000 description 5
- 230000009471 action Effects 0.000 description 4
- 230000005684 electric field Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000005405 multipole Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
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
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 infinity∞And 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:
wherein equivalent current on the closed sideElectromagnetic currentG (r, r ') is a free space Green function, and r' represents a field vector of any point; define integration operator L, K:
then, a matrix relationship between the electromagnetic field and the equivalent surface electromagnetic current can be obtained:
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:
(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:
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:
wherein the subscript S represents an equivalent plane,andfor an equivalent incident current and an equivalent incident magnetic current,andfor equivalent scattering currents and equivalent incident magnetic currents,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:
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:
wherein the superscript h is the equivalent plane.
Step 2.2, the S operator can be used for processingConverting into equivalent radiation electromagnetic current to be solvedObtaining a specific equation for solving equivalent radiation electromagnetic currents on two equivalent surfaces:
wherein I is an identity matrix, and I is an identity matrix,for equivalently radiating electromagnetic currents For equivalent incident electromagnetic flowS11Is 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:
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:
wherein the subscript S represents an equivalent plane,andfor an equivalent incident current and an equivalent incident magnetic current,andare equivalent scattering currents and equivalent scattering magnetic currents,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:
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:
wherein the content of the first and second substances,representing the equivalent radiation current and magnetic current of the array element 1;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 withConverting into equivalent radiation electromagnetic current to be solvedObtaining a specific equation for solving equivalent radiation electromagnetic currents on two equivalent surfaces:
wherein, I is an identity matrix,for an equivalent radiated electromagnetic current on the equivalent plane 1,for an equivalent radiated electromagnetic current on the equivalent plane 2,for an equivalent incident electromagnetic flow on the equivalent plane 1,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:
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.
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115935671A (en) * | 2022-12-20 | 2023-04-07 | 安徽大学 | Regional decomposition electromagnetic simulation method and system |
-
2020
- 2020-11-09 CN CN202011241448.4A patent/CN112668213A/en active Pending
Non-Patent Citations (1)
Title |
---|
邵汉儒: "等效原理算法在三维电磁问题中的分析及其应用", 《中国优秀博硕士学位论文全文数据库(博士) 基础科学辑》 * |
Cited By (2)
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 |