Disclosure of Invention
In order to solve the above-mentioned problems of the prior art, an object of the present invention is to provide a low-complexity method for estimating the packing rate of a super nyquist system, so as to estimate the packing rate of a random super nyquist system and reduce the implementation complexity.
In order to achieve the above object, the present invention provides a low complexity super nyquist system packing rate estimation method, including:
The method comprises the steps of obtaining a transmission pilot frequency block and mapped symbols in a super Nyquist system, constructing a transmission frame based on the transmission pilot frequency block and mapped symbols, up-sampling the transmission frame and forming a baseband to generate a transmission signal, and carrying out noise adding and filtering on the transmission signal to generate a receiving signal;
Acquiring an initial assumed downsampling factor, carrying out iterative updating on the initial assumed downsampling factor, respectively carrying out downsampling on a received signal based on each assumed downsampling factor in the iterative updating process to obtain a downsampled pilot block, and carrying out coherent correlation processing on the downsampled pilot block and a transmission pilot block until the assumed downsampling factor updated last time is larger than an upsampling parameter in the upsampling process to obtain a judgment value corresponding to each assumed downsampling factor;
and extracting an assumed downsampling factor corresponding to the maximum decision value, and generating the packing rate of the super Nyquist system according to the extracted assumed downsampling factor.
Optionally, the construction process of the transmission frame includes:
wherein b k denotes the kth transmission frame, b n denotes the nth symbol in the corresponding transmission frame, p is the transmission pilot block, p= [ p 0,p1,p2,…,pL-1 ], L denotes the length of the transmission pilot block, Representing the length of the mapped symbol, a n is the mapped nth symbol.
Optionally, the generating process of the transmission signal includes:
Where s (T) represents a transmission signal, E s is an average power of a signal in a transmission frame, N is a symbol length in the transmission frame, b n represents an nth symbol in a corresponding transmission frame, N is a reference number of the transmission frame, c (T) is a root-raised cosine shaped pulse of unit energy corresponding to time T, and T is a symbol period of the root-raised cosine shaped pulse.
Optionally, the generating process of the received signal includes:
Where r k represents the kth received symbol in the received signal, g (T) represents the impulse response of the nyquist system transfer function corresponding to time T, E s is the average power of the signal in the transmission frame, N is the symbol length in the transmission frame, b n represents the nth symbol in the corresponding transmission frame, N is the label of the transmission frame, M represents the up-sampling multiple of the nyquist system, T is the symbol period of the root-raised cosine shaped pulse, Q is the up-sampling parameter of the root-raised cosine shaped pulse, and η k is the colored noise.
Optionally, the process of downsampling the received signal includes:
Where y k represents the kth downsampled symbol, g (T) represents the impulse response of the nyquist system transfer function corresponding to time T, E s is the average power of the signal in the transmission frame, N is the symbol length in the transmission frame, b n represents the nth symbol in the corresponding transmission frame, N is the label of the transmission frame, M represents the upsampling multiple of the nyquist system, T is the symbol period of the root-raised cosine shaped pulse, Q is the upsampling parameter of the root-raised cosine shaped pulse, The representation assumes that the downsampling factor is, Is the colored noise after downsampling.
Optionally, the process of performing coherent correlation processing on the pilot block after downsampling and the transmission pilot block includes:
Wherein, The i symbol of the kth pilot block after downsampling is represented, p i is the i symbol of the transmission pilot block, L represents the length of the transmission pilot block, and i represents the symbol index in the transmission pilot block.
Optionally, the generating process of the packing rate of the super nyquist system includes:
where τ is the packing fraction, Q is the upsampling parameter of the root-raised cosine shaped pulse, Is the hypothesized downsampling factor corresponding to the maximum decision value.
Compared with the prior art, the invention has the following advantages and technical effects:
The invention provides a low-complexity super Nyquist system packing rate estimation method by using a pilot frequency-based super Nyquist transmission model, which can convert the packing rate estimation problem into a downsampling factor estimation problem. And the method comprises the steps of carrying out downsampling on the signals after matched filtering by using different assumed downsampling factors, so as to obtain a downsampled pilot block corresponding to each assumed downsampling factor, carrying out coherent correlation operation on the transmission pilot block and the downsampled pilot block, obtaining a judgment value corresponding to the assumed downsampling factors, finding out the assumed downsampling factor corresponding to the maximum judgment value, calculating the packing rate of the corresponding super Nyquist system, and adopting the corresponding packing rate to design a transmission scheme of the corresponding random super Nyquist system, thereby ensuring safe and reliable implementation of signal detection and physical layer and low implementation complexity.
Detailed Description
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer executable instructions, and that although a logical order is illustrated in the flowcharts, in some cases the steps illustrated or described may be performed in an order other than that illustrated herein.
In order to solve the above-mentioned problems in the prior art, an object of the present invention is to provide a low-complexity super nyquist system packing rate estimation method, so as to estimate a random super nyquist system packing rate and reduce implementation complexity.
The invention provides a method for estimating the packing rate of a low-complexity super Nyquist system, which comprises the steps of framing mapped symbols, constructing a transmission frame of the super Nyquist based on pilot frequency, using the super Nyquist transmission model to convert the packing rate estimation problem into a downsampling factor estimation problem, carrying out upsampling and baseband forming on the transmission frame to obtain a transmission signal, carrying out additive Gaussian white noise channel and matched filtering on the transmission signal to describe the transmission process of the transmission signal to obtain a receiving signal, carrying out downsampling on the matched and filtered receiving signal by using different postulated downsampling factors to generate downsampled symbols, extracting the acquired pilot frequency blocks from the downsampled symbols, carrying out coherent correlation operation on the transmission pilot frequency blocks and the downsampled pilot frequency blocks to obtain judgment values corresponding to different postulated downsampling factors, finding the postulated downsampling factors corresponding to the maximum judgment values, and calculating the packing rate of the super Nyquist system according to the found postulated downsampling factors to obtain the estimated packing rate.
As some embodiments, the constructing of the k-th transmission frame based on the pilot frequency comprises constructing the transmission frame based on the transmission pilot frequency block and the mapped symbol:
wherein b k denotes the kth transmission frame, k is the corresponding reference number of the transmission frame, b n denotes the nth symbol in the corresponding transmission frame, p= [ p 0,p1,p2,…,pL-1 ] is the transmission pilot block, L denotes the length of the transmission pilot block, Representing the mapped symbol block length, a n is the mapped nth symbol.
As some embodiments, up-sampling and baseband shaping are performed on the transmission frame and transmitting, and a transmission signal is obtained:
Where s (T) represents a transmission signal, E s is an average power of a transmission frame signal, N is a symbol corresponding to a symbol in the transmission frame, N is a length of the symbol in the transmission frame, c (T) is a root-raised cosine shaped pulse of unit energy corresponding to time T, and T is a symbol period of the root-raised cosine shaped pulse.
As some embodiments, the transmitted signal is subjected to an additive white gaussian noise channel and matched filtering to obtain a received signal:
Where r k represents the kth received symbol in the received signal, g (T) represents the impulse response of the nyquist system transfer function corresponding to time T, E s is the average power of the signal in the transmission frame, N is the symbol length in the transmission frame, b n represents the nth symbol in the corresponding transmission frame, N is the label of the transmission frame, M represents the up-sampling multiple of the nyquist system, T is the symbol period of the root-raised cosine shaped pulse, Q is the up-sampling parameter of the root-raised cosine shaped pulse, and η k is the colored noise.
The specific acquisition process of the received signal is as follows:
wherein r (t) represents the received signal, Representing a linear convolution of the data,Representing the impulse response of the super Nyquist system transfer function corresponding to time t, n (t) representing white noiseRepresenting the matched filtered colored noise.
When the assumed downsampling factor is not greater than Q, i.eDownsampling the received signal using a hypothetical downsampling factor to obtain downsampled symbols:
Where y k represents the kth downsampled symbol, M represents the upsampling multiple of the Nyquist system, Q is the upsampling parameter of the root-raised cosine shaped pulse, Representing an assumed downsampling factor, which takes a value not greater than the upsampling parameter Q of the root-raised cosine shaped pulse,Is the colored noise after downsampling.
Extracting pilot blocks from the corresponding positions in the downsampled symbols to obtain downsampled pilot blocks
Performing coherent correlation operation on the transmission pilot block and the pilot block after downsampling to obtain a decision value corresponding to the assumed downsampling factor:
where Λ represents the decision value corresponding to the hypothesized downsampling factor, Representing the ith symbol of the kth pilot block after downsampling, p i is the ith symbol of the transmission pilot block.
And in a certain range, the numerical value of the hypothesized downsampling factor is continuously adjusted through iterative updating, the operations of downsampling the matched and filtered received signal, acquiring a downsampled pilot block and related operations are sequentially iterated based on each hypothesized downsampling factor adjusted in the iterative updating process, and the generated judgment value in the iterative process is stored until the hypothesized downsampling factor is larger than the upsampling parameter of the root raised cosine shaping pulse, and the iteration is stopped.
As some embodiments, among the stored decision values corresponding to different assumed downsampling factors, finding the assumed downsampling factor corresponding to the maximum decision value, and calculating the packing rate τ of the ultranyquist system:
Wherein, Is the hypothesized downsampling factor corresponding to the maximum decision value.
The invention discloses the technical effects that a pilot frequency-based super Nyquist transmission model is used, and the invention provides a low-complexity super Nyquist system packing rate estimation method which can convert the packing rate estimation problem into a downsampling factor estimation problem. And carrying out downsampling on the signals after matching and filtering by using different assumed downsampling factors, so as to obtain a downsampled pilot block corresponding to each assumed downsampling factor, then carrying out coherent correlation operation on the transmission pilot block and the downsampled pilot block, obtaining a judgment value corresponding to the assumed downsampling factor, finding out the assumed downsampling factor corresponding to the maximum judgment value, thus calculating the packing rate of the super Nyquist system, ensuring safe and reliable implementation of signal detection and physical layer, and having low implementation complexity.
The foregoing is described in detail with reference to the accompanying drawings:
the invention aims to overcome the defects of the prior art, and provides a low-complexity method for estimating the packing rate of a super Nyquist system, so as to obtain the packing rate of the super Nyquist and reduce the implementation complexity.
In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.
The pilot-based super nyquist transmission system adopted by the invention with reference to fig. 1-2 mainly comprises a data source module, a constellation mapping module, a framing module, an up-sampling module, a baseband shaping module, a channel module, a matched filtering module and a packing rate estimation module, wherein:
The data source module is used for generating bit data required to be transmitted by the super Nyquist system and transmitting the bit data to the constellation mapping module;
The constellation mapping module maps the bit data into symbols according to constellation mapping rules and transmits the mapped symbols to the framing module;
The framing module is used for inserting a known transmission pilot frequency block into each mapped symbol block, forming a transmission frame so as to construct the transmission frame based on the transmission pilot frequency block and the mapped symbols, and transmitting the transmission frame to the up-sampling module;
the up-sampling module performs zero value interpolation on the transmission frame to perform up-sampling on the transmission frame and transmits the transmission frame subjected to zero value interpolation to the baseband shaping module;
The baseband forming module is used for performing super Nyquist forming on the up-sampled transmission frame to perform baseband forming on the transmission frame, transmitting the transmission frame after the baseband forming, namely a transmission signal, and transmitting the transmission frame to the channel module;
The channel module is used for adding Gaussian white noise to the transmission frame after the baseband shaping and transmitting the transmission frame after the Gaussian white noise addition to the matched filtering module;
The matched filtering module carries out matched filtering on a received signal which is a transmission frame added with Gaussian white noise, and transmits the filtered received signal to the packing rate estimation module;
And the packing rate estimation module is used for downsampling the filtered received signals by using different assumed downsampling factors to obtain a downsampled pilot block, performing coherent correlation operation on the transmission pilot block and the downsampled pilot block to obtain a judgment value corresponding to the assumed downsampling factors, sequentially iterating the downsampling operation and the correlation operation, finding the assumed downsampling factor corresponding to the maximum judgment value, and calculating the packing rate of the super Nyquist system.
Referring to fig. 3, the steps of the present invention for packet rate estimation using the above-mentioned nyquist system are as follows:
step 1, the construction of the k transmission frame of the super nyquist based on the pilot frequency comprises the steps of transmitting a pilot frequency block and mapped symbols:
wherein b k denotes the kth transmission frame, k is the corresponding reference number of the transmission frame, b n denotes the nth symbol in the corresponding transmission frame, p= [ p 0,p1,p2,…,pL-1 ] is the transmission pilot block, L denotes the length of the transmission pilot block, Representing the mapped symbol block length, a n is the mapped nth symbol.
Step 2, up-sampling and baseband shaping are carried out on the transmission frame, and transmitting signals are obtained:
Where s (T) represents a transmission signal, E s is an average power of a transmission frame signal, N is a symbol corresponding to a symbol in the transmission frame, N is a length of the symbol in the transmission frame, c (T) is a root-raised cosine shaped pulse of unit energy corresponding to time T, and T is a symbol period of the root-raised cosine shaped pulse.
Step 3, the transmitting signal is subjected to an additive Gaussian white noise channel and matched filtering to obtain a receiving signal:
Where r k represents the kth received symbol in the received signal, g (T) represents the impulse response of the nyquist system transfer function corresponding to time T, E s is the average power of the signal in the transmission frame, N is the symbol length in the transmission frame, b n represents the nth symbol in the corresponding transmission frame, N is the label of the transmission frame, M represents the up-sampling multiple of the nyquist system, T is the symbol period of the root-raised cosine shaped pulse, Q is the up-sampling parameter of the root-raised cosine shaped pulse, and η k is the colored noise.
The specific acquisition process of the received signal is as follows:
wherein r (t) represents the received signal, Representing a linear convolution of the data,Representing the impulse response of the super nyquist system transfer function corresponding to time t, n (t) represents white noise,Representing the matched filtered colored noise. When the assumed downsampling factor is not greater than Q, downsampling the matched filtered received signal using the assumed downsampling factor to obtain downsampled symbols:
Where y k represents the kth downsampled symbol, M represents the upsampling multiple of the Nyquist system, Q is the upsampling parameter of the root-raised cosine shaped pulse, Representing an assumed downsampling factor, which takes a value not greater than the upsampling parameter Q of the root-raised cosine shaped pulse,Is the colored noise after downsampling.
Extracting pilot blocks from the corresponding positions in the downsampled symbols to obtain downsampled pilot blocks
Performing coherent correlation operation on the transmission pilot block and the pilot block after downsampling to obtain a decision value corresponding to the assumed downsampling factor:
where Λ represents the decision value corresponding to the hypothesized downsampling factor, Representing the ith symbol of the kth pilot block after downsampling, p i is the ith symbol of the transmission pilot block.
And in a certain range, the numerical value of the hypothesized downsampling factor is continuously adjusted through iterative updating, the operations of downsampling the matched and filtered received signal, acquiring the downsampled pilot block and related operations are sequentially iterated based on each adjusted hypothesized downsampling factor in the iterative updating process, and the generated judgment value in the iterative process is stored until the hypothesized downsampling factor is larger than the upsampling parameter of the root raised cosine shaping pulse, and the iteration is stopped.
It should be noted that, the above-mentioned transmission signal, the received signal and the symbol after downsampling are all simulation signal models for signal transmission, where the simulation signal model corresponding to the symbol after downsampling only includes "assumed downsampling factor" as a variable, under the condition of parameter determination of the nyquist system and transmission frame determination, different assumed downsampling factors can generate the symbol after downsampling correspondingly, where the value of the assumed downsampling factor obtained by each iteration is substituted into the calculation formula of the symbol after downsampling, so as to extract the pilot block after downsampling, further calculate the decision value through correlation operation, obtain several groups of different downsampling factors through several groups of iterations, further obtain several groups of different decision values, select the maximum decision value from them, and finally calculate the packing rate by using the downsampling factor corresponding to the maximum decision value.
In the iterative updating process of the assumed downsampling factor, the numerical value of the assumed downsampling factor of the previous iteration is increased by a fixed step length to perform iterative updating, wherein the initial value and the fixed step length of the initial assumed downsampling factor are set according to the upsampling parameters of the root raised cosine shaping pulse, the initial value of the downsampling factor is smaller than the upsampling parameters of the root raised cosine shaping pulse, and the fixed step length is 1.
Step 4, finding an assumed downsampling factor corresponding to the maximum decision value, and calculating the packing rate tau of the super Nyquist system:
Wherein, Is the hypothesized downsampling factor corresponding to the maximum decision value.
In this embodiment, the effect of this embodiment is further described in connection with a simulation experiment;
1. simulation conditions
The simulation experiments of this example were performed under MATLAB2022B software. In the simulation experiment of this embodiment, BPSK is used as the modulation scheme of the pilot block, and the length of a single pilot block is 32.
The roll-off factor of the root-raised cosine shaped pulse and the matched filter is set to 0.3.
2. Simulation content and result analysis
Simulation 1, under the above conditions, when the slow change of the packing rate is considered, the packing rate is estimated by the present invention, and the result is shown in fig. 4.
Simulation 2, under the above conditions, when the packing rate is considered to be randomly and uniformly distributed in [0.6,0.7,0.8,0.9,1], the packing rate is estimated by using the method, and the result is shown in fig. 5.
The horizontal axis in fig. 4 and 5 represents the bit signal-to-noise ratio of the nyquist system in dB (decibel) and the vertical axis represents the false alarm probability PFA (ProbabilityofFalseAlarm).
As can be seen from fig. 4, when the slow change of the packing rate is considered, the invention has the highest estimation accuracy for the packing rate of 1 at a low bit signal-to-noise ratio, and has the highest estimation accuracy for the packing rate of 0.8 as the bit signal-to-noise ratio increases, and has the lowest estimation accuracy for the packing rate of 0.6 regardless of the bit signal-to-noise ratio. As can be seen from fig. 5, under the same bit signal-to-noise ratio, the false alarm probability when the packing rate is randomly and uniformly distributed is higher than the false alarm probability when the packing rate is slowly changed, and the estimation accuracy is reduced compared with the case when the packing rate is slowly changed. The method can estimate the packing rate of the random super Nyquist system, and ensure the safe and reliable implementation of signal detection and physical layer.
The present application is not limited to the above-mentioned embodiments, and any changes or substitutions that can be easily understood by those skilled in the art within the technical scope of the present application are intended to be included in the scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.