A symbol rate calculation method and device for a phase amplitude modulation mode

Abstract

The application relates to a symbol rate calculation method for a phase amplitude modulation mode, which uses a density clustering and statistical method to calculate the symbol rate according to the self characteristics of a radio phase amplitude modulation mode signal; specifically, a signal forward difference sequence of a phase is subjected to statistical analysis to find a suspected symbol position threshold in a distance sequence, suspected symbols whose positions are greater than the suspected symbol position threshold are found, the suspected symbols are counted and an interval sequence is reconstructed, then clustering calculation is adopted, and finally the symbol rate is obtained. A device for the symbol rate calculation method for the phase amplitude modulation mode and a storage medium are also disclosed. The application has the beneficial effects that even if noise exists in the signal, the calculation accuracy of the symbol rate is not affected, and the application has strong practicability and robustness.

Description

Technical Field

[0001] This invention relates to the field of communication technology, and in particular to a method and apparatus for calculating symbol rate for phase amplitude modulation, as well as a corresponding storage medium. Background Technology

[0002] Since its invention, radio communication has played a vital role in numerous fields and is now inextricably linked to human society. An increasing number of scientists and engineers are dedicating themselves to the research and optimization of radio communication technology. Due to the inherent characteristics of radio waves, transmission is affected by factors such as weather, obstacles, and electromagnetic fields, causing signal attenuation or interference. These factors severely impact radio signal analysis. During radio signal demodulation, the symbol rate needs to be calculated with high precision, and noise directly affects the accuracy of this calculation. Therefore, in practical computing scenarios, a high signal-to-noise ratio (SNR) is generally required for radio signals.

[0003] In recent years, machine learning technology has been applied to more and more industries, achieving excellent results; combining machine learning with radio signal technology has become a trend. In radio signal technology, symbol rate calculation plays a crucial role in demodulation, and how to improve the accuracy of symbol rate and how to accurately calculate the symbol number under noisy conditions has become an increasingly popular research topic. Summary of the Invention

[0004] The purpose of this invention is to overcome the shortcomings of the prior art and provide a symbol rate calculation method for phase amplitude modulation, which will not affect the accuracy of symbol rate calculation even if there is noise in the signal, and has strong practicality and robustness.

[0005] The objective of this invention is achieved through the following technical solution: a method for calculating the symbol rate of a phase-amplitude modulation scheme, comprising:

[0006] a. After processing the input signal IQ data, a forward differential sequence of the signal with respect to phase is obtained;

[0007] b. Locate the phase positions in the forward differential sequence of the signal that are most concentrated;

[0008] Calculate the distance between each value in the forward difference sequence of the signal and the point where the distribution is most concentrated, and form a distance sequence;

[0009] Calculate the threshold for the location of suspected symbols in the distance sequence;

[0010] c. Find out the suspected symbols whose positions are greater than the suspected symbol position threshold value from the signal forward difference sequence, count the total number of the found suspected symbols, and mark it as L;

[0011] d. Divide the number L into multiple levels with different quantities;

[0012] Sort the data of each level, calculate the interval distance between the positions of two adjacent suspected symbols in each level, and put all of them into an interval sequence;

[0013] e. Obtain the reference interval by performing reference interval analysis on the interval sequence;

[0014] f. Calculate the sampling points occupied by a single symbol according to the suspected symbol position and the reference interval;

[0015] g. The ratio of the original frequency sampling rate to the number of sampling points occupied by a single symbol is the symbol rate of the phase amplitude modulation mode; the original frequency sampling rate refers to the frequency sampling rate of the input signal in a.

[0016] It should be noted that the signal IQ data refers to providing the collected signal to the user in the form of signal IQ, and the user obtains a series of IQ sample points of the signal.

[0017] In b, the suspected symbol position distance threshold value is estimated by statistical analysis of the distance sequence. There are various estimation methods; it can be the estimation method of wave peak and wave trough; it can be the estimation analysis of probability. Analysis and estimation are prior art in this field, and will not be described in detail here.

[0018] Here, the so-called suspected symbol position distance threshold value is the full name of the suspected symbol position distance phase difference threshold value.

[0019] Understanding of c: It can be compared to that there are many points in the X, Y coordinate system, draw a line in the coordinate system, find out the points above the line, and then count the number of found points. Here, the point is the suspected symbol, and the drawn line is the suspected symbol position threshold value; find out the suspected symbol whose suspected symbol position is greater than the suspected symbol position threshold value, that is, find out the point above the suspected symbol position threshold value, and then count it.

[0020] Understanding of d: Essentially, it is to divide the number L into levels with different quantities, then process each level, and finally put it into an interval sequence. When multiple levels are put into an interval sequence, there is overlapping data between levels; when put into an interval sequence, the data base is expanded. The repeated data is more accurate data; by increasing the number of accurate data, the influence of inaccurate data is reduced.

[0021] Further, in a, first, the input signal IQ data is subjected to phase analysis calculation, and then a phase data sequence is obtained; second, the signal phase sequence is subjected to forward difference processing, and finally a signal forward difference sequence is obtained.

[0022] Further preferably, in a, the signal forward difference sequence is analyzed by using density clustering to obtain a phase position where the signal phase forward difference distribution is most concentrated.

[0023] Further, in d, the following steps are performed:

[0024] The number L is divided into three levels of L / 4, L / 2 and L;

[0025] From the total number of L suspected symbols, the top L / 4, L / 2 and L data in phase difference ranking are found out, and three groups of data are constructed.

[0026] Further preferably, the finding out of the top L / 4, L / 2 and L data in phase difference ranking from the total number of L suspected symbols means:

[0027] From the total number of L suspected symbols, the L / 4, L / 2 and L data with the largest distance between the suspected symbol position and the distance threshold value are selected.

[0028] Further, in e, the following steps are performed:

[0029] e1, performing minimum length cleaning on the interval sequence, and removing all interval values with a length less than 4 from the sequence;

[0030] e2, performing statistics on the cleaned interval sequence, and putting intervals with a frequency greater than N times into a suspected interval set;

[0031] e3, performing similar multiple interval merging and statistics on the values in the suspected interval set, and when the following formula is satisfied, the suspected interval frequency is incremented by one;

[0032]

[0033] Wherein, gap is the suspected interval value, x is the value in the cleaned interval sequence, and round represents rounding to the nearest integer;

[0034] Finally, the most frequent suspected interval is found, which is the reference interval, and gap1 is assigned as the reference interval;

[0035] e4, filtering intervals that satisfy the above formula from the basic interval to form a calculation interval.

[0036] It should be noted that the number of N is determined according to the required accuracy, and is not limited, for example, it can be 3, it can be 4 and other.

[0037] Further preferably, in the f, the following conditions are met:

[0038] According to the calculation interval, the reference gap, the sampling points occupied by a single symbol are calculated by the following calculation formula

[0039]

[0040] Wherein, g n is the interval value in the calculation interval set, m is the number of elements in the interval set, gap1 is the reference interval, f(g n ) is the multiple of the reference interval gap1 occupied by the calculation interval.

[0041] A symbol rate calculation device for a phase amplitude modulation mode, comprising:

[0042] A memory for storing a computer program;

[0043] A processor for executing the steps of a symbol rate calculation method for a phase amplitude modulation mode when the computer program is executed.

[0044] A computer readable storage medium:

[0045] Including a program, the program can be executed by the processor to support the completion of a symbol rate calculation method for a phase amplitude modulation mode.

[0046] The present application has the following advantages:

[0047] (1) The symbol rate calculation method of the present application has high accuracy;

[0048] The traditional phase amplitude modulation mode also analyzes and processes the phase forward difference sequence of the signal, but the sequence contains noise, which often leads to the accuracy of the symbol rate calculation cannot be guaranteed;

[0049] According to the signal characteristics of the phase amplitude modulation mode, the signal parameters are calculated by the overall distribution law of the signal through hierarchical statistics and clustering, even if there is noise in the signal, the accuracy of the signal can be guaranteed;

[0050] (2) The symbol rate calculation method of the present application has high accuracy;

[0051] The traditional phase amplitude modulation mode will have a precision problem when calculating the symbol rate due to the introduction of noise;

[0052] The present application is in the processing of data, by dividing into a plurality of different levels in d and constructing three groups of data, and putting the constructed three groups of data in an interval sequence, because three groups of data in an interval sequence, let the correct, useful data is repeated, and the noise data is put in only once, thereby increasing the amount of correct data, and then after analyzing and processing the data, the noise information is not involved in the final analysis, ensuring the final calculation accuracy of the symbol rate;

[0053] (3) The essence of the present application in the process of signal analysis and processing, i.e. e, f, g, is to use signal data priori and clustering method; The effect is that even if there is noise and partial interference in the signal, the symbol rate can still be calculated correctly, the longer the duration of the signal, the higher the accuracy of the symbol rate calculation, and it has strong robustness. BRIEF DESCRIPTION OF DRAWINGS

[0054] Figure 1 The statistical method of the present application is shown in the whole flow chart

[0055] Figure 2 The statistical method of the present application is shown in the whole flow chart

[0056] Figure 3 The signal IQ data forward difference result effect diagram of the embodiment of the present application is shown in the figure

[0057] Figure 4 The signal difference sequence of the embodiment of the present application is shown in the effect diagram

[0058] Figure 5 The threshold value screening position extraction effect diagram of the embodiment of the present application is shown in the figure

[0059] Figure 6 The basic interval analysis flow chart of the embodiment of the present application is shown in the figure DETAILED DESCRIPTION

[0060] The present application will be further described below in conjunction with the drawings, but the protection scope of the present application is not limited to the following description.

[0061] A symbol rate calculation method for phase amplitude modulation mode, the whole flow chart is shown in Figure 1 、 Figure 6 As shown, it comprises:

[0062] a. First, input signal IQ data for phase analysis and calculation, get phase data sequence;

[0063] Then the signal phase sequence is processed by forward difference, and the signal forward difference sequence is obtained as shown in Figure 3 ;

[0064] b. Analyzing the signal forward difference sequence in a way of using density clustering to obtain the phase position where the signal phase forward difference distribution is most concentrated;

[0065] Obtaining the distance value of each value in the signal forward difference sequence from the most concentrated position of the distribution, and forming a distance sequence;

[0066] After statistical analysis of the distance sequence, as shown in Figure 3 ;

[0067] Finally, estimate the suspected symbol position distance threshold;

[0068] Among them,

[0069] There are many ways to estimate, such as peak and trough estimation; It can be a probabilistic estimation analysis; Analysis and estimation are prior art in this field, and will not be described in detail here;

[0070] It should be noted that,

[0071] c. From the signal forward difference sequence, find the suspected symbol whose position is greater than the suspected symbol position threshold, and the result after finding is shown in Figure 5 ; Count the total number of suspected symbols found and record it as L;

[0072] Analogous to,

[0073] There are many points in the X, Y coordinate system, draw a line in the coordinate system, find the points above the line, and then count the number of points found. Here the point is the suspected symbol, and the line drawn is the suspected symbol position threshold; Find the suspected symbol whose position is greater than the suspected symbol position threshold, that is, find the points above the suspected symbol position threshold, and then count them;

[0074] d. Divide the number L into L / 4, L / 2, and L levels;

[0075] From the total number of suspected symbols L, find the top L / 4, L / 2, and L data in the phase difference, and construct three groups of data;

[0076] Sort the data in each level, calculate the interval distance between the positions of the adjacent two suspected symbols in each level, and put them all into an interval sequence;

[0077] Specifically, it means,

[0078] Sort the total number of suspected symbols L, select the L / 4, L / 2, and L data with the largest distance between the suspected symbol position and the suspected symbol position distance threshold, and divide them into three groups;

[0079] For each group of data, the interval distance between two adjacent suspected symbol positions in each group of data is calculated, and interval distances corresponding to the three groups of data are put into an interval sequence;

[0080] e. The reference interval is obtained by performing reference interval analysis on the interval sequence;

[0081] The reference interval analysis method is shown in the following formula: Figure 6

[0082] e1. The minimum length cleaning is performed on the interval sequence, and all interval values with a length less than 4 are removed from the sequence;

[0083] e2. The cleaned interval sequence is counted, and intervals with a frequency greater than N times are put into a suspected interval set;

[0084] e3. The values in the suspected interval set are combined and counted according to similar multiple intervals, and when the following formula is satisfied, the frequency of the suspected interval is incremented by one,

[0085]

[0086] Wherein, gap is the suspected interval value, x is the value in the cleaned interval sequence, and round represents rounding to the nearest integer;

[0087] Finally, the most frequent suspected interval is found, which is the reference interval, and gap1 is assigned as the reference interval;

[0088] e4. The intervals that satisfy the above formula are selected to form a calculation interval sequence;

[0089] f. According to the calculation interval sequence and the reference interval, the sampling points occupied by a single symbol are calculated;

[0090] Specifically, the calculation of the sampling points occupied by a single symbol is calculated according to the following formula,

[0091]

[0092] Wherein, g n is the interval value in the calculation interval sequence, m is the number of elements in the interval set, gap1 is the reference interval, and f(g n ) is the multiple of the reference interval gap1 occupied by the calculation interval;

[0093] g. The ratio of the original frequency sampling rate to the number of sampling points occupied by a single symbol is the symbol rate of the phase amplitude modulation mode. The original frequency sampling rate refers to the frequency sampling rate of the input signal in a.

[0094] Optionally, the application also provides a symbol rate calculation device for a phase amplitude modulation mode.​

[0095] Specifically, the application comprises a memory for storing a computer program; and a processor for executing the computer program to implement steps of a symbol rate calculation method for a phase amplitude modulation mode.

[0096] Optionally, the application also provides a computer readable storage medium.

[0097] Specifically, the application comprises a program, which can be executed by a processor to support completion of a symbol rate calculation method for a phase amplitude modulation mode.

[0098] The above embodiments only express the more preferred embodiments, which are described in more detail and in more detail, but cannot be understood as limiting the patent scope of the application. It should be noted that for those skilled in the art, without departing from the concept of the present application, a number of modifications and improvements can be made, which are within the scope of protection of the present application.

Claims

1. A method for calculating the symbol rate for a phase amplitude modulation scheme, characterized in that: The method comprises the following steps: a. After processing the input signal IQ data, a signal forward difference sequence about phase is obtained; b. The most concentrated phase position in the signal forward difference sequence is found; The distance value of each value in the signal forward difference sequence from the most concentrated position is calculated, and a distance sequence is formed; The threshold of the suspected symbol position in the distance sequence is calculated; it can be analogized as that a plurality of points are distributed in an X-Y coordinate system, a line is drawn in the coordinate system, and the points above the line are found; the points in the coordinate system are the suspected symbols, and the drawn line is the threshold of the suspected symbol position; c. The suspected symbols whose positions are greater than the threshold of the suspected symbol position are found from the signal forward difference sequence, the total number of the found suspected symbols is counted, and the total number is recorded as L; the suspected symbols whose positions are greater than the threshold of the suspected symbol position are found, that is, the points above the threshold of the suspected symbol position are found, and then counted; d. The number L is divided into a plurality of levels with different numbers; The data of each level is sorted, the interval distance between two adjacent suspected symbol positions in each level is calculated, and all the interval distances are put into an interval sequence; e. The reference interval is obtained by performing reference interval analysis on the interval sequence; In e: e1. The interval sequence is cleaned by the minimum length, and all interval values with a length less than 4 are removed from the sequence; e2. The cleaned interval sequence is counted, and the intervals with an occurrence frequency greater than N times are put into a suspected interval set; e3. The values in the suspected interval set are combined and counted according to similar multiple intervals, and when the following formula is satisfied, the occurrence frequency of the suspected interval is increased by one, Wherein, gap is the suspected interval value, x is the value in the cleaned interval sequence, and round represents rounding off; Finally, the most frequent suspected interval is found, which is the reference interval, and gap1 is assigned as the reference interval; e4. The intervals satisfying the above formula are screened out, and a calculation interval sequence is formed; f. The sampling points occupied by a single symbol are calculated according to the suspected symbol position and the reference interval; The ratio of the original frequency sampling rate to the number of sampling points occupied by a single symbol is the symbol rate of the phase amplitude modulation mode; the original frequency sampling rate refers to the frequency sampling rate of the input signal in a.

2. The method according to claim 1, wherein the method is for a phase amplitude modulation scheme. In a: The phase data sequence is obtained by performing phase analysis and calculation on the input signal IQ data; The signal forward difference sequence is obtained by performing forward difference processing on the signal phase sequence.

3. The method of claim 2, wherein the method is for a phase amplitude modulation scheme. In b: the most concentrated phase position of the signal phase forward difference is obtained by using density clustering to analyze the signal forward difference sequence.

4. The method of claim 1, wherein the method is used for a phase amplitude modulation scheme. In d: The number L is divided into three levels of L / 4, L / 2 and L; From the total number L of suspected symbols, the data with the top L / 4, L / 2 and L in the phase difference ranking are found, and three groups of data are constructed.

5. The symbol rate calculation method for phase amplitude modulation according to claim 4, characterized in that: The data with the top L / 4, L / 2 and L in the phase difference ranking are found from the total number L of suspected symbols, which means that From the total number of L suspected symbols, select the suspected symbol position, and the distance between the suspected symbol position and the distance threshold value, the maximum L / 4 data, L / 2 data, L data.

6. The method of claim 1, wherein the method is characterized by: The f is: Single symbol occupies sampling point calculation is calculated by the formula as shown below wherein, to calculate the gap value in the gap sequence, m is the number of elements in the gap set, gap1 is the reference gap, to calculate the multiple of the basic gap gap1 occupied by the gap.

7. A device for calculating a symbol rate for a phase amplitude modulation scheme, characterized by: Including: Memory, for storing computer programs; Processor, for executing computer programs to implement the steps of the symbol rate calculation method for phase amplitude modulation mode according to any one of claims 1 to 6.

8. A computer readable storage medium, characterized in that: Including programs, the programs can be executed by the processor to support the completion of the symbol rate calculation method for phase amplitude modulation mode according to any one of claims 1 to 6.

Images