CN110138427A - Extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection - Google Patents

Extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection Download PDF

Info

Publication number
CN110138427A
CN110138427A CN201910400887.6A CN201910400887A CN110138427A CN 110138427 A CN110138427 A CN 110138427A CN 201910400887 A CN201910400887 A CN 201910400887A CN 110138427 A CN110138427 A CN 110138427A
Authority
CN
China
Prior art keywords
matrix
coding matrix
optimal
simulation pre
error
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
Application number
CN201910400887.6A
Other languages
Chinese (zh)
Other versions
CN110138427B (en
Inventor
庞立华
吴文捷
赵恒�
牛晓娟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GUANGZHOU ITS COMMUNICATION EQUIPMENT Co.,Ltd.
Original Assignee
Xian University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Science and Technology filed Critical Xian University of Science and Technology
Priority to CN201910400887.6A priority Critical patent/CN110138427B/en
Publication of CN110138427A publication Critical patent/CN110138427A/en
Application granted granted Critical
Publication of CN110138427B publication Critical patent/CN110138427B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0686Hybrid systems, i.e. switching and simultaneous transmission
    • H04B7/0695Hybrid systems, i.e. switching and simultaneous transmission using beam selection
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0408Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas using two or more beams, i.e. beam diversity
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radio Transmission System (AREA)

Abstract

The invention discloses a kind of extensive multiple-input and multiple-output mixed-beam forming algorithms based on part connection, given part connects the optimal without constraint precoder, primary simulation pre-coding matrix and convergence tolerance of architecture system, and initial number pre-coding matrix is calculated according to without constraint precoder and primary simulation pre-coding matrix, and then calculate initial error;When initial error is less than or equal to convergence tolerance, the optimal simulation pre-coding matrix and optimal digital precode matrix of architecture system are connected using primary simulation pre-coding matrix and initial number pre-coding matrix as part;When initial error is greater than convergence tolerance, to minimize optimal simulation pre-coding matrix and optimal digital precode matrix after error is iterated calculating as target;Part connection architecture system generates and issues mixed-beam according to optimal simulation pre-coding matrix and optimal digital precode matrix;The present invention improves the performance of algorithm, is close to the performance of digital beam-forming system.

Description

Extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection
[technical field]
The invention belongs to field of communication technology more particularly to a kind of extensive multiple-input and multiple-output based on part connection are mixed Close beamforming algorithm.
[background technique]
Massive MIMO technology can be obviously improved the performance and capacity of 5G wireless communication system.Massive mimo system Dual-mode antenna it is in a large number, in systems use digital beam forming technology, will lead to system cost of implementation and energy consumption be very high, Obstruction is brought for the application of Massive MIMO beamforming technique.
In order to overcome this problem, there is Massive MIMO mixed-beam figuration scheme (i.e. part connection framework), Compared to digital beam forming technology, system cost of implementation and energy consumption can be reduced using mixed-beam figuration technology, i.e., in part It connecting in framework, each rf chain is only connected with some antennas, although system sacrifices the gain of part wave beam forming, It is the implementation complexity for greatly reducing hardware, still, in part connection framework, some algorithms based on encoder design, Performance loss is than more serious.
As shown in Figure 1, being Massive mimo system downlink communication mode under partially connected architecture.At this In system, transmitting terminal is equipped with NtRoot antenna, receiving end are equipped with single antenna, send N from transmitting terminalsData flows to receiving end.For The transmission of realization multiple data stream, transmitting terminal are equipped withRf chain connects L on every rf chaintRoot antenna, receiving end are equipped with Single rf chain.
Assuming that array antenna model is ULA, the spacing of antenna is the half of wavelength.Define dimension vector fi, and dimension Vector fiIn each element meetSo entire simulation pre-coding matrix of transmitting terminal is
In order to simplify following description, index of definition collection Ω1, by FRFThe position of neutral element forms, then as (i, j) ∈ Ω1 When,Wherein, i and j respectively indicates the line number and columns of simulation pre-coding matrix,Indicate simulation precoding square The element of i-th row jth column in battle array, in transmitting end structure, digital precode matrix FBBDimension beTransmitting terminal Total power constraint passes through normalization FBBRealize, meet | | FRFFBB| |=Ns
Consider a narrow band blocks flat fading channel, y=HF can be expressed as by receiving signalRFFBBS+n, y are that dimension is Nr× 1 reception vector, H are that dimension is Nr×NtChannel matrix, NrIndicate that the antenna amount of receiving end, s are to send symbol arrow Amount, meets the mathematics phase Expression dimension is Ns×NsUnit vector, n is noise vector, is obeyed only The vertical Gaussian Profile with distribution, mean value 0, variance σ2
In above-mentioned part connection framework, since analog domain pre-coding matrix is permanent mould, so optimization problem right and wrong Convex, therefore, the optimization problem for solving transmitting terminal is extremely complex, needs to consume a large amount of time and can just solve the complete of optimization problem The optimal pre-coding matrix of office, this results in increasing the energy loss of existing part connection architecture system, and greatly improves Network delay.
[summary of the invention]
The object of the present invention is to provide a kind of extensive multiple-input and multiple-output mixed-beam figuration calculations based on part connection Method improves the performance of algorithm, is close to the performance of digital beam-forming system under the premise of cost and energy consumption control, Reduce the energy loss of part connection architecture system.
The invention adopts the following technical scheme: a kind of extensive multiple-input and multiple-output mixed-beam based on part connection is assigned Shape algorithm, comprising the following steps:
Given part connects the optimal public without constraint precoder, primary simulation pre-coding matrix and convergence of architecture system Difference, and initial number pre-coding matrix is calculated according to without constraint precoder and primary simulation pre-coding matrix, and then calculate Initial error out;
When initial error is less than or equal to convergence tolerance, by primary simulation pre-coding matrix and initial number precoding square Battle array connects the optimal simulation pre-coding matrix and optimal digital precode matrix of architecture system as part;
When initial error is greater than convergence tolerance, it is iterated calculating as target to minimize error, until error is less than Or it is equal to convergence tolerance, or complete maximum number of iterations, corresponding simulation pre-coding matrix and digital pre-coding matrix are made For optimal simulation pre-coding matrix and optimal digital precode matrix;
Part connection architecture system is generated and issued according to optimal simulation pre-coding matrix and optimal digital precode matrix Mixed-beam.
Further, optimal to be obtained without constraint precoder by following steps:
The channel state information that architecture system is connected according to part, obtains channel matrix
Wherein, NcIndicate group variety number, NpIndicate the number of path in every cluster, αilIndicate the l articles biography in the i-th cluster reflector The gain factor in defeated path obeys the multiple Gauss distribution that mean value is 0, variance is 1, for (i, l) strip diameter, θr,ilWith Respectively indicate the azimuth for leaving angle and pitch angle, θt,ilWithAzimuth and the pitch angle of angle of arrival are respectively indicated,WithRespectively indicate azimuth angle thetar,ilWithPitching angle thetat,ilWithCorresponding receiving array Response and emission array response;
Pass through H=U Σ VHSingular value decomposition is carried out to channel matrix H, unit matrix V is obtained, is proposed from unit matrix V Its preceding NsColumn obtain matrix V1;Pass throughDiagonal matrix Γ is calculated, and passes through F*=V1Γ obtains optimal pre- without constraint Encoder F*;Wherein, U is unit matrix, and Σ is diagonal matrix.
Further, pass throughObtain initial number pre-coding matrixAnd pass throughObtain initial error ε0
Further, the detailed process of iterative calculation are as follows:
Set objective functionConversion | | F*-FRFFBB||F It can obtainWherein, TnFor F*In n-th of matrix-block, fnIndicate composition FRFIn N-th of matrix-block,Indicate FBBLine n vector;
Then to objective function optimize for
Wherein, tn,mIt indicatesM column,It indicatesM-th of element,The F obtained for kth time iterationRF In n-th of matrix-block,To indicateM-th of element phase increment, It is F after k iterationBBLine n vector;It indicatesThe average value of phase increment;
By to obtaining simulation pre-coding matrix after the objective function optimizationAnd from which further follow that digital precode square Battle arrayAnd error εk
By error εkIt is compared with convergence tolerance, whenWhen, by corresponding simulation pre-coding matrixAnd number Pre-coding matrixAs optimal simulation pre-coding matrix and optimal digital precode matrix;WhenWhen, repeat on Step is stated, until error is less than or equal to convergence tolerance, or maximum number of iterations is completed, by corresponding simulation pre-coding matrixWith digital pre-coding matrixAs optimal simulation pre-coding matrix and optimal digital precode matrix.
The beneficial effects of the present invention are: the present invention is based on Massive MIMO mixed-beam shaping systems, to promote frequency spectrum Efficiency is target, mixing pre-coding matrix is designed using a kind of alternative optimization algorithm based on matrix decomposition, to channel matrix Singular value decomposition is carried out, design is optimal to design final precoding by alternative optimization without constraint digital precode device and synthesizer Device improves the performance of algorithm, is close to the property of digital beam-forming system under the premise of controlling cost and energy consumption as far as possible Can, solve the problems, such as performance loss, and effective verifying has been carried out to the performance of algorithm.
[Detailed description of the invention]
Fig. 1 is Massive mimo system downlink communication mode figure under partially connected architecture in the prior art;
Fig. 2 is algorithm flow chart of the invention;
Fig. 3 is that the spectrum efficiency of many algorithms in verifying embodiment of the invention compares figure;
Fig. 4 is different data streams mesh in verifying embodiment of the invention to the influence diagram of algorithm spectrum efficiency.
[specific embodiment]
The following describes the present invention in detail with reference to the accompanying drawings and specific embodiments.
A kind of extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection is disclosed in the present invention, is obtained Take base station to the channel state information between all user terminalsWherein, HkIndicate base station to use Fading information between the k of family, K are the total number of users of cell.Optimal digital precode matrix based on singular value decomposition design base station end Digital precode matrix and simulation pre-coding matrix are designed by alternative optimization with synthesizer matrix.Information data transmission is opened Begin, sends signal and first pass around a digital precoder FBB, then pass throughA radio frequency link, usingA simulation Precoder, later LtSignal is sent to wireless channel simultaneously by root antenna.In receiving end, a simulation precoder receives channel In signal, by a rf chain, using a digital precoder WBBReceive signal.
As shown in Fig. 2, method of the invention specifically:
Given part connects the optimal public without constraint precoder, primary simulation pre-coding matrix and convergence of architecture system Difference, and initial number pre-coding matrix is calculated according to without constraint precoder and primary simulation pre-coding matrix, and then calculate Initial error out.
It is optimal to be obtained without constraint precoder by following steps:
The channel state information that architecture system is connected according to part, obtains channel matrix
Wherein, NcIndicate group variety number, NpIndicate the number of path in every cluster, αilIndicate the l articles biography in the i-th cluster reflector The gain factor in defeated path obeys the multiple Gauss distribution that mean value is 0, variance is 1, for (i, l) strip diameter, θr,ilWith Respectively indicate the azimuth for leaving angle and pitch angle, θt,ilWithAzimuth and the pitch angle of angle of arrival are respectively indicated,WithRespectively indicate azimuth angle thetar,ilWithPitching angle thetat,ilWithCorresponding receiving array Response and emission array response;
Pass through H=U Σ VHSingular value decomposition is carried out to channel matrix H, unit matrix V is obtained, is proposed from unit matrix V Its preceding NsColumn obtain matrix V1;Pass throughDiagonal matrix Γ is calculated, and passes through F*=V1Γ obtains optimal pre- without constraint Encoder F*;Wherein, U is unit matrix, and Σ is diagonal matrix.
Pass throughObtain initial number pre-coding matrixAnd pass throughObtain initial error ε0
When initial error is less than or equal to convergence tolerance, by primary simulation pre-coding matrix and initial number precoding square Battle array connects the optimal simulation pre-coding matrix and optimal digital precode matrix of architecture system as part.
When initial error is greater than convergence tolerance, it is iterated calculating as target to minimize error, until error is less than Or it is equal to convergence tolerance, or complete maximum number of iterations, corresponding simulation pre-coding matrix and digital pre-coding matrix are made For optimal simulation pre-coding matrix and optimal digital precode matrix.
The detailed process of iterative calculation are as follows:
Set objective functionIteration optimization two Parameter, as simulation pre-coding matrix and digital pre-coding matrix.
Conversion | | F*-FRFFBB||FIt can obtainWherein, TnFor F*In N-th of matrix-block, fnIndicate composition FRFIn n-th of matrix-block,Indicate FBBLine n vector, it is maximum in the present invention The number of iterations can be set as Ku=100, convergence tolerance is set as
When then being optimized to objective function, due to that should meet in the systemSo in optimization process In, temporarily remove normalization constraint
Simplified optimization problem are as follows:
The kth time iteration of definition mixing precoding isAssuming that known initialThenClosed solutions beConversely, when knownWhen, it can update
It is updatingWhen, following derivation is first done, F is worked asBBTo timing, definitionI-th column vector be vi, then can be with It obtainsIn view of FRFStructure, then derivation process are as follows:
Wherein, TnIt is F*N-th of block matrix, dimension Lt×Ns
Then above-mentioned optimization problem can be converted toA subproblem, wherein n-th of subproblem be
In order to solve optimization problem, it is assumed thatA small range search pairIt is updated, definesM-th yuan Element phase beIt can be expressed asSo,It can be expressed as
Wherein,It indicatesThe increment of m-th of element phase, whenWhen very little, Taylor expansion approximation is done
Wherein,It is a vector, symbolIt is defined as Hadamard product, optimization problem can be reconstructed into
Above formula optimization problem is a convex quadratic objective function, and has been contemplated that permanent modular constraint, and above formula is based on approximationCome what is rebuild, so, only whenIt is just able to satisfy when very little, thereforeConstraint must add, So optimization problem is converted into
It is a convex optimization problem at this time, once it obtainsIt is obtained with
Further definitionM-th of element beM be classified as tn,m, above formula can derive are as follows:
Optimization problem, which can be converted into, solves LtA subproblem, wherein m-th of subproblem is
By to obtaining simulation pre-coding matrix after the objective function optimizationAnd from which further follow that digital precode square Battle arrayAnd error εk
By error εkIt is compared with convergence tolerance, whenWhen, by corresponding simulation pre-coding matrixAnd number Pre-coding matrixAs optimal simulation pre-coding matrix and optimal digital precode matrix;WhenWhen, repeat on Step is stated, until error is less than or equal to convergence tolerance, or maximum number of iterations is completed, by corresponding simulation pre-coding matrixWith digital pre-coding matrixAs optimal simulation pre-coding matrix and optimal digital precode matrix.
Part connection architecture system is generated and issued according to optimal simulation pre-coding matrix and optimal digital precode matrix Mixed-beam.
Verify embodiment one:
Architecture system is connected for part of the invention, corresponding spectrum efficiency is
Wherein,It is signal-to-noise ratio, P indicates the average transmission power between base station and mobile station, Σ1It represents to angular moment First dimension of battle array Σ is Ns×NsBlock, be defined as
When what is sent in the channel is Gauss symbol, then spectrum efficiency is
Wherein,It is Ns×NsUnit vector, P is average transmission power, NsIt is number of data streams, RnIt is noise association side Poor matrix, Rn2FRFFBB, σ expression variance, HkIt is channel matrix;
Mixing Precoding Design can be obtained by the optimal solution of the optimization problem of one transmitting terminal of solution, this optimization The target of problem is that spectrum efficiency maximizes, i.e.,
Wherein,Be simulate pre-coding matrix in the i-th row, jth column element,It is digital precode matrix In the i-th row, jth column element.
System uses the narrowband the Saleh-Valenzuela group variety channel model of geometry.The aerial array of base station and user terminal It is all made of uniform straight line array, its array response vector is
If the receiving array response in the channel matrix of above formula is replaced with emission array response using uniform straight line array αULA(θ), wherein N indicates that array number in line array, λ indicate that carrier wavelength, d indicate the spacing between antenna, and angle/arrive is left in θ expression Up to the attitude angle at angle.
In emulation part, millimeter wave cluster channel model is not only considered, it is also contemplated that Rayleigh channel model.Believe in Rayleigh fading In road model, normalize channel matrix H in each element obey independent identically distributed mean value be 0, variance σ2Multiple height This distribution.
Numerical result shows that the performance of mixing method for precoding of the invention can be close to be assigned in high-dimensional digital wave beam The performance of shape system, and it is much higher than the performance of analog beam shaping system.The present invention simulates the gains such as Precoding Design is based on Transmission, digital precode matrix design are based on singular value decomposition.Final mixing precoder design is based on alternative optimization, connects it Nearly optimal no constraint matrix, so that performance of the system performance close to digital beamforming algorithm.
Verify embodiment two:
Verified in the present embodiment by Matlab emulation the spectrum efficiency of the alternative optimization algorithm based on matrix decomposition with The advantage of other algorithm performances.
The antenna number that sending and receiving end is equipped in the embodiment is 128 and 1, and the number of rf chain is all 4, number of data streams It is 2, using the millimeter wave channel model of ULA Array Model, each parameter occurrence is as described in Table 1 in the model.
Parameter assignment in 1 embodiment of table
As shown in figure 3, which show the spectrum efficiency comparison results of many algorithms.The performance of analog beam shaping system is most The performance of difference, digital beam forming system is best.The performance of algorithm proposed by the invention is substantially better than analog beam figuration Performance, and better than with algorithm performance in the prior art, closest to the performance of digital beam-forming system.
As shown in Figure 4, it is shown that influence of the number of data streams to algorithm spectrum efficiency.Compare be number of data streams be 1, 2, the spectrum efficiency of 4 and 8 lower proposed algorithms, as number of data streams increases, the spectrum efficiency of system increases, system performance It becomes better and better.

Claims (4)

1. it is a kind of based on part connection extensive multiple-input and multiple-output mixed-beam forming algorithm, which is characterized in that including with Lower step:
Given part connects the optimal without constraint precoder, primary simulation pre-coding matrix and convergence tolerance of architecture system, and Initial number pre-coding matrix is calculated without constraint precoder and primary simulation pre-coding matrix according to described, and then is calculated Initial error;
When the initial error is less than or equal to the convergence tolerance, by the primary simulation pre-coding matrix and initial number Optimal simulation pre-coding matrix and optimal digital precode matrix of the pre-coding matrix as part connection architecture system;
When the initial error is greater than the convergence tolerance, it is iterated calculating as target to minimize error, until error Less than or equal to the convergence tolerance, or maximum number of iterations is completed, corresponding simulation pre-coding matrix and number are prelisted Code matrix is as optimal simulation pre-coding matrix and optimal digital precode matrix;
The part connection architecture system generates simultaneously according to the optimal simulation pre-coding matrix and optimal digital precode matrix Issue mixed-beam.
2. a kind of extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection as described in claim 1, It optimal is obtained without constraint precoder by following steps it is characterized in that, described:
The channel state information that architecture system is connected according to the part, obtains channel matrix
Wherein, NcIndicate group variety number, NpIndicate the number of path in every cluster, αilIndicate the l articles transmission road in the i-th cluster reflector The gain factor of diameter obeys the multiple Gauss distribution that mean value is 0, variance is 1, for (i, l) strip diameter, θr,ilWithRespectively Azimuth and the pitch angle at angle, θ are left in expressiont,ilWithAzimuth and the pitch angle of angle of arrival are respectively indicated,WithRespectively indicate azimuth angle thetar,ilWithPitching angle thetat,ilWithCorresponding receiving array Response and emission array response;
Pass through H=U Σ VHSingular value decomposition is carried out to the channel matrix H, unit matrix V is obtained, from the unit matrix V It is proposed its preceding NsColumn obtain matrix V1;Pass throughDiagonal matrix Γ is calculated, and passes through F*=V1Γ obtains optimal without about Beam precoder F*;Wherein, U is unit matrix, and Σ is diagonal matrix.
3. a kind of extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection as claimed in claim 2, It is characterized in that, passing throughObtain initial number pre-coding matrixAnd pass throughObtain initial error ε0
4. a kind of extensive multiple-input and multiple-output mixed-beam figuration based on part connection is calculated as claimed in claim 2 or claim 3 Method, which is characterized in that the detailed process of the iterative calculation are as follows:
Set objective functionConversion | | F*-FRFFBB||FIt can obtainWherein, TnFor F*In n-th of matrix-block, fnIndicate composition FRFIn n-th of square Battle array block,Indicate FBBLine n vector;
Then to objective function optimize for
Wherein, tn,mIt indicatesM column,It indicatesM-th of element,The F obtained for kth time iterationRFIn N matrix-block,To indicateM-th of element phase increment,
For F after kth time iterationBBLine n vector;It indicatesPhase increment is averaged Value;
By to obtaining simulation pre-coding matrix after the objective function optimizationAnd from which further follow that digital precode matrix And error εk
By error εkIt is compared with convergence tolerance, whenWhen, by corresponding simulation pre-coding matrixIt prelists with number Code matrixAs optimal simulation pre-coding matrix and optimal digital precode matrix;WhenWhen, repeat above-mentioned step Suddenly, until error is less than or equal to the convergence tolerance, or maximum number of iterations is completed, by corresponding simulation pre-coding matrixWith digital pre-coding matrixAs optimal simulation pre-coding matrix and optimal digital precode matrix.
CN201910400887.6A 2019-05-15 2019-05-15 Large-scale multi-input multi-output hybrid beam forming algorithm based on partial connection Active CN110138427B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910400887.6A CN110138427B (en) 2019-05-15 2019-05-15 Large-scale multi-input multi-output hybrid beam forming algorithm based on partial connection

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910400887.6A CN110138427B (en) 2019-05-15 2019-05-15 Large-scale multi-input multi-output hybrid beam forming algorithm based on partial connection

Publications (2)

Publication Number Publication Date
CN110138427A true CN110138427A (en) 2019-08-16
CN110138427B CN110138427B (en) 2020-09-08

Family

ID=67574168

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910400887.6A Active CN110138427B (en) 2019-05-15 2019-05-15 Large-scale multi-input multi-output hybrid beam forming algorithm based on partial connection

Country Status (1)

Country Link
CN (1) CN110138427B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111277311A (en) * 2020-02-10 2020-06-12 电子科技大学 Active and passive combined beam forming design method for millimeter wave symbiotic communication system
CN111726144A (en) * 2020-06-24 2020-09-29 中南大学 Hybrid precoding design method, device, medium and equipment based on initial value optimization
CN112039565A (en) * 2020-09-11 2020-12-04 成都大学 Large-scale MIMO mixed pre-coding method based on distributed part connection
CN112910521A (en) * 2021-02-27 2021-06-04 中电万维信息技术有限责任公司 Deep learning-based MIMO mixed beam forming method
CN114598363A (en) * 2022-02-09 2022-06-07 国网电力科学研究院有限公司 Large-scale MIMO (multiple input multiple output) safety precoding method for 3D (three-dimensional) space wireless channel

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107135024A (en) * 2017-05-12 2017-09-05 西安交通大学 A kind of mixed-beam figuration Iterative Design method of low complex degree
CN108599825A (en) * 2018-02-12 2018-09-28 大连理工大学 A kind of hybrid coding method based on MIMO-OFDM millimeters of wave structures
CN108736943A (en) * 2018-05-22 2018-11-02 湘潭大学 A kind of mixing method for precoding suitable for extensive mimo system
EP3404843A1 (en) * 2017-05-17 2018-11-21 Mitsubishi Electric R&D Centre Europe B.V. Method for enabling both analog and digital beamforming
CN109450501A (en) * 2018-12-14 2019-03-08 东南大学 A kind of extensive mimo system mixing Precoding Design method in broadband

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107135024A (en) * 2017-05-12 2017-09-05 西安交通大学 A kind of mixed-beam figuration Iterative Design method of low complex degree
EP3404843A1 (en) * 2017-05-17 2018-11-21 Mitsubishi Electric R&D Centre Europe B.V. Method for enabling both analog and digital beamforming
CN108599825A (en) * 2018-02-12 2018-09-28 大连理工大学 A kind of hybrid coding method based on MIMO-OFDM millimeters of wave structures
CN108736943A (en) * 2018-05-22 2018-11-02 湘潭大学 A kind of mixing method for precoding suitable for extensive mimo system
CN109450501A (en) * 2018-12-14 2019-03-08 东南大学 A kind of extensive mimo system mixing Precoding Design method in broadband

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
ANNA VIZZIELLO 等: "A Kalman Based Hybrid Precoding for Multi-User Millimeter Wave MIMO Systems", 《IEEE ACCESS ( VOLUME: 6 )》 *
CMCC: "Hybrid Beamforming for Massive MIMO", 《3GPP TSG RAN WG1 #85 R1-164893》 *
YANG ZHANG 等: "On Precoder Design for Amplify-and-Forward MIMO Relay Systems", 《2011 IEEE VEHICULAR TECHNOLOGY CONFERENCE (VTC FALL)》 *
丁青锋 等: "毫米波MIMO系统离散化混合预编码设计", 《计算机应用研究》 *
向建伟 等: "毫米波大规模MIMO系统中低复杂度混合预编码方法", 《电信科学 2016年第9期》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111277311A (en) * 2020-02-10 2020-06-12 电子科技大学 Active and passive combined beam forming design method for millimeter wave symbiotic communication system
CN111277311B (en) * 2020-02-10 2022-03-25 电子科技大学 Active and passive combined beam forming design method for millimeter wave symbiotic communication system
CN111726144A (en) * 2020-06-24 2020-09-29 中南大学 Hybrid precoding design method, device, medium and equipment based on initial value optimization
CN112039565A (en) * 2020-09-11 2020-12-04 成都大学 Large-scale MIMO mixed pre-coding method based on distributed part connection
CN112910521A (en) * 2021-02-27 2021-06-04 中电万维信息技术有限责任公司 Deep learning-based MIMO mixed beam forming method
CN112910521B (en) * 2021-02-27 2022-04-05 中电万维信息技术有限责任公司 Deep learning-based MIMO mixed beam forming method
CN114598363A (en) * 2022-02-09 2022-06-07 国网电力科学研究院有限公司 Large-scale MIMO (multiple input multiple output) safety precoding method for 3D (three-dimensional) space wireless channel

Also Published As

Publication number Publication date
CN110138427B (en) 2020-09-08

Similar Documents

Publication Publication Date Title
CN110138427A (en) Extensive multiple-input and multiple-output mixed-beam forming algorithm based on part connection
CN111294095B (en) IRS (inter-range instrumentation Standard) assisted large-scale MIMO (multiple input multiple output) wireless transmission method based on statistical CSI (channel State information)
CN111181619B (en) Millimeter wave hybrid beam forming design method based on deep reinforcement learning
CN107359921B (en) Mixed precoding method of large-scale MIMO system based on standard orthogonalization
CN107809274B (en) Hybrid precoding method based on novel phase-shifting switch network
CN107135024A (en) A kind of mixed-beam figuration Iterative Design method of low complex degree
CN109714091B (en) Iterative hybrid precoding method based on hierarchical design in millimeter wave MIMO system
CN110166103A (en) A kind of novel mixed-beam shaped structure and setting method of millimeter wave MU-MISO system
CN110943768B (en) Mixed precoding codebook joint design method of millimeter wave large-scale MIMO system
CN110365388B (en) Low-complexity millimeter wave multicast beam forming method
CN106031068A (en) Beamforming-based communication method and apparatus
CN105429686A (en) Transmission device and method for separated type asymmetric hybrid beam formation
CN113114317B (en) IRS-assisted phase shift optimization method for downlink multi-user communication system
CN113193893B (en) Millimeter wave large-scale MIMO intelligent hybrid beam forming design method
Zhang et al. Near-optimal design for hybrid beamforming in mmWave massive multi-user MIMO systems
CN110011712A (en) A kind of mixing method for precoding towards the extensive multiple-input and multiple-output of millimeter wave
CN113098574B (en) Irregular subarray efficient solving method for large-scale MIMO system
Ren et al. Machine learning-based hybrid precoding with robust error for UAV mmWave massive MIMO
CN109981154A (en) Low complex degree array antenna multi-input multi-output system mixing precoding algorithms
CN114465643B (en) Mixed precoding method of millimeter wave large-scale MIMO antenna system based on gradient descent method
CN112910520A (en) Convolutional neural network-based MIMO system beam training method
CN109067446B (en) Mixed precoding method for multi-antenna multi-user large-scale antenna
CN112398513A (en) Beam forming method of massive MIMO system
CN107104719A (en) A kind of millimeter wave digital analog mixed Precoding Design method based on geometrical construction
CN104821840A (en) Anti-interference method of massive multiple input multiple output downlink system

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210125

Address after: 510000 3rd floor, building B2, No. 11, Kaiyuan Avenue, Science City, Guangzhou high tech Industrial Development Zone, Guangdong Province

Patentee after: Guangzhou aeters communication equipment Co.,Ltd.

Address before: 710054 No. 58 Yanta Middle Road, Beilin District, Xi'an City, Shaanxi Province

Patentee before: XI'AN University OF SCIENCE AND TECHNOLOGY

CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 510000 3rd floor, building B2, No. 11, Kaiyuan Avenue, Science City, Guangzhou high tech Industrial Development Zone, Guangdong Province

Patentee after: GUANGZHOU ITS COMMUNICATION EQUIPMENT Co.,Ltd.

Address before: 510000 3rd floor, building B2, No. 11, Kaiyuan Avenue, Science City, Guangzhou high tech Industrial Development Zone, Guangdong Province

Patentee before: Guangzhou aeters communication equipment Co.,Ltd.