CN114614927A - Nakagami fading complex channel simulation method - Google Patents

Nakagami fading complex channel simulation method Download PDF

Info

Publication number
CN114614927A
CN114614927A CN202210274362.4A CN202210274362A CN114614927A CN 114614927 A CN114614927 A CN 114614927A CN 202210274362 A CN202210274362 A CN 202210274362A CN 114614927 A CN114614927 A CN 114614927A
Authority
CN
China
Prior art keywords
distribution
nakagami
complex channel
envelope
simulation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210274362.4A
Other languages
Chinese (zh)
Inventor
孙成刚
岳红霞
张剑锋
周武林
唐庆生
吴翠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu Zhongxiangtiandi Network Technology Co ltd
Original Assignee
Chengdu Zhongxiangtiandi Network Technology Co ltd
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 Chengdu Zhongxiangtiandi Network Technology Co ltd filed Critical Chengdu Zhongxiangtiandi Network Technology Co ltd
Priority to CN202210274362.4A priority Critical patent/CN114614927A/en
Publication of CN114614927A publication Critical patent/CN114614927A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3911Fading models or fading generators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3912Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

The invention discloses a Nakagami fading complex channel simulation method, which comprises the following steps: s1: obtaining a time domain signal sequence of Gamma distribution with parameters of (m, omega/m) by calling a Gamma function, and squaring the sequence element by element to obtain an envelope r so that the envelope r follows the Nakagami-m distribution of the parameters of (m, omega); s2: the corresponding phase phi is obtained by a rejection method to generate uniform distribution
Figure 100004_DEST_PATH_IMAGE001
Generates a random number x obeying uniform distribution
Figure 887977DEST_PATH_IMAGE002
Random number of
Figure 100004_DEST_PATH_IMAGE003
B is
Figure 725483DEST_PATH_IMAGE004
Is at a maximum value of
Figure 100004_DEST_PATH_IMAGE005
Then, then
Figure 465163DEST_PATH_IMAGE006
Otherwise, repeating the step; s3: combining envelopes r and
Figure 100004_DEST_PATH_IMAGE007
phase, establishing a corresponding complex channel of

Description

Nakagami fading complex channel simulation method
Technical Field
The invention relates to a simulation method, in particular to a Nakagami fading complex channel simulation method.
Background
The wireless fading channel modeling is to establish a mathematical model based on channel characteristic analysis and perform simulation similar to an actual channel in a laboratory environment. Compared with the traditional field actual measurement, the method can greatly reduce the difficulty and the cost of system test and is widely applied.
When the radio signal encounters an undulating terrain during its propagation, it causes absorption and penetration of energy and reflection, scattering and diffraction of radio waves, and the signal arriving at the mobile station antenna is not a single path but a composite of many reflected waves from many paths. Since the distance of the radio wave passing through each path is different, the arrival time of the reflected wave from each path is different, and the phase is different. Multiple signals of different phases are superimposed at the receiving end, sometimes with in-phase superimposition to enhance them, and sometimes with reverse superimposition to attenuate them. Therefore, the amplitude of the received signal will change sharply, i.e., fading is generated. Research has shown that the envelope of a wireless fading received signal generally follows a rayleigh distribution, a rice distribution, or a Nakagami-m distribution. The Nakagami-m distribution (Nakagami distribution for short) is more general, different fading conditions such as severe, moderate, slight and no fading can be described by changing the value of the fading factor m, and the actual measurement result shows that the distribution is closer to the actual situation than the distribution such as Rayleigh and Rice.
At present, the modeling and simulation aiming at the Nakagami fading channel which are common at home and abroad can be divided into three categories: a channel decomposition method, an inverse transform method, and a discard method. Although the channel decomposition method is simpler, the performance is poorer when m is not an integer multiple of 0.5; because the cumulative probability density function of Nakagami distribution has no inverse function, the inverse transformation method can only adopt polynomial approximation, and a large number of numerical value approximations are needed for different m values to obtain a coefficient table, which is very complex; the rejection method is a general method suitable for arbitrary fading distribution, and the difficulty is to find an efficient hat function. According to related reports, the efficiency of the hat function which is common at present is different from 60% to 65%.
Disclosure of Invention
The technical problem to be solved by the invention is that the verification result is simply compared with a theoretical value and a simulation value, and a correlation method for scientifically verifying the validity and the accuracy of the model is not mentioned, and the invention aims to provide a Nakagami fading complex channel simulation method and solve the problems.
The invention is realized by the following technical scheme:
a method for simulating a Nakagami fading complex channel, the method comprising the steps of:
s1: obtaining a time domain signal sequence of Gamma distribution with parameters of (m, omega/m) by calling a Gamma function, and squaring the sequence element by element to obtain an envelope r so that the envelope r follows the Nakagami-m distribution of the parameters of (m, omega);
s2: the corresponding phase phi is obtained by a rejection method to generate uniform distribution
Figure 702793DEST_PATH_IMAGE001
Generates a random number x obeying uniform distribution
Figure DEST_PATH_IMAGE002
Random number of
Figure 377488DEST_PATH_IMAGE003
B is
Figure DEST_PATH_IMAGE004
Is at a maximum value of
Figure 838556DEST_PATH_IMAGE005
Then, then
Figure DEST_PATH_IMAGE006
Otherwise, repeating the step;
s3: binding envelope
Figure 704619DEST_PATH_IMAGE007
And
Figure DEST_PATH_IMAGE008
phase, establishing a corresponding complex channel of
Figure 883927DEST_PATH_IMAGE009
In wireless communication, a transmission signal is often influenced by occlusion, absorption, reflection, refraction and diffraction caused by various objects in the environment during the propagation process, and multiple path signal components are formed to reach a receiver. The signal components of different paths have different propagation delays, phases and amplitudes and are added with channel noise, and their superposition causes the composite signals to cancel or enhance each other, resulting in severe fading. Such fading can reduce the available power of the desired signal and increase the effects of interference, resulting in distortion, broadening, overlapping and distorting the received signal at the receiver, and even causing substantial errors in the demodulator output of the communication system, resulting in complete failure to communicate. Therefore, in order to improve the quality of wireless communication, intensive research and analysis on the transmission characteristics of wireless channels are required.
Wireless channel modeling refers to establishing a mathematical model to analyze the characteristics of a channel and simulate the actual conditions of the channel. Compared with the traditional field actual measurement, the method can be widely applied due to the fact that the test workload and difficulty can be greatly reduced, and the cost can be saved.
For the simulation of the flat fading channel, a Clarke model based on scattering waves is often adopted, and the model indicates that in a channel environment with uniform scattering and no direct path, the signal envelope conforms to Rayleigh distribution; in a channel environment with a direct path, the signal envelope conforms to a rice distribution. However, the description of the model on the fast fading in the long-distance channel is quite rough, the phenomenon is firstly observed by Nakagami, and the density function based on the variable parameter gamma distribution is established to fit the obtained experimental data to obtain the approximate distribution. For the simulation method of the Nakagami channel, Brute force method, inverse transformation method and truncation method are commonly used at present.
Further, b in the step S2 is
Figure DEST_PATH_IMAGE010
Is at a maximum value of
Figure 362313DEST_PATH_IMAGE011
Maximum value of
Figure DEST_PATH_IMAGE012
(ii) a If it is
Figure 943467DEST_PATH_IMAGE013
Maximum value
Figure DEST_PATH_IMAGE014
Further, in the actual simulation, at
Figure 983361DEST_PATH_IMAGE015
Approach to 0 time
Figure DEST_PATH_IMAGE016
Get it
Figure 649966DEST_PATH_IMAGE017
And (4) inputting simulation values to realize actual simulation.
Further, the air conditioner is provided with a fan,
Figure DEST_PATH_IMAGE018
taking a value within a function, wherein
Figure 666463DEST_PATH_IMAGE019
(ii) a The requirements are met,
Figure DEST_PATH_IMAGE020
i.e. gamma distribution.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the invention discloses a Nakagami fading complex channel simulation method, which generates a Nakagami-m complex channel by calling a Gamma function to obtain a parameter, quickly generates a random number obeying Nakagami-m distribution by using a simple algorithm, thereby knowing that phase information of the Nakagami-m distribution is related to a fading parameter m value, and then according to an envelope
Figure 102124DEST_PATH_IMAGE021
And
Figure 608191DEST_PATH_IMAGE015
phase, establishing corresponding complex channel.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not used as limitations of the present invention.
Examples
The invention relates to a Nakagami fading complex channel simulation method, which comprises the following steps:
s1: obtaining a time domain signal sequence of Gamma distribution with the parameter of (m, omega/m) by calling a Gamma function, and squaring the sequence element by element to obtain an envelope
Figure 27671DEST_PATH_IMAGE021
So as to envelop
Figure 877553DEST_PATH_IMAGE021
A Nakagami-m distribution obeying the parameters (m, Ω);
s2: corresponding phase
Figure 698879DEST_PATH_IMAGE015
By obtaining with a discardment method, a uniform distribution of obedience is generated
Figure DEST_PATH_IMAGE022
Generates a random number x obeying uniform distribution
Figure 579110DEST_PATH_IMAGE023
Random number of
Figure DEST_PATH_IMAGE024
B is the maximum value of
Figure 220307DEST_PATH_IMAGE025
If, if
Figure DEST_PATH_IMAGE026
Then, then
Figure 578607DEST_PATH_IMAGE027
Otherwise, repeating the step;
s3: binding envelope
Figure 785598DEST_PATH_IMAGE021
And with
Figure 869354DEST_PATH_IMAGE015
Phase, establishing a corresponding complex channel of
Figure DEST_PATH_IMAGE028
In the step S2, b is
Figure 732267DEST_PATH_IMAGE029
Is at a maximum value of
Figure DEST_PATH_IMAGE030
Maximum value
Figure 894258DEST_PATH_IMAGE031
(ii) a If it is
Figure DEST_PATH_IMAGE032
Maximum value
Figure 893439DEST_PATH_IMAGE033
In the actual simulation, at
Figure 646631DEST_PATH_IMAGE015
Approach to 0 time
Figure DEST_PATH_IMAGE034
Get it
Figure 495376DEST_PATH_IMAGE035
And (4) inputting simulation values to realize actual simulation.
Figure 992216DEST_PATH_IMAGE025
Taking a value within a function, wherein
Figure DEST_PATH_IMAGE036
(ii) a The requirements are met,
Figure 111482DEST_PATH_IMAGE037
i.e. gamma distribution.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (4)

1. A Nakagami fading complex channel simulation method is characterized by comprising the following steps:
s1: obtaining a time domain signal sequence of Gamma distribution with the parameter of (m, omega/m) by calling a Gamma function, and squaring the sequence element by element to obtain an envelope r so that the envelope r follows the Nakagami-m distribution of the parameter (m, omega);
s2: the corresponding phase phi is obtained by a rejection method to generate uniform distribution
Figure DEST_PATH_IMAGE001
Generates a random number x obeying uniform distribution
Figure 433865DEST_PATH_IMAGE002
Random number of
Figure DEST_PATH_IMAGE003
B is
Figure 288688DEST_PATH_IMAGE004
Is at a maximum value of
Figure DEST_PATH_IMAGE005
Then, then
Figure 6109DEST_PATH_IMAGE006
Otherwise, repeating the step;
s3: binding envelope
Figure DEST_PATH_IMAGE007
And
Figure 73422DEST_PATH_IMAGE008
phase, establishing a corresponding complex channel of
Figure DEST_PATH_IMAGE009
2. The method as claimed in claim 1, wherein b in step S2 is
Figure 324012DEST_PATH_IMAGE010
Is at a maximum value of
Figure DEST_PATH_IMAGE011
Maximum value
Figure 84158DEST_PATH_IMAGE012
(ii) a If it is
Figure DEST_PATH_IMAGE013
Maximum value
Figure 23295DEST_PATH_IMAGE014
3. The method as claimed in claim 2, wherein the simulation is performed in real time at
Figure DEST_PATH_IMAGE015
Approach to 0 time
Figure 159878DEST_PATH_IMAGE016
Taking out
Figure DEST_PATH_IMAGE017
And (4) inputting simulation values to realize actual simulation.
4. The method of claim 2, wherein the Nakagami fading complex channel simulation method,
Figure 267905DEST_PATH_IMAGE018
taking a value within a function, wherein
Figure DEST_PATH_IMAGE019
(ii) a Satisfy the requirement of
Figure 198952DEST_PATH_IMAGE020
I.e. gamma distribution.
CN202210274362.4A 2022-03-21 2022-03-21 Nakagami fading complex channel simulation method Pending CN114614927A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210274362.4A CN114614927A (en) 2022-03-21 2022-03-21 Nakagami fading complex channel simulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210274362.4A CN114614927A (en) 2022-03-21 2022-03-21 Nakagami fading complex channel simulation method

Publications (1)

Publication Number Publication Date
CN114614927A true CN114614927A (en) 2022-06-10

Family

ID=81864661

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210274362.4A Pending CN114614927A (en) 2022-03-21 2022-03-21 Nakagami fading complex channel simulation method

Country Status (1)

Country Link
CN (1) CN114614927A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102130734A (en) * 2011-04-22 2011-07-20 南京航空航天大学 Method for modelling and simulating Nakagami fading channel
US20120129466A1 (en) * 2010-11-23 2012-05-24 I Shou University Evaluation Device and Method for Providing a Transceiver System with Performance Information Thereof
CN104052557A (en) * 2014-06-24 2014-09-17 西安电子科技大学 Method for modeling Nakagami repeated fading channel
CN105846926A (en) * 2016-04-15 2016-08-10 西安电子科技大学 Time domain self-correlation Nakagami-m fading complex channel simulation method
CN110138481A (en) * 2019-04-30 2019-08-16 河海大学 A kind of Nakagami based on inverse transformation method answers random digit generation method
CN112511241A (en) * 2020-11-10 2021-03-16 河海大学 Composite fading channel random number generation method based on lognormal distribution approximation

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120129466A1 (en) * 2010-11-23 2012-05-24 I Shou University Evaluation Device and Method for Providing a Transceiver System with Performance Information Thereof
CN102130734A (en) * 2011-04-22 2011-07-20 南京航空航天大学 Method for modelling and simulating Nakagami fading channel
CN104052557A (en) * 2014-06-24 2014-09-17 西安电子科技大学 Method for modeling Nakagami repeated fading channel
CN105846926A (en) * 2016-04-15 2016-08-10 西安电子科技大学 Time domain self-correlation Nakagami-m fading complex channel simulation method
CN110138481A (en) * 2019-04-30 2019-08-16 河海大学 A kind of Nakagami based on inverse transformation method answers random digit generation method
CN112511241A (en) * 2020-11-10 2021-03-16 河海大学 Composite fading channel random number generation method based on lognormal distribution approximation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MICHEL DAOUD YACOUB ET AL.: ""Nakagami-m phase-envelope joint distribution: An improved model"", 《2009 SBMO/IEEE MTT-S INTERNATIONAL MICROWAVE AND OPTOELECTRONICS CONFERENCE》 *
周涛 等: ""Nakagami-m复衰落信道的仿真模型研究"", 《电子测量技术》 *

Similar Documents

Publication Publication Date Title
CN103532644B (en) Multi-path shadow compound fading channel simulation device and work method thereof
CN105049142A (en) Dual-path static baseband channel simulating device and method
US9288696B2 (en) Method and apparatus for wireless device performance testing
CN110515038B (en) Self-adaptive passive positioning device based on unmanned aerial vehicle-array and implementation method
CN109714120B (en) Method for simulating propagation characteristics of coupling multi-antenna indoor space fading channel
CN102130734B (en) Method for modelling and simulating Nakagami fading channel
Raju et al. Evaluation of BER for AWGN, Rayleigh fading channels under M-QAM modulation scheme
Fard et al. An FPGA-based simulator for high path count Rayleigh and Rician fading
Qureshi et al. SEP performance of triangular QAM with MRC spatial diversity over fading channels
CN102201876A (en) Rapid Nakagami channel model simulation method
US20190109611A1 (en) Mitigating passive intermodulation interference in a mimo wireless network
CN114614927A (en) Nakagami fading complex channel simulation method
CN107276657B (en) Millimeter wave mixed beam forming method based on parallel algorithm
Skima et al. A multi-criteria comparative analysis of different Rayleigh fading channel simulators
Al-Omary Investigating the visibility of 28 and 73 GHz frequency bands for outdoor MIMO channel
Ruiz-García et al. SDR-based channel emulator for vehicular communications
CN102263599B (en) Intelligent antenna array simulation method and apparatus thereof
Papazafeiropoulos et al. An extended generalized Rice model for wireless communications
Shuli et al. A real-time simulation design of multi-path fading channel based on SOS method
Shen et al. A simple simulation method for Nakagami fading channel
Komatsu et al. Development of tool for evaluation of automotive conformity of FM receivers using two-stage method
JP3500428B2 (en) Antenna directivity control device
CN101557250A (en) Intelligent antenna weight value generating method based on direction prediction and device thereof
CN114268398A (en) Simulation algorithm of suzuki channel
Emmer et al. Measurements of base station two-branch space and polarization diversity reception and a comparison of the diversity gain based on the CDF of signal level and simulations of BER in a GSM 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
RJ01 Rejection of invention patent application after publication

Application publication date: 20220610