CN103560864B - Blind channel self-adaptation method and device - Google Patents

The invention discloses a blind channel self-adaptation method and device. The blind channel self-adaptation method comprises the following steps of estimating phase noise of input signals, estimating amplitude attenuation of the input signals, detecting pulse interference of the input signals, weighing a phase noise estimation value and an amplitude attenuation estimation value to obtain a channel quality estimation value, carrying out threshold comparison to obtain the self-adaptation information of a channel, and by utilizing the channel self-adaptation information, finishing a channel self-adaptation process in a self-adaptation sending mode through a sending end. The blind channel self-adaptation device comprises a phase angle calculation and quantization module, a demapping module, a phase noise estimation module, an amplitude calculation and quantization module, an amplitude attenuation and estimation module, an interference detection module, a weighing module, a channel self-adaptation module and a self-adaptation sending module. Through coefficient-variable weighing processing on phase noise and amplitude attenuation, the blind channel self-adaptation method and device improve the accuracy of channel estimation, are high in applicability, and are high in transmission efficiency. The sending end adjusts the parameters of sending signals in a self-adaptation mode according to the channel quality, and makes full use of the channel resources.

A kind of blind channel self-adaptation method and device

Technical field

The present invention relates to digital communicating field is and in particular to a kind of blind Channel self adaptation being applied to multi-carrier reception system Method and its device.

Background technology

At present, in general communication system, channel circumstance all compares badly, except signal itself is with communication distance Beyond increase decay is gradually increased, also there are various noise jamming, such as Gaussian noise, impulsive noise etc., in addition, Multipath Transmission The receptivity of signal can be produced serious influence.So, want to obtain in channel so complicated and changeable and reliably lead to It is necessary to take multiple channel interference protection measures, wherein, channel self-adapting is commonly used one kind in existing communication system to letter Improve the technology of systematic function.Channel self-adapting technology is the estimated value according to present channel quality, and transmitter adaptively changes Become modulation system, signal bandwidth, coded system, transmission power, power adaptive distribution etc., to realize ensureing receipt signal matter On the premise of amount, substantially utilize channel resource.

The premise of channel self-adapting is the estimation of channel quality, and its method can be divided mainly into two classes, and a class is by calculating Bit error rate（BER）, the bit error rate（SER）Or packet loss（PLR）To obtain, but error rate be not only relevant with channel quality, also with Modulation system, coded system etc. are relevant, so current channel quality condition can not be reflected exactly, in addition, mistake of statistics Rate needs to send mass data.Another kind of be by estimate receipt signal Signal Interference and Noise Ratio（SINR, Signal to Interference plus Noise Ratio）Obtain, at present, what most systems used is all the method.Lot of domestic and international The method of the estimation to SINR for the experts and scholars has also done substantial amounts of work, and its method is broadly divided into two kinds, and one kind is to be based on Data-aided non-blind estimating method, the method needs the reliability that substantial amounts of pilot signal or training sequence guarantee are estimated Property, a kind of is blind estimating method based on unbound nucleus, and the method adopts special algorithm directly receipt signal to be processed.This Two methods respectively have pluses and minuses, and non-blind estimating method estimates that performance is relatively preferable, but reduces efficiency of transmission, and the availability of frequency spectrum is relatively Low, in today of frequency spectrum resource growing tension, it is not very high occasion typically without the method to estimation performance requirement.Blind estimate The characteristic that method makes full use of signal itself is estimated the assistance data it is not necessary to transmitting terminal, and efficiency of transmission is high, has wide Application space, but directly the algorithm complex of blind estimate SINR is high at present, realizes difficult, practical application is not strong, performance Have much room for improvement.

Content of the invention

It is an object of the invention to provide a kind of complexity is little, the availability of frequency spectrum is high, realization is simple and accuracy is high Blind channel self-adaptation method and device.

The technical scheme realizing above-mentioned purpose is：

A kind of blind channel self-adaptation method of one of the present invention, comprises the following steps：

Step S0, phase noise is estimated, according to input signal in the phase offset of planisphere, unites in time domain and frequency domain respectively Meter phase noise；

Step S1, amplitude fading is estimated, according to the amplitude of input signal, in time domain and frequency domain, statistics amplitude declines respectively Subtract；

Step S2, Interference Detection, based on default interference threshold, detects the impulse disturbances of input signal；

Step S3, weighting, according to default weight coefficient, weighting processes phase noise and amplitude fading, obtains channel matter The estimated value of amount；

Step S4, channel self-adapting, according to the estimated value of channel quality, based on preset channel thresholding, obtain channel adaptive Answer information；

Step S5, self adaptation sends, according to the channel self-adapting information of receiving end feedback, the transmission ginseng of adjustment data of making a start Number.

Above-mentioned blind channel self-adaptation method, wherein, described step S0 includes：

Calculate the phase value amount of input signal；

Phase quantization, based on Interference Detection, is melted into some bits the phase value amount even amount obtaining；

According to the phase quantization result obtaining, the calculating phase offset on planisphere, and respectively in time domain and frequency domain Carry out statistical average.

Above-mentioned blind channel self-adaptation method, wherein, the input signal of described step S0, for differential modulation mode, be Refer to the signal after differential ference spiral.

Above-mentioned blind channel self-adaptation method, wherein, described phase quantization result, also by non-uniform quantizing, completes to solve Mapping process, exports soft bit, to follow-up soft decoding process.

Above-mentioned blind channel self-adaptation method, wherein, described step S1 includes：

Real part according to input signal and imaginary part, calculate amplitude；

Based on Interference Detection, if the uniform amplitude obtaining is quantized into kilo byte；

According to amplitude quantization result, statistical average is carried out on time domain and frequency domain respectively.

Above-mentioned blind channel self-adaptation method, wherein, the input signal of described step S1, the planisphere made a start has Identical amplitude.

Above-mentioned blind channel self-adaptation method, wherein, described amplitude, its standard value is real part and the quadratic sum of imaginary part Arithmetic square root.

Above-mentioned blind channel self-adaptation method, wherein, described step S2 refers specifically to：By the amplitude of comparator input signal with The size of default interference threshold obtains disturbing indication signal.

Above-mentioned blind channel self-adaptation method, wherein, described step S3 specifically includes：

Step S30, deducts current phase noise with the maximum that phase noise occurs and obtains a new phase noise table Reach value；

Step S31, amplitude fading is normalized into and has identical bit wide with new phase noise expression values；

Step S32, the phase noise expression values that step S30 is obtained are multiplied by weight coefficient α, the width that step S31 is obtained Weight coefficient β is multiplied by degree decay, and acquired results carry out being added the estimated value obtaining channel quality.

Above-mentioned blind channel self-adaptation method, wherein, the span of described weight coefficient α and β is [0,1], and full Sufficient alpha+beta=1.

Above-mentioned blind channel self-adaptation method, wherein, described step S4 includes：

Step S40, compares the size of symbol channel quality estimates and predetermined symbol channel thresholding and carries out symbol quality Count；

Step S41, compares the size of carrier channel quality evaluation value and default carrier channel thresholding and carries out carrier-quality Count；

Step S42, calculates the overall average of the quality being all worth to channel of carrier channel quality evaluation value；

Step S43, the count value based on step S40, S41 and the average of step S42 rule out modulation system；

Step S44, selects effective subcarrier number, estimates sub-band decay, and calculates transmission power.

Above-mentioned blind channel self-adaptation method, wherein, described step S5 includes：

Step S50, based on the modulation system selecting, is mapped to data in corresponding constellation point；

Step S51, selects only secured transmission of payload data on effective subcarrier or on sub-band；

Step S52, according to sub-band attenuation, on the premise of keeping frequency domain total transmission power constant, in quality relatively Reduce transmission power on good sub-band, transmission power is increased on second-rate sub-band；

Step S53, raises whole transmission power in time domain.

A kind of blind Channel self-reacting device based on the described blind channel self-adaptation method of one of present invention of the two of the present invention, Including calculating phase angle quantization modules, De-mapping module, phase noise estimation module, calculate amplitude quantization modules, amplitude Decay behavior module, interference detection module, weighting block, channel self-adapting module and self adaptation sending module；Described calculating phase Parallactic angle quantization modules connection De-mapping module and phase noise estimation module；Described phase noise estimation module connects described adding Power module；Described calculating amplitude the quantization modules described interference detection module of connection and amplitude fading estimation module；Described amplitude Decay behavior module connects weighting block；Described interference detection module connects phase noise estimation module and amplitude fading estimates mould Block；Described weighting block connects described channel self-adapting module；Described channel self-adapting module connects described self adaptation and sends mould Block, wherein：

Calculate phase angle quantization modules, for calculating the phase angle of input signal, even amount passes after being melted into some bits To described De-mapping module and phase noise estimation module；

De-mapping module, according to the phase angle quantized value receiving, realizes the inverse process of mapping, after giving by non-uniform quantizing Continuous decoder module；

Phase noise estimation module, based on phase offset on planisphere for the phase angle quantized value, on symbol and carrier wave Upper statistical average phase noise, result passes to described weighting block；

Calculate amplitude quantization modules, for calculating the amplitude of input signal, even amount passes to institute after being melted into some bits State amplitude fading estimation module and interference detection module；

Amplitude fading estimation module, based on amplitude quantizing value, statistical average amplitude fading, result on symbol and on carrier wave Pass to described weighting block；

Interference detection module, by the size of the amplitude of comparison signal and default interference threshold value, during more than thresholding, output Interference indication signal gives described phase noise estimation module and amplitude fading estimation module；

Weighting block, according to weight coefficient, weighting processes phase noise and amplitude attenuation, the channel quality estimation obtaining Value passes to described channel self-adapting module；

Channel self-adapting module, according to the channel quality estimates of input, based on preset channel thresholding, is compared, counts Number and judgement obtain channel self-adapting information, and result is passed to described self adaptation sending module；

Self adaptation sending module, according to the channel self-adapting information of input, by adjusting the transmission parameter of data, realizes letter Road self adaptation.

Above-mentioned blind Channel self-reacting device, wherein, described calculating phase angle quantization modules include phase angle calculate single Unit and phase quantization unit, wherein：

Phase angle computing unit, is used for realizing rectangular coordinate to polar conversion, obtains the phase angle of input signal Value；

Phase quantization unit, for the value even amount of phase angle is melted into some bits, passes to described De-mapping module With phase noise estimation module.

Above-mentioned blind Channel self-reacting device, wherein, described phase quantization unit includes the first multiplier, the first divider With round unit, wherein：

Input phase angle and 2 realized by first multiplierⁿMultiplication, n be quantization bit；

First divider is used for the output of the first multiplier divided by 2 π, makes output normalize to [0,2ⁿ]；

Round unit and take out integer-bit output, output valve is [0,2^n-1] between integer.

Above-mentioned blind Channel self-reacting device, wherein, described De-mapping module includes form memory element and lookup unit, Wherein：

Form memory element, for storing phase angle to the index of mapped bits；

Lookup unit, for realizing phase angle to the Index process of mapped bits.

Above-mentioned blind Channel self-reacting device, wherein, described phase noise estimation module includes phase deviation computing unit With phase noise statistic unit, wherein：

Phase deviation computing unit, based on Interference Detection indication signal, for calculating deviation on planisphere for the phase place, by The first adder being sequentially connected in series, unit, the first subtractor and the unit that takes absolute value composition more than mould, wherein：

First adder realizes the signal and 2 of its inputⁿ¹Addition, n1 be quantization bit；

More than mould, unit realizes taking the remainder of the signal of its input；

The signal that first subtractor realizes its input deducts 2ⁿ¹；

The signal normalization that the unit that takes absolute value inputs it is to 2ⁿ¹In the range of；

Phase noise statistic unit, for statistical average phase noise on symbol and on carrier wave respectively, obtains symbol phase Position noise and Carrier Phase Noise, make an uproar including the first accumulator, the second accumulator, Carrier Phase Noise memorizer, symbol phase The acoustic memory, the second divider and the 3rd divider, wherein：

First accumulator, for adding up, current sign Carrier Phase Noise is made an uproar with the carrier phase of corresponding symbol before Sound；

Second accumulator, for all Carrier Phase Noises of the current sign that adds up；

Carrier Phase Noise memorizer, for store the cumulative of each symbol carrier phase noise and, length is equal to carrier wave Number；

Symbol phase noise memory, for storing the cumulative of each symbols phase noise and length equal symbol Number；

Second divider, realizes the cumulative of Carrier Phase Noise and the division with symbol numbers, obtains Carrier Phase Noise Average；

3rd divider, realizes the cumulative of symbol phase noise and the division with variable number, obtains symbol phase noise Average.

Above-mentioned blind Channel self-reacting device, wherein, described calculating amplitude and quantization modules include amplitude computing unit and Amplitude quantization unit, wherein：

Amplitude computing unit, for calculating the amplitude of input signal；

Amplitude quantization unit, if for amplitude quantization is become kilo byte, estimate to described interference detection module and amplitude fading Meter module.

Above-mentioned blind Channel self-reacting device, wherein, described amplitude fading estimation module is single with described phase noise statistics The composition of unit is identical.

Above-mentioned blind Channel self-reacting device, wherein, described weighting block includes phase noise processing unit, amplitude fading Processing unit, the second multiplier, the 3rd multiplier and second adder, wherein,

Phase noise processing unit, is individually subtracted symbol phase noise and the carrier wave phase of input with the maximum of phase noise Position noise, obtains phase noise processing costs；

Amplitude fading processing unit, for keeping amplitude fading and phase noise processing costs to have identical bit wide, obtains Amplitude fading processing costs；

Second multiplier, is multiplied by the processing costs of phase noise for realizing weight coefficient α；

3rd multiplier, is used for realizing weight coefficient β being multiplied by amplitude fading processing costs；

Second adder, for realizing the weighting of the phase noise processing costs after multiplication process and amplitude fading processing costs, Obtain symbol quality evaluation value and carrier channel quality evaluation value.

Above-mentioned blind Channel self-reacting device, wherein, described channel self-adapting module includes comparator, enumerator and judgement Device, wherein：

Comparator, is used for comparing symbol channel quality estimates and carrier channel quality evaluation value and corresponding predetermined symbol Channel thresholding and the size of default carrier channel thresholding；

Enumerator, for counting less than the symbol numbers of predetermined symbol channel thresholding with more than default carrier channel thresholding Variable number；

Mean module, for calculating the average of carrier channel quality evaluation value；

Modulation system determining unit, for determining modulation system based on count value and average；

Effectively subcarrier determining unit, based on a determination that modulation system, determine corresponding effective subcarrier sequence number；

Calculate sub-band gain coefficient elements, for carrier channel quality evaluation value is converted into corresponding sub-band gain Coefficient；

Calculate transmission power coefficient elements, based on a determination that modulation system and channel quality estimation average, calculate transmitting work( The coefficient of rate.

Above-mentioned blind Channel self-reacting device, wherein, described self adaptation sending module includes modulation system select unit, has Effect sub-carrier allocation unit, subband power allocation unit and transmission power adjustment unit, wherein：

Modulation system select unit, the modulation system being come according to receiving end feedback, for being mapped to corresponding star data On seat point；

Effectively sub-carrier allocation unit, feeds back, according to receiving end, effective subcarrier sequence number that comes, in effective subcarrier Upper placement transmission data；

Subband power allocation unit, the sub-band fading gain coefficient being come according to receiving end feedback, input sub-band letter Number it is multiplied by the adjustment that corresponding gain coefficient realizes subband power；

Transmission power adjustment unit, according to the receiving end transmission power coefficient that comes of feedback it is ensured that less than emission maximum On the premise of power, input signal is multiplied by the adjustment that transmission power coefficient realizes transmission power.

The invention has the beneficial effects as follows：The blind channel self-adaptation method of the present invention, because not requiring to send assistance data, thus Improve transfer rate, the availability of frequency spectrum of system is high；Meanwhile, for different systems, by using different weight coefficients Process phase noise and amplitude fading is estimated to obtain channel quality so that accuracy of estimation is high, and simple；In addition, this Interference Detection is passed through in invention, decreases impulse disturbances well.

Brief description

Fig. 1 is the flow chart of the blind channel self-adaptation method of one of the present invention；

Fig. 2 is the schematic flow sheet that phase noise provided in an embodiment of the present invention is estimated；

Fig. 3 is the process schematic of mapping provided in an embodiment of the present invention and demapping；

Fig. 4 is the process schematic that phase deviation provided in an embodiment of the present invention calculates；

Fig. 5 is the schematic diagram calculating phase noise time domain and frequency domain statistical average provided in an embodiment of the present invention；

Fig. 6 is the schematic flow sheet that amplitude fading provided in an embodiment of the present invention is estimated；

Fig. 7 is weighting procedure schematic diagram provided in an embodiment of the present invention；

Fig. 8 is self adaptation transmission flow schematic diagram of making a start provided in an embodiment of the present invention；

Fig. 9 is the structured flowchart of the two blind Channel self-reacting device of the present invention；

Figure 10 is the structured flowchart of calculating phase angle provided in an embodiment of the present invention quantization modules；

Figure 11 is the structured flowchart of phase noise estimation module provided in an embodiment of the present invention；

Figure 12 is the structured flowchart of weighting block provided in an embodiment of the present invention；

Figure 13（a）Amplitude-versus-frequency curve for multipath channel in the embodiment of the present invention；

Figure 13（b）For receipt signal amplitude-versus-frequency curve under the multi-path channel environment of emulation in the embodiment of the present invention；

Figure 13（c）For phase of received signal noise characteristic curve under the multi-path channel environment of emulation in the embodiment of the present invention；

Figure 14 be the embodiment of the present invention in emulation Gaussian channel environment under receipt signal average amplitude attenuation curve and Average phase noise curve；

Figure 15（a）For the time domain waveform of the receipt signal of emulation and spectrum waveform in the embodiment of the present invention；

Figure 15（b）Average amplitude decay waveform for emulation in the embodiment of the present invention；

Figure 15（c）Average phase noise pattern for emulation in the embodiment of the present invention.

Specific embodiment

Below in conjunction with accompanying drawing, the invention will be further described.

The application scenarios of the present invention are half-duplex, symmetric channel, and receiving-transmitting sides adopt constant amplitude phase-modulation, i.e. constellation mapping Point is on same circle.

Refer to Fig. 1, the blind channel self-adaptation method of one of present invention, comprise the following steps：

Step S0, phase noise estimate, according to input signal planisphere phase offset, respectively in time domain（Based on symbol Number）And frequency domain（Based on carrier wave）Statistics phase noise；

Step S1, amplitude fading is estimated, according to the amplitude of input signal, in time domain and frequency domain, statistics amplitude declines respectively Subtract；

Step S2, Interference Detection, based on default interference threshold, detects the impulse disturbances of input signal；

Step S3, weighting, according to default weight coefficient, weighting processes phase noise and amplitude fading, obtains channel The estimated value of quality；

Step S4, channel self-adapting, according to the estimated value of channel quality, based on preset channel thresholding, obtain channel adaptive Answer information；

Step S5, self adaptation sends, according to the channel self-adapting information of receiving end feedback, the transmission ginseng of adjustment data of making a start Number.

Wherein, as shown in Fig. 2 step S0 specifically includes：

Step S00, calculates input signal（It is made up of real part I and imaginary part Q）Phase angle, realize input signal from right angle sit Mark polar conversion；Refer specifically to calculate arctan (Q/I) value, span be [0,2 π), wherein, if input signal exists Make a start and carried out differential modulation, here with the value after differential ference spiral；

Step S01, phase value amount even amount that step S00 is obtained is melted into some bits, compromise consider resource consumption and Performance selects suitable quantization digit, uses 7 bit quantizations, concrete quantization method is in embodimentWherein θ is input Phase angle, [] is rounding operation；

Step S02, based on the phase quantization value of step S01 output, by non-uniform quantizing, completes de-mapping process, defeated Go out soft bit, facilitate follow-up soft decoding；3 bit quantizations are selected, as shown in figure 3, being the BPSK that embodiment is given in embodiment （Binary phase shift keying）And QPSK（Quaternary phase-shift keying (PSK)）Mapping and demapping process schematic；

For BPSK modulation system, in Fig. 3（a）It is the signal constellation and mapping made a start, in Fig. 3（b）It is receiving end step S01 In carry out the planisphere after 7 bit quantizations, in Fig. 3（c）It is the planisphere of further non-uniform quantizing demapping, non-homogeneous amount Change mode is as shown in table 1 below.Wherein, the first hurdle corresponds to the 7bit quantized value of step S01, and the second hurdle is corresponding soft bit, This quantification manner has carried out big step size quantization in the high region of confidence level, and the low region of confidence level has carried out little step size quantization, Thus improve the bit error rate of demapping.

Table 1

For QPSK modulation system, in Fig. 3（d）It is the signal constellation and mapping made a start, using in Fig. 3（b）Quantization star Seat figure, carries out non-uniform quantizing respectively to first bit and second bit, similar with above-mentioned BPSK quantization method, in Fig. 3 's（e）It is the planisphere of first bit demapping, corresponding to table 2（a）, in Fig. 3（f）It is the star of second bit demapping Seat figure, corresponding to table 2（b）.

Table 2（a）

Table 2（b）

Above-mentioned demapping mode is a preferred embodiment, and this non-uniform quantizing mode can improve the property of demapping Energy；Can also higher order modulation system, or other demapping mode.

Step S03, based on the phase quantization value of input sample point, calculates phase deviation, as shown in figure 4, being to combine step The schematic diagram of S02 embodiment.

For BPSK modulation system, phase deviation be pi/2 to the maximum, the quantized value corresponding to step S01 is 32, and phase place is inclined Difference calculating process such as (a) in Fig. 4～（E) shown in, 7bit quantized result [0,127], add 32 first, coordinate axess rotate pi/2, obtain Arrive [0,127], then more than mould（MOD）64, obtain [0,63], after deducting 32 [- 32,31], obtaining after taking absolute value [0, 32], as phase noise estimated value.

For QPSK modulation system, phase deviation be π/4 to the maximum, the quantized value corresponding to step S01 is 16, and phase place is inclined Difference calculating process such as (f) in Fig. 4～（J) shown in, quantized result [0,127], first+16, coordinate axess rotate π/4, obtain [0, 127], then more than mould（MOD）32, obtain [0,31], after deducting 32 [- 16,15], obtain [0,16], that is, after taking absolute value For phase noise estimated value.

If in addition, indicating that this sampled point is big impulse disturbances, phase noise is provided directly as deviation maximum, keeps away Exempt to disturb the phase-detection error causing.For other modulation systems, the process calculating phase noise is similar to.

Step S04, obtains phase noise according to step S03, counts phase noise in time domain and frequency domain respectively, process is illustrated Figure is as shown in Figure 5.

In the present embodiment, the symbol numbers of blind estimate are M, and effective subcarrier number is N, each being obtained by step S03 The phase noise of sampled point is P_noise (i, j), and wherein i represents symbol sequence number, and j represents subcarrier sequence number, can obtain：

Wherein, 0≤i≤M-1,0≤j≤N-1.

Wherein, as shown in fig. 6, step S1 specifically includes：

Step S10, calculates the amplitude of input signal, exact value is real part（I）With imaginary part（Q）Quadratic sum arithmetic square RootIn concrete system, it is simplified operation, possible treating excess syndrome portion absolute value or imaginary part absolute value or both coefficients Weighting, for example, if the point after constellation mapping concentrates on real axis or the imaginary axis, can a treating excess syndrome portion absolute value or imaginary part absolute value As amplitude；

Step S11, according to the amplitude range of input signal, uniform quantization amplitude, quantization method is similar with step S01, this It is quantified as 7bit in embodiment, if this sampled point has big interference, be directly quantized into 0, amplitude range can be according to mean power Obtain；

Step S12, based on quantization amplitude, respectively in time domain and frequency domain statistical average amplitude fading, method and step S04 class Seemingly, frequency domain amplitude fading average Mag_c and temporal amplitude decay average Mag_s are obtained.

Wherein, step S2, is to detect whether this sampled point is to there are big impulse disturbances, the power of sampled point is dry with default Disturb thresholding to compare, if exporting an interference indication signal more than thresholding.The programmable setting of default interference threshold is it is preferable that can It is set to the 10dB of average energy value.

Wherein, as shown in fig. 7, step S3 specifically includes：

Step S30, phase noise is processed, and deducts current phase noise with the maximum that phase noise occurs and obtains one newly Phase noise expression values；Phase noise P_noise_c/P_noise_s value is bigger, and channel quality is more bad, and Mag_c/ Mag_s value is bigger, and channel quality is better, represents the value of phase noise with Phase_c/Phase_s, and this value is bigger, channel quality Better, then：

Phase_c=（max(P_noise_c)-P_noise_c）,

Phase_s=（max(P_noise_s)-P_noise_s）.

In the present embodiment, for BPSK modulation system, Phase_c=（32-P_noise_c）, Phase_s=（32-P_ noise_s）；For QPSK modulation system, Phase_c=（16-P_noise_c）, Phase_s=（16-P_noise_s）.

Step S31, amplitude fading is processed, and this step is optional, in order to keep Phase_c/Phase_s and Mag_c/ Mag_s has identical bit wide, carries out cut position or increases position process, amplitude fading is normalized into and the expression of new phase noise Value has identical bit wide；

Step S32, weighting, the phase noise expression values that step S30 is obtained are multiplied by weight coefficient α, step S31 are obtained Amplitude fading be multiplied by weight coefficient β, acquired results carry out being added the estimated value obtaining channel quality；Weighting expression formula is as follows：

Channel_c(j)=α*Phase_c(j)+β*Mag_c(j)；

Channel_s(i)=α*Phase_s(i)+β*Mag_s(i)；

Wherein, 0≤i≤M-1,0≤j≤N-1；

Thus obtaining the estimated value of channel quality（Channel_c/Channel_s）, the span of weight coefficient α and β For [0,1], and meet alpha+beta=1, in concrete system, different degrees of to receipt signal according to phase noise and amplitude fading Impact, weight coefficient is adjusted flexibly, such as, in differential system, receipt signal is mainly affected by phase noise, a kind of simple Processing mode be, arrange α=1, β=0.

Wherein, step S4, according to the estimated value of channel quality（Channel_c/Channel_s）, based on preset channel door Limit value, counts subcarrier quality, and the frequency that time-domain symbol pulse occurs is estimated overall channel quality average, and ruled out modulation Mode, selects effective subcarrier number, estimates sub-band decay and calculates transmission power, these channel self-adaptings letter Cease to making a start.

There are three kinds of optional modulation systems in embodiment system, be robust-BPSK respectively（Data message repeated before mapping N times）, BPSK, and QPSK, preset two carrier channels thresholding Th_c1, Th_c2, and Th_c1<Th_c2, value is less, and phase place is made an uproar Sound is bigger, characterizes phase noise grade with this two values；Two symbol thresholding Th_s1, Th_s2, and Th_s1<Th_s2, Value is less, and amplitude fading is bigger, characterizes amplitude fading grade with this two values；Step S4 specifically includes：

Step S40, by Channel_s (i) respectively with Th_s1, Th_s2 compares, if Channel_s (i) is less than Th_s1, Cnt_s1 adds 1, if Channel_s (i) is less than Th_s2, Cnt_s2 adds 1；Wherein, Cnt_s is used for counting different amplitude fadings etc. The symbol numbers of level, Cnt_s1 counts the symbol numbers less than thresholding Th_s1, and Cnt_s2 counts the symbol less than thresholding Th_s2 Number.

Step S41, by Channel_c (j) respectively with Th_c1, Th_c2 compares, if Channel_c (j) is more than Th_c1, Cnt_c1 adds 1, carrier_map1（j）Put 1, otherwise set to 0；If Channel_c (j) is more than Th_c2, Cnt_c2 adds 1, carrier_map2（j）Put 1, otherwise set to 0;Wherein, Cnt_c is used for counting the variable number of out of phase noise grade, Cnt_ C1 counts the symbol numbers more than thresholding Th_c1, and Cnt_c2 counts the symbol numbers more than thresholding Th_c2;Carrier_map table Lattice are used for recording the effective or invalid of the subcarrier of out of phase noise grade, and carrier_map1 corresponds under BPSK pattern Effectively subcarrier maps table, carrier_map2 corresponds to the effective subcarrier maps form under QPSK mode, puts this son of 1 expression Carrier wave is effectively usable, set to 0 expression this subcarrier invalid without；

Step S43, reference table 3 selects to meet condition modulation system, if being unsatisfactory for the requirement of form, modulation system selects Shandong Rod-BPSK, if meet BPSK and QPSK modulation simultaneously, compares Cnt_c1, the size of Cnt_c2, if meet Cnt_c1>2*Cnt_ C2, selects BPSK modulation, otherwise selects QPSK.

Table 3

In condition stub in above table 3, the number of parameter and big I are adjusted according to specific system, here only Give one group of enforcement parameter.

Step S44 is it is determined that after modulation system, the carrier_map corresponding to step S41 is effective subcarrier sequence number number Group, if being chosen as robust-BPSK, carrier_map is default setting；If being chosen as BPSK, effective subcarrier sequence number array is Carrier_map1, if being chosen as QPSK, effective subcarrier sequence number array is carrier_map2；Channel_c corresponds to son The channel quality of carrier wave, be worth more big corresponding frequency band quality better it is also possible to according to this value calculate subcarrier packet matter Amount, i.e. the channel quality of sub-band, subband power partition coefficient can be calculated based on this value；

Certainly other channel self-adapting information can also be calculated based on Channel_s and Channel_c, such as encode Mode, interleaving mode etc., are not limited to the present embodiment.

Wherein, as shown in figure 8, step S5 specifically includes：

Step S50, based on the modulation system selecting, is mapped to data in corresponding constellation point, if channel quality compares Good, a constellation point adopts many bit maps, and on the contrary, if channel quality is bad, a constellation point adopts few bit map, On the premise of ensureing transmission quality, realize the maximization of channel resource；

Step S51, according to effective subcarrier sequence number array carrier_map information, only in the good subcarrier of channel quality（Or Sub-band）Upper placement valid data, improve the reliability of transmission quality；

Step S52, according to the channel quality Channel_c information of subcarrier, frequency domain increases sub-band（Or subcarrier）Work( Rate, the power on sub-carrier carries out reasonable distribution, on the premise of keeping frequency domain total transmission power constant, reduces quality preferable Sub-band on power, increase the power on second-rate sub-band, realize the performance of system by reasonable distribution power Raising；

Step S53, according to Channel_avg average and selected modulation system, time domain adjusts transmission power, is protecting In the case of card performance, select little transmission power as far as possible, reduce power consumption.

Refer to Fig. 9, the present invention two the blind Channel self-reacting device based on blind channel self-adaptation method, including：Meter Calculate phase angle quantization modules 01, De-mapping module 02, phase noise estimation module 03, calculating amplitude quantization modules 04, width Degree decay behavior module 05, interference detection module 06, weighting block 07, channel self-adapting module 08 and self adaptation sending module 09；Calculate phase angle and quantization modules 01 connect De-mapping module 02 and phase noise estimation module 03；Phase noise estimates mould Block 03 connects weighting block 07；Calculate amplitude and quantization modules 04 connect interference detection module 06 and amplitude fading estimation module 05；Amplitude fading estimation module 05 connects weighting block 07；Interference detection module 06 connects phase noise estimation module 03 and width Degree decay behavior module 05；Weighting block 07 connecting channel adaptation module 08；Channel self-adapting module 08 connects self adaptation and sends out Send module 09, wherein：

Calculate phase angle and quantization modules 01 be used for calculating the phase angle of input signal, and even amount is melted into some bits, Pass to De-mapping module 02 and phase noise estimation module 03；As shown in Figure 10, calculate phase angle and quantization modules 01 include：Phase Parallactic angle computing unit 011 and phase quantization unit 012, wherein：

Phase angle computing unit 011 is used for realizing rectangular coordinate to polar conversion, uses CORDIC（Coordinate Rotation Digital Computational methods）Algorithm realize, obtain the value of the phase angle of input signal, span be [0,2 π)；

The value even amount of 012 phase angle of phase quantization unit is melted into some bits, by the first multiplier 0121, first Divider 0122 and round unit 0123 and form, wherein：

Input phase angle and 2 realized by first multiplier 0121ⁿ（N is quantization bit）Multiplication, also can be by moving to right n-bit Realize multiplication；Corresponding to embodiment of the method step S01, n=7；

First divider 0122 is used for the output of the first multiplier 0121 divided by 2 π, makes output normalize to [0,2ⁿ]；

Round unit 0123 and take out integer-bit output, output valve is [0,2^n-1] between integer.

De-mapping module 02, according to phase angle quantized value, realizes the inverse process of mapping, to follow-up by non-uniform quantizing Decoder module, De-mapping module 02 includes form memory element and lookup unit, wherein：

Form memory element, for storing phase angle to the index of mapped bits；In conjunction with above-mentioned embodiment of the method, if adjusting Mode processed is robust-BPSK or BPSK, look-up table 1, if modulation system QPSK, look-up table 2（a）With table 2（b）；In concrete system, Typically there is plurality of optional modulation system, often will store multiple mapping tables；

Lookup unit, for realizing phase angle to the Index process of mapped bits.

Phase noise estimation module 03, based on phase offset on planisphere for the phase angle quantized value, on symbol and carry Statistical average phase noise on ripple, as shown in figure 11, phase noise estimation module 03 includes：

Phase deviation computing unit 031, for calculating deviation on planisphere for the phase place, by the first addition being sequentially connected in series More than device 0311, mould, unit 0312, the first subtractor 0313 and the unit 0314 that takes absolute value form, wherein：

First adder 0311 realizes the signal and 2 of its inputⁿ¹Addition, n1 be quantization bit, associated methods embodiment, For robust-BPSK or BPSK, n1=5, for QPSK, n1=4；

More than mould, unit 0312 realizes taking the remainder of the signal of its input；

The signal that first subtractor 0313 realizes its input deducts 2ⁿ¹；

Take absolute value its signal normalization of inputting of unit 0314 to 2ⁿ¹In the range of.

Phase deviation computing unit 031 can associated methods embodiment step S03 implement.

Phase noise statistic unit 032, for statistical average phase noise on symbol and on carrier wave respectively, obtains symbol Phase noise and Carrier Phase Noise, as shown in figure 11, are made an uproar by the first accumulator 0321, the second accumulator 0322, carrier phase The acoustic memory 0323, symbol phase noise memory 0324, the second divider 0325 and the 3rd divider 0326 form, wherein：

First accumulator 0321 is used for the carrier phase of current sign Carrier Phase Noise and corresponding symbol before that adds up Noise；

Second accumulator 0322 is used for all Carrier Phase Noises of current sign that add up；

Carrier Phase Noise memorizer 0323 is used for storing the cumulative of each symbol carrier phase noise and length is equal to and carries Ripple number；

Symbol phase noise memory 0324 is used for storing the cumulative of each symbols phase noise and length is equal to symbol Number number；

Second divider 0325 realizes the cumulative of Carrier Phase Noise and the division with symbol numbers M, obtains carrier phase The average of noise；

3rd divider 0326 realizes the cumulative of symbol phase noise and the division with variable number N, obtains symbol phase The average of noise.

Phase noise statistic unit 032 can associated methods embodiment step S04 implement.

Calculate amplitude quantization modules 04, for calculating the amplitude of input signal, and even amount is melted into some bits, bag Include：

Calculate amplitude unit, for calculating the amplitude of input signal, the exact value of amplitude is real part（I）With imaginary part（Q）'s The arithmetic square root of quadratic sumThen by squarer, adder and extraction of square root unit form, for simplifying fortune this unit Calculate, the possible treating excess syndrome portion absolute value of amplitude or imaginary part absolute value or both coefficients weight, and now this unit includes taking absolute value Unit, multiplier and adder.

Amplitude quantization unit, if for uniform amplitude is quantized into kilo byte, with the phase quantization unit 012 in Figure 10 Composition is similar, including multiplier, divider with round unit.

Amplitude fading estimation module 05, based on amplitude quantizing value, statistical average amplitude fading on symbol and on carrier wave, with The composition of the phase noise statistic unit 032 in Figure 11 is identical.

Interference detection module 06, the amplitude of comparison signal and the size of default interference threshold value, during more than thresholding, export dry Disturb indication signal, including a comparator, comparator is a subtractor, realizes the subtraction of interference threshold and amplitude.

Weighting block 07, according to weight coefficient, weighting processes phase noise and amplitude attenuation, obtains channel quality estimation Value, as shown in figure 12, including：

Phase noise processing unit 071, is individually subtracted symbol phase noise and carrier phase with the maximum of phase noise Noise, obtains phase noise processing costs, realizes this module by a subtractor；

Amplitude fading processing unit 072, cut position or increasing position process amplitude attenuation value so as to have with phase noise processing costs Identical bit wide, obtains amplitude fading processing costs；

Second multiplier 073, realizes the multiplication of weight coefficient α and the processing costs of phase noise；

3rd multiplier 074, realizes the multiplication of weight coefficient β and amplitude fading processing costs；

Second adder 075, realizes the addition of the phase noise processing costs after multiplication process and amplitude fading processing costs, obtains To symbol quality evaluation value and carrier channel quality evaluation value；

Weighting block 07 can associated methods embodiment step S3 implement.

Channel self-adapting module 08, the estimated value of the channel quality according to input, based on preset channel thresholding, is compared Relatively, count and judgement obtains channel self-adapting information, obtaining channel self-adapting information to self adaptation sending module, channel is adaptive Module 08 is answered to include：

Comparator, is used for comparing symbol channel quality estimates and carrier channel quality evaluation value and corresponding predetermined symbol Channel thresholding and the size of default carrier channel thresholding；

Enumerator, for metering less than the symbol numbers of predetermined symbol channel thresholding with more than default carrier channel thresholding Variable number；

Mean module, for calculating the average of carrier channel quality evaluation value；

Modulation system determining unit, for determining modulation system based on count value and average；

Effectively subcarrier determining unit, based on a determination that modulation system, determine corresponding effective subcarrier sequence number；

Calculate sub-band gain coefficient elements, for carrier channel quality evaluation value is converted into corresponding sub-band gain Coefficient；

Calculate transmission power coefficient elements, based on a determination that modulation system and channel quality estimation average, calculate transmitting work( The coefficient of rate.

Step S4 in associated methods embodiment, channel self-adapting module 08 includes two comparator Cmp_1/Cmp_2, four Enumerator Cnt_s1/Cnt_s2/Cnt_c1/Cnt_c2, two subcarrier mapping table carrier_map1/carrier_map2 with And mean module.

Comparator Cmp_1 is used for comparing the size of Channel_s (i) and Th_s1 and Th_s2；If Channel_s (i) is little In Th_s1, enumerator Cnt_s1 adds 1, if Channel_s (i) is less than Th_s2, enumerator Cnt_s2 adds 1；

Comparator Cmp_2 is used for comparing Channel_c（J) and Th_c1, and the size of Th_c2, if Channel_c (j) is big In Th_c1, enumerator Cnt_c1 adds 1, carrier_map1（j）Put 1, otherwise set to 0；If Channel_c (j) is more than Th_c2, meter Number device Cnt_c2 adds 1, carrier_map2（j）Put 1, otherwise set to 0；

What mean module was used for calculating Channel_c (j) is all worth to Channel_avg.

Modulation system determining unit is used for selecting suitable modulation system, reference table 3, if being unsatisfactory for the requirement of form, adjusts Mode processed selects robust-BPSK, if meet BPSK and QPSK modulation simultaneously, compares Cnt_c1, the size of Cnt_c2, if meeting Cnt_c1>2*Cnt_c2, selects BPSK modulation, otherwise selects QPSK.

Effectively subcarrier determining unit, based on the determination of modulation system, selects corresponding subcarrier maps table carrier_ Map1/carrier_map2 or acquiescence carrier_map.

Calculate sub-band gain coefficient elements, the value of Channel_c is directly proportional to the channel quality of subcarrier, calculate The inverse of Channel_c, is encoded into sub-band gain coefficient according to gain step size；

Calculate transmission power coefficient elements, the value of Channel_avg is proportional with total channel quality, calculate The inverse of Channel_avg, and it is encoded into transmission power gain coefficient according to certain step-length.

Self adaptation sending module 09, for the channel self-adapting information according to input, by adjusting the transmission parameter of data, Realize channel self-adapting, including：

Modulation system select unit, the modulation system being come according to receiving end feedback, select corresponding modulation module, data It is mapped in corresponding constellation point；

Effectively sub-carrier allocation unit, feeds back, according to receiving end, effective subcarrier sequence number that comes（carrier_map）, only having The subcarrier transmitting data of effect, other carrier waves place dummy data or zero setting, in addition, this unit can also be a son frequency Band select unit；

Subband power allocation unit, the sub-band fading gain coefficient being come according to receiving end feedback, input sub-band letter Number it is multiplied by the adjustment that corresponding gain coefficient realizes subband power, realized with multiplier；

Transmission power adjustment unit, according to the receiving end transmission power coefficient that comes of feedback it is ensured that less than emission maximum On the premise of power, input signal is multiplied by the adjustment that transmission power coefficient realizes transmission power, is realized with multiplier.

Using MATLAB instrument, the present embodiment is emulated, Figure 13（a）It is 15 multipath channels adopting in embodiment Amplitude-versus-frequency curve；In Figure 13（a）Channel circumstance under, amplitude-versus-frequency curve such as Figure 13 of receipt signal（b）Shown.

Phase noise characteristic curve such as Figure 13 of receipt signal（c）Shown, from Figure 13（b）And Figure 13（c）It can be seen that connecing The amplitude of the collection of letters number can reflect the change of multipath channel well, and phase noise reflection is very rough；Wherein, frequency corresponds to Subcarrier sequence number in diagram.

Figure 14 is the average amplitude attenuation curve of receipt signal under Gaussian channel environment in embodiment and average phase is made an uproar Acoustic curve it can be seen that increase with signal to noise ratio, average amplitude decay value changes less, the change of average phase noise substantially, So, average phase decay can reflect the size of signal to noise ratio well.

As Figure 15（a）Shown, it is the time domain waveform of the receipt signal of emulation and the spectrum waveform in embodiment, simulated environment Based on Figure 13（a）Multi-path channel environment, the white Gaussian noise of superposition 0dB, superposition two big time domain impulse interference（impus_ T1 and impus_t2）With superposition frequency-domain impulse interference on fixing subcarrier（impus_f）, concrete such as Figure 15（a）Shown, letter Road environment is very severe.

Based on Figure 15（a）, Figure 15（b）It is the Symbol average amplitude fading of emulation in embodiment（Mag_s in above-mentioned）'s Waveform and the decay of carrier wave average amplitude（Mag_c in above-mentioned）Waveform it can be seen that Mag_s can reflect two big time domains Interference on symbol for the pulse（Detect_t1 and detect_t2）, Mag_c can reflect the variation tendency of multipath channel, be based on Interference Detection, can reflect the subcarrier sequence number more than interference threshold（At detect_f）.

Based on Figure 15（a）, Figure 15（c）It is the symbol phase noise of emulation in embodiment（P_noise_s in above-mentioned）'s Waveform and Carrier Phase Noise（P_noise_c in above-mentioned）Figure it can be seen that value on each symbol for the P_noise_s Basically identical it is impossible to reflect interference on symbol for the big time domain impulse, this is because sending signal has carried out time-domain difference, P_noise_c can clearly reflect the subcarrier sequence number more than interference threshold（At detect_f）.

In sum, count phase noise respectively and amplitude fading can reflect exactly using on symbol and on carrier wave The change of multipath and the change of signal to noise ratio, can recognize that the big pulse of time domain and the big pulse of frequency domain, based on the weighting of both, can be very Realize channel self-adapting well.

Above example is used for illustrative purposes only, rather than limitation of the present invention, about the technology people of technical field Member, without departing from the spirit and scope of the present invention, can also make various conversion or modification, therefore all equivalent Technical scheme also should belong to scope of the invention, should be limited by each claim.