CN106059685A - Large scale MIMO channel simulation apparatus under time evolution and simulation method thereof - Google Patents
Large scale MIMO channel simulation apparatus under time evolution and simulation method thereof Download PDFInfo
- Publication number
- CN106059685A CN106059685A CN201610327570.0A CN201610327570A CN106059685A CN 106059685 A CN106059685 A CN 106059685A CN 201610327570 A CN201610327570 A CN 201610327570A CN 106059685 A CN106059685 A CN 106059685A
- Authority
- CN
- China
- Prior art keywords
- channel
- unit
- signal
- subsystem
- fading
- 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
- 238000004088 simulation Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims description 22
- 238000005562 fading Methods 0.000 claims abstract description 38
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 238000006243 chemical reaction Methods 0.000 claims abstract description 31
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 21
- 230000006854 communication Effects 0.000 claims abstract description 16
- 238000004891 communication Methods 0.000 claims abstract description 15
- 238000003786 synthesis reaction Methods 0.000 claims description 30
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000007423 decrease Effects 0.000 claims description 10
- 230000003993 interaction Effects 0.000 claims description 10
- 238000010606 normalization Methods 0.000 claims description 7
- 102100040862 Dual specificity protein kinase CLK1 Human genes 0.000 claims description 6
- 238000009826 distribution Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 4
- 238000013139 quantization Methods 0.000 claims description 4
- 102100040844 Dual specificity protein kinase CLK2 Human genes 0.000 claims description 3
- 102100040856 Dual specificity protein kinase CLK3 Human genes 0.000 claims description 3
- 101000749294 Homo sapiens Dual specificity protein kinase CLK1 Proteins 0.000 claims description 3
- 101000749291 Homo sapiens Dual specificity protein kinase CLK2 Proteins 0.000 claims description 3
- 101000749304 Homo sapiens Dual specificity protein kinase CLK3 Proteins 0.000 claims description 3
- 238000005070 sampling Methods 0.000 abstract description 4
- 230000002194 synthesizing effect Effects 0.000 abstract 1
- 230000005540 biological transmission Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000010295 mobile communication Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000005039 memory span Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/0082—Monitoring; Testing using service channels; using auxiliary channels
- H04B17/0087—Monitoring; Testing using service channels; using auxiliary channels using auxiliary channels or channel simulators
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3911—Fading models or fading generators
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Radio Transmission System (AREA)
Abstract
The invention discloses a large scale MIMO channel simulation apparatus under time evolution. According to the invention, a user inputs communication scenario parameters to a master control PC subsystem, and a parameter computing unit computes fading factors of each MIMO sub-channel. A channel fading matrix is separately transmitted to a channel fading RAM of each channel simulation unit. A channel input signal undergoes a signal decomposition unit of a signal decomposition sub-system to obtain a plurality group of simulation signals. The simulation signals undergo an AD conversion unit of a signal processing sub-system and convert to digital signals, and the digital signals are transmitted to the channel simulation units, and at the same time a filter RAM is read through different sampling rates to obtain MIMO channel fading under time evolution. The channel simulation units conduct multiplication, adding and outputting on the digital signals and the channel fading. Each component of a channel output signal undergoes the AD conversion unit of the signal processing sub-system, is connected to a signal synthesizing sub-system, and then undergoes a noise superposing unit to obtain a final channel output signal.
Description
Technical field:
The present invention relates to a kind of time evolution extensive mimo channel analog and analogy method thereof, belong to wireless messages
Transmission field, propagates the extensive mimo channel analogy method under scene particular for dynamic time evolution.
Background technology:
Radio communication has become as a part indispensable in people's daily life, no matter and wireless communication system is setting
The meter development phase, or in the planning of communication network and deployment, the impact of wireless propagation environment will be considered.Such as, radio
Ripple signal is affected by landform, landforms and transmission range in communication process exists path loss, due to reflection, scattering and diffraction
Multipath fading also can be produced etc. factor.Meanwhile, the movement of sending and receiving end makes electric wave signal produce Doppler frequency shift, additionally, by week
Enclose building and impact that landforms block, it is also possible to there is shadow fading.In order to effectively assess and validation of wireless communication system
Scheme and performance, reduce R&D costs, shorten exploitation test period, the electromagnetism needing to simulate in laboratory environments complexity passes
Broadcast environment, i.e. develop the hardware simulator possessing real-time channel simulator function.
Multiple-input and multiple-output (Multiple-input Multiple-out, MIMO) technology can improve the channel of system
Capacity, the rate of information throughput and reliability, it has also become one of core technology of forth generation mobile communication system (4G).But, with
The universal of intelligent terminal and move new business demand sustainable growth, wireless transmission rate demand is exponentially increased, it is contemplated that
The wireless communication transmissions rate requirement of the year two thousand twenty will be thousand times of currently operation system.To this end, next-generation mobile communications system
System (5G) considers the number by dramatically increasing base station side configuration antenna, i.e. uses many at present with large scale array antenna replacement
Antenna, thus form extensive mimo wireless communication environment, with the Radio Resource of degree of depth digging utilization Spatial Dimension, solve future
The spectrum efficiency of mobile communication and power efficiency problem.Meanwhile, sharply increasing along with number of antennas, traditional MIMO
Channel simulator will be difficult to completion system test checking work.
On the other hand, along with movement-movement (Mobile-to-Mobile, M2M) and high-speed railway (High Speed
Train, HST) etc. the going deep into of high-speed mobile scene lower channel measurement research.A large amount of measured datas show, transceiver communication end is all located
Under the communication scenes of mobile status, the statistical property such as the multipath fading distribution of radio propagation channel, doppler spectral and delay profile
It is respectively provided with time-varying characteristics, belongs to time-varying non-stationary fading channel.Therefore the master pattern for following 5G channel adopts non-stationary
The thought of Channel Modeling, this also causes the mimo channel simulator for the steady fading channel of tradition will be the most applicable.
Therefore, developing the hardware simulator of non stationary channel under mobile time-varying scene is to solve following extensive MIMO communication
System algorithm checking and the important step of Performance Evaluation.But, when number of antennas is the most, need to solve MT×MRIndividual MIMO
Subchannel process resource limit problem in hardware realizes.On the other hand, the data interaction between each subchannel, will cause hard
There is more serious interconnection resources and limit in part.
Summary of the invention:
The present invention provide a kind of for the time evolution extensive mimo channel analog under the mobile context of sending and receiving end and
Its analogy method, the method can be with the mimo channel decline of any number of antennas under real-time Simulation reproduction dynamic scene.
The present invention adopts the following technical scheme that a kind of time evolution extensive mimo channel analog, including master control PC
Subsystem, signal decomposition subsystem, signal processing subsystem and signal syntheses subsystem；Described master control PC subsystem includes user
Interactive unit and parameter calculation unit；Described signal decomposition subsystem includes signal decomposition unit；Described signal processing subsystem
Including AD conversion unit, channel simulation unit and D/A conversion unit；Described signal syntheses subsystem includes signal syntheses list
Unit and noise superpositing unit；The input interface of the output interface of described parameter calculation unit and described channel simulation unit is with PCIE
Bus is connected；The described output interface of signal decomposition unit is connected with the input interface of AD conversion unit；Described analog digital conversion
The output interface of unit is connected with the input interface of channel simulation unit；The output interface of described channel simulation unit turns with digital-to-analogue
The input interface changing unit is connected；The output interface of described D/A conversion unit is connected with the input interface of signal synthesis unit；
The output interface of described signal synthesis unit is connected with noise superpositing unit.
The present invention also adopts the following technical scheme that the simulation side of a kind of time evolution extensive mimo channel analog
Method, comprises the steps:
The first step, user is on master control PC subsystem, and by user interaction unit input communication scenario parameters, parameter calculates
Unit, according to input parameter, calculates MIMO each subchannel fading factor；
Second step, is respectively transmitted to each channel simulation list by PCIE high speed data bus interface by channel fading matrix
In the channel fading RAM of unit；
3rd step, channel input signal, through the signal decomposition unit of signal decomposition subsystem, obtains some groups of simulation letters
Number；
4th step, analogue signal becomes digital signal by the AD conversion unit on signal processing subsystem, is transferred to
Channel simulation unit, reads wave filter RAM by different sample rates simultaneously, it is thus achieved that the mimo channel decline of time evolution；
5th step, digital signal is multiplied accumulating with time evolution channel fading, obtains channel by channel simulation unit mutually
Each digital component of output signal；
6th step, each component of channel output signal, through the D/A conversion unit of signal processing subsystem, obtains letter
Each analogue component of road output signal；
7th step, is connected to a signal syntheses subsystem, Jing Guoxin by important for the institute of same channel output signal
The synthesis of number synthesis unit, then obtain final channel output signal through noise superpositing unit.
Further, specifically to produce step as follows for the first step:
1) calculate path loss α (t) of each subchannel, and carry out fixed point quantization；
2) the multiple gauss potential function G of independent same distribution of power normalization is produced；
3) the cross-correlation matrix R of the acquisition each subchannel of MIMO is calculatedH；
4) if RHFor positively definite matrix, utilize Cholesky to decompose and obtain lower triangle coefficient matrix L, otherwise use eigenvalue to divide
Solution method solves L；
5) utilizeAnd after fixed point, obtain normalization fading matrix；
6) 1 is repeated)-5) produce L sub-channels fading matrix.
Further, to implement step as follows for the 4th step
1) sample rate inputted by user interaction unit, parameter calculation unit obtains first order interpolation multiple, the second level
Interpolation multiple, third level interpolation multiple, it is transferred to wave filter RAM；
2) producing 4 any clock sources by high accuracy DDS, CLK1 is the clock reading RAM, and CLK2 is first order interpolation
Clock, CLK3 is interpolated clock in the second level, and CLK 4 is interpolated clock in the third level, and the phase place of 4 clock sources is identical；
3) first order interpolation uses high performance multiphase filter structure to realize 1～100 interpolations；The wave filter of multiphase filter
Coefficient is calculated by parameter calculation unit, and is stored in wave filter RAM；
4) second level uses half-band filter to realize 1～2 times of interpolation；
5) third level uses cic filter to realize 1～5 times of interpolation；Employ 5 depositor 5 data of storage, by the 3rd
Level interpolation multiple controls cumulative register number, i.e. controls filter order；
6) by changing sample rate, 1 is repeated)-5) i.e. can obtain time evolution and decline continuously.
There is advantages that
(1) extensive mimo channel is decomposed into fundamental matrix computing unit, and two inputs two using FPGA to be core
Output primary channel analogue unit, builds a kind of simulator being applicable to extensive mimo channel, and this simulator is prone to hardware and expands
Exhibition, can be used for the channel hardware simulation of any number of antennas；
(2) have employed the technology that signal decomposes at analog domain and merges, solve the simulation of large-scale F PGA primary channel single
Signal stationary problem between unit；
(3) propose a kind of scheme reading channel fading based on variable sampling rate, thus realize channel output decline
Time evolution effect, had both improve system processing capability in real time, had solved again the problem that memory span is limited.
Accompanying drawing illustrates:
Fig. 1 is mimo system illustraton of model.
Fig. 2 is the hardware implementations figure that apparatus of the present invention use.
Fig. 3 is the production method figure of present invention time evolution channel fading.
Detailed description of the invention:
Consider a configuration MTIndividual transmitting antenna and MRThe mimo system (as shown in Figure 1) of individual reception antenna, channel input-
Output relation is
Y (t)=H (τ, t) * x (t)+n (t) (1)
Wherein,For launching signal phasor；For
Interchannel noise vector；For receiving signal phasor；The multipath number of clusters mesh assuming each subchannel is
L, then mimo channel model is represented by
Wherein, hi,j(τ t) represents that jth launches the impulse response between antenna and i-th reception antenna, τlRepresent l
The time delay in footpath, L represents discretization multipath number of clusters mesh.
For above-mentioned mimo channel model, conventional simulator hardware implementation mode includes two kinds: prestoring type and real-time formula, its
The generation algorithm of middle prestoring type, amount of storage is big, it is difficult to support time-varying dynamic scene；The generation algorithm of formula in real time, operand is big, difficult
To support extensive multichannel occasion.
The present invention relates to a kind of time evolution extensive mimo channel analog, the method initially with matrix decomposition changes
Write-channel model representation form, order
And
Formula (2) is rewritten as
I.e.
For formula (6) model form, present invention time evolution extensive mimo channel analog includes master control PC subsystem
System, signal decomposition subsystem, signal processing subsystem and signal syntheses subsystem (as shown in Figure 2), master control PC subsystem includes
User interaction unit 1-1 and parameter calculation unit 1-2；Signal decomposition subsystem includes signal decomposition unit 1-3；Signal processing
System includes AD conversion unit 1-4, channel simulation unit 1-5 and D/A conversion unit 1-6；Signal syntheses subsystem includes letter
Number synthesis unit 1-7 and noise superpositing unit 1-8；The output interface of parameter calculation unit 1-2 is defeated with channel simulation unit 1-5's
Incoming interface is connected with PCIE bus；The output interface of signal decomposition unit 1-3 and the input interface phase of AD conversion unit 1-4
Even；The output interface of AD conversion unit 1-4 is connected with the input interface of channel simulation unit 1-5；Channel simulation unit 1-5's
Output interface is connected with the input interface of D/A conversion unit 1-6；The output interface of D/A conversion unit 1-6 and signal syntheses list
The input interface of unit 1-7 is connected；The output interface of signal synthesis unit 1-7 is connected with noise superpositing unit 1-8.
The handling process of present invention time evolution extensive mimo channel analog internal signal is as follows: channel input letter
Number (xn-1,xn) be decomposed into through signal decomposition subsystem and to organize analogue signal more；In signal processing subsystem, each group simulation letter
Number being input to channel simulation unit 1-5 through AD conversion unit 1-4, the channel transmitted with parameter calculation unit 1-2 declines
Fall superposition, and finally by D/A conversion unit 1-6, output comprises (xn-1,xn) informationComponent；By same output
The institute of signal is important is connected to a signal syntheses subsystem, is synthesized by signal synthesis unit 1-7, then by noise superposition
Unit 1-8, finally gives channel output signal
The analogy method of present invention time evolution extensive mimo channel analog, comprises the steps:
The first step, user is on master control PC subsystem, by user interaction unit 1-1 input communication scenario parameters, channel
Parameter calculation unit 1-2, according to input parameter, calculates MIMO each subchannel fading factor, and concrete generation step is as follows:
(1) calculate path loss α (t) of each subchannel, and carry out fixed point quantization；
(2) the multiple gauss potential function G of independent same distribution of power normalization is produced；
(3) the cross-correlation matrix R of the acquisition each subchannel of MIMO is calculatedH；
(4) if RHFor positively definite matrix, utilize Cholesky to decompose and obtain lower triangle coefficient matrix L, otherwise use eigenvalue
Decomposition method solves L；
(5) utilizeAnd after fixed point, obtain normalization fading matrix；
(6) 1 is repeated)-5) produce L sub-channels fading matrix；
Second step, is respectively transmitted to each channel simulation list by PCIE high speed data bus interface by channel fading matrix
In the channel fading RAM of unit 1-5；
3rd step, channel input signal, through the signal decomposition unit 1-3 of signal decomposition subsystem, obtains some groups of simulations
Signal；
4th step, analogue signal becomes digital signal by AD conversion unit 1-4 on signal processing subsystem, transmission
To channel simulation unit 1-5, read wave filter RAM by different sample rates, it is thus achieved that the mimo channel decline of time evolution simultaneously,
Implement step following (as shown in Figure 3):
1) by user interaction unit 1-1 input sampling rate, parameter calculation unit 1-2 obtains first order interpolation multiple, the
Two grades of interpolation multiples, third level interpolation multiple, it is transferred to wave filter RAM；
2) producing 4 any clock sources, CLK1 is the clock reading RAM, and CLK2 is interpolated clock in the first order, and CLK3 is the
Interpolated clock in two grades, CLK 4 is interpolated clock in the third level, and the phase place of 4 clock sources is identical；
3) first order interpolation uses high performance multiphase filter structure to realize 1～100 interpolations；The wave filter of multiphase filter
Coefficient is calculated by parameter calculation unit 1-2, and is stored in wave filter RAM；
4) second level uses half-band filter to realize 1～2 times of interpolation；
5) third level uses cic filter to realize 1～5 times of interpolation；Employ 5 depositor 5 data of storage, by the 3rd
Level interpolation multiple controls cumulative register number, i.e. controls filter order；
6) by changing sample rate, 1 is repeated)-5) i.e. can obtain time evolution and decline continuously；
5th step, channel simulation unit 1-5 is by digital signal and being multiplied shown in time evolution channel fading progressive form (6)
Cumulative, obtain each digital component of channel output signal；
6th step, each component of channel output signal, through the D/A conversion unit 1-6 of signal processing subsystem, obtains
Each analogue component of channel output signal；
7th step, is connected to a signal syntheses subsystem, Jing Guoxin by important for the institute of same channel output signal
Number synthesis unit 1-7 synthesis, then obtain final channel output signal through noise superpositing unit 1-8.
For making the object, technical solutions and advantages of the present invention clearer, below as a example by the 3 footpath mimo channels of 8 × 8 also
In conjunction with the accompanying drawing of the present invention, technical scheme is carried out clear, complete description.
The method utilizing matrix decomposition of the present invention, the 3 footpath mimo channels of 8 × 8 can be expressed as
I.e.
Visible, the mimo channel simulation of 8 × 8 needs 4 signal decomposition subsystems, 16 signal processing subsystems and 4
Signal syntheses subsystem；Wherein 8 path channels input signals are resolved into 32 tunnel analogue signals by 4 signal decomposition subsystems, and pass
It is passed to 16 signal processing subsystems；The digital signal that process is obtained by AD conversion unit 1-4 on signal processing subsystem
Being input to channel simulation unit 1-5, channel simulation unit 1-5 and process the matrix operations in formula (7), the digital signal obtained passes through
D/A conversion unit 1-6 obtains 32 tunnel analogue signals, and the component of same channel output signal is connected to a signal respectively
In synthon system, synthesize through signal synthesis unit 1-7, then by noise laminating module 1-8, the certain letter of final output
Make an uproar than y1～y8Signal.
It is embodied as step as follows:
The first step, user is on master control PC, by user interaction unit 1-1 input communication scenario parameters, channel parameter
Computing unit 1-2, according to input parameter, calculates MIMO each subchannel fading factor, produces step as follows:
1) calculate path loss α (t) according to formula (9), and carry out fixed point quantization；
α (t)=32.44+20lg (fMHz)+20lg(dkm) (9)
In formula, fMHz,dkmRepresent that frequency and communication distance and unit use MHz, km respectively；
2) with reference to non-patent literature 1 [PATZOLD M, Mobile fading channel [M] .New York:
Wiley.2002:51-81.] the harmonic superposition method in produces the independent same distribution channel fading matrix G of power normalization；
Wherein,N represents scattering number of branches, is taken as 32, f hereindTable
Show Doppler frequency, αn,ij,φnRefer to the random angle of incidence of each scattering branch road and initial phase；
3) transmitting antenna related matrix R that user inputs is utilizedTxWith reception antenna correlation matrix RRxKronecker amass
Calculate the space cross-correlation matrix R of mimo channelH；
4) if RHFor positively definite matrix, utilize Cholesky to decompose and obtain lower triangle coefficient matrix L, otherwise use feature
Value decomposition method solves L；
5) vectorization of usage factor matrix L and channel matrix G is calculated MIMO correlated fading channels matrix H, i.e.
Vec (H)=L vec (G)=RH 1/2·[g11…g18…g81…g88]T (11)
Second step, is respectively transmitted to each channel simulation list by PCIE high speed data bus interface by channel fading matrix
In the channel fading RAM of unit 1-5；
3rd step, channel input signal, through the signal decomposition unit 1-3 of signal decomposition subsystem, obtains some groups of simulations
Signal；
4th step, analogue signal becomes digital signal by AD conversion unit 1-4 on signal processing subsystem, transmission
To channel simulation unit 1-5, read wave filter RAM by different sample rates, it is thus achieved that the mimo channel decline of time evolution simultaneously,
Implement step as follows:
1) by User Interface 1-1 input sampling rate 100KHz, calculating first order interpolation multiple is 100 times, the
Two grades of interpolation multiples are 2 times, and third level interpolation multiple is 5 times, are transferred to channel simulation unit 1-5；
2) producing any clock source, clock source 1 is the clock 100KHz reading RAM, and clock source 2 is interpolated clock in the first order
10MHz, clock source 3 is interpolated clock 20MHz in the second level, and clock source 4 is interpolated clock 100MHz in the third level；
3) first order interpolation uses high performance multiphase filter structure to realize 100 times of interpolations；The wave filter of multiphase filter
Coefficient is calculated by parameter calculation unit 1-2, and is stored in wave filter RAM；
4) second level uses half-band filter to realize 2 times of interpolations；
5) third level uses cic filter, uses 5 depositor 5 data of storage, exports, i.e. obtain in 5 times after adding up
The effect inserted, the final channel fading exporting 1000 times of interpolations；
6) by changing sample rate, 1 is repeated)-5) i.e. can obtain the channel fading of time evolution；
5th step, digital signal is multiplied accumulating mutually by channel simulation unit 1-5 with the decline of time evolution continuous channel,
Each digital component to channel output signal；
6th step, each component of channel output signal, through the D/A conversion unit 1-6 of signal processing subsystem, obtains
Each analogue component of channel output signal；
7th step, is connected to a signal syntheses subsystem, Jing Guoxin by important for the institute of same channel output signal
Number synthesis unit 1-7 synthesis, then obtain final channel output signal through noise superpositing unit 1-8.
The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For Yuan, can also make some improvement under the premise without departing from the principles of the invention, these improvement also should be regarded as the present invention's
Protection domain.
Claims (4)
1. a time evolution extensive mimo channel analog, it is characterised in that: include master control PC subsystem, signal decomposition
Subsystem, signal processing subsystem and signal syntheses subsystem；Described master control PC subsystem include user interaction unit (1-1) and
Parameter calculation unit (1-2)；Described signal decomposition subsystem includes signal decomposition unit (1-3)；Described signal processing subsystem
Including AD conversion unit (1-4), channel simulation unit (1-5) and D/A conversion unit (1-6)；Described signal syntheses subsystem
Including signal synthesis unit (1-7) and noise superpositing unit (1-8)；The output interface of described parameter calculation unit (1-2) and institute
The input interface stating channel simulation unit (1-5) is connected with PCIE bus；The output interface of described signal decomposition unit (1-3) with
The input interface of AD conversion unit (1-4) is connected；The output interface of described AD conversion unit (1-4) and channel simulation unit
(1-5) input interface is connected；The output interface of described channel simulation unit (1-5) and the input of D/A conversion unit (1-6)
Interface is connected；The output interface of described D/A conversion unit (1-6) is connected with the input interface of signal synthesis unit (1-7)；Institute
The output interface stating signal synthesis unit (1-7) is connected with noise superpositing unit (1-8).
2. the analogy method of a time evolution extensive mimo channel analog, it is characterised in that: comprise the steps
The first step, user is on master control PC subsystem, by user interaction unit (1-1) input communication scenario parameters, parameter meter
Calculate unit (1-2) according to input parameter, calculating MIMO each subchannel fading factor；
Second step, is respectively transmitted to each channel simulation unit by PCIE high speed data bus interface by channel fading matrix
(1-5) in channel fading RAM；
3rd step, channel input signal, through the signal decomposition unit (1-3) of signal decomposition subsystem, obtains some groups of simulation letters
Number；
4th step, analogue signal becomes digital signal by the AD conversion unit (1-4) on signal processing subsystem, is transferred to
Channel simulation unit (1-5), reads wave filter RAM by different sample rates simultaneously, it is thus achieved that the mimo channel decline of time evolution；
5th step, digital signal is multiplied accumulating with time evolution channel fading, obtains channel by channel simulation unit (1-5) mutually
Each digital component of output signal；
6th step, each component of channel output signal, through the D/A conversion unit (1-6) of signal processing subsystem, obtains letter
Each analogue component of road output signal；
7th step, is connected to a signal syntheses subsystem by important for the institute of same channel output signal, closes through signal
Become unit (1-7) synthesis, then obtain final channel output signal through noise superpositing unit (1-8).
3. the analogy method of time evolution extensive mimo channel analog as claimed in claim 2, it is characterised in that: the
It is as follows that one step specifically produces step:
1) calculate path loss α (t) of each subchannel, and carry out fixed point quantization；
2) the multiple gauss potential function G of independent same distribution of power normalization is produced；
3) the cross-correlation matrix R of the acquisition each subchannel of MIMO is calculatedH；
4) if RHFor positively definite matrix, utilize Cholesky to decompose and obtain lower triangle coefficient matrix L, otherwise use Eigenvalues Decomposition side
Method solves L；
5) utilizeAnd after fixed point, obtain normalization fading matrix；
6) 1 is repeated)-5) produce L sub-channels fading matrix.
4. the analogy method of time evolution extensive mimo channel analog as claimed in claim 3, it is characterised in that: the
It is as follows that four steps implement step
1) sample rate inputted by user interaction unit (1-1), parameter calculation unit (1-2) obtains first order interpolation multiple,
Second level interpolation multiple, third level interpolation multiple, it is transferred to wave filter RAM；
2) producing 4 any clock sources by high accuracy DDS, CLK1 is the clock reading RAM, when CLK2 is first order interpolation
Clock, CLK3 is interpolated clock in the second level, and CLK 4 is interpolated clock in the third level, and the phase place of 4 clock sources is identical；
3) first order interpolation uses high performance multiphase filter structure to realize 1～100 interpolations；The filter coefficient of multiphase filter
Calculated by parameter calculation unit (1-2), and be stored in wave filter RAM；
4) second level uses half-band filter to realize 1～2 times of interpolation；
5) third level uses cic filter to realize 1～5 times of interpolation；Employ 5 depositor 5 data of storage, by the third level
Insert multiple and control cumulative register number, i.e. control filter order；
6) by changing sample rate, 1 is repeated)-5) i.e. can obtain time evolution and decline continuously.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610327570.0A CN106059685B (en) | 2016-05-17 | 2016-05-17 | The extensive mimo channel simulator of time evolution and its analogy method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610327570.0A CN106059685B (en) | 2016-05-17 | 2016-05-17 | The extensive mimo channel simulator of time evolution and its analogy method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106059685A true CN106059685A (en) | 2016-10-26 |
| CN106059685B CN106059685B (en) | 2019-02-26 |
Family
ID=57176995
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610327570.0A Expired - Fee Related CN106059685B (en) | 2016-05-17 | 2016-05-17 | The extensive mimo channel simulator of time evolution and its analogy method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106059685B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109067482A (en) * | 2018-08-14 | 2018-12-21 | 南京航空航天大学 | Reconfigurable network channel simulation method and device towards car networking communication |
| CN112655154A (en) * | 2018-08-30 | 2021-04-13 | 华为技术加拿大有限公司 | Method and system for linear signal processing using signal decomposition |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1988424A (en) * | 2005-12-23 | 2007-06-27 | 中兴通讯股份有限公司 | Channel simulator |
| CN101217324A (en) * | 2008-01-04 | 2008-07-09 | 中兴通讯股份有限公司 | A test system and method of WIMAX upstream cooperation MIMO |
| CN102801482A (en) * | 2011-05-26 | 2012-11-28 | 中兴通讯股份有限公司 | Device, method and system for dynamic range adjustment for channel simulation system |
| CN105049142A (en) * | 2015-07-16 | 2015-11-11 | 中国电子科技集团公司第四十一研究所 | Dual-path static baseband channel simulating device and method |
| CN105187232A (en) * | 2015-07-16 | 2015-12-23 | 中国电子科技集团公司第四十一研究所 | Dynamic channel simulation device and method |
-
2016
- 2016-05-17 CN CN201610327570.0A patent/CN106059685B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1988424A (en) * | 2005-12-23 | 2007-06-27 | 中兴通讯股份有限公司 | Channel simulator |
| CN101217324A (en) * | 2008-01-04 | 2008-07-09 | 中兴通讯股份有限公司 | A test system and method of WIMAX upstream cooperation MIMO |
| CN102801482A (en) * | 2011-05-26 | 2012-11-28 | 中兴通讯股份有限公司 | Device, method and system for dynamic range adjustment for channel simulation system |
| CN105049142A (en) * | 2015-07-16 | 2015-11-11 | 中国电子科技集团公司第四十一研究所 | Dual-path static baseband channel simulating device and method |
| CN105187232A (en) * | 2015-07-16 | 2015-12-23 | 中国电子科技集团公司第四十一研究所 | Dynamic channel simulation device and method |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109067482A (en) * | 2018-08-14 | 2018-12-21 | 南京航空航天大学 | Reconfigurable network channel simulation method and device towards car networking communication |
| CN109067482B (en) * | 2018-08-14 | 2021-02-26 | 南京航空航天大学 | Reconfigurable network channel simulation method and device for Internet of vehicles communication |
| CN112655154A (en) * | 2018-08-30 | 2021-04-13 | 华为技术加拿大有限公司 | Method and system for linear signal processing using signal decomposition |
| CN112655154B (en) * | 2018-08-30 | 2022-08-26 | 华为技术加拿大有限公司 | Method and system for linear signal processing using signal decomposition |
| US11601157B2 (en) | 2018-08-30 | 2023-03-07 | Huawei Technologies Canada Co., Ltd. | Method and system for linear signal processing with signal decomposition |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106059685B (en) | 2019-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhu et al. | A novel 3D non-stationary wireless MIMO channel simulator and hardware emulator | |
| CN103532644B (en) | Multi-path shadow compound fading channel simulation device and work method thereof | |
| CN1702986B (en) | Baseband multi-path fading channel simulator | |
| CN111787558A (en) | MIMOME system and safety rate optimization method thereof | |
| CN109067482A (en) | Reconfigurable network channel simulation method and device towards car networking communication | |
| CN102130734B (en) | Method for modelling and simulating Nakagami fading channel | |
| CN109039508A (en) | Wireless multipath fading channel simulator system and method | |
| CN105306150A (en) | Device and method for channel simulation | |
| CN105187232A (en) | Dynamic channel simulation device and method | |
| CN111817806B (en) | Automobile air interface communication performance test method and system | |
| CN106059685B (en) | The extensive mimo channel simulator of time evolution and its analogy method | |
| CN105099588A (en) | Device and method for simulating interference propagation channel of aviation communication system | |
| CN112865895A (en) | Dual non-stationary channel model of Internet of vehicles | |
| Gherm et al. | HF propagation in a wideband ionospheric fluctuating reflection channel: Physically based software simulator of the channel | |
| CN101951297A (en) | Measurement model for fading channel of broadband wireless communication system and realization method | |
| CN102122996B (en) | Analog radio frequency signal generation method and system | |
| CN201018512Y (en) | Radio communication analog simulation system | |
| CN102752250B (en) | Intermediate frequency processing method and device for uplink and downlink, and radio remote unit (RRU) | |
| Picol et al. | Hardware simulator for MIMO radio channels: Design and features of the digital block | |
| CN106533593A (en) | Synchronous-random-access-memory-based dynamic multi-path time delay simulation apparatus and method | |
| CN101304275B (en) | Wireless channel modeling and simulation method for mobile communication network | |
| CN105577352B (en) | The long-range phase synchronization method of digitlization and system realized based on phase field | |
| CN103888196A (en) | Two-way mobile communication environment simulation system | |
| CN103684640B (en) | A kind of emulation mode of large-scale complex UWB channel | |
| CN103906252A (en) | Virtual multiple input and multiple output user pairing method and system and base station |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190226 |