Time-Varying Communication Channel High Altitude Platform Station Link Budget and Channel Modeling

: Because of the high BER (Bit Error Rate), low time delay and low channel transmission efﬁciency of HAPS (High Altitude Platform Station) in the near space. The link budget of HAPS and channel model are proposed in this paper. According to the channel characteristic, the channel model is set up, combined with different CNR (Carrier Noise Ratio), elevation angle, coding situations of wireless communication link by using Hamming code, PSK (Pulse Shift Keying) and Golay code respectively, then the situations of link quality and BER are analyzed. The simulation results show that the established model of the link budget and channel are suitable for the theoretical analysis results. The elevation of the HAPS communication link is smaller while the BER is increasing. The case of channel in the coding is better than in the un-coded situation. When every bit power and thermal noise power spectral density is larger, the BER of the HAPS communication link is becoming smaller.


Introduction
As an airspace that humans have not systematically developed and utilized, near space is an important operational space for the future air and space weaponry system, and it is also a new "strategic place" for the great powers [1,2].Near space often refers to an airspace with a height of 20 to 100 km above the flying altitude and below the satellite's orbit.The high-altitude platform station (HAPS) is a platform that is stuck in a space of 20 to 50 km in height and is stationary relative to a specific location on the earth [3,4].
HAPS can be seen as a communication system between a terrestrial communication system and a satellite communication system.The aim is to develop the potential benefits of a high space between land and space, improve communication capacity and spectrum utilization, and reduce system equipment cost and complexity.HAPS is smaller than the path fading of terrestrial communication systems.And has a smaller delay than the satellite communication system [5,6].To improve observation and coverage, the construction of surveillance and communication systems on near space platform has received further attention.HAPS has a similar development history to satellite communications in the 1960s [7].The international telecommunication union (ITU) describes HAPS as "a new technology with a promising application prospect that can subvert the development of the telecommunications industry."It is the foundation of next-generation wireless communications and has the ability to take full advantage of wireless spectrum resources, the system has large user capacity, good communication quality, and low operational risk.And you can upgrade the communication load operation [8][9][10] at any time.
Arnau et al. [11] studied the inter-frequency forwarding of space satellite communications.The signal-to-noise ratio of the signal quality of the spatial communication link is quantitatively analyzed under the conditions of multipath fading and rain attenuation.Ammari et al. [12] analyzed the HATS IMT-2000 system.Mathematical modeling and simulation analysis of ultra-high frequency (49 GHz) link characteristics were carried out.Jiang et al. [13] analyzed the uplink capacity limitation of near space and the method of expanding the uplink and downlink capacity.Wei [14] analyzed the relationship between Adaptive Modulation and Coding (AMC) and Automatic Repeat Request (ARQ) and Quality of Service (QoS) in the time domain.The time domain parameters such as modulation mode and coding mode can be dynamically adjusted according to the channel state.Perlaza et al. [15] combines spatial channel transmission characteristics.The AMC and the joint physical layer are studied based on the cross-layer optimization design of spatial information based on link adaptation.ARQ technology at the data link layer [16][17][18].
The authors of [19][20][21][22][23], for high-altitude platform CDMA communication systems, established a wireless link model in quasi-stationary state.The main impact of the swing on the communication system is analyzed.Literatures [24][25][26] studied the effect of horizontal movement of the platform on HAPS.The probability of handover between terminals in the cell caused by the horizontal movement of the platform is analyzed.The fixed beam-width coverage cell model is given in [27][28][29][30].And analyzes the impact of platform rotation on the probability of switching between terminal cells.The rotation has the greatest influence on the outer honeycomb.The authors of [31][32][33] studied the state of platform motion (horizontal and vertical motion), changes in the ground coverage area.
Alotaibi et al. [34] proposed a dynamic downlink transmission power of femtocells.Each femtocell adjusts its transmission power autonomously based on measured cost function unit.The transmission power level of femtocell is constrained by the rate of interference that femtocell produced to adjacent femtocells.Eirini et al. [35] addressed the problem of efficient power allocation in the uplink of two-tier closed-access femtocell.
Networks and the operational effectiveness of the proposed approach is evaluated through modeling and simulation.Sheikhzadeh et al. [36] considered a downlink resource allocation in an Orthogonal Frequency Division Multiple Access (OFDMA) based heterogonous cellular network and proposed Non-Exclusive Spectrum Allocation scheme (NESA).Tsiropoulou et al. [37] tackled the problem of joint users' uplink transmission power and data rate allocation in multi-service two-tier femtocell networks.Khamidehi et al. [38] investigated the optimal power allocation problem for a femtocell network where the utilized air interface is single carrier frequency division multiple access (SC-FDMA) and derived the optimal power allocation when we protect the data rate of macro users by imposing a temperature limit on the interference arising from femto users to the macro base station.Tsiropoulou et al. [39] proposed a novel utility-based game theoretic frame work towards simultaneously allocating users' uplink transmission power and rate in multi-service two-tier open-access femtocell networks.These literatures did not consider the Hamming, PSK and Golay code mode.This paper introduced Hamming, PSK and Golay code into the HAPS wireless communication channel.This paper mainly combines spatial communication links under different carrier-to-noise ratio, elevation angle, and coding.The HAPS channel model of its near space is established and simulated.
The remainder of this paper is organized as follows: In Section 2 we propose the system model.HAPS link budget and channel modeling are researched in Section 3. We discuss our findings by simulation analysis in Section 4, and conclude in Section 5.

System Model
The HAPS-based air-ground integrated aeronautical telecommunication network architecture is shown in Figure 1.The HAPS network includes a set of HAPS communication nodes that can perform routing, forwarding, and traffic management, and microwave communication can be implemented between each HAPS node.The link problem between HAPS is a technical problem, and the geometric relationship of the communication system is shown in Figure 1.Use a stable communication platform as a microwave relay station in the near-Earth space.a communication system with ground control equipment, entry equipment, and a variety of wireless users.The high-altitude platform can be integrated with the satellite ground or separately from the ground.Consider the communication platform as a synchronous communication satellite, that is, the communication platform is synchronized with the Earth's rotation.Can stay in the air for a long time.High-altitude platform communication utilizes good radio wave transmission characteristics, realizing communication connections between ground users, between platforms or between platforms and satellites through platforms.It has the advantages of flexible layout, wide application, low cost, safety and reliability.The schematic diagram of the HAPS space communication network is shown in Figure 2. Use a stable communication platform as a microwave relay station in the near-Earth space.a communication system with ground control equipment, entry equipment, and a variety of wireless users.The high-altitude platform can be integrated with the satellite ground or separately from the ground.Consider the communication platform as a synchronous communication satellite, that is, the communication platform is synchronized with the Earth's rotation.Can stay in the air for a long time.High-altitude platform communication utilizes good radio wave transmission characteristics, realizing communication connections between ground users, between platforms or between platforms and satellites through platforms.It has the advantages of flexible layout, wide application, low cost, safety and reliability.The schematic diagram of the HAPS space communication network is shown in Figure 2. Use a stable communication platform as a microwave relay station in the near-Earth space.a communication system with ground control equipment, entry equipment, and a variety of wireless users.The high-altitude platform can be integrated with the satellite ground or separately from the ground.Consider the communication platform as a synchronous communication satellite, that is, the communication platform is synchronized with the Earth's rotation.Can stay in the air for a long time.High-altitude platform communication utilizes good radio wave transmission characteristics, realizing communication connections between ground users, between platforms or between platforms and satellites through platforms.It has the advantages of flexible layout, wide application, low cost, safety and reliability.The schematic diagram of the HAPS space communication network is shown in Figure 2. It can be seen from Figure 2 that HAPS is a highly dynamic space platform.The communication between them is related to the atmospheric environment, multipath fading, space flicker and other factors.

HAPS Network
According to the geometric relationship of Figure 2, it can be calculated using the following equation: min min sin( 90) sin( 90) where E R is the radius of earth.

H h is the altitude of HAPS,
S h is the altitude from HAPS to earth ground.min


is the pitch angle of the HAPS,  is the central angle of the HAPS, which can be seen in Figure 2b.According to the geometric relationship, the central angle  can be obtained through the following Equation (2).
In Equation ( 2), the relevant coefficients have been explained in Equation ( 1).The central angle  of the HAPS can be calculated by the following equation: where SH is the distance between HAPS and point H .

E R is the radius of earth.
H h is the altitude of HAPS.The pitch angle SH  of the HAPS can be described by the Equation (4).

= arctan cos sin
As can be seen from Equation ( 4), the pitch angle where i is the -th i slot time,  denotes the time step.It can be seen from Figure 2 that HAPS is a highly dynamic space platform.The communication between them is related to the atmospheric environment, multipath fading, space flicker and other factors.
According to the geometric relationship of Figure 2, it can be calculated using the following equation: sin(90 where R E is the radius of earth.h H is the altitude of HAPS, h S is the altitude from HAPS to earth ground.ε min is the pitch angle of the HAPS, δ is the central angle of the HAPS, which can be seen in Figure 2b.According to the geometric relationship, the central angle δ can be obtained through the following Equation (2).
In Equation ( 2), the relevant coefficients have been explained in Equation (1).
The central angle β of the HAPS can be calculated by the following equation: where SH is the distance between HAPS and point H. R E is the radius of earth.h H is the altitude of HAPS.The pitch angle ε SH of the HAPS can be described by the Equation (4).
As can be seen from Equation ( 4), the pitch angle ε SH is related to the central angle β, radius of earth R E , HAPS altitude h H and h S (the distance from HAPS to earth ground), it has a nonlinear function relationship.x(t) is the distance at time t.
where i is the i-th slot time, ∆ denotes the time step.

HAPS Link Budget
The main goal of the link budget is to ensure the availability of the HAPS communication link.Considering the cost of space segments and ground station equipment, the HAPS link needs to be carefully designed.To fully optimize and save all available resources.The impairment of the HAPS link performance evaluation can be represented by the link diagram shown in Figure 3.
Information 2018, 10, x FOR PEER REVIEW 5 of 13

HAPS Link Budget
The main goal of the link budget is to ensure the availability of the HAPS communication link.Considering the cost of space segments and ground station equipment, the HAPS link needs to be carefully designed.To fully optimize and save all available resources.The impairment of the HAPS link performance evaluation can be represented by the link diagram shown in Figure 3. N is the thermal noise power, k is Boltzmann constant.B is IF (Intermediate Frequency) effective bandwidth.As can be seen from Figure 4, for different users ,, i j k , etc., there is interference between them.The smaller the elevation angle is, the longer the distance from the HAPS to its users.
The energy per bit b E of the near space HAPS can be expressed as: where C is the carrier wave average power.
b T indicates the time it takes to accept 1 bit(s).The signal to noise ratio can be expressed as: Where, P EIRP indicates the equivalent isotropic radiated power of the transmitter, AWGN is Additive white Gaussian noise.l FS is the loss of free space, G AT indicates the gain of acceptance antenna.s i (t) expresses the input signal and s o (t) denotes the output signal.N is the thermal noise power, k is Boltzmann constant.B is IF (Intermediate Frequency) effective bandwidth.T AR indicates antenna noise temperature, T e indicates the effective input noise temperature of the overall acceptance system, P TX is the transmit power, G A denotes the amplifier gain.The uplink interference of the near space HAPS is shown in Figure 4.
Information 2018, 10, x FOR PEER REVIEW 5 of 13

HAPS Link Budget
The main goal of the link budget is to ensure the availability of the HAPS communication link.Considering the cost of space segments and ground station equipment, the HAPS link needs to be carefully designed.To fully optimize and save all available resources.The impairment of the HAPS link performance evaluation can be represented by the link diagram shown in Figure 3. N is the thermal noise power, k is Boltzmann constant.B is IF (Intermediate Frequency) effective bandwidth.As can be seen from Figure 4, for different users ,, i j k , etc., there is interference between them.The smaller the elevation angle is, the longer the distance from the HAPS to its users.
The energy per bit b E of the near space HAPS can be expressed as: where C is the carrier wave average power.
b T indicates the time it takes to accept 1 bit(s).The signal to noise ratio can be expressed as: As can be seen from Figure 4, for different users i, j, k, etc., there is interference between them.The smaller the elevation angle is, the longer the distance from the HAPS to its users.
The energy per bit E b of the near space HAPS can be expressed as: where C is the carrier wave average power.T b indicates the time it takes to accept 1 bit(s).The signal to noise ratio can be expressed as: where R = 1/T b is the bit rate (b/s).T b indicates the spent time when transmit 1-bit energy.N 0 denotes the thermal noise power spectral density, then E b /N 0 can be expressed as: where P EIRP represents the equivalent isotropic radiated power of the transmitter.l FS denotes the loss of free space.g AR indicates the antenna gain.k is the Boltzmann constant.T s is the effective temperature in HAPS system.R E is the radius of earth.
As for 2PSK or QPSK, bit error rate P e and E b /N 0 have the following relationship: where er f c denotes the error function.For example, er f c( x e −t 2 dt.When the bit error rate P e = 10 −4 is required, the normalized ideal threshold SNR can be described by the following equation: where k is the Boltzmann constant.P b is the power when receive 1 bit.
The elevation angle of HAPS can be expressed as where δ is the central angle of the HAPS.φ 2 , φ 1 are the elevation angle of two different stations to the HAPS platform, respectively.In the HAPS wireless communication space, the change relationship of the elevation angle φ should be considered.The error rate P e in the un-coded case can be expressed as: If PSK (8 phase shift keying) and Hamming code (15,11) and Golay code (23,12) are used to encode near space HAPS channels, the bit error rate P e can be expressed as: where P e is the bit error rate, the value of P e is different based on different code mode.φ is the elevation angle of HAPS.

Channel Modeling
For the given time slot (k, l), the corresponding sub-channel h r,q [k, l] can be arranged in a where w[k, l] is M R dimension noise vector for every time slot (k, l).s[k, l] is the input vector with M T dimension; r[k, l] is the output vector with M R dimension.
Limited average power and peak power the upper and lower boundaries of the MIMO channel capacity of Equation ( 15) can be expressed as: where T denotes the entire transmission time period.E{•} is the mean value, P expresses the transmission power.K indicates the maximum number of slot time.s[k, l] is the input vector: where ζ is the adjustment factor.L is the maximum value of the slot time with M T dimension.Under the assumption that the input signal satisfies the average power limited Equation ( 16) and the peak power limited Equation ( 17), the non-coherent capacity of channel Equation ( 17) has an upper boundary of C ≤ U mimo 1 , and wireless communication channel capacity U mimo 1 can be denoted by the following equation: where B represents the channel bandwidth.sup indicates the minimum upper bound of one set.For example, sup{x} is the minimum upper bound of the set x. ℘ denotes the jitter factor of HAPS platform.℘ max expresses the jitter maximum value.Without loss of generality, ψ(B) can be expressed as: where C denotes the wireless communication signal power.N is the noise power.
If related parameters are obtained, the value of ψ r (B) could be calculated.

Analysis in the Femtocell
In the process of practical application, clustering has also proved to be an effective method to reduce interference in the femtocell.This method is widely used in ad hoc wireless networks.However, if this method is applied to a femtocell network, it is first necessary to know the number of clusters.The shared registered spectrum band is then divided into m sub-bands, each of which uses a sub-band.Since there are different numbers of femtocells in each cluster, this means different rate requirements.Therefore, the bandwidth allocation of each sub-band will directly affect the performance of the femtocell network.If the number of clusters is large, the bandwidth allocated by each cluster is relatively small.This will not meet the minimum requirements of femtocells.However, when the bandwidth allocated by each cluster is relatively large, the interference between the femtocells cannot be effectively reduced.Therefore, a new method is needed to determine the number of clusters and the bandwidth allocated for each cluster.This will be the next step in this paper.

Impact on the Resource Management in HAPS Wireless Communication System
The goal of wireless communication resource management is to provide service quality assurance for wireless user terminals in the network under limited bandwidth conditions.The basic starting point is that the network traffic is unevenly distributed, and the channel characteristics are fluctuating due to channel weakness and interference.It is necessary for flexible allocation and dynamic adjustment of the available resources of the wireless transmission, maximize wireless spectrum utilization and prevent network congestion and keep the signaling load as small as possible.
One of the main features of the HAPS communication system is the existence of a large number of non-real-time packet data services.Because different users have different rates.The sum of all user rates in a base station tends to exceed the channel capacity that can be transmitted by the base station's own frequency band.Therefore, it is necessary to have a scheduler in the base station according to the user QoS requirements.The channel resources will be allocated to different users according to the type of service.

Discussion
The noise spectral density is normalized by P(1W/Hz) = 3.5 × 10 6 s −1 .Transmit power is 0.5 Mw.The level of receiver thermal noise is 200 dBm/Hz.Free space path loss is 10 m.The receiver noise coefficient is 15 dB, β is set by 1.
The relationship between the transmit power of the HAPS and the antenna gain at different carrier-to-noise ratios is shown in Figure 5. can be seen from Figure 5, as the carrier-to-noise ratio C/N increases, the transmit power of the HAPS also increases proportionally.At the same time, the antenna gain is also increased with C/N.The duration simulation of the HAPS wireless communication link at different altitudes is shown in Figure 6.As can be seen from Figure 6, the higher the HAPS space base station, the longer its average link duration.Here the HAPS height is increased from 20 km to 30 km.But the average link duration does not increase linearly, as  increases, the average link duration also increases with a non-linear growth trend.It can be seen from Figure 7 that the higher the HAPS space base station, the higher its elevation angle, and as  increases, the elevation angle also increases accordingly.Under the same  , the higher the HAPS, the higher the elevation angle.
Whether the channel is encoded or not has a large influence on the bit error rate (BER), and its bit error rate e P at different  As can be seen from Figure 6, the higher the HAPS space base station, the longer its average link Here the HAPS height is increased from 20 km to 30 km.But the average link duration does not increase linearly, as ∆ increases, the average link duration also increases with a non-linear growth trend.
It can be seen from Figure 7 that the higher the HAPS space base station, the higher its elevation angle, and as ∆ increases, the elevation angle also increases accordingly.Under the same ∆, the higher the HAPS, the higher the elevation angle.As can be seen from Figure 6, the higher the HAPS space base station, the longer its average link duration.Here the HAPS height is increased from 20 km to 30 km.But the average link duration does not increase linearly, as  increases, the average link duration also increases with a non-linear growth trend.It can be seen from Figure 7 that the higher the HAPS space base station, the higher its elevation angle, and as  increases, the elevation angle also increases accordingly.Under the same  , the higher the HAPS, the higher the elevation angle.
Whether the channel is encoded or not has a large influence on the bit error rate (BER), and its bit error rate e P at different  Whether the channel is encoded or not has a large influence on the bit error rate (BER), and its bit error rate P e at different E b /N 0 can be shown in Figure 8.It can be seen from Figure 8 that when the value of 0 / b EN is small, the un-coded performance is better than the coded (8 PSK) channel case; when 0 / b EN is greater than 3.5, the coded channel has a bit error rate BER smaller than that of the un-coded channel.Overall, the bit error rate e P in the coding case is lower than the un-coded.
The influence of the elevation angle between the HAPS platform and the user in the near space on the transmission performance of the HAPS cannot be ignored.This is because if the elevation angle is small, the inclination path from the platform to the user will be large.This will increase free space loss and rain attenuation.Therefore, it will affect the HAPS link communication performance.Figure 9 simulates the error rate at different elevation angles  .As can be seen from Figure 9, in the case of the same 0 / b EN value, the larger the elevation angle  , the smaller the bit error rate.When the elevation angle C is 20°, PSK (8 phase shift keying) is used respectively, coding with the Golay code (23,12) and the Hamming code (15,11), the simulation results are shown in Figure 10.It can be seen from Figure 8 that when the value of E b /N 0 is small, the un-coded performance is better than the coded (8 PSK) channel case; when E b /N 0 is greater than 3.5, the coded channel has a bit error rate BER smaller than that of the un-coded channel.Overall, the bit error rate P e in the coding case is lower than the un-coded.
The influence of the elevation angle between the HAPS platform and the user in the near space on the transmission performance of the HAPS cannot be ignored.This is because if the elevation angle is small, the inclination path from the platform to the user will be large.This will increase free space loss and rain attenuation.Therefore, it will affect the HAPS link communication performance.Figure 9 simulates the error rate at different elevation angles φ.It can be seen from Figure 8 that when the value of 0 / b EN is small, the un-coded performance is better than the coded (8 PSK) channel case; when 0 / b EN is greater than 3.5, the coded channel has a bit error rate BER smaller than that of the un-coded channel.Overall, the bit error rate e P in the coding case is lower than the un-coded.
The influence of the elevation angle between the HAPS platform and the user in the near space on the transmission performance of the HAPS cannot be ignored.This is because if the elevation angle is small, the inclination path from the platform to the user will be large.This will increase free space loss and rain attenuation.Therefore, it will affect the HAPS link communication performance.Figure 9 simulates the error rate at different elevation angles  .As can be seen from Figure 9, in the case of the same 0 / b EN value, the larger the elevation angle  , the smaller the bit error rate.When the elevation angle C is 20°, PSK (8 phase shift keying) is used respectively, coding with the Golay code (23,12) and the Hamming code (15,11), the simulation results are shown in Figure 10.As can be seen from Figure 9, in the case of the same E b /N 0 value, the larger the elevation angle φ, the smaller the bit error rate.When the elevation angle C is 20 • , PSK (8 phase shift keying) is used respectively, coding with the Golay code (23,12) and the Hamming code (15,11), the simulation results are shown in Figure 10.

Conclusions
The characteristics of the near space HAPS are analyzed under the fast time-varying channel.The HAPS link budget mathematical model and wireless communication channel model are proposed respectively.The bit error rate under different HAPS elevation angles, signal to noise ratios are considered in this paper.The simulation results show that the greater the signal-to-noise ratio, the greater the gain and transmission power of the HAPS link.The bit error rate of the link channel in the coding case is better than that of the un-coded one.The key technical issue to be studied in the next step is to describe the time-varying network in terms of duration.Multi-source multi-suspension capacity model of HAPS link and super-radius overlapping coverage scene of near space interference domain capacity model.It can be seen from Figure 10 that when E b /N 0 ∈ (0, 4), the bit error rate P e of the PSK is lower than the Golay code and the Hamming code; when E b /N 0 ∈ (4, 10), the bit error rate of the Golay code is lower than the PSK and Hamming codes.

Conclusions
The characteristics of the near space HAPS are analyzed under the fast time-varying channel.The HAPS link budget mathematical model and wireless communication channel model are proposed respectively.The bit error rate under different HAPS elevation angles, signal to noise ratios are considered in this paper.The simulation results show that the greater the signal-to-noise ratio, the greater the gain and transmission power of the HAPS link.The bit error rate of the link channel in the coding case is better than that of the un-coded one.The key technical issue to be studied in the next step is to describe the time-varying network in terms of duration.Multi-source multi-suspension capacity model of HAPS link and super-radius overlapping coverage scene of near space interference domain capacity model.

Figure 1 .
Figure 1.System frame diagram of the High Altitude Platform Station (HAPS).

Figure 1 .
Figure 1.System frame diagram of the High Altitude Platform Station (HAPS).

Figure 1 .
Figure 1.System frame diagram of the High Altitude Platform Station (HAPS).

SH
is related to the central angle  , radius of earth E R , HAPS altitude H h and S h (the distance from HAPS to earth ground), it has a nonlinear function relationship.() xt is the distance at time t .

Figure 2 .
Figure 2. Space network diagram and link geometry of HAPS.(a) Network diagram; (b) Geometry diagram.

Figure 3 .
Figure 3. Sketch map of HAPS link budget.

Figure 3 .
Figure 3. Sketch map of HAPS link budget.

Figure 3 .
Figure 3. Sketch map of HAPS link budget.

Figure 5 .
Figure 5.The transmitted power and antenna gain under different carrier noise ratio.

Figure 7 .
Figure 7. Relationship between elevation angles at different HAPS altitudes.

EN
can be shown in Figure8.

Figure 7 .
Figure 7. Relationship between elevation angles at different HAPS altitudes.

EN
can be shown in Figure8.

Figure 7 .
Figure 7. Relationship between elevation angles at different HAPS altitudes.

Figure 8 .
Figure 8.The BER (Bit Error Rate) under different encoding.

Figure 9 .
Figure 9.The BER (Bit Error Rate) under different elevation angle.

Figure 8 .
Figure 8.The BER (Bit Error Rate) under different encoding.

Figure 8 .
Figure 8.The BER (Bit Error Rate) under different encoding.

Figure 9 .
Figure 9.The BER (Bit Error Rate) under different elevation angle.

Figure 9 .
Figure 9.The BER (Bit Error Rate) under different elevation angle.

Author Contributions:
The authors performed the experiments and analyzed the results together.Introduction, methodology, cosmetic detection algorithms and proposed model were written by X.L. and H.L.; while data collection, experimental results and conclusion sections were written by C.L. and Y.L. Funding: This research was funded by National Planning Office of Philosophy and Social Science of China (No. 17XXW004), Young Fund Project of Humanities and Social Sciences Research of Ministry of Education of China (No. 16YJC860010), Social Science of Humanity of Chongqing Municipal Education Commission (No. 17SKG144), Science and Technology Research Program of Chongqing Municipal Education Commission (No. KJ1600923, No. KJ17092060, No. KJ1600928), 2018 Chongqing Science and Technology Commission Technology Innovation and Application Demonstration Project (No. cstc2018jscx-msybX0049). Natural Science Foundation of China (No. 61503052, No. 61571069, No. 61501065, No. 61502064).Open Fund Project of Chongqing Technology and Business University, Research Center of Chongqing University Network Public Opinion and Ideological Dynamic (No. KFJJ2017024).China Postdoctoral Science Foundation (No. 2017M612911), Research Foundation of the Natural Foundation of Chongqing City (No. cstc2016jcyjA0076).The author Xiaoyang Liu also thanks for the financial support from CSC (China Scholarship Council, No. 201608505142).