CN110516336A - Directional diagram reconstructable pixel antenna optimization method based on built-in multiport algorithm - Google Patents
Directional diagram reconstructable pixel antenna optimization method based on built-in multiport algorithm Download PDFInfo
- Publication number
- CN110516336A CN110516336A CN201910766927.9A CN201910766927A CN110516336A CN 110516336 A CN110516336 A CN 110516336A CN 201910766927 A CN201910766927 A CN 201910766927A CN 110516336 A CN110516336 A CN 110516336A
- Authority
- CN
- China
- Prior art keywords
- value
- ant
- cross
- directional diagram
- matrix
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention proposes a kind of directional diagram reconstructable pixel antenna optimization method based on built-in multiport algorithm, for solving the poor technical problem of the non-reconstruct segmental stability of reconfigurable antenna existing in the prior art.Realize step are as follows: determine the structure and optimization object of antenna to be optimized, the radiation characteristic of antenna to be optimized is calculated using built-in multiport algorithm, using the structural parameters of particle swarm algorithm optimization antenna, using binary strings genetic algorithm optimization duplexer state distribution, the directional diagram reconstructable pixel antenna that optimization is completed is obtained.This method improves the stability of the non-reconstruct part of reconfigurable antenna, has higher engineering practical value by increasing the optimization process of structural parameters on the basis of guaranteeing optimization efficiency same as the prior art.Present invention can apply to the optimization designs of directional diagram reconstructable pixel antenna.
Description
Technical field
The invention belongs to antenna technical fields, are related to a kind of antenna optimization method, and in particular to one kind is based on built-in multiterminal
The directional diagram reconstructable pixel antenna optimization method of mental arithmetic method, can be applied to the optimization of directional diagram reconstructable pixel antenna, and can
To improve the stability of the non-reconstruct part of antenna.
Background technique
The fast development of wireless communication makes wireless system increasingly complicated, in face of more complicated demand, generally require be
System can work under more scenes, and reconfigurable antenna can be under the premise of not changing overall structure by adjusting the resonance of antenna
Length and current path realize different radiation characteristics, can be competent at the need that current wireless system works for more scenes well
It asks.
The radiation characteristic of antenna mainly includes the reflection coefficient under center of antenna frequency, the polarization mode of antenna and antenna
Directional diagram scanning angle.Different for restructural antenna radiation characteristics are participated in, reconfigurable antenna can be divided into frequency to be weighed
Structure, polarize restructural and directional diagram reconstructable aerial.
Pixel antenna is that a kind of pixel patch using numerous smaller sizes forms array as the restructural day of irradiation structure
Line, to switch connection between each pixel patch, designer can be distributed by changing the on off operating mode of duplexer fast and effective
Ground changes its current path, resonance length, to obtain excellent reconfigurability energy；But restructural pixel antenna is due to using
Fairly large switch proposes huge test to design optimization process, how sufficiently to mention to realize changeable radiation characteristic
Height optimization cost, saving optimization time become for urgent problem.
Reconfigurable antenna it is non-reconstruct part stability refer to except restructural antenna radiation characteristics, to be not involved in weight
The characterization of other radiation characteristic variation degree of the antenna of structure, variation degree is lower, and stability is higher；The reason of influencing the stability
Mainly during adjusting antenna resonance length and current path, shadow is produced to the radiation characteristic of non-restructural part
It rings, the state before having caused a deviation from, or even becomes entirely different or even opposite state；It is non-heavy how reconfigurable antenna is improved
The stability of structure part is of great significance in practice in engineering.
2 months 2014 Sichao Song and Ross D.Murch are published in " IEEE TRANSACTIONS ON
Entitled " the An Efficient Approach for Optimizing of ANTENNAS AND PROPAGATION "
The paper of Frequency Reconfigurable Pixel Antennas Using Genetic Algorithms " proposes
A kind of method using built-in multiport algorithm optimization frequency reconfigurable pixel antenna, implementation method are to extract antenna to be optimized
Port parameter calculate reflection coefficient of the antenna to be optimized under assigned frequency using built-in multiport algorithm, and calculated by heredity
The switch state distribution of the restructural pixel antenna of method optimization frequency obtains the optimized switching state distribution under assigned frequency, realizes frequency
Rate reconfigurability energy.This method replaces full-wave simulation to calculate the radiation characteristic of antenna using built-in multiport algorithm, greatly improves
Optimization efficiency saves the optimization time, but the deficiency of this method is to only considered antenna reflection coefficient in optimization process
It calculates, the polarization of antenna is not analyzed with pattern characteristics, ignore switch state distribution to the shadow of non-reconstruct part
It rings, causes the stability of the non-reconstruct part of antenna poor.
Summary of the invention
It is an object of the invention to overcome the problems of the above-mentioned prior art, propose a kind of based on the mental arithmetic of built-in multiterminal
The directional diagram reconstructable pixel antenna optimization method of method, for solving the stabilization of the non-reconstruct part of antenna existing in the prior art
The poor technical problem of property.
To achieve the goals above, the technical solution adopted by the present invention includes the following steps:
(1) structure and optimization object of directional diagram reconstructable pixel antenna to be optimized are determined:
(1a) directional diagram reconstructable pixel antenna Ant to be optimized, including parasitic agent plate, the driving dielectric-slab being sequentially arranged
With feed dielectric plate, it is provided with air layer between parasitic agent plate and driving dielectric-slab, drives dielectric-slab and feed dielectric plate phase
Mutually stacking, the upper surface of parasitic agent plate is printed with the m of periodic arrangementp×npA shape is the pixel patch of rectangle, mp>=2,
np>=2, it is switched between adjacent pixel patch and is connected by rectangle, the quantity of rectangle switch is N；Drive the upper surface print of dielectric-slab
It is formed with the microband paste that shape is rectangle, microband paste is provided centrally with rectangular channel；Feed dielectric plate upper surface is printed with feedback
Electric microstrip line, lower surface are printed with metal floor, and the incline of the feed dielectric plate is equipped with the feed with feeding microstrip line connection
Port；
The optimization object of (1b) directional diagram reconstructable pixel antenna to be optimized includes structural parameters and rectangle switch on and off shape
State distribution, the quantity of structural parameters are d, length and the wide, length of rectangular channel and width, the length of pixel patch including microband paste and
The height of width and air layer；
(2) experiment vector set U needed for calculating Ant reflection coefficient and directional diagram is established:
By number 0 to number 2N- 1 is converted to corresponding 2NA N binary code, and according to sequence from small to large by 2N
A N binary code forms a line, and obtains the binary matrix U, M=2 that size is M × NN, calculating is used as after U is pressed row piecemeal
Experiment vector set needed for Ant reflection coefficient and directional diagram:
(3) impedance control vector set Z needed for calculating Ant reflection coefficient and directional diagram is established:
(3a) establishes the square of characteristic impedance when size matches for the value of 2 × 1 and each element equal to feed port in Ant
Battle array Z1；Establish equivalent impedance when size is equal to rectangle switch conduction in Ant for the value of 2 × N and the first row element, the second row member
The value of element is equal to the matrix Z of equivalent impedance when rectangle switch disconnects in Ant2；Establish the value that size is N × 1 and each element
Equal to the matrix Z of equivalent impedance when rectangle switch disconnects in Ant3；The value that size is established as N × N and diagonal entry is equal to
Equivalent impedance in Ant when rectangle switch conduction, the value of remaining all elements are equal to equivalent resistance when rectangle switch disconnects in Ant
Anti- matrix Z4；
(3b) is to matrix Z1、Z2、Z3And Z4It is combined, obtains the matrix Z that size is (N+2) × (N+1), and by Z by row
It carries out again after piecemeal as impedance control vector set needed for calculating Ant reflection coefficient and directional diagram:
(4) the directional diagram reconstructable pixel antenna Ant' to be optimized of load lump port is obtained:
The N number of rectangle connected between adjacent pixel patch in Ant switch is replaced by N number of lump port, obtains adjacent picture
The directional diagram reconstructable pixel antenna Ant' to be optimized connected by lump port is connected between mourning card piece；
(5) script function of Ant' is write:
Foundation includes the function f of input, output and main body in MATLAB software, and using f as the script letter of Ant'
Number, in which:
The input of function f is s and z, and wherein s indicates that length is equal to d, element includes the length for being arranged successively microband paste, micro-
The knot of the length of width, rectangular channel, the width of rectangular channel, the length of pixel patch, the height of the width of pixel patch and air layer with patch
Structure dominant vector；The resistance that z expression length is equal to N+1, element includes the feed port included in Ant' and being sequentially arranged
Anti-, first lump port impedance ..., the impedance control vector of the impedance of n-th lump port；
The output of function f is S, Etotal(θ)、Eco(θ) and Ecross(θ), wherein S indicates that size is equal to (N+1) × (N+1)
Full port scattering parameter matrix；Etotal(θ)、Eco(θ) and EcrossIt is θ that (θ), which respectively indicates scanning angle, and length is equal to 181
Total intensity directional diagram vector, main polarization field strength pattern vector sum cross polarization field strength pattern vector, -90≤θ≤90；
The main body of function f includes modeling code, calculation code and return code, wherein modeling code, for reading f's
Input, and the modeling functions in the control script MATLAB-HFSS-API of Calling MATLAB software, carry out abstract to Ant' and build
Mould；Calculation code carries out analysis meter to the abstract model of foundation for the analytic function in Calling MATLAB-HFSS-API
It calculates；Return code, for exporting calculation code result calculated；
(6) the optimum structure parameter of Ant' is calculated based on particle swarm optimization algorithm:
(6a) initializes the parameter of particle swarm optimization algorithm, including inertial factor w ∈ (0,1), the first Studying factors c1∈
(1,3), the second Studying factors c2∈ (1,3), the dimension D=d of particle, the quantity P ∈ (50,500) of particle, search space range
[value_min, value_max], velocity interval [v_min, v_max], the number of iterations tp, maximum number of iterations Tp >=50, most
The length of excellent fitness list F and optimal particle list R, F and R are equal to Tp, and the tp element F (tp) in F represents tp
The adaptive optimal control angle value that secondary iteration obtains, the tp element R (tp) in R represent the optimal particle that the tp times iteration obtains, and enable
Tp=1；
(6b) generates the particle mass matrix X and rate matrices V that size is P × D at random, wherein the i-th row X [i] in X
I-th of particle is represented, the i-th row V [i] in V represents the speed of i-th of particle, and element is from left to right distinguished in X [i] and V [i]
Represent the length of microband paste in Ant', the width of microband paste, the length of rectangular channel, the width of rectangular channel, the length of pixel patch, pixel patch
The width of piece and the height of air layer, and the value of each element is greater than value_min and is less than value_max in X, it is every in V
The value of one element is greater than the opposite number of value_min and is less than value_min；
(6c) enables s=X [i], z=zs, and by s and z input f, obtain the corresponding full port scattering parameter matrix of X [i]
Stp,i；
(6d) enables s=X [i], z=zE0,zE1,...,zEN, and by s and z input f, obtain the corresponding main polarization field X [i]
Strong direction set of graphs { Eco}iAnd cross polarization field strength pattern set { Ecross}i:
{Eco}i={ Eco(θ)tp,0,...,Eco(θ)tp,n,...,Eco(θ)tp,N}i,
{Ecross}i={ Ecross(θ)tp,0,...,Ecross(θ)tp,n,...,Ecross(θ)tp,N}i；
(6e) uses built-in multiport algorithm, passes through Stp,i、{Eco}i{ Ecross}iIt is corresponding to calculate each particle X [i] in X
Reflection coefficient set { SLtp}i, main polarization radiation direction set of graphs { AFco}iAnd cross polarization radiations direction set of graphs
{AFcross}i:
{SLtp}i={ SLtp,1,...,SLtp,m,...,SLtp,M}i,
{AFco}i={ AFco(θ)tp,1,...,AFco(θ)tp,m,...,AFco(θ)tp,M}i,
{AFcross}i={ AFcross(θ)tp,1,...,AFcross(θ)tp,m,...,AFcross(θ)tp,M}i；
(6f) passes through { SLtp}i、{AFco}i{ AFcross}i, calculate the ideal adaptation angle value of each particle X [i] in X
fitPSO(i), and the smallest F (tp) progress in adaptive optimal control degree list F in P calculated individual fitness value is chosen
It updates, obtains updated adaptive optimal control degree list F', choose the corresponding particle of minimum fitness value to optimal particle list R's
R (tp) is updated, and obtains updated optimal particle list:
In formula, SLtp,mFor { SLtp}iIn m-th of reflection coefficient, SL0For the reflection coefficient reference value of setting, Cptp,mIt is logical
Cross { AFco}i{ AFcross}iCorresponding m-th of the cross polarization of the X [i] of calculating inhibits ratio, Cp0Inhibit for the cross polarization of setting
Than reference value, max { } representative takes biggish element, ∑ in { } to represent sum operation:
AF in formulaco(θ)tp,mFor { AFco}iIn m-th of main polarization antenna pattern, AFcross(θ)tp,mFor { AFcross}iIn
M-th of cross polarization radiations directional diagram, θco_maxFor AFco(θ)tp,mThe corresponding scanning angle of middle maximum value；
(6g) is updated every a line of V and every a line of X respectively, obtains updated rate matrices V' and population
Matrix X', and the value for being less than value_min in X' is replaced with value_min, the value value_max greater than value_max
Replacement, obtains replaced particle mass matrix X ", wherein the calculation formula of V' and X' is respectively as follows:
V'[i]=w × V [i]+c1×random()×(R'(tp)-X[i])+c2×random()×(R'(ind)-X
[i]),
X'[i]=X [i]+V'[i],
Wherein, random () indicates to generate the random number between one 0 to 1, and ind is indicated in adaptive optimal control degree list F' most
Position corresponding to small value；
(6h) enables X=X ", enables V=V', enables F=F', enables R=R', and whether judge tp=Tp true, if so, executing step
Suddenly (6i) otherwise enables tp=tp+1, and executes step (6c)；
(6i) determines to whether there is 0 element in F, if so, using the particle in the corresponding R of 0 element as the optimal knot of Ant'
Structure parameter pbest, otherwise, using the particle in the corresponding R of least member in F as the optimum structure parameter pbest of Ant'；
(7) the corresponding full port scattering parameter matrix S of optimum structure parameter pbest is obtainedbestWith total intensity direction atlas
Close { Etotal}best:
(7a) enables s=pbest, z=zs, and by s and z input f, obtain the corresponding full port scattering parameter square of pbest
Battle array Sbest；
(7b) enables s=pbest, z=zE0,zE1,...,zEN, and by s and z input f, obtain the corresponding total intensity of pbest
Direction set of graphs { Etotal}best:
{Etotal}best={ Etotal(θ)0,...,Etotal(θ)n,...,Etotal(θ)N}best；
(8) the optimized switching on off operating mode distribution of Ant' is calculated based on binary strings genetic algorithm:
(8a) set the crossover probability Pc ∈ (0.4,0.99) of binary strings genetic algorithm, mutation probability Pm ∈ (0.0001,
0.001), the length L=N of chromosome, the quantity Q of chromosome, the number of iterations tg, maximum number of iterations Tg >=50, adaptive optimal control
The length for spending list G and optimal chromosome list C, G and C is equal to Tg, and the tg element G (tg) in G represents the tg times repeatedly
For obtained adaptive optimal control angle value, the tg Elements C (tg) in C represents the optimal chromosome that the tg times iteration obtains, and enables tg
=1；
It is row that (8b), which is generated at random with Q, is that the jth row Y [j] in binary chromosome the matrix Y, Y of column represents jth with L
The element of a chromosome, Y [j] from left to right respectively represents the resistance for the feed port included in Ant' and being sequentially arranged
Anti-, first lump port impedance ..., the impedance of n-th lump port；
(8c) uses built-in multiport algorithm, passes through Sbest{ Etotal}bestCalculate Y [j] corresponding global radiation directional diagram
AFtotal(θ)tg,j；
(8d) passes through AFtotal(θ)tg,jCalculate the ideal adaptation angle value fit of Y [j]GA(j), Q individual fitness value is obtained,
And choose minimum fitness value therein and G (tg) in adaptive optimal control degree list G is updated, it obtains updated optimal suitable
Response list G', while choosing the corresponding chromosome of minimum fitness value and the C (tg) of optimal chromosome list C is updated,
Obtain updated optimal particle list C':
fitGA(j)=| θgoal-θtotal_max|
θ in formulagoalRepresent the directional diagram target orientation angle of setting, θtotal_maxRepresent AFtotal(θ)tg,jMaximum value institute is right
The directional diagram scanning angle answered, | | represent the operation that takes absolute value；
(8e) by Y [j] carry out crossover operation intersected after dyeing volume matrix Y', to Y'[j] carry out variation behaviour
Dyeing volume matrix Y " after being made a variation；
(8f) enables Y=Y ", enables G=G', enables C=C', and judges whether tg=Tg is true, if so, step (8g) is executed, it is no
Then, tg=tg+1 is enabled, and executes step (8c)；
(8g) determines to whether there is 0 element in G, if so, using the chromosome in the corresponding C of 0 element as the optimal of Ant'
Switch on and off state is distributed gbest, and otherwise, the chromosome in the corresponding C of fitness value minimum in G is opened as the optimal of Ant'
It closes on off operating mode and is distributed gbest；
(9) optimum results of Ant structural parameters and the distribution of switch on and off state are obtained:
The optimum results being distributed using pbest and gbest as Ant structural parameters and switch on and off state.
Compared with prior art, the present invention having the advantage that
The present invention increases on the basis of existing built-in multiport algorithm calculates antenna reflection coefficient to antenna direction
The computing function of figure, and in the optimization process of the structural parameters to directional diagram reconstructable pixel antenna to be optimized using improve after
Built-in multiport algorithm analyze the frequency stability of directional diagram reconstructable pixel antenna to be optimized and the stability of polarized state,
It has obtained optimization rear center's frequency and polarization mode stablizes constant antenna to be optimized, and passed through optimization directional diagram on this basis
The switch on and off state of restructural pixel antenna is distributed, and obtains the switch on and off state point of antenna under different directions figure scanning angle
Cloth realizes the performance of directional diagram reconstructable.Therefore compared with prior art, optimization method proposed by the present invention ensure that and show
On the basis of having technology same efficiency, solves the poor skill of the non-reconstruct segmental stability of reconfigurable antenna of the existing technology
Art problem has higher engineering practical value.
Detailed description of the invention
Fig. 1 is implementation flow chart of the present invention；
Fig. 2 is the structural schematic diagram for the directional diagram reconstructable pixel antenna to be optimized that the present invention uses；
Specific embodiment
In the following with reference to the drawings and specific embodiments, present invention is further described in detail:
With reference to Fig. 1, the present invention includes the following steps:
Step 1) determines the structure and optimization object of directional diagram reconstructable pixel antenna to be optimized；
The structure of directional diagram reconstructable pixel antenna Ant to be optimized is as shown in Fig. 2, include the parasitic agent being sequentially arranged
Plate, driving dielectric-slab and feed dielectric plate are provided with air layer between parasitic agent plate and driving dielectric-slab, and the present embodiment is hollow
Gas-bearing formation is replaced using the relative dielectric constant cystosepiment equal with air, and by parasitic agent plate, cystosepiment, driving dielectric-slab and
Feed dielectric plate is layered on top of each other by sequence from top to bottom, and four laminate material length and width are 50mm, and parasitic agent plate dielectric constant is
2.2, it is highly 1mm；Cystosepiment dielectric constant is 1, is highly ha；The dielectric-slab is driven with feed dielectric version dielectric constant to be
3.55, it is highly 1mm；The upper surface of parasitic agent plate is printed with the pixel that 3 × 3 shapes of periodic arrangement are square
Patch, side length wp, spacing is 1.5mm between adjacent pixel patch, is switched and is connected by rectangle, and the quantity of rectangle switch is N
=12, having a size of 1.5 × 0.8mm2；The upper surface of driving dielectric-slab is printed with the microband paste that shape is rectangle, having a size of ld
× wd, microband paste is provided centrally with rectangular channel, having a size of ls × ws；Feed dielectric plate upper surface is printed with fed microstrip
Line, having a size of 25 × 2.5mm2, lower surface is printed with metal floor, and metal floor size is identical as feed dielectric plate lower surface,
The incline of the feed dielectric plate is equipped with the feed port with feeding microstrip line connection；
Determining that optimization object includes structural parameters and the distribution of rectangle switch on and off state, the quantity of structural parameters is d=6,
Including ld, wd, ls, ws, wp and ha；
Step 2), which is established, calculates Ant reflection coefficient and experiment vector set U needed for directional diagram:
Number 0 is converted into corresponding 4096 12 binary codes to number 4095, and according to sequence from small to large
4096 12 binary codes are formed a line, binary matrix U, the M=4096 that size is M × N are obtained, U is pressed into row piecemeal
Afterwards as experiment vector set needed for calculating Ant reflection coefficient and directional diagram:
Impedance control vector set Z needed for step 3) establishes calculating Ant reflection coefficient and directional diagram:
Establish the matrix of characteristic impedance when size matches for the value of 2 × 1 and each element equal to feed port in Ant
Z1；Establish size be 2 × 12 and first row element value be equal to Ant in rectangle switch conduction when equivalent impedance, the second row member
The value of element is equal to the matrix Z of equivalent impedance when rectangle switch disconnects in Ant2；Establish the value that size is 12 × 1 and each element
Equal to the matrix Z of equivalent impedance when rectangle switch disconnects in Ant3；Establish the value etc. that size is 12 × 12 and diagonal entry
Equivalent impedance in Ant when rectangle switch conduction, the value of remaining all elements are equal to equivalent when rectangle switch disconnects in Ant
The matrix Z of impedance4；
To matrix Z1、Z2、Z3And Z4It is combined, obtains the matrix Z that size is 14 × 13, and Z is divided again by row
As impedance control vector set needed for calculating Ant reflection coefficient and directional diagram after block:
Characteristic impedance when feed port matches in the present embodiment is 50 Ω, and equivalent impedance when switch conduction is 10-9
Ω, equivalent impedance when switch disconnects are 109Ω；
Step 4) obtains the directional diagram reconstructable pixel antenna Ant' to be optimized of load lump port:
12 rectangles connected between adjacent pixel patch in Ant switch is replaced by 12 lump ports, is obtained adjacent
The directional diagram reconstructable pixel antenna Ant' to be optimized connected by lump port is connected between pixel patch；
Step 5) writes the script function of Ant':
Foundation includes the function f of input, output and main body in MATLAB software, and using f as the script letter of Ant'
Number, in which:
The input of function f is s and z, and it includes ld, wd, ls, ws, the wp being arranged successively that wherein s, which indicates that length is equal to 6, element,
With the structure control vector of ha；It includes the feed end included in Ant' and being sequentially arranged that z, which indicates that length is equal to 13, element,
Mouthful impedance, first lump port impedance ..., the impedance control vector of the impedance of the 12nd lump port；
The output of function f is S, Etotal(θ)、Eco(θ) and Ecross(θ), wherein S indicates the full end that size is equal to 13 × 13
Mouth scattering parameter matrix；Etotal(θ)、Eco(θ) and EcrossIt is θ that (θ), which respectively indicates scanning angle, and length is equal to 181 total intensity
Directional diagram vector, main polarization field strength pattern vector sum cross polarization field strength pattern vector, -90≤θ≤90；
The main body of function f includes modeling code, calculation code and return code, wherein modeling code, for reading f's
Input, and the modeling functions in the control script MATLAB-HFSS-API of Calling MATLAB software, carry out abstract to Ant' and build
Mould；Calculation code carries out analysis meter to the abstract model of foundation for the analytic function in Calling MATLAB-HFSS-API
It calculates；Return code, for exporting calculation code result calculated；
Step 6) calculates the optimum structure parameter of Ant' based on particle swarm optimization algorithm:
Step 6a) initialization particle swarm optimization algorithm parameter, including inertial factor w=0.5, the first Studying factors c1=
2, the second Studying factors c2=2, the dimension D=6 of particle, the quantity P=100 of particle, search space range [value_min,
Value_max], velocity interval [- 1,1] the number of iterations tp, maximum number of iterations Tp=50, adaptive optimal control degree list F and most
The length of excellent particle list R, F and R are equal to Tp, and the tp element F (tp) in F represents the adaptive optimal control that the tp times iteration obtains
Angle value, the tp element R (tp) in R represent the optimal particle that the tp times iteration obtains, enable tp=1, enable value_min=
[10,20,5,0.1,10,5], value_max=[20,30,10,0.5,10,15]；
Step 6b) particle mass matrix X and rate matrices V that size is P × D are generated at random, wherein the i-th row X in X
[i] represents i-th of particle, and the i-th row V [i] in V represents the speed of i-th of particle, and element from left to right divides in X [i] and V [i]
Ld, wd, ls, ws, wp and ha in Ant' are not represented, and the value of each element is greater than value_min and is less than value_ in X
The value of each of max, V element is greater than the opposite number of value_min and is less than value_min；
Step 6c) enable s=X [i], z=zs, and by s and z input f, obtain the corresponding full port scattering parameter square of X [i]
Battle array Stp,i；
Step 6d) enable s=X [i], z=zE0,zE1,...,zEN, and by s and z input f, obtain the corresponding main pole of X [i]
Change field strength pattern set { Eco}iAnd cross polarization field strength pattern set { Ecross}i:
{Eco}i={ Eco(θ)tp,0,...,Eco(θ)tp,n,...,Eco(θ)tp,N}i,
{Ecross}i={ Ecross(θ)tp,0,...,Ecross(θ)tp,n,...,Ecross(θ)tp,N}i；
Step 6e) built-in multiport algorithm is used, pass through Stp,i、{Eco}i{ Ecross}iCalculate each particle X [i] in X
Corresponding reflection coefficient set { SLtp}i, main polarization radiation direction set of graphs { AFco}iAnd cross polarization radiations direction set of graphs
{AFcross}i:
{SLtp}i={ SLtp,1,...,SLtp,m,...,SLtp,M}i,
{AFco}i={ AFco(θ)tp,1,...,AFco(θ)tp,m,...,AFco(θ)tp,M}i,
{AFcross}i={ AFcross(θ)tp,1,...,AFcross(θ)tp,m,...,AFcross(θ)tp,M}i；
The built-in multiport algorithm realizes step are as follows:
By full port scattering parameter matrix S piecemeal as follows:
In formula, S=Stp,i；
Construct the binary vector u that length is N:
U=[u (1) ... u (n) ... u (N)]
In formula, u=um；
Construct the transmission matrix Γ that size is N × NBB, ΓBBOff diagonal element be 0, diagonal entry is respectively
[τ1,...,τn,...,τN] and:
Calculate reflection coefficient:
SL=SAA+SAB(I-ΓBBSBB)-1ΓBBSBA
I is the unit matrix that size is N × N in formula；
The full port incidence wave column vector a that length is N+1 is constructed, and piecemeal is carried out to a by following equation:
Wherein aA=1, aB=(I- ΓBBSBB)-1ΓBBSBAaA；
By being weighted to field strength pattern set, calculating antenna pattern:
In formula, E (θ)nFor field strength pattern set { E (θ)1,...,E(θ)n,...,E(θ)NIn n-th of field strength direction
Figure calculates { AFco}iWhen, { E (θ)1,...,E(θ)n,...,E(θ)N}={ Eco(θ)tp,0,...,Eco(θ)tp,n,...,Eco
(θ)tp,N}i, calculate { AFcross}iWhen, { E (θ)1,...,E(θ)n,...,E(θ)N}={ Ecross(θ)tp,0,...,Ecross
(θ)tp,n,...,Ecross(θ)tp,N}i；
Step 6f) pass through { SLtp}i、{AFco}i{ AFcross}i, calculate the ideal adaptation angle value of each particle X [i] in X
fitPSO(i), and the smallest F (tp) progress in adaptive optimal control degree list F in P calculated individual fitness value is chosen
It updates, obtains updated adaptive optimal control degree list F', choose the corresponding particle of minimum fitness value to optimal particle list R's
R (tp) is updated, and obtains updated optimal particle list:
In formula, SLtp,mFor { SLtp}iIn m-th of reflection coefficient, SL0For the reflection coefficient reference value of setting, SL0=10,
Cptp,mTo pass through { AFco}i{ AFcross}iCorresponding m-th of the cross polarization of the X [i] of calculating inhibits ratio, Cp0For the intersection of setting
Polarization inhibits than reference value, Cp0=15, max { } representative take biggish element, ∑ in { } to represent sum operation:
AF in formulaco(θ)tp,mFor { AFco}iIn m-th of main polarization antenna pattern, AFcross(θ)tp,mFor { AFcross}iIn
M-th of cross polarization radiations directional diagram, θco_maxFor AFco(θ)tp,mThe corresponding scanning angle of middle maximum value；
Step 6g) every a line of V and every a line of X are updated respectively, obtain updated rate matrices V' and grain
Subgroup matrix X', and the value for being less than value_min in X' is replaced with value_min, the value value_ greater than value_max
Max replacement, obtains replaced particle mass matrix X ", wherein the calculation formula of V' and X' is respectively as follows:
V'[i]=w × V [i]+c1×random()×(R'(tp)-X[i])+c2×random()×(R'(ind)-X
[i]),
X'[i]=X [i]+V'[i],
Wherein, random () indicates to generate the random number between one 0 to 1, and ind is indicated in adaptive optimal control degree list F' most
Position corresponding to small value；
Step 6h) X=X " is enabled, V=V' is enabled, F=F' is enabled, enables R=R', and judges whether tp=Tp is true, if so, executing
Step 6i), otherwise, tp=tp+1 is enabled, and execute step 6c)；
Step 6i) determine to whether there is 0 element in F, if so, using the particle in the corresponding R of 0 element as the optimal of Ant'
Structural parameters pbest, otherwise, using the particle in the corresponding R of least member in F as the optimum structure parameter pbest of Ant'；
Step 7) enables s=pbest, z=zs, and by s and z input f, obtain the corresponding full port scattering parameter of pbest
Matrix Sbest；
Enable s=pbest, z=zE0,zE1,...,zEN, and by s and z input f, obtain the corresponding total intensity direction pbest
Set of graphs { Etotal}best:
{Etotal}best={ Etotal(θ)0,...,Etotal(θ)n,...,Etotal(θ)N}best；
Step 8) calculates the optimized switching on off operating mode distribution of Ant' based on binary strings genetic algorithm:
Step 8a) set the crossover probability Pc=0.5 of binary strings genetic algorithm, mutation probability Pm=0.0005, chromosome
Length L=N=12, the quantity Q=100 of chromosome, the number of iterations tg, maximum number of iterations Tg=50, adaptive optimal control degree list G
And the length of optimal chromosome list C, G and C are equal to Tg, the tg element G (tg) in G represents what the tg times iteration obtained
Adaptive optimal control angle value, the tg Elements C (tg) in C represent the optimal chromosome that the tg times iteration obtains, enable tg=1；
Step 8b) it is generated at random with Q as row, it is jth row Y [j] representative in binary chromosome the matrix Y, Y arranged with L
The element of j-th of chromosome, Y [j] from left to right respectively represents the feed port included in Ant' and being sequentially arranged
Impedance, the impedance of first lump port ..., the impedance of n-th lump port；
Step 8c) built-in multiport algorithm is used, pass through Sbest{ Etotal}bestCalculate the corresponding global radiation direction Y [j]
Scheme AFtotal(θ)tg,j；
The built-in multiport algorithm, realizes that steps are as follows:
By full port scattering parameter matrix S piecemeal as follows:
In formula, S=Sbest；
Construct the binary vector u that length is N:
U=[u (1) ... u (n) ... u (N)]
In formula, u=Y [j]；
Construct the transmission matrix Γ that size is N × NBB, ΓBBOff diagonal element be 0, diagonal entry is respectively
[τ1,...,τn,...,τN] and:
Calculate reflection coefficient:
SL=SAA+SAB(I-ΓBBSBB)-1ΓBBSBA
I is the unit matrix that size is N × N in formula；
The full port incidence wave column vector a that length is N+1 is constructed, and piecemeal is carried out to a by following equation:
Wherein aA=1, aB=(I- ΓBBSBB)-1ΓBBSBAaA；
By being weighted to field strength pattern set, calculating antenna pattern:
In formula, E (θ)nFor field strength pattern set { E (θ)1,...,E(θ)n,...,E(θ)NIn n-th of field strength direction
Figure, { E (θ)1,...,E(θ)n,...,E(θ)N}={ Etotal(θ)0,...,Etotal(θ)n,...,Etotal(θ)N}best；
Step 8d) pass through AFtotal(θ)tg,jCalculate the ideal adaptation angle value fit of Y [j]GA(j), Q individual adaptation degree is obtained
Value, and choose minimum fitness value therein and G (tg) in adaptive optimal control degree list G is updated, obtain it is updated most
Excellent fitness list G', while choosing the corresponding chromosome of minimum fitness value and the C (tg) of optimal chromosome list C is carried out more
Newly, updated optimal particle list C' is obtained:
fitGA(j)=| θgoal-θtotal_max|
θ in formulagoalRepresent the directional diagram target orientation angle of setting, θ in the present embodimentgoal- 15,0,15 are taken respectively,
θtotal_maxRepresent AFtotal(θ)tg,jDirectional diagram scanning angle corresponding to maximum value, | | represent the operation that takes absolute value；
Step 8e) by carrying out the dyeing volume matrix Y' after crossover operation is intersected to Y [j], to Y'[j] make a variation
Operate the dyeing volume matrix Y " after being made a variation；
Step 8f) Y=Y " is enabled, enable G=G', enable C=C', and judge whether tg=Tg true, if so, executing step 8g),
Otherwise, tg=tg+1 is enabled, and executes step 8c)；
Step 8g) determine to whether there is 0 element in G, if so, most using the chromosome in the corresponding C of 0 element as Ant'
Excellent switch on and off state is distributed gbest, otherwise, using the chromosome in the corresponding C of fitness value minimum in G as the optimal of Ant'
Switch on and off state is distributed gbest；
The optimum results that step 9) is distributed using pbest and gbest as Ant structural parameters and switch on and off state.
It is that any limitation of the invention is not constituted, it is clear that for ability to a specific embodiment of the invention above
It, all may be in the feelings without departing substantially from the principle of the invention and step for the professional in domain, then after understanding the content of present invention and principle
Under condition, various modifications and variations in form and details are carried out, but these modifications and variations based on inventive concept still exist
In claim and protection scope of the invention.
Claims (2)
1. a kind of directional diagram reconstructable pixel antenna optimization method based on built-in multiport algorithm, which is characterized in that including such as
Lower step:
(1) structure and optimization object of directional diagram reconstructable pixel antenna to be optimized are determined:
(1a) directional diagram reconstructable pixel antenna Ant to be optimized, including the parasitic agent plate, driving dielectric-slab and feedback being sequentially arranged
Dielectric plate is provided with air layer between parasitic agent plate and driving dielectric-slab, drives dielectric-slab and feed dielectric plate phase alternating layers
Folded, the upper surface of parasitic agent plate is printed with the m of periodic arrangementp×npA shape is the pixel patch of rectangle, mp>=2, np≥
2, it is switched between adjacent pixel patch and is connected by rectangle, the quantity of rectangle switch is N；The upper surface of driving dielectric-slab is printed with
Shape is the microband paste of rectangle, and microband paste is provided centrally with rectangular channel；It is micro- that feed dielectric plate upper surface is printed with feed
Band line, lower surface are printed with metal floor, and the incline of the feed dielectric plate is equipped with the feed port with feeding microstrip line connection；
The optimization object of (1b) directional diagram reconstructable pixel antenna to be optimized includes structural parameters and rectangle switch on and off state point
Cloth, the quantity of structural parameters are d, length and the wide, length of rectangular channel and width, the length of pixel patch and width including microband paste, with
And the height of air layer；
(2) experiment vector set U needed for calculating Ant reflection coefficient and directional diagram is established:
By number 0 to number 2N- 1 is converted to corresponding 2NA N binary code, and according to sequence from small to large by 2NIt is N a
Binary code forms a line, and obtains the binary matrix U, M=2 that size is M × NN, by U by anti-as Ant is calculated after row piecemeal
Experiment vector set needed for penetrating coefficient and directional diagram:
(3) impedance control vector set Z needed for calculating Ant reflection coefficient and directional diagram is established:
(3a) establishes the matrix Z of characteristic impedance when size matches for the value of 2 × 1 and each element equal to feed port in Ant1；
Equivalent impedance when size is equal to rectangle switch conduction in Ant for the value of 2 × N and the first row element is established, the second row element
Value is equal to the matrix Z of equivalent impedance when rectangle switch disconnects in Ant2；The value that size is established as N × 1 and each element is equal to
The matrix Z of equivalent impedance when rectangle switch disconnects in Ant3；Size is established to be equal in Ant for the value of N × N and diagonal entry
Equivalent impedance when rectangle switch conduction, the value of remaining all elements are equal to equivalent impedance when rectangle switch disconnects in Ant
Matrix Z4；
(3b) is to matrix Z1、Z2、Z3And Z4It is combined, obtains the matrix Z that size is (N+2) × (N+1), and Z is carried out by row
Again as impedance control vector set needed for calculating Ant reflection coefficient and directional diagram after piecemeal:
(4) the directional diagram reconstructable pixel antenna Ant' to be optimized of load lump port is obtained:
The N number of rectangle connected between adjacent pixel patch in Ant switch is replaced by N number of lump port, obtains adjacent pixel patch
The directional diagram reconstructable pixel antenna Ant' to be optimized connected by lump port is connected between piece；
(5) script function of Ant' is write:
Established in MATLAB software include input, output and main body function f, and using f as the script function of Ant',
In:
The input of function f is s and z, and wherein s indicates that length is equal to d, element includes the length for being arranged successively microband paste, micro-strip patch
The structure control of the width of piece, the length of rectangular channel, the width of rectangular channel, the length of pixel patch, the height of the width of pixel patch and air layer
Vector processed；Z indicates that length is equal to the impedance that N+1, element include the feed port included in Ant' and being sequentially arranged, the
The impedance of one lump port ..., the impedance control vector of the impedance of n-th lump port；
The output of function f is S, Etotal(θ)、Eco(θ) and Ecross(θ), wherein S indicates that size is complete equal to (N+1) × (N+1)
Port scattering parameter matrix；Etotal(θ)、Eco(θ) and EcrossIt is θ that (θ), which respectively indicates scanning angle, and length is equal to 181 resultant field
Strong directional diagram vector, main polarization field strength pattern vector sum cross polarization field strength pattern vector, -90≤θ≤90；
The main body of function f includes modeling code, calculation code and return code, wherein modeling code, for reading the defeated of f
Enter, and the modeling functions in the control script MATLAB-HFSS-API of Calling MATLAB software, abstract is carried out to Ant' and is built
Mould；Calculation code carries out analysis meter to the abstract model of foundation for the analytic function in Calling MATLAB-HFSS-API
It calculates；Return code, for exporting calculation code result calculated；
(6) the optimum structure parameter of Ant' is calculated based on particle swarm optimization algorithm:
(6a) initializes the parameter of particle swarm optimization algorithm, including inertial factor w ∈ (0,1), the first Studying factors c1∈(1,3)、
Second Studying factors c2∈ (1,3), the dimension D=d of particle, the quantity P ∈ (50,500) of particle, search space range
[value_min, value_max], velocity interval [v_min, v_max], the number of iterations tp, maximum number of iterations Tp >=50, most
The length of excellent fitness list F and optimal particle list R, F and R are equal to Tp, and the tp element F (tp) in F represents tp
The adaptive optimal control angle value that secondary iteration obtains, the tp element R (tp) in R represent the optimal particle that the tp times iteration obtains, and enable
Tp=1；
(6b) generates the particle mass matrix X and rate matrices V that size is P × D at random, wherein the i-th row X [i] in X is represented
I-th of particle, the i-th row V [i] in V represent the speed of i-th of particle, and element from left to right respectively represents in X [i] and V [i]
The length of microband paste in Ant', the width of microband paste, the length of rectangular channel, the width of rectangular channel, the length of pixel patch, pixel patch
Wide and air layer height, and the value of each element is greater than value_min and is less than value_max, each of V in X
The value of element is greater than the opposite number of value_min and is less than value_min；
(6c) enables s=X [i], z=zs, and by s and z input f, obtain the corresponding full port scattering parameter matrix S of X [i]tp,i；
(6d) enables s=X [i], z=zE0,zE1,...,zEN, and by s and z input f, obtain the corresponding main polarization field strength side X [i]
To set of graphs { Eco}iAnd cross polarization field strength pattern set { Ecross}i:
{Eco}i={ Eco(θ)tp,0,...,Eco(θ)tp,n,...,Eco(θ)tp,N}i,
{Ecross}i={ Ecross(θ)tp,0,...,Ecross(θ)tp,n,...,Ecross(θ)tp,N}i；
(6e) uses built-in multiport algorithm, passes through Stp,i、{Eco}i{ Ecross}iIt is corresponding anti-to calculate each particle X [i] in X
Penetrate coefficient sets { SLtp}i, main polarization radiation direction set of graphs { AFco}iAnd cross polarization radiations direction set of graphs
{AFcross}i:
{SLtp}i={ SLtp,1,...,SLtp,m,...,SLtp,M}i,
{AFco}i={ AFco(θ)tp,1,...,AFco(θ)tp,m,...,AFco(θ)tp,M}i,
{AFcross}i={ AFcross(θ)tp,1,...,AFcross(θ)tp,m,...,AFcross(θ)tp,M}i；
(6f) passes through { SLtp}i、{AFco}i{ AFcross}i, calculate the ideal adaptation angle value fit of each particle X [i] in XPSO
(i), it and chooses the smallest F (tp) in adaptive optimal control degree list F in P calculated individual fitness value and is updated,
Updated adaptive optimal control degree list F' is obtained, chooses the corresponding particle of minimum fitness value to the R (tp) of optimal particle list R
It is updated, obtains updated optimal particle list:
In formula, SLtp,mFor { SLtp}iIn m-th of reflection coefficient, SL0For the reflection coefficient reference value of setting, Cptp,mTo pass through
{AFco}i{ AFcross}iCorresponding m-th of the cross polarization of the X [i] of calculating inhibits ratio, Cp0Inhibit ratio for the cross polarization of setting
Reference value, max { } representative take biggish element, ∑ in { } to represent sum operation:
AF in formulaco(θ)tp,mFor { AFco}iIn m-th of main polarization antenna pattern, AFcross(θ)tp,mFor { AFcross}iIn m-th
Cross polarization radiations directional diagram, θco_maxFor AFco(θ)tp,mThe corresponding scanning angle of middle maximum value；
(6g) is updated every a line of V and every a line of X respectively, obtains updated rate matrices V' and particle mass matrix
X', and the value for being less than value_min in X' is replaced with value_min, the value greater than value_max is replaced with value_max,
Replaced particle mass matrix X " is obtained, wherein the calculation formula of V' and X' is respectively as follows:
V'[i]=w × V [i]+c1×random()×(R'(tp)-X[i])+c2× random () × (R'(ind)-X [i]),
X'[i]=X [i]+V'[i],
Wherein, random () indicates to generate the random number between one 0 to 1, and ind indicates minimum value in adaptive optimal control degree list F'
Corresponding position；
(6h) enables X=X ", enables V=V', enables F=F', enables R=R', and judges whether tp=Tp is true, if so, executing step
(6i) otherwise enables tp=tp+1, and executes step (6c)；
(6i) determines to whether there is 0 element in F, if so, joining the particle in the corresponding R of 0 element as the optimum structure of Ant'
Number pbest, otherwise, using the particle in the corresponding R of least member in F as the optimum structure parameter pbest of Ant'；
(7) the corresponding full port scattering parameter matrix S of optimum structure parameter pbest is obtainedbestWith total intensity direction set of graphs
{Etotal}best:
(7a) enables s=pbest, z=zs, and by s and z input f, obtain the corresponding full port scattering parameter matrix of pbest
Sbest；
(7b) enables s=pbest, z=zE0,zE1,...,zEN, and by s and z input f, obtain the corresponding total intensity direction pbest
Set of graphs { Etotal}best:
{Etotal}best={ Etotal(θ)0,...,Etotal(θ)n,...,Etotal(θ)N}best；
(8) the optimized switching on off operating mode distribution of Ant' is calculated based on binary strings genetic algorithm:
(8a) set the crossover probability Pc ∈ (0.4,0.99) of binary strings genetic algorithm, mutation probability Pm ∈ (0.0001,0.001),
The length L=N of chromosome, the quantity Q of chromosome, the number of iterations tg, maximum number of iterations Tg >=50, adaptive optimal control degree list G
And the length of optimal chromosome list C, G and C are equal to Tg, the tg element G (tg) in G represents what the tg times iteration obtained
Adaptive optimal control angle value, the tg Elements C (tg) in C represent the optimal chromosome that the tg times iteration obtains, enable tg=1；
It is row that (8b), which is generated at random with Q, is that the jth row Y [j] in binary chromosome the matrix Y, Y of column represents j-th of dye with L
Colour solid, the element of Y [j] from left to right respectively represent the impedance for the feed port included in Ant' and being sequentially arranged,
The impedance of one lump port ..., the impedance of n-th lump port；
(8c) uses built-in multiport algorithm, passes through Sbest{ Etotal}bestCalculate Y [j] corresponding global radiation directional diagram AFtotal
(θ)tg,j；
(8d) passes through AFtotal(θ)tg,jCalculate the ideal adaptation angle value fit of Y [j]GA(j), Q individual fitness value is obtained, and is selected
It takes minimum fitness value therein to be updated the G (tg) in adaptive optimal control degree list G, obtains updated adaptive optimal control degree
List G', while choosing the corresponding chromosome of minimum fitness value and the C (tg) of optimal chromosome list C is updated, it obtains
Updated optimal particle list C':
fitGA(j)=| θgoal-θtotal_max|
θ in formulagoalRepresent the directional diagram target orientation angle of setting, θtotal_maxRepresent AFtotal(θ)tg,jCorresponding to maximum value
Directional diagram scanning angle, | | represent the operation that takes absolute value；
(8e) by Y [j] carry out crossover operation is intersected after dyeing volume matrix Y', to Y'[j] progress mutation operation obtain
Dyeing volume matrix Y " after to variation；
(8f) enables Y=Y ", enables G=G', enables C=C', and judges whether tg=Tg is true, if so, step (8g) is executed, otherwise,
Tg=tg+1 is enabled, and executes step (8c)；
(8g) determines to whether there is 0 element in G, if so, using the chromosome in the corresponding C of 0 element as the optimized switching of Ant'
On off operating mode is distributed gbest, and otherwise, the chromosome in the corresponding C of fitness value minimum in G is led to as the optimized switching of Ant'
Disconnected state is distributed gbest；
(9) optimum results of Ant structural parameters and the distribution of switch on and off state are obtained:
The optimum results being distributed using pbest and gbest as Ant structural parameters and switch on and off state.
2. the directional diagram reconstructable pixel antenna optimization method according to claim 1 based on built-in multiport algorithm,
It is characterized in that, built-in multiport algorithm described in step (6d) and step (8c), realizes step are as follows:
By full port scattering parameter matrix S piecemeal as follows:
For step (6d), S=Stp,i, for step (8c), S=Sbest；
Construct the binary vector u that length is N:
U=[u (1) ... u (n) ... u (N)]
For step (6d), u=um, for step (8c), u=Y [j]；
Construct the transmission matrix Γ that size is N × NBB, ΓBBOff diagonal element be 0, diagonal entry be respectively [τ1,...,
τn,...,τN] and:
Calculate reflection coefficient:
SL=SAA+SAB(I-ΓBBSBB)-1ΓBBSBA
I is the unit matrix that size is N × N in formula；
The full port incidence wave column vector a that length is N+1 is constructed, and piecemeal is carried out to a by following equation:
Wherein aA=1, aB=(I- ΓBBSBB)-1ΓBBSBAaA；
By being weighted to field strength pattern set, calculating antenna pattern:
In formula, E (θ)nFor field strength pattern set { E (θ)1,...,E(θ)n,...,E(θ)NIn n-th of field strength pattern, for
Step (6d) calculates { AFco}iWhen, { E (θ)1,...,E(θ)n,...,E(θ)N}={ Eco(θ)tp,0,...,Eco(θ)tp,n,...,
Eco(θ)tp,N}i, calculate { AFcross}iWhen, { E (θ)1,...,E(θ)n,...,E(θ)N}={ Ecross(θ)tp,0,...,Ecross
(θ)tp,n,...,Ecross(θ)tp,N}i, for step (8c), { E (θ)1,...,E(θ)n,...,E(θ)N}={ Etotal
(θ)0,...,Etotal(θ)n,...,Etotal(θ)N}best。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910766927.9A CN110516336B (en) | 2019-08-20 | 2019-08-20 | Directional diagram reconfigurable pixel antenna optimization method based on built-in multi-port algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910766927.9A CN110516336B (en) | 2019-08-20 | 2019-08-20 | Directional diagram reconfigurable pixel antenna optimization method based on built-in multi-port algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110516336A true CN110516336A (en) | 2019-11-29 |
| CN110516336B CN110516336B (en) | 2023-01-06 |
Family
ID=68626719
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910766927.9A Active CN110516336B (en) | 2019-08-20 | 2019-08-20 | Directional diagram reconfigurable pixel antenna optimization method based on built-in multi-port algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110516336B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022257194A1 (en) * | 2021-06-10 | 2022-12-15 | 深圳技术大学 | Miniaturized pattern reconfigurable pixel antenna and application method |
| AT526407A1 (en) * | 2022-08-08 | 2024-02-15 | Univ Innsbruck | PIXELATED ANTENNAS |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014172383A2 (en) * | 2013-04-16 | 2014-10-23 | Paneratech, Inc. | Antenna and method for optimizing the design thereof |
| CN106935973A (en) * | 2017-02-16 | 2017-07-07 | 北京科技大学 | A kind of method for designing of coaxial feed antenna |
| CN107423529A (en) * | 2017-08-30 | 2017-12-01 | 同济大学 | Metamaterial Precise spraying method |
-
2019
- 2019-08-20 CN CN201910766927.9A patent/CN110516336B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014172383A2 (en) * | 2013-04-16 | 2014-10-23 | Paneratech, Inc. | Antenna and method for optimizing the design thereof |
| CN106935973A (en) * | 2017-02-16 | 2017-07-07 | 北京科技大学 | A kind of method for designing of coaxial feed antenna |
| CN107423529A (en) * | 2017-08-30 | 2017-12-01 | 同济大学 | Metamaterial Precise spraying method |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2022257194A1 (en) * | 2021-06-10 | 2022-12-15 | 深圳技术大学 | Miniaturized pattern reconfigurable pixel antenna and application method |
| AT526407A1 (en) * | 2022-08-08 | 2024-02-15 | Univ Innsbruck | PIXELATED ANTENNAS |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110516336B (en) | 2023-01-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dong et al. | Fast multi-objective optimization of multi-parameter antenna structures based on improved BPNN surrogate model | |
| Jiao et al. | A dynamic inertia weight particle swarm optimization algorithm | |
| CN106355245A (en) | Method for integrating array antenna directional images on basis of neural network algorithms | |
| CN110516336A (en) | Directional diagram reconstructable pixel antenna optimization method based on built-in multiport algorithm | |
| CN103646144B (en) | aperiodic array antenna design method | |
| CN110504550A (en) | It is a kind of to radiate and scatter integrated information metamaterial surface and its application | |
| CN109670140A (en) | A kind of array antenna Electrical Analysis method based on array element feature mode | |
| CN104899374B (en) | Based on small echo variation wind Drive Optimization algorithm collinear array Pattern Synthesis method | |
| CN101872894A (en) | Medium resonator antenna with reconfigurable directional diagram and phased array thereof | |
| CN110148838A (en) | A kind of directional diagram reconstructable planar array antenna and its control method based on digital coding characterization | |
| Li et al. | Tri-objective compact log-periodic dipole array antenna design using MOEA/D-GPSO | |
| CN106777601A (en) | Based on the planar array antenna Pattern Synthesis method that MIFT is combined with QP algorithms | |
| CN112100811A (en) | Antenna array directional diagram synthesis method based on adaptive wind-driven optimization algorithm | |
| CN105975746B (en) | MIMO radar square matrix based on PSO class algorithm is structured the formation optimization method | |
| Panduro et al. | Design of electronically steerable linear arrays with evolutionary algorithms | |
| CN116205143A (en) | Design method for realizing antenna pattern based on physical information neural network | |
| CN116562143A (en) | Antenna topology and parameter mixing optimization method based on normalized Gaussian network | |
| Zhang et al. | Optimal design of fragment-type antenna structure based on PSO-CNN | |
| CN103928763A (en) | Multi-beam antenna based on transformation optics | |
| Gupta et al. | Bandwidth Enhancement in multilayer microstrip proximity coupled array | |
| Patnaik et al. | Modeling frequency reconfigurable antenna array using neural networks | |
| Zheng | Application of Swarm Intelligence Optimization Algorithm in Antenna Design | |
| Ram et al. | IPSO based performance comparison of optimum uniform and non-uniform spacing of circular antenna arrays | |
| Lu et al. | Development of ESPAR antenna array using numerical modelling techniques | |
| Bian et al. | Large thinned array design based on multi-objective cross entropy algorithm |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |