• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Adaptive protograph-based BICM-ID relying on the RJ-MCMC algorithm: a reliable and efficient transmission solution for plasma sheath channels

    2022-05-05 01:49:06YanmingLIU劉彥明HailiangWEI魏海亮LeiSHI石磊andBoYAO姚博
    Plasma Science and Technology 2022年4期
    關(guān)鍵詞:石磊

    Yanming LIU (劉彥明), Hailiang WEI (魏海亮), Lei SHI (石磊),?and Bo YAO (姚博)

    1 Key Laboratory of Information and Structure Efficiency in Extreme Environment, Ministry of Education of China, Xi’an 710071, People’s Republic of China

    2 School of Aerospace Science and Technology, Xidian University, Xi’an 710071, People’s Republic of China

    Abstract For reentry communication, owing to the influence of the highly dynamic plasma sheath (PS),the parasitic modulation effect can occur and the received phase shift keying (PSK)signal constellation can be severely rotated, leading to unacceptable demodulation performance degradation.In this work, an adaptive non-coherent bit-interleaved coded modulation with iterative decoding(BICM-ID)system with binary PSK(BPSK)modulation and protograph lowdensity parity-check under the PS channel is proposed.The proposed protograph-based BICMID(P-BICM-ID)system can achieve joint processing of demodulation and decoding,where the soft information is adaptively estimated by reversible-jump Markov chain Monte Carlo (RJMCMC)algorithms.Simulation results indicate that compared to existing algorithms, the proposed system can adapt well to the dynamic characteristics of the PS channel and can obtain a 5 dB performance improvement at a bit error rate of 10-6.

    Keywords: plasma sheath, BICM-ID, Markov model, joint demodulation and decoding

    1.Introduction

    Due to the wide application prospects of hypersonic vehicles in civil and military fields, research on the communication technology of hypersonic aircraft has attracted widespread attention.However,when reentering or moving in the atmosphere at a speed exceeding Mach 8,the vehicle will be surrounded by a plasma layer (also called the plasma sheath)[1, 2], which can attenuate electromagnetic(EM)waves.The plasma sheath(PS)channel presents highly dynamic and non-stationary characteristics [3, 4] with parasitic modulation effects [5].These effects seriously affect the quality of reentry communication and can even cause radio blackouts [6-9].The problem of radio blackouts has not been completely solved yet.

    Many methods have been developed to mitigate communication blackouts by lowering the PS electron density [10-12],such as aerodynamic shaping [13], liquid quenchant injection[14], and electromagnetic windowing [15, 16].However, the practical shortcomings and engineering obstacles of these methods have not been satisfactorily resolved.A probably feasible solution,aiming at adapting to the high dynamics and large-scale fading characteristics of the plasma sheath channel, is to design adaptive communication methods to alleviate communication blackout,such as increasing transmitted power or communication frequency [17], employing low-frequency communication [18],channel prediction relied on the reflected waves [19], and usingan adaptive data rate control method based on the antenna voltage standing wave ratio(VSWR)[20].Regarding modulation and demodulation methods over the PS channel, [21] and [22]investigated the performance of traditional modulation methods over the PS channel from theoretical and experimental perspectives, respectively.Research has shown that non-coherence frequency shift keying (NC-FSK)is better than phase-shift keying(PSK)and amplitude shift keying(ASK)modulations.Owing to the parasitic modulation effect, PSK modulation over the PS channel presents a great challenge.To address this problem,Yang et al proposed a novel demodulation method based on spectral clustering for PSK signals over the PS channel[23].Liu et al developed a joint demodulation and estimation algorithm based on a Bayesian learning method for extracting curves on the signal constellation and identifying the transmitted symbols[24].

    Table 1.Typical values of Ne ( z , t).

    However, these existing methods rarely involve coding schemes, nor do they exploit joint processing techniques for signal detection and decoding.Due to the non-stationary nature of PS channels, the performance improvement of these methods is limited.In fact, in many communication scenarios, the joint processing technique for signal detection and decoding can yield joint gains in addition to coding gains,and the LDPC-BICM-ID system[25-28]has been proved to be very effective and popular.Bit-interleaved coded modulation and iterative decoding(BICMID), which practically utilize iterative information exchange between the demodulator and the channel decoder while obtaining knowledge of other bit values, can achieve highly reliable wireless communication.The application of an interleaver in BICM-ID can also address sudden changes in the channel.Therefore,the LDPC-BICM-ID system can adapt to deep fading and dynamic channels.Note that the PS channel is just a kind of deep fading and highly dynamic channel.PSK and quadrature amplitude modulation (QAM)are widely used in LDPC-BICMID systems.However,high-order PSK modulations with poor bit error rate(BER)performance are difficult to apply to PS channels.Therefore,to achieve highly reliable and efficient communication performance under the PS channel, the LDPC-BICM-ID system based on binary PSK (BPSK)may be a feasible solution.

    Low-density parity-check (LDPC)codes [29, 30] and protograph codes[31,32]have been extensively researched and combined with BICM-ID systems over many channels,such as fast-fading,block-fading,and satellite communication channels.However, designing an adaptive BICM-ID system over the PS channel remains a challenge.This is because PS channels possess a parasitic modulation effect and non-stationary characteristics, and the prior probability density function of the received signal is unknown.Thus, the soft information for maximum a posteriori (MAP)detector cannot be accurately calculated.

    Against this background,an adaptive modified protographbased BICM-ID (P-BICM-ID)framework with BPSK modulation for PS channels is proposed.The main contribution of this work is that an adaptive soft information estimation method based on reversible-jump Markov chain Monte Carlo (RJMCMC)for our P-BICM-ID system is developed to calculate the soft information.Considering the uncertainty of the prior channel information of PS channels,the PDF of the channel soft information for BICM-ID is adaptively estimated using the RJMCMC algorithm, and then the soft information is calculated and used to decode iteratively.

    The remainder of this paper is organized as follows.In section 2, the dynamic PS channel is modeled, and the parasitic modulation effect on the BPSK signal is analyzed.Section 3 proposes an adaptive P-BICM-ID system over PS channels.Then, numerical examples and analyses under different channel conditions are presented in section 4.Finally,this work is summarized in section 5.

    2.Dynamic PS channel modeling and its impact on BPSK signals

    2.1.Electron density modeling

    During reentry, electron density fluctuates violently and varies in the spatial and temporal domains.The spatiotemporal electron densityNe(z,t)[3] can be described as follows:

    where(z,t)is the average distribution of electron density from the surface of the vehicle outward, θ is the angle of attack(AOA),and A has a typical value of 0.1 per degree[3].δ(t)is the Gaussian jitter [33] with zero mean and standard deviation σmax.Suppose that θ follows a sinusoidal function with amplitude Aθand frequency fθ.In addition,presents large-scale changes in electron density and can be represented by [34]:

    whereNepeakis the peak value of the electron density,zpeakis the position ofNepeak,zmaxis the sheath thickness, c1and c2are parameters,and B is a proportion coefficient with a typical value of 1/14 per degree [35].The typical values ofNe(z,t)are listed in table 1.

    In our spatiotemporal electron density model, the largescale change of electron density by flight speed and height,mesoscale change by the angle of attack, and small-scale disturbance are considered.Therefore, the model is suitable for all deep fading, non-stationary, highly dynamic PS channels.Reference [36] also constructs a simple model, which only considers the large-scale and small-scale factors.Our model is more reasonable than this model.What is more,[37]derived the outage probability of NC-MFSK(non-coherent M-ary frequency shift keying)modulation based on our channel model.

    2.2.Time-varying plasma sheath channel

    To ensure the accuracy and efficiency of calculation,a layered PS model is used to approximate the real continuous model.The PS is modeled as M neighboring uniform thin PS slabs with the complex dielectric constant ε(zm, t)(m = 1, …, M),as shown in figure 1.The layer thickness is dm= zmax/M.The plasma sheath channel h(t)can be precisely calculated by using the transmission line analogy method [33].

    ε(zm, t)can be determined by

    wherefp(z m,t)is the plasma frequency of the mth layer,v is the electron collision frequency, ω = 2πfEM, and fEMis the carrier frequency.The propagation constantk(z m,t)can be described as

    where μ0is vacuum permeability.

    The wave impedances of the mth layer is

    The transmission line matrix can be calculated by

    Finally, the overall transmission matrix can be obtained as

    The PS channel h(t)can be calculated by

    whereZ(zM+1,t)donates the wave impedance of the(M + 1)th layer andZ(z0,t)presents the wave impedances of the vacuum layer.

    2.3.Impact of the dynamic PS channel on BPSK signals

    For BPSK modulation with the initial phase±π/2, the signal transmitted in one symbol period is given by

    where Esis the symbol energy and Tsis the symbol period.The received signal is expressed as

    where n(t)is assumed to be additive white circularly complex Gaussian noise with zero mean and variance σ2.The bit signal-to-noise-ratio (SNR)γbis defined as Es/σ2.

    To comprehensively understand the effect of PS channels on BPSK signals, two typical simulations are performed[5, 38].Simulation case 1:Nepeak= 4 × 1018m-3, σmax=0.05, Aθ= 4°, v = 1 GHz, fEM= 32 GHz, and γb= 20 dB;simulation case 2:Nepeak= 1 × 1019m-3, σmax= 0.1, Aθ= 0°, v = 1 GHz, fEM= 32 GHz, and γb= 20 dB.

    As shown in figure 2, there are two scattered sets representing bits 0 and 1 in the constellation diagram of the BPSK signals and the constellations rotate dynamically due to the random jitter of the electron density and AOA change,indicating that the dynamics and parasitic modulation effect of PS channels are severe.The statistical characteristics of the signal phases are shown in figure 3.Under case 1, the phase conditional probability density function(CPDF)for bit 0 or 1 presents multipeak features because the change in AOA makes the channel exhibit non-stationary characteristics.Under case 2,since the AOA is equal to 0°,each phase CPDF exhibits one dominating peak.However, two scatter sets representing bits 0 and 1 overlap in constellations whenNepeakincreases to 1 × 1019m-3.This is because the increase inNepeakleads to a more dynamic signal, thereby deteriorating the communication performance.Moreover, unexpected phase shifts occur in the CPDFs, and the phase centers move away from ±π/2 due to the parasitic modulation effect,which brings great difficulties for correct demodulation.

    From figures 2 and 3, the PS channel exhibits serious signal attenuation and phase shift effects.Thus,the prior PDF of the received signal is unknown,which makes the traditional soft detector of BICM-ID for PS channels no longer effective.Therefore, the design of a P-BICM-ID system over the PS channel is a great challenge.

    3.System model and adaptive scheme

    To overcome the above challenge, this work proposes a P-BICM-ID system in which the RJ-MCMC algorithm adaptively estimates soft information.At the receiver, the phase detector is used to calculate and correct the signal phase,and the soft detector is used to adaptively estimate the CPDF of the signal and calculate the soft information.The CPDF of the signal fluctuates randomly with the change in the non-stationary PS channel,and it can be expressed as a mixed Gaussian function.The RJ-MCMC algorithm can be used to adaptively estimate the state parameters of the multi-Gaussian process.The soft information is then calculated based on the estimated PDF of the PSK signal and is used to decode iteratively.

    3.1.System model

    The proposed adaptive P-BICM-ID system relying on the RJMCMC algorithm is illustrated in figure 4.

    The signal to be transmitted is denoted by s = [s1,s2,…,sL], a vector of L BPSK symbols.After the signal travels through the time-varying PS channel, r = [r1, r2, …, rL] is received at the receiver as shown in figure 4.At the transmitter, the information bit stream b is encoded by an LDPC encoder [28] producing the bitstream c, which is then fed to the random bit interleaver π1.The generated stream u is then sent to the unity-rate code (URC)encoder [32], where the produced coded bitstream c′ is interleaved by the interleaver π2, generating the bitstream u′.The URC encoder is used to avoid the BER floor [39].Then, the modulator maps the bits to symbols s.At time instant n, the BPSK signal sntravels through PS channel hn, and the received signal rnis obtained by the receiver.

    At the receiver, to achieve better joint demodulation and decoding performance, a phase detector, soft detector, URC decoder, and belief propagation (BP)decoder are employed,as shown in figure 4.The phase detector is used to accurately calculate and correct the phase φ of the received signal r.The URC decoder processes the logarithmic-likelihood ratios(LLRs)L(u′)produced by the soft detector, deinterleaved by π2,the generated L(c′)is passed to the iterative URC and BP decoders, and the a posteriori LLRs Lch,o(c)is computed.Iterative detection is conducted by exchanging extrinsic information between the outer BP decoder and the inner URC decoder.Finally, the estimated bitstreamb?is judged and output.

    For LDPC codes, the protograph LDPC codes retain the advantages (approximate to the capacity limit)of traditional LDPC codes and have the advantage of simple structures.In this study, we adopt the promising protograph LDPC code family from the 3GPP TS 38.212 protocol [40], which has a low decoding threshold and a minimum distance that grows linearly with increasing block-length.The code rate Rc=L/N, where L is the information bit number per codeword of length N.The transmitted SNR is defined as γ = γb+ 10log10(Rc).

    3.2.Phase detector

    Owing to the serious parasitic modulation effect, unexpectedly high phase shifts occur in CPDFs, which brings great difficulties for correct demodulation.To demodulate the information correctly,the signal phase is first amended by the phase detector before estimating the phase CPDFs.

    Two problems must be solved when calculating the signal phase.As shown in figure 2, there are obviously twophase clusters corresponding to bits 0 and 1, but the corresponding relationship cannot be determined.In addition,under the effect of parasitic modulation the position of the phase CPDF is shifted and the center phases move away from±π/2, as shown in figure 3.The first problem can be solved by inserting a small number of pilot symbols to distinguish the clusters.For the second problem, an appropriate phase correction method should be executed so that the two dominating peaks are close to the two fixed values ±π/2.In this work, ideal timing synchronization and carrier frequency synchronization is assumed.In actual systems, many synchronization methods can be transmitted to ensure timing synchronization and carrier synchronization, such as Zadoffchu pilots with given frame headers [41, 42], and phaselocked loop technologies [43, 44].

    The phase correction is performed as follows:

    (1)The phase sequence{φn}of r(t)is first calculated.Next,we divide{φn}into two subsets with a boundary 0,and then acquire the mean value of each subset, p0and p1.

    (2)A phase shift |p1+ p0|/2 is performed so that the two dominating peak phases are shifted to approximately ±π/2.

    (3)Step(1)to Step(2)are repeated tmaxtimes,where tmaxis the maximum iteration number.Finally, the corrected {φn} is divided into two-phase subsets with a boundary 0.

    Letχ0andχ1denote bit subsets representing bits 0 and 1,respectively; Φ0and Φ1can denote the corresponding phase subsets.

    3.3.Adaptive soft detector based on RJ-MCMC algorithm

    Considering that the CPDF of the signal phase randomly fluctuates with the change in the non-stationary PS channel,prior channel information is difficult to achieve at the receiving end.Therefore,the RJ-MCMC algorithm is used to track the phase CPDF.

    Since the corrected phase is disturbed around ±π/2 under the non-stationary PS channel, it can be regarded as a multistate quasi-stationary stochastic variable.The multistate numbers of the two-phase sequences corresponding to bits 0 and 1 are the same because the two sequences experience the same channel.To describe the multistate statistical characteristics, the two received phase sequences can be characterized by mixed Gaussian processes [45] as follows:

    where K is the state number for each mixed Gaussian process,and ωk,μk,andσk2(k = 1,…,2K)are the weight,mean,and variance of the kth sub-distribution, respectively.Furthermore, the RJ-MCMC algorithm [46] is adopted to adaptively distinguish the state number and obtain the distribution parameters.The estimation process is elaborated as follows.

    A two-dimensional sequence (φn, zn)is adopted to represent the multistate time-varying process of the phase sequence φnthrough the PS channel, where zn∈{1, …, 2K}indicates the sequence of the corresponding hidden states.

    With respect to the RJ-MCMC algorithm,the parameters to be searched areλ=(μ,σ,z,ω),whereμ= (μ),σ = (σk), and ω = (ωk).Their prior distributions are set as follows:

    (1)The state number K follows a Poisson distribution with[1, 2..., Kmax], where Kmaxis the largest acceptable value.The initial value K is set as 1.

    (2)The weight vector ω obeys the Dirichlet distribution:ω ~D(δ1,δ2,...,δ2K),with the initial valueofδ1,δ2,...,δ2Kequal to1.

    (3)μk~N(ξ,κ-1)andσ2k~G(α,β),where the parametersβ~G(α,β),α = 2,κ=?-2,ξ= φmin+?/2,g = 0.2,h = 10/?2,?= φmax- φmin, φmax= max{φn}, φmin=min{φn}.G(·)denotes the Gamma distribution.

    The complete joint PDF, as the distribution of interest[47],is described as

    whereη= (α,β,ξ,κ).The parameter spacex=(η,K,ω,z,μ,σ,φ)containing all distribution parameters at iteration N is updated to a new parameter space, i.e.x′ = (η′ ,K′ ,ω′ ,z′ ,μ′ ,σ′ ,φ′)at iteration N + 1.Each iteration from x to x′is categorized in five consecutive moves as follows:

    (1)Updateω:

    where is the number of states k during iteration N.

    (2)Update (μ,σ):

    (3)Update z:

    (4)Update β:

    (5)Split or combine states.

    For the combined states,two channel states are randomly chosen and combined into a new one.Correspondingly, the hidden states are also combined.Conversely, for the split step, a randomly selected state is separated into two.The choice between the split and combine steps depends on the probability dkwhich is a random number following a uniform distribution over (0, 1).

    Figure 5 compares the statistical curves and mixed Gaussian curves of the corrected phase CPDFs under two typical cases, where the Gaussian fitting results obtained by the RJ-MCMC algorithm are in red dashed lines and the statistical results are in black solid lines.The two dominating phases corresponding to bits 1 and 0 are shifted to near±π/2,and the PDF shapes exhibit quasi-symmetric characteristics,which also proves the high calculation accuracy of the phase detector.The six-state Gaussian distribution(2K = 6)fits the phase PDF in case 1 and the two-state Gaussian distribution(2K = 2)fits the phase PDF in case 2 by our algorithm.Observe that the mixed Gaussian curves can fit the phase CPDFs well in the two cases,meaning that the phase CPDFs fitted by the RJ-MCMC algorithm can be used in BICM-ID systems.Additionally, notice that the mixed Gaussian curves have some fitting discrepancies near the extreme point of the curves in figure 5(a).The one reason for the fitting discrepancy is the presence of noise, and another reason is thatthe fitted Gaussian state number is small.The more the number of states, the better the fitting effect, but the higher the computational complexity.However,only the intersection of two-phase PDF curves will affect communication performance, so the fitting discrepancy near the extreme point of the curves does not affect the method performance.The larger the intersection area,the worse the performance.Therefore, the fitting result of the curve edges determines the accuracy of the method.Our fitting method can guarantee this with a small state number.

    Figure 1.(a)Plasma sheath and(b)layered model of plasma sheath.

    Figure 2.Received signals under (a)case 1 and (b)case 2.

    Figure 3.Phase CPDFs under (a)case 1 and (b)case 2.

    Figure 4.Adaptive P-BICM-ID model over plasma sheath channels.

    Figure 5.Results of the statistical model and mixed Gaussian model after phase correction in (a)case 1 and (b)case 2.(Statistical curve displayed in black solid lines and mixed Gaussian curve displayed in red dashed lines.)

    Figure 6.BER comparison of our adaptive P-BICM-ID system with those of [23] and [24] under condition 1.

    Table 2.States parameters of the mixed Gaussian distributions for case 1.

    Table 3.States parameters of the mixed Gaussian distributions for case 2.

    Considering the two signal phase sequences corresponding to bits 0 and 1 experience the same PS channel, the multistate numbers of the two sequences are the same and can be determined after the total state number 2K is adaptively estimated by RJ-MCMC algorithm.All state parameters of the mixed Gaussian distribution on two typical cases are listed in tables 2 and 3.It can be seen that the state number for each mixed Gaussian process is 3 for case 1, and the one is 1 for case 2.

    After all space parameters are achieved by the RJMCMC, the phase CPDFs for bits 0 and 1,fφ∣χ0(φ∣χ0)andfφ∣χ1(φ∣χ1),can be written as

    and

    respectively.

    At the receiving end, the soft decoder produces the soft phase information for the coded bitstream u′, which is denoted by the LLRL(un′)for the nth bit.Owing to the agreeable Gaussian fitting results, the output of the soft decoder for the nth bit,L(un′),can be represented as

    3.4.Iterative URC-BP-decoder

    To improve the system performance, a URC detector and a BP decoder with an iterative feedback loop are employed.The key problem lies in designing the interfaces from the soft detector and BP decoder to the URC decoder input port.

    For the non-stationary channel, LM( ′un), as the output of the soft decoder, is given by

    The a posteriori probability P(un|φn)may be computed using thea prioriinformation L(un′)and the informationLa,i(un)gleaned from the outer decoder, which is defined as

    Explicitly, according to Bayes’ rule, we have

    for b∈{0,1}.Based on equations(20)-(22),the URC decoder generates the extrinsic soft information

    being passed back into the BP decoder.Then, the iterative decoding is performed.

    Generally, compared with the traditional BICM-ID system, the program running time of our proposed adaptive system increases by 11%.However, the highly dynamic and non-stationary characteristics and parasitic modulation effects of the PS channel make the traditional BICM-ID system no longer effective.Our solution can use a small computational cost to solve these problems, and the simulations will be introduced in section 4.

    4.Simulation results and discussion

    This section describes examples of several simulations that were performed to validate the high performance of the new P-BICM-ID system for joint demodulation and decoding under the PS channel.The proposed scheme was compared with the existing methods in terms of BERs under different SNRs.The effects of PS channel parameters and protographs with different code rates and information bit numbers on the new system were also investigated.The simulations prove the superiority of the new method over the other two methods.The proposed adaptive P-BICM-ID system exhibits better performance and has high robustness and adaptability to the highly dynamic PS channel.The plasma sheath parameters in our simulations are described in section 2.Protographs with Rc= {1/2,2/3,3/4,5/6}and L = {528,1056,2112,4224}bits are considered.In addition,we denote the inner and outer iteration numbers by Iiand Io, respectively.

    4.1.Method comparison

    The adaptive BICM-ID system with different iteration numbers Iois first compared with the methods in[23]and[24]for demodulation and decoding under case 1, as shown in figure 6.For the protograph code, Rc= 1/2, L = 528, and Ii= 10.The performance of the methods in [23] and [24]under case 1 is relatively similar but greatly inferior to that of our adaptive P-BICM-ID system.The adaptive P-BICM-ID system has the best performance with a steep-waterfall region since corrected signal phase CPDFs are used to estimate the soft information and the joint optimization of the LDPC iterative soft decoding and P-BICM-ID system brings both high coding gain and excellent error performance.The performance of the adaptive new system improves as Ioincreases and increases by more than 11 dB at BER = 10-6compared to those of the other two methods.In addition, when Ii=10 and Io=10,the system performance is already very good,and the performance improvement space is limited.Therefore,both Iiand Iocan be chosen as 10 in the analysis below.

    Furthermore, we analyze the superiority of the adaptive P-BICM-ID system accurately.It is compared with the methods in[23] and[24]with LDPC codes under case 1 and case 2, as shown in figure 7.The proposed system has a greater than 5 dB improvement over the methods with LDPC codes in [23] and [24] at BER = 10-6under case 1 in figure 7(a).As seen in figure 7(b), the proposed system can reach a BER of 10-6at 22 dB.However, for the other two systems,performance floors appear around BER=10-3as γbincreases.AsNepeakand σmaxincrease, the plasma sheath channel worsens, and the phase jitter becomes stronger.As a result,the schemes with LDPC codes in[23]and[24]cannot cope with this channel situation.From the scheme comparisons, the superiority of the adaptive P-BICM-ID system is proven.

    4.2.Robustness analysis

    Figures 8 and 9 depict the BER curves of the proposed method under the influence ofNepeakand σmaxor Aθ, where Rc= 1/2 and L = 4224.From the figures, the larger Aθ,σmax,andNepeakare,the stronger the channel non-stationarity,and the worse the system performance.In telemetry communication, 10-6is generally set to be a BER threshold to judge whether there is reliable communication.WhenNepeak=1 × 1018m-3,the BERs for different σmaxor Aθare less than 10-6at γb≥4.5 dB, and the BER curves almost overlap.This means that at a lowNepeak,the impact of σmaxor Aθis trivial.However,whenNepeak=1 × 1019m-3,Nepeakis relatively high, the performance worsens, and the impact of σmaxor Aθis obvious.Therefore, to achieve a BER below 10-6, a larger γbis required.The required γbvalues are 8.2 and 11.6 dB when σmax= 0.05 and 0.1, respectively, as shown in figure 8.The required γbvalues are 8.2 dB,11.0 dB,and 20.9 dB at Aθ= 0°,2°,and 4°,respectively, as shown in figure 9.In general, with increasing σmax, Aθ, andNepeak,the system performance decreases.However, although Aθ, σmax,andNepeakreach 4°, 0.1, and 1 × 1019m-3, respectively, the γbrequired to meet BER = 10-6is still within the acceptable range.Thus, the new scheme is considerably robust to the effects of changes in the PS channel on the performance of demodulation and decoding.

    Take the RAM-C reentry flight test [48] as an example.When the reentry flight altitude h = 71 km,Nepeak= 1 × 1017m-3, and Aθ≤4°; when h = 47 km,Nepeak= 1 × 1018m-3, and Aθ≤ 4°; when h = 30 km,Nepeakreaches a maximum of 1 × 1019m-3, and Aθ≤ 4°.The existing communication scheme adopts S- or X- band communication, and the RAM-C vehicle is almost in radio blackout during the flight from 71 km to 30 km [21].However,if the proposed adaptive P-BICM-ID system is adopted,reliable communication during the whole reentry can be ensured when γb≥ 21 dB; reliable communication from 71 km to 47 km can be ensured when γb≥ 4 dB.Therefore,the new system can greatly alleviate radio blackouts.

    4.3.Impact of protographs with different Rc and L

    Furthermore, we analyze the impact of Rcand L on the decoding performance of the proposed scheme.Figure 10 illustrates the BER variation with γbatNepeak=1 × 1018m-3, and figure 11 depicts the BER variation withNepeakat γb= 8 dB,where protographs with different Rcand L are used; σmax= 0.1 and Aθ= 0°.The BERs decrease significantly with increasing γband increase rapidly with increasingNepeak.This is because with a decrease in γbor an increase inNepeak,the relative intensity of the received signal decreases, which consequently increases the BER rapidly.Moreover, with decreasing Rcor increasing L, considerable code gain occurs and brings a significant reduction in the BER.When a protograph with Rc= 5/6 and L = 528 switches to a protograph with Rc= 1/2 and L = 4224,the γbrequired to reduce the BER to 10-5is lowered from 7.8 dB to 4.3 dB, and theNepeaktolerance limit increases to 8.4 × 1018m-3from 1 × 1018m-3for the same target.

    Figure 7.BER curves of our adaptive P-BICM-ID system, [23] with LDPC, and [24] with LDPC under (a)case 1 and (b)case 2.

    Figure 8.BER curves with Ne peak = 1×{1018 m-3, 1019 m-3},σmax = 0.05/0.1, Aθ = 0°.

    Figure 9.BER curves with Ne peak = 1×{1018 m-3, 1019 m-3}, Aθ= 0°-4°, σmax = 0.05.

    Figure 10.BER versus γb, protographs with different Rc and L,Ne peak = 1 × 1018 m-3, Aθ = 0°, and σmax = 0.1.

    Figure 11.BER versus Ne peak ,protographs with different Rc and L,γb = 8 dB, Aθ = 0°, and σmax = 0.1.

    IfNepeakcontinues to increase or γbcontinues to decrease, our adaptive BICM-ID method with a lower Rcor longer L can be adopted to maintain the data link.

    5.Conclusions

    In this work, we conceive an adaptive P-BICM-ID system based on the RJ-MCMC algorithm for the problem of amplitude fading and phase shift caused by the PS channel.This scheme couples demodulation and decoding by combining protograph LDPC iterative decoding and the BICM-ID framework to make a soft decision over the PS channel, where the phase PDF is adaptively estimated by the RJ-MCMC algorithm and then soft information is calculated.Protograph LDPC codes with high coding gain and low encoding complexity are used to design the coding scheme.Compared to the existing algorithms, our adaptive P-BICMID system has a steep-waterfall performance,obtains a greater than 5 dB performance improvement at 10-6BER, and can cope with a worse plasma sheath environment.Simulation results show that even in the worst case for the plasma sheath during the reentry process, i.e.Nepeak,Aθ, and σmaxreach 1 × 1019m-3, 4°, and 0.1, respectively, the proposed adaptive P-BICM-ID system with an acceptable γbcan ensure reliable communication.In addition, protographs with lower Rcor longer L can be used to further improve the system performance.

    Acknowledgments

    The authors gratefully acknowledge financial support from National Natural Science Foundation of China (Nos.61871302, 62101406 and 62001340)and the Fundamental Research Funds for the Central Universities (No.JB211311).

    ORCID iDs

    猜你喜歡
    石磊
    “巨嬰”老公總讓我收拾爛攤子,忍無可忍我決定分居
    婦女生活(2025年2期)2025-02-20 00:00:00
    Nitrogen-tailored quasiparticle energy gaps of polyynes
    Momentum-space polarization fields in two-dimensional photonic-crystal slabs: Physics and applications
    Inverse synthetic aperture radar range profile compensation of plasma-sheathenveloped reentry object
    PERIODIC AND ALMOST PERIODIC SOLUTIONS FOR A NON-AUTONOMOUS RESPIRATORY DISEASE MODEL WITH A LAG EFFECT*
    Energy dissipation and power deposition of electromagnetic waves in the plasma sheath
    闕 題
    親愛的,你現(xiàn)在可以求婚了
    伴侶(2019年11期)2019-08-09 08:47:31
    Probabilistic Teleportation of an Arbitrary Two-Qubit State via Positive Operator-Valued Measurement with Multi Parties?
    昔日創(chuàng)客的美麗蝶變
    国产成人影院久久av| 91久久精品国产一区二区成人 | 亚洲美女视频黄频| 亚洲av成人一区二区三| 亚洲成人久久爱视频| 俄罗斯特黄特色一大片| www.自偷自拍.com| 精品国内亚洲2022精品成人| 亚洲国产色片| 欧美激情在线99| 91麻豆精品激情在线观看国产| 超碰成人久久| 美女被艹到高潮喷水动态| 国产激情久久老熟女| 精品熟女少妇八av免费久了| 噜噜噜噜噜久久久久久91| 最新中文字幕久久久久 | 欧美日韩综合久久久久久 | 久久这里只有精品中国| 午夜两性在线视频| 久久国产精品人妻蜜桃| 特大巨黑吊av在线直播| 女警被强在线播放| 天天躁日日操中文字幕| 黄色视频,在线免费观看| 在线免费观看的www视频| 在线永久观看黄色视频| 五月伊人婷婷丁香| 久久欧美精品欧美久久欧美| 国产精品99久久99久久久不卡| xxx96com| 成熟少妇高潮喷水视频| 又爽又黄无遮挡网站| 成人性生交大片免费视频hd| 亚洲五月天丁香| 噜噜噜噜噜久久久久久91| 长腿黑丝高跟| 久久精品国产清高在天天线| 午夜福利在线在线| 12—13女人毛片做爰片一| 欧美日韩瑟瑟在线播放| 亚洲av成人一区二区三| 黑人欧美特级aaaaaa片| 伦理电影免费视频| 国产精品久久电影中文字幕| 亚洲人与动物交配视频| 久久久久九九精品影院| 精品国产乱码久久久久久男人| 国产野战对白在线观看| 女同久久另类99精品国产91| 国产黄片美女视频| 国产精品99久久久久久久久| 日日夜夜操网爽| 久久午夜综合久久蜜桃| 伦理电影免费视频| 美女扒开内裤让男人捅视频| 夜夜爽天天搞| 国产97色在线日韩免费| 国产av不卡久久| 最近最新中文字幕大全电影3| 国产精品av视频在线免费观看| 国产 一区 欧美 日韩| 午夜福利成人在线免费观看| 亚洲avbb在线观看| 两人在一起打扑克的视频| 少妇裸体淫交视频免费看高清| 美女 人体艺术 gogo| 天天躁日日操中文字幕| 色精品久久人妻99蜜桃| 一个人免费在线观看电影 | 久久午夜亚洲精品久久| 麻豆av在线久日| 18美女黄网站色大片免费观看| 九九热线精品视视频播放| 久久香蕉国产精品| 亚洲国产欧美人成| 午夜福利免费观看在线| 国产精品自产拍在线观看55亚洲| 国产免费av片在线观看野外av| 一a级毛片在线观看| 色噜噜av男人的天堂激情| 男女午夜视频在线观看| 久久久久久人人人人人| 国产三级中文精品| 国产男靠女视频免费网站| 伊人久久大香线蕉亚洲五| 此物有八面人人有两片| 国产精品久久久久久亚洲av鲁大| 亚洲欧美精品综合一区二区三区| 成人精品一区二区免费| 怎么达到女性高潮| 国产又黄又爽又无遮挡在线| 国产精品久久电影中文字幕| 欧美一级毛片孕妇| 成熟少妇高潮喷水视频| 亚洲精品乱码久久久v下载方式 | 久久久久国内视频| 国产精华一区二区三区| 久久精品91无色码中文字幕| 性色avwww在线观看| 亚洲精品久久国产高清桃花| 男女床上黄色一级片免费看| 在线观看一区二区三区| 亚洲精品在线美女| 国产午夜福利久久久久久| 曰老女人黄片| 在线a可以看的网站| h日本视频在线播放| 亚洲色图av天堂| 又黄又爽又免费观看的视频| 一个人看视频在线观看www免费 | 一本一本综合久久| 噜噜噜噜噜久久久久久91| 成人特级av手机在线观看| 免费观看人在逋| 国产精品av久久久久免费| 亚洲精品国产精品久久久不卡| 十八禁网站免费在线| 最近最新中文字幕大全免费视频| 免费人成视频x8x8入口观看| 最近视频中文字幕2019在线8| 久久精品夜夜夜夜夜久久蜜豆| av福利片在线观看| 一二三四在线观看免费中文在| 国产成人精品久久二区二区免费| 国产成人欧美在线观看| 精品熟女少妇八av免费久了| 午夜免费激情av| 国产欧美日韩一区二区精品| 麻豆成人午夜福利视频| 禁无遮挡网站| 男人和女人高潮做爰伦理| 国产欧美日韩一区二区精品| 国产麻豆成人av免费视频| 最近在线观看免费完整版| 最近在线观看免费完整版| 国产成人一区二区三区免费视频网站| 国产欧美日韩精品亚洲av| 亚洲精品色激情综合| 1024香蕉在线观看| 99热只有精品国产| 日韩三级视频一区二区三区| 国产一区二区三区在线臀色熟女| 午夜a级毛片| 久久国产乱子伦精品免费另类| 日韩中文字幕欧美一区二区| 欧美日韩黄片免| 国产精品久久久久久亚洲av鲁大| 国产精品亚洲一级av第二区| 两个人的视频大全免费| 伊人久久大香线蕉亚洲五| 国产成人系列免费观看| 久久这里只有精品中国| 国产av不卡久久| 精品久久久久久久毛片微露脸| 亚洲精品美女久久av网站| 亚洲av日韩精品久久久久久密| 亚洲精品美女久久久久99蜜臀| 久久这里只有精品中国| 国产视频一区二区在线看| 国产精品电影一区二区三区| 国产精品电影一区二区三区| netflix在线观看网站| 又黄又爽又免费观看的视频| 国产精品亚洲一级av第二区| 每晚都被弄得嗷嗷叫到高潮| 久久精品人妻少妇| 欧美+亚洲+日韩+国产| 国产精品综合久久久久久久免费| 久久久久国产精品人妻aⅴ院| 国产乱人伦免费视频| 一级黄色大片毛片| 国产99白浆流出| 老熟妇仑乱视频hdxx| 色在线成人网| 国产又色又爽无遮挡免费看| 国产精品一及| 成在线人永久免费视频| 18禁美女被吸乳视频| 这个男人来自地球电影免费观看| 岛国在线免费视频观看| 嫁个100分男人电影在线观看| 国产精品av久久久久免费| 国产黄色小视频在线观看| 中亚洲国语对白在线视频| 国产美女午夜福利| 免费在线观看亚洲国产| 熟女少妇亚洲综合色aaa.| 国产综合懂色| 亚洲最大成人中文| 黄色 视频免费看| 18禁黄网站禁片免费观看直播| 免费看a级黄色片| 一区福利在线观看| 久久精品国产清高在天天线| 免费观看人在逋| 欧美成人性av电影在线观看| 久9热在线精品视频| 蜜桃久久精品国产亚洲av| 欧美最黄视频在线播放免费| 国产精品99久久久久久久久| 欧美激情久久久久久爽电影| 国内少妇人妻偷人精品xxx网站 | 午夜福利免费观看在线| 亚洲七黄色美女视频| 性色avwww在线观看| 99久久精品国产亚洲精品| 人妻久久中文字幕网| 人妻丰满熟妇av一区二区三区| 国产久久久一区二区三区| 免费在线观看日本一区| 不卡av一区二区三区| 在线观看免费视频日本深夜| 男人舔女人的私密视频| 久久99热这里只有精品18| 变态另类丝袜制服| 99国产精品一区二区蜜桃av| 成年人黄色毛片网站| 99久国产av精品| 99国产精品一区二区三区| 午夜影院日韩av| 嫩草影院精品99| 久久天躁狠狠躁夜夜2o2o| 99久久国产精品久久久| 国产av一区在线观看免费| 欧美高清成人免费视频www| 中文亚洲av片在线观看爽| 成人三级黄色视频| 黄片小视频在线播放| 国产亚洲欧美在线一区二区| 亚洲精华国产精华精| 99热6这里只有精品| 国产 一区 欧美 日韩| 天天一区二区日本电影三级| 日韩免费av在线播放| 美女黄网站色视频| 欧美性猛交╳xxx乱大交人| 一本久久中文字幕| 91麻豆av在线| 亚洲无线观看免费| 久久久久久人人人人人| 久久这里只有精品中国| а√天堂www在线а√下载| 国产成人aa在线观看| 久久久久久久久免费视频了| 亚洲午夜精品一区,二区,三区| www日本在线高清视频| 亚洲电影在线观看av| 久久这里只有精品中国| 一个人免费在线观看电影 | 99国产精品99久久久久| 国产亚洲精品久久久久久毛片| 日本免费一区二区三区高清不卡| 村上凉子中文字幕在线| 成年女人永久免费观看视频| 日本 欧美在线| 在线看三级毛片| 九色国产91popny在线| 婷婷精品国产亚洲av| 搡老岳熟女国产| 国产精品九九99| 桃红色精品国产亚洲av| 欧美中文综合在线视频| 亚洲国产精品合色在线| 精品久久久久久,| 欧美日韩中文字幕国产精品一区二区三区| 国产一区二区激情短视频| a在线观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 又黄又粗又硬又大视频| 久久人人精品亚洲av| 久久亚洲精品不卡| 亚洲精品色激情综合| 一个人观看的视频www高清免费观看 | 久久精品91蜜桃| 久久国产乱子伦精品免费另类| 午夜a级毛片| 搡老妇女老女人老熟妇| 免费在线观看成人毛片| 亚洲中文日韩欧美视频| 日本黄大片高清| 巨乳人妻的诱惑在线观看| 嫩草影院精品99| 最近在线观看免费完整版| 色尼玛亚洲综合影院| 国产亚洲精品综合一区在线观看| 婷婷亚洲欧美| 男女下面进入的视频免费午夜| 亚洲精品国产精品久久久不卡| 亚洲性夜色夜夜综合| 在线免费观看不下载黄p国产 | 成年女人毛片免费观看观看9| 婷婷丁香在线五月| 欧美高清成人免费视频www| 1000部很黄的大片| 亚洲国产欧洲综合997久久,| 国产成人影院久久av| av黄色大香蕉| 亚洲最大成人中文| 精品电影一区二区在线| 91麻豆精品激情在线观看国产| 国产一区二区在线av高清观看| 天天躁日日操中文字幕| 精品99又大又爽又粗少妇毛片 | 免费看a级黄色片| 亚洲精品在线观看二区| 亚洲精品美女久久久久99蜜臀| 国产黄色小视频在线观看| 又爽又黄无遮挡网站| 国产熟女xx| 亚洲美女黄片视频| 听说在线观看完整版免费高清| 精品欧美国产一区二区三| 免费在线观看日本一区| 欧美一区二区精品小视频在线| www日本黄色视频网| 18禁美女被吸乳视频| 亚洲欧美精品综合久久99| 91在线观看av| 国产三级在线视频| 琪琪午夜伦伦电影理论片6080| 亚洲性夜色夜夜综合| 亚洲国产精品成人综合色| 91麻豆av在线| 免费在线观看日本一区| 国产乱人伦免费视频| 亚洲自偷自拍图片 自拍| 国语自产精品视频在线第100页| 精品国产亚洲在线| 2021天堂中文幕一二区在线观| 村上凉子中文字幕在线| 日韩三级视频一区二区三区| 精品国产乱码久久久久久男人| 最近视频中文字幕2019在线8| 日日干狠狠操夜夜爽| 一本一本综合久久| 嫩草影院入口| 午夜精品久久久久久毛片777| 国产单亲对白刺激| 日韩欧美国产在线观看| 狂野欧美激情性xxxx| 国产精品 欧美亚洲| 国产精品亚洲美女久久久| 搡老岳熟女国产| 国产亚洲欧美在线一区二区| 91久久精品国产一区二区成人 | 精品人妻1区二区| 国产私拍福利视频在线观看| 国产精品精品国产色婷婷| 女人被狂操c到高潮| 伦理电影免费视频| 亚洲最大成人中文| 亚洲美女视频黄频| 亚洲av电影在线进入| 听说在线观看完整版免费高清| 久久天堂一区二区三区四区| 欧美成人一区二区免费高清观看 | 99精品欧美一区二区三区四区| 天堂√8在线中文| 久久亚洲真实| 日本撒尿小便嘘嘘汇集6| 久久精品亚洲精品国产色婷小说| 久久精品亚洲精品国产色婷小说| 黄色成人免费大全| 国内揄拍国产精品人妻在线| 国产一区二区在线av高清观看| 亚洲无线在线观看| 俺也久久电影网| 久久久国产成人精品二区| 国产激情欧美一区二区| 亚洲电影在线观看av| av欧美777| 欧美成人一区二区免费高清观看 | 亚洲国产欧美人成| 亚洲欧美日韩高清专用| 国产精品久久电影中文字幕| 99国产精品99久久久久| 亚洲精品在线观看二区| 国产一区二区激情短视频| 亚洲欧美激情综合另类| 啦啦啦观看免费观看视频高清| 国产视频一区二区在线看| 国产三级在线视频| 宅男免费午夜| 日韩av在线大香蕉| 久久婷婷人人爽人人干人人爱| 国产激情欧美一区二区| 久久性视频一级片| 精华霜和精华液先用哪个| 黑人巨大精品欧美一区二区mp4| 免费在线观看影片大全网站| 草草在线视频免费看| av在线天堂中文字幕| 99热这里只有精品一区 | 亚洲成av人片免费观看| 啦啦啦韩国在线观看视频| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站| 偷拍熟女少妇极品色| 免费搜索国产男女视频| 久久伊人香网站| 国产精品一区二区三区四区免费观看 | 网址你懂的国产日韩在线| 国产精品香港三级国产av潘金莲| 一进一出好大好爽视频| 在线观看午夜福利视频| 久久久久精品国产欧美久久久| 欧美乱色亚洲激情| 黑人巨大精品欧美一区二区mp4| 好男人电影高清在线观看| 国产91精品成人一区二区三区| 这个男人来自地球电影免费观看| 麻豆成人午夜福利视频| 99久国产av精品| 国产精品亚洲美女久久久| 男人舔奶头视频| 五月伊人婷婷丁香| 变态另类丝袜制服| 男人的好看免费观看在线视频| 亚洲第一欧美日韩一区二区三区| 美女被艹到高潮喷水动态| 国产三级在线视频| 嫩草影院精品99| 久久久色成人| 一个人看视频在线观看www免费 | 99riav亚洲国产免费| 色尼玛亚洲综合影院| 我要搜黄色片| 美女大奶头视频| 精品国产乱码久久久久久男人| 国产亚洲av嫩草精品影院| 制服人妻中文乱码| h日本视频在线播放| 午夜福利欧美成人| 成年人黄色毛片网站| 丁香六月欧美| 久久精品91蜜桃| 国产精品女同一区二区软件 | 国产亚洲欧美在线一区二区| 91久久精品国产一区二区成人 | 国产伦精品一区二区三区四那| 少妇熟女aⅴ在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品成人综合色| 日本三级黄在线观看| 免费在线观看影片大全网站| 国产三级黄色录像| 特大巨黑吊av在线直播| 岛国在线观看网站| 日韩三级视频一区二区三区| 女警被强在线播放| 精品福利观看| 免费av毛片视频| cao死你这个sao货| 亚洲国产精品久久男人天堂| 手机成人av网站| 国产一区二区在线av高清观看| 一本久久中文字幕| 亚洲无线观看免费| 少妇人妻一区二区三区视频| 久久久精品欧美日韩精品| www.熟女人妻精品国产| 成年版毛片免费区| 少妇人妻一区二区三区视频| 毛片女人毛片| 99精品在免费线老司机午夜| avwww免费| 两性夫妻黄色片| 国产伦一二天堂av在线观看| www日本黄色视频网| 国产成人av教育| 国产精品一区二区三区四区免费观看 | 亚洲欧美精品综合一区二区三区| 最好的美女福利视频网| 岛国在线免费视频观看| 国产成人av激情在线播放| 黄色 视频免费看| 99久久无色码亚洲精品果冻| 亚洲第一电影网av| www.自偷自拍.com| 俄罗斯特黄特色一大片| 亚洲一区高清亚洲精品| 欧美乱色亚洲激情| 一a级毛片在线观看| 99riav亚洲国产免费| 亚洲va日本ⅴa欧美va伊人久久| 91字幕亚洲| 国产成人福利小说| 伦理电影免费视频| 一级黄色大片毛片| 波多野结衣高清作品| 午夜精品在线福利| 国产高清三级在线| 国产精品综合久久久久久久免费| 久久精品影院6| 在线播放国产精品三级| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 成人欧美大片| 久久精品国产清高在天天线| 熟女人妻精品中文字幕| 日韩国内少妇激情av| 露出奶头的视频| 国产精品 欧美亚洲| 2021天堂中文幕一二区在线观| 91av网站免费观看| 免费在线观看日本一区| 亚洲美女黄片视频| 狂野欧美白嫩少妇大欣赏| 熟妇人妻久久中文字幕3abv| 欧美在线黄色| 国产麻豆成人av免费视频| 日韩精品中文字幕看吧| 欧美性猛交黑人性爽| 成熟少妇高潮喷水视频| 午夜免费成人在线视频| 国产精品精品国产色婷婷| 成人三级做爰电影| 免费看a级黄色片| 亚洲精华国产精华精| 日韩精品中文字幕看吧| 日韩欧美三级三区| 国产精品 欧美亚洲| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 欧美不卡视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 99久久精品国产亚洲精品| 老汉色∧v一级毛片| 日韩欧美免费精品| 成人亚洲精品av一区二区| 亚洲熟女毛片儿| 88av欧美| 久久久久精品国产欧美久久久| 搞女人的毛片| 床上黄色一级片| 最新在线观看一区二区三区| 国内精品久久久久精免费| av中文乱码字幕在线| 长腿黑丝高跟| 99久久精品国产亚洲精品| 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 男人的好看免费观看在线视频| 国产欧美日韩一区二区精品| 免费看光身美女| 搡老妇女老女人老熟妇| 国产高清三级在线| 1024香蕉在线观看| 狂野欧美白嫩少妇大欣赏| 日日摸夜夜添夜夜添小说| www.999成人在线观看| 国产成人aa在线观看| 精品电影一区二区在线| 免费大片18禁| 国产精品亚洲美女久久久| 精品乱码久久久久久99久播| 日本在线视频免费播放| 不卡一级毛片| 毛片女人毛片| 成人永久免费在线观看视频| 最近最新中文字幕大全免费视频| 欧美丝袜亚洲另类 | 国产亚洲av高清不卡| 免费大片18禁| 一个人看视频在线观看www免费 | 99热精品在线国产| 成年女人永久免费观看视频| 啦啦啦观看免费观看视频高清| 最好的美女福利视频网| 亚洲精品中文字幕一二三四区| 国产精品电影一区二区三区| 久久精品aⅴ一区二区三区四区| 高清毛片免费观看视频网站| 日日夜夜操网爽| 久久精品人妻少妇| 中文字幕精品亚洲无线码一区| 亚洲精品色激情综合| 精品国产乱子伦一区二区三区| 999久久久国产精品视频| av天堂中文字幕网| 国产一区二区在线av高清观看| 欧美乱妇无乱码| 一级黄色大片毛片| 国产精品爽爽va在线观看网站| 黄色丝袜av网址大全| 村上凉子中文字幕在线| 日本撒尿小便嘘嘘汇集6| 亚洲人成网站在线播放欧美日韩| 精品免费久久久久久久清纯| 久久热在线av| 国产精品久久久久久精品电影| av女优亚洲男人天堂 | 人人妻,人人澡人人爽秒播| 99国产精品一区二区蜜桃av| 手机成人av网站| 国产黄a三级三级三级人| 国产1区2区3区精品| 国产91精品成人一区二区三区| 中出人妻视频一区二区| 成人av一区二区三区在线看| 午夜激情欧美在线| 亚洲国产精品合色在线| 久久亚洲真实| 国产探花在线观看一区二区| 国产三级在线视频| 国产亚洲av嫩草精品影院| 欧美精品啪啪一区二区三区| 国产午夜精品久久久久久| 欧美国产日韩亚洲一区| 91麻豆精品激情在线观看国产| 精品熟女少妇八av免费久了| 国产毛片a区久久久久| 免费人成视频x8x8入口观看| 黑人操中国人逼视频|