High-speed aircraft channel detection communication method and device based on combined ZC sequence

Abstract

The invention discloses a high-speed aircraft channel detection communication method and device based on a combined ZC sequence, and the method comprises the following steps: S1, selecting E ZC sequences which have the sequence length of Nzc and meet a non-periodic self-correlation characteristic, and mapping information bits into a root value sequence Q; s2, generating a combined ZC sequence of multiple root values; s3, performing sliding correlation on the SEQ sequence and the E ZC sequences in the S1, recording the position and amplitude phase of a correlation peak, obtaining a position sequence, and extracting multipath time delay; the invention discloses a multi-path time delay detection communication method, which comprises the following steps of S1, obtaining position sequences of which the position difference is integer multiples of Nzc, S2, obtaining position sequences of which the position difference is integer multiples of Nzc, S4, combining the position sequences of which the position difference is integer multiples of Nzc to obtain a root value sequence Q'of the SEQ sequence, and S5, obtaining phases and amplitude sequences of each path according to amplitude phase information, calculating frequency offset and amplitude, and completing detection communication.Under the condition that the signal bandwidth is not reduced, the detection range of multi-path time delay is expanded, and the detection contradiction between fast time variation of a channel and multi-path time delay difference is solved.

Description

Technical Field

[0001] This invention relates to the field of high-speed aircraft telemetry, tracking, and communication, and specifically to a high-speed aircraft channel detection and communication method and apparatus based on combined ZC sequences. Background Technology

[0002] The basic idea behind wireless channel measurement is to send a probe signal of a certain form at the transmitting end, then process the received signal accordingly, and finally obtain channel state information through channel estimation. Wireless channel measurement can be divided into two categories: time-domain measurement and frequency-domain measurement.

[0003] Frequency domain measurement involves the transmitter progressively changing the frequency of a probe signal within a certain range, and then using a vector network analyzer to process the signals at the transmitting and receiving ends to obtain the frequency response at discrete frequency points. This response is the frequency domain expression of the channel impulse response. Frequency sweep measurement requires precise time positioning and synchronization, and the vector network analyzer requires a physical connection; therefore, it is mostly used for indoor channel measurements.

[0004] Common time-domain measurement methods include direct pulse measurement and spread-spectrum sliding correlation. Direct pulse measurement involves repeatedly inputting an extremely narrow pulse signal into the channel, allowing the channel measurement receiver to directly obtain the time-domain channel impulse response without complex signal processing. However, ideal pulse signals are difficult to achieve; their high peak-to-average power ratio and interference sensitivity limit the application of direct pulse measurement, making it popular only in the early days.

[0005] Spread-spectrum sliding correlation (SCR) methods utilize the correlation characteristics of pseudo-noise (PN) sequences or constant-envelope zero-autocorrelation (CAZAC) sequences to obtain the channel impulse response. The transmitted probe signal, composed of the PN or CAZAC sequence, carries information characterizing the channel after propagation through the channel. Using the same local sequence as the probe signal, a sliding correlation is performed with the received signal. When the received sequence and the transmitted sequence are aligned, a correlation peak is generated. Analysis of the correlation peak allows for the extraction of parameters such as Doppler and multipath delay. Spread-spectrum sliding correlation measurement offers processing gain and effectively eliminates interference, making it the most commonly used channel measurement method currently.

[0006] In channel measurement methods, some have proposed using ZC sequences with a single root value to achieve multipath separation detection of the integrated channel. However, this method does not resolve the contradiction between the rapid time-varying nature of the channel and the detection of multiple path delay differences. Other methods have proposed a spread-spectrum sliding correlation detection scheme based on a combination of ZC and Frank sequences, which uses dual sequences to address the contradiction between the rapid time-varying nature of the channel and the detection of multiple path delay differences. However, this scheme does not achieve communication functionality and requires a large bandwidth for the detection signal. Therefore, this paper proposes a channel detection and communication method and device for high-speed aircraft based on combined ZC sequences. Summary of the Invention

[0007] The purpose of this invention is to provide a high-speed aircraft channel detection communication method and apparatus based on combined ZC sequences, so as to overcome the shortcomings of existing technologies in detecting the contradiction between fast time-varying channels and multipath delay differences.

[0008] To achieve the above objectives, the present invention provides the following technical solutions: In a first aspect, the present invention provides a high-speed aircraft channel sounding communication method based on combined ZC sequences, comprising the following steps: S1, Select a sequence with length N zc E ZC sequences that satisfy the aperiodic autocorrelation property are mapped to the root value sequence Q. S2, Generate a ZC sequence of multiple root values ​​based on the root value sequence Q and name it the SEQ sequence; S3, the SEQ sequence is slide correlated with the ZC sequence in step S1, and the position of the correlation peak and the amplitude and phase information of the correlation peak are recorded after the slide correlation. The correlation peak position sequence is obtained based on the correlation peak position and the amplitude and phase information of the correlation peak, and the multipath delay of the correlation peak position sequence is extracted. S4, the merge position difference is the sequence length N. zc The sequence of related peak positions is an integer multiple of the sequence of related peak positions. Based on the root value information of the first path correlation peak positions, the root value sequence Q' of the SEQ sequence in step S3 is obtained. S5. Based on the amplitude and phase information of each relevant peak recorded in step S3, obtain the phase sequence and amplitude sequence of each path in step S4, calculate the frequency offset and amplitude of each path, and complete the detection communication.

[0009] Furthermore, the specific process of S1 is as follows: Select a sequence length of N zc The E ZC sequences that satisfy the aperiodic autocorrelation property are E ZC sequences whose main peak is much larger than the secondary peak of the aperiodic autocorrelation function, and the recorded root values ​​are q1, q2, ..., q. E ; The information bits required for communication are mapped to a root value sequence Q=[Q1,Q2,…,Q]. M ], where Q m=q1, q2, ..., q E m = 1, 2, ..., M, where M is the number of roots in the root sequence Q.

[0010] Furthermore, the specific process of S2 is as follows: Based on the root value sequence Q in step S1, generate a combination ZC sequence of multiple root values ​​and name it the SEQ sequence, SEQ=[ZC1(n),ZC2(n),…,ZC…]. M [(n)], where Q m =q1, q2, ..., q E m = 1, 2, ...; M, where M is the number of roots in the root sequence Q; Send the SEQ sequence to the receiving end.

[0011] Furthermore, the specific process of S3 is as follows: The receiving end receives the SEQ sequence; The E different implanted ZC sequences selected in step S1 are used as local sequences and slide correlated with the SEQ sequences respectively, and the positions and amplitude and phase information of all correlation peaks are recorded. The correlation peak positions obtained by sliding correlation of the local ZC sequence for each root value are merged in ascending order to form the correlation peak position sequence D. qe D qe =[D qe (1), D qe (2),…, D qe (T e )] indicates that the root value is q e The sequence of related peak positions; where e = 1, 2, ..., E; T e D represents the number of relevant peaks in the sequence. qe (t e ) represents the t-th element in the sequence. e The positions of the relevant peaks, t e =1,2,…,T e ;q e It is a variable, q E It is q e The maximum value, D qE The root value is q E The sequence of related peak positions; Extract the multipath delay of the relevant peak position sequence.

[0012] Furthermore, the step of extracting multipath delay based on the relevant peak position sequence specifically includes: S31, note D qe (1) is D qe The position of the first relevant peak in D qe Extract all elements related to Dqe (1) The position difference is the sequence length N zc Integer multiples of D qe (1) Forming position sequence D qe 1 D qe 1 =[D qe 1 (1), D qe 1 (2), ..., D qe 1 (B)], where B represents D qe 1 The number of relevant peak positions in D qe 1 (b) indicates D qe 1 The position of the b-th related peak in the data, where b = 1, 2, ..., B; S32, find D qe 1 The positions of the two adjacent related peaks with the largest position difference are D qe 1 (b') and D qe 1 (b'+1), where b'=1, 2, ..., B; D qe 1 (b') and D qe 1 (b'+1) represents D qe 1 The positions of the two adjacent related peaks with the largest position difference; S33, Extract D qe Located in D qe 1 (b') and D qe 1 The sequence of related peak positions between (b'+1) belonging to the same path, D qe r =D qe 2 , ...D qe R , where D qe r The root value is q e In the relevant peak position sequence, it belongs to the r-th path and is located in D qe 1 (b') and D qe 1 For any sequence of related peak positions between (b'+1), r=2,…,R; S34, extract D qe 2,…, D qe R Position of the relevant peak at the end of the middle D qe r (end) =D qe 2 (end), ..., D qe R (end), D qe r (end) indicates D qe r The position of the relevant peak at the end of the middle is used to calculate D. qe r (end)-D qe 1 (b') can then be used to obtain the time delay difference between the r-th path and the first path; S35, calculate the path frequency offset using the phase sequence of each path correlation peak, using the following formula: ; Among them, F r Let W represent the frequency offset of the r-th path, and W be the bandwidth of the transmitted signal. Indicates to Perform a floor operation; Among them, P r =[P r (1), P r (2),…P r (m),P r [M] represents the phase sequence of the correlation peak of the r-th path, where r = 1, 2, ..., R.

[0013] Furthermore, the specific process of S4 is as follows: The D recorded in step S3 is merged q1 D q2 ,…,D qE The positional difference is the sequence length N. zc Multiples of D1, D2, ..., D2 are obtained. R D r =[D r (1), D r (2),…,D r [M] represents the sequence of correlation peak positions for the r-th path, r = 1, 2, ..., R; D r (m)=D qe (t e ) indicates that the position of the m-th correlation peak in this sequence originates from D. qe The tth e There are several related peak positions, m = 1, 2, ..., M; The root value information of each relevant peak position in the first path merging result D1 yields the root value sequence Q'=[Q1',Q2',…,Q] of the ZC sequence in the received signal. M ').

[0014] Furthermore, the specific process of S5 is as follows: Based on the amplitude and phase information of each related peak recorded in step S3, the phase sequence P of related peaks belonging to the same path is merged in step S4. r =P1, P2, ..., P R and amplitude sequence A r =A1, A2, ..., A R ;where P r =[P r (1), P r (2),…,P r [(M)] represents the correlation peak phase sequence of the r-th (r=1,2,…,R) path, where P r =[P r (1), P r (2),…P r (m),P r [M] represents the phase sequence of the correlation peak of the r-th path, r=1,2,…,R, P r (m) represents the phase of the m-th correlation peak in the sequence, m = 1, 2, ..., M; A r =[A r (1), A r (2), …A r (m), A r [M] represents the correlation peak amplitude sequence of the r-th path, r = 1, 2, ..., R; A r (m) represents the amplitude of the m-th related peak in the sequence, m=1,2,…,M; R represents that there are R different paths in total; The frequency offset and amplitude of each path are calculated using the phase sequence and amplitude sequence of the correlation peaks, thus completing the detection and communication.

[0015] Secondly, the present invention provides a high-speed vehicle channel sounding communication system based on combined ZC sequences, which applies the aforementioned high-speed vehicle channel sounding communication method based on combined ZC sequences, including: Mapping module: Select sequence length N zc E ZC sequences that satisfy the aperiodic autocorrelation property are mapped to the root value sequence Q. Sending module: Generates a ZC sequence of multiple root values ​​based on the root value sequence Q and names it as the SEQ sequence; Extraction module: The SEQ sequence is correlated with the ZC sequence in step S1, and the position and amplitude and phase information of the correlation peak after the correlation are recorded. The correlation peak position sequence is obtained based on the correlation peak position and amplitude and phase information of the correlation peak, and the multipath delay of the correlation peak position sequence is extracted. Merge module: Merge position difference is sequence length N zc The sequence of related peak positions that are integer multiples of the first path correlation peak position is used to obtain the root value sequence Q' of the SEQ sequence in the extraction module. Calculation module: Based on the amplitude and phase information of each relevant peak recorded in the extraction module, the phase sequence and amplitude sequence of each path in the merging module are obtained, the frequency offset and amplitude of each path are calculated, and the detection communication is completed.

[0016] Thirdly, the present invention provides a computer device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the method described above.

[0017] Fourthly, a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the method described above.

[0018] Compared with the prior art, the present invention has the following beneficial technical effects: This invention provides a high-speed aircraft channel sounding communication method based on combined ZC sequences. Utilizing the good correlation characteristics between ZC sequences with different root values, it expands the detection range of multipath delay without reducing signal bandwidth, thus meeting the detection requirements for fast-changing channels and resolving the contradiction between fast-changing channels and multipath delay differences. By using the same combined ZC sequence to simultaneously complete information transmission and estimation of multipath delay, frequency offset, and amplitude, it eliminates the need to send additional pilot sequences, saving spectrum resources. By utilizing the correlation peak position characteristics of ZC sequences with different root values ​​and merging correlation peak position sequences of the same path, multipath delay information can be accurately extracted. In particular, by distinguishing correlation peaks whose position difference is an integer multiple of the sequence length, it effectively avoids mutual confusion between multipaths and improves the accuracy of delay estimation in complex multipath environments. Attached Figure Description

[0019] Figure 1 This is a schematic diagram of a high-speed aircraft channel sounding communication method based on a combined ZC sequence in an embodiment of the present invention.

[0020] Figure 2 This is a schematic diagram of the sliding correlation results of the local ZC sequence with a root value of 1 in an embodiment of the present invention.

[0021] Figure 3This is a schematic diagram of the sliding correlation results of the local ZC sequence with a root value of 119 in an embodiment of the present invention.

[0022] Figure 4 This is a schematic diagram of the sliding correlation results of the local ZC sequence with a root value of 120 in an embodiment of the present invention.

[0023] Figure 5 This is a schematic diagram of the sliding correlation results of the local ZC sequence with a root value of 238 in an embodiment of the present invention.

[0024] Figure 6 This is a schematic diagram illustrating the merging of correlation peaks belonging to path one in the sliding correlation results of all local sequences in an embodiment of the present invention.

[0025] Figure 7 This is a schematic diagram illustrating the merging of correlation peaks belonging to path two in the sliding correlation results of all local sequences in an embodiment of the present invention.

[0026] Figure 8 This is a schematic diagram of the integrated detection and communication method based on the ZC sequence of multiple root values ​​in an embodiment of the present invention. Detailed Implementation

[0027] In the following description, only certain exemplary embodiments are briefly described. As those skilled in the art will recognize, the described embodiments can be modified in various ways without departing from the spirit or scope of the invention. Therefore, the drawings and description are considered to be exemplary in nature and not restrictive.

[0028] In the description of this invention, it should be understood that the terms "center," "longitudinal," "lateral," "length," "width," "thickness," "upper," "lower," "front," "rear," "left," "right," "vertical," "horizontal," "top," "bottom," "inner," "outer," "clockwise," "counterclockwise," "axial," "radial," and "circumferential" indicate the orientation or positional relationship based on the orientation or positional relationship shown in the accompanying drawings. They are used only for the convenience of describing this invention and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, or be constructed and operated in a specific orientation. Therefore, they should not be construed as limitations on this invention.

[0029] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of this invention, "a plurality of" means two or more, unless otherwise explicitly specified.

[0030] In this invention, unless otherwise explicitly specified and limited, the terms "installation," "connection," "linking," and "fixing," etc., should be interpreted broadly. For example, they can refer to a fixed connection, a detachable connection, or an integral part; they can refer to a mechanical connection, an electrical connection, or a communication connection; they can refer to a direct connection or an indirect connection through an intermediate medium; they can refer to the internal communication of two components or the interaction between two components. Those skilled in the art can understand the specific meaning of the above terms in this invention according to the specific circumstances.

[0031] In this invention, unless otherwise explicitly specified and limited, "above" or "below" the second feature can include direct contact between the first and second features, or contact between the first and second features through another feature between them. Furthermore, "above," "over," and "on top" of the second feature includes the first feature being directly above or diagonally above the second feature, or simply indicates that the first feature is at a higher horizontal level than the second feature. "Below," "below," and "under" the second feature includes the first feature being directly above or diagonally above the second feature, or simply indicates that the first feature is at a lower horizontal level than the second feature.

[0032] It should be understood that, when used in this specification and the appended claims, the terms "comprising" and "including" indicate the presence of the described features, integrals, steps, operations, elements and / or components, but do not exclude the presence or addition of one or more other features, integrals, steps, operations, elements, components and / or collections thereof.

[0033] It should also be understood that the terminology used in this specification is for the purpose of describing particular embodiments only and is not intended to limit the invention. As used in this specification and the appended claims, the singular forms “a,” “an,” and “the” are intended to include the plural forms unless the context clearly indicates otherwise.

[0034] It should also be further understood that the term "and / or" as used in this specification and the appended claims refers to any combination of one or more of the associated listed items and all possible combinations, and includes such combinations.

[0035] The accompanying drawings illustrate various structural schematic diagrams according to embodiments disclosed in this invention. These drawings are not to scale, and some details have been enlarged for clarity, and some details may have been omitted. The shapes of the various regions and layers shown in the drawings, as well as their relative sizes and positional relationships, are merely exemplary and may deviate from reality due to manufacturing tolerances or technical limitations. Furthermore, those skilled in the art can design regions / layers with different shapes, sizes, and relative positions as needed.

[0036] Example 1 See Figure 1 and Figure 8 This invention provides a high-speed aircraft channel sounding communication method based on combined ZC sequences, comprising the following steps: S1, Select a sequence with length N zc E ZC sequences that satisfy the aperiodic autocorrelation property are mapped to the root value sequence Q. In one specific embodiment, the sequence length is selected as N. zc The E ZC sequences that satisfy the aperiodic autocorrelation property are E ZC sequences whose main peak is much larger than the secondary peak of the aperiodic autocorrelation function, and the recorded root values ​​are q1, q2, ..., q. E ; The information bits required for communication are mapped to a root value sequence Q=[Q1,Q2,…,Q]. M ], where Q m =q1, q2, ..., q E m = 1, 2, ..., M, where M is the number of roots in the root sequence Q.

[0037] S2, Generate a ZC sequence of multiple root values ​​based on the root value sequence Q and name it the SEQ sequence, and send the SEQ sequence to the receiving end; In one specific embodiment, a combination ZC sequence of multiple root values ​​is generated based on the root value sequence Q in step S1 and named the SEQ sequence, where SEQ = [ZC1(n), ZC2(n), ..., ZC...]. M [(n)], where Q m =q1, q2, ..., q E m = 1, 2, ...; M, where M is the number of roots in the root sequence Q; Send the SEQ sequence to the receiving end.

[0038] S3, the SEQ sequence is slide correlated with the ZC sequence in step S1, and the position of the correlation peak and the amplitude and phase information of the correlation peak are recorded after the slide correlation. The correlation peak position sequence is obtained based on the correlation peak position and the amplitude and phase information of the correlation peak, and the multipath delay of the correlation peak position sequence is extracted. In one optional embodiment, the receiving end receives the SEQ sequence; The E different implanted ZC sequences selected in step S1 are used as local sequences and slide correlated with the SEQ sequences respectively, and the positions and amplitude and phase information of all correlation peaks are recorded. The correlation peak positions obtained by sliding correlation of the local ZC sequence for each root value are merged in ascending order to form the correlation peak position sequence D. qe D qe =[Dqe (1), D qe (2),…, D qe (T e )] indicates that the root value is q e The sequence of related peak positions; where e = 1, 2, ..., E; T e D represents the number of relevant peaks in the sequence. qe (t e ) represents the t-th element in the sequence. e The positions of the relevant peaks, t e =1,2,…,T e D qE D qe The maximum value that can be obtained; Extract the multipath delay of the relevant peak position sequence.

[0039] In one specific embodiment, the step of extracting multipath delay based on the relevant peak position sequence is as follows: S31, note D qe (1) is D qe The position of the first relevant peak in D qe Extract all elements related to D qe (1) The position difference is the sequence length N zc Integer multiples of D qe (1) Forming position sequence D qe 1 D qe 1 =[D qe 1 (1), D qe 1 (2), ..., D qe 1 (B)], where B represents D qe 1 The number of relevant peak positions in D qe 1 (b) indicates D qe 1 The position of the b-th related peak in the data, where b = 1, 2, ..., B; S32, find D qe 1 The positions of the two adjacent related peaks with the largest position difference are D qe 1 (b') and D qe 1 (b'+1), where b'=1, 2, ..., B; D qe 1 (b') and D qe 1 (b'+1) represents Dqe 1 The positions of the two adjacent related peaks with the largest position difference; S33, Extract D qe Located in D qe 1 (b') and D qe 1 The sequence of related peak positions between (b'+1) belonging to the same path, D qe r =D qe 2 ...D qe R , where D qe r The root value is q e In the relevant peak position sequence, it belongs to the r-th path and is located in D qe 1 (b') and D qe 1 For any sequence of related peak positions between (b'+1), r=2,…,R; S34, extract D qe 2 ,…, D qe R Position of the relevant peak at the end of the middle D qe r (end) =D qe 2 (end), ..., D qe R (end), D qe r (end) indicates D qe r The position of the relevant peak at the end of the middle is used to calculate D. qe r (end)-D qe 1 (b') can then be used to obtain the time delay difference between the r-th path and the first path; S35, calculate the path frequency offset using the phase sequence of each path correlation peak, using the following formula: ; Among them, F r Let W represent the frequency offset of the r-th path, and W be the bandwidth of the transmitted signal. Indicates to Perform a floor operation; Among them, P r =[P r (1), P r (2),…P r (m),P r[M] represents the phase sequence of the correlation peak of the r-th path, where r = 1, 2, ..., R.

[0040] S4, the merge position difference is the sequence length N. zc The sequence of related peak positions is an integer multiple of the sequence of related peak positions. Based on the root value information of the first path correlation peak positions, the root value sequence Q' of the SEQ sequence in step S3 is obtained. In one specific embodiment, the D recorded in step S3 is merged. q1 D q2 ,…,D qE The positional difference is the sequence length N. zc Multiples of D1, D2, ..., D2 are obtained. R D r =[D r (1), D r (2),…,D r [M] represents the sequence of correlation peak positions for the r-th path, r = 1, 2, ..., R; D r (m)=D qe (t e ) indicates that the position of the m-th correlation peak in this sequence originates from D. qe The tth e There are several related peak positions, m = 1, 2, ..., M; The root value information of each relevant peak position in the first path merging result D1 yields the root value sequence Q'=[Q1',Q2',…,Q] of the ZC sequence in the received signal. M ').

[0041] S5. Based on the amplitude and phase information of each relevant peak recorded in step S3, obtain the phase sequence and amplitude sequence of each path in step S4, calculate the frequency offset and amplitude of each path, and complete the detection communication.

[0042] In one specific embodiment, based on the amplitude and phase information of each related peak recorded in step S3, the phase sequence P of related peaks belonging to the same path merged in step S4 is... r =P1, P2, ..., P R and amplitude sequence A r =A1, A2, ..., A R ;where P r =[P r (1), P r (2),…P r (m),P r [M] represents the phase sequence of the correlation peak of the r-th path, r=1,2,…,R; P r (m) represents the phase of the m-th correlation peak in the sequence, m = 1, 2, ..., M; A r =[A r(1), A r (2), …A r (m), A r [M] represents the correlation peak amplitude sequence of the r-th path, r = 1, 2, ..., R; A r (m) represents the amplitude of the m-th related peak in the sequence, m=1,2,…,M; R represents that there are R different paths in total; The frequency offset and amplitude of each path are calculated using the phase sequence and amplitude sequence of the correlation peaks, thus completing the detection and communication.

[0043] Example 2 See Figures 1 to 8 The present invention provides an embodiment of the method, which is mainly implemented through the following steps: S1: Select 4 sequences with length N zc The ZC sequence with a main peak much larger than the secondary peak of the aperiodic autocorrelation function is 239, with root values ​​q1=1, q2=119, q3=120, and q4=238. The information bits to be transmitted are mapped to the four ZC sequences with different root values ​​selected in S1, resulting in a root value sequence Q. The root value sequence of the combined ZC sequences is: Q=[1,238,120,1,238,120,119,119,238,1,120,119,1,238,119,1,238,120,1,119]; S2: Generate a ZC sequence of multiple root values ​​based on the root value sequence Q and name it the SEQ sequence. Send the SEQ sequence to the receiving end. In this embodiment, the transmit signal bandwidth is set to 23.9MHz, the channel is configured as a two-path channel, and the delay difference is 840×23.9e6. -1 The first path frequency offset is 20kHz, the second path frequency offset is 10kHz, the first path amplitude is 1, the second path amplitude is 0.5, and the fixed signal-to-noise ratio at the receiver is 10dB.

[0044] S3: Four different rooted ZC sequences were selected as local sequences and subjected to sliding correlation with the SEQ sequences. The root values ​​of the four different rooted ZC sequences were 1, 119, 120, and 238, respectively. The results are as follows: Figures 2 to 5 As shown, the correlation peaks extracted by q1 are q1(1)~q1(12), the correlation peaks extracted by q2 are q2(1)~q2(10), the correlation peaks extracted by q3 are q3(1)~q3(8), and the correlation peaks extracted by q4 are q4(1)~q4(10). The position information and amplitude and phase information of all correlation peaks are recorded. The correlation peak positions obtained by sliding correlation of the local ZC sequence of each root value are merged into a correlation peak position sequence in ascending order to obtain D. q1=[239,956,1079,1796,2390,3107,3230,3824,3947,4541,4664,5381]; D q2 =[1673,1912,2513,2752,2868,3585,3708,4425,4780,5620]; D q3 =[717,1434,1557,2274,2629,3469,4302,5142];D q4 =[478,1195,1318,2035,2151,2991,3346,4063,4186,4903]; D q3 The relevant peak position sequence belonging to the first path, D q3 1 D in =[717,1434,2629,4302] q3 1 (3) and D q3 1 (4) The positional difference is greatest at D. q3 Located in D q3 1 (3) and D q3 1 (4) The correlation peak position is 3469. The time delay difference between the second path and the first path is calculated as (3469-2629)×23.9e6. -1 s = 840 × 23.9e6 -1 s.

[0045] S4: Combine the peak positions of each path, such as Figure 6 and Figure 7 As shown, from Figure 4 First path merging result: D1=[D q1 (1),D q4 (1),D q3 (1),D q1 (2),D q4 (2),D q3 (2),D q2 (1),D q2 (2),D q4 (5),D q1 (5),D q3 (5),D q2 (5),D q1 (6),D q4 (7),D q2 (6),D q1 (8),D q4 (8),Dq3 (7),D q1 (10),D q2 (9)]; The root value information of the SEQ sequence is obtained as follows: Q'=[1,238,120,1,238,120,119,119,238,1,120,119,1,238,119,1,238,120,1,119]; Communication is complete.

[0046] S5: Obtain the phase sequence of the first path based on the amplitude and phase information corresponding to the position information of each relevant peak recorded in S3: P1=[0.61,1.89,3.11,-1.89,-0.62,0.66,1.9,3.08,-1.9,-0.67,0.58,1.87,-3.13,-1.85,-0.65,0.66,1.87,-3.12,-1.89,-0.65]; And the amplitude sequence of the first path: A1=[0.94,0.95,0.93,0.95,0.97,0.92,0.93,0.95,0.94,0.91,0.95,0.93,0.92,0.91,0.94,0.92,0.95,0.99,0.95,0.91]; Phase sequence of the second path: P2=[0.23,0.85,1.58,2.27,2.81,-2.75,-2.11,-1.38,-1,-0.39,0.28,0.83,1.54,2.15,2.8,-2.78,-2.09,-1.54,-0.91,-0.36]; And the amplitude sequence of the second path: A2 = [0.51, 0.47, 0.5, 0.5, 0.52, 0.53, 0.46, 0.47, 0.42, 0.48, 0.46, 0.48, 0.54, 0.55, 0.61, 0.43, 0.5, 0.51, 0.48, 0.49]. The first path frequency offset is calculated to be 19.997 kHz, the second path frequency offset to be 10.032 kHz, the first path amplitude to be 1.003, and the second path amplitude to be 0.504, thus completing the integrated detection and communication.

[0047] Example 3 This invention provides a high-speed aircraft channel sounding communication system based on combined ZC sequences, applying the aforementioned high-speed aircraft channel sounding communication method based on combined ZC sequences, including: Mapping module: Select sequence length Nzc E ZC sequences that satisfy the aperiodic autocorrelation property are mapped to the root value sequence Q. Sending module: Generates a ZC sequence of multiple root values ​​based on the root value sequence Q and names it as the SEQ sequence; Extraction module: The SEQ sequence is correlated with the ZC sequence in step S1, and the position and amplitude and phase information of the correlation peak after the correlation are recorded. The correlation peak position sequence is obtained based on the correlation peak position and amplitude and phase information of the correlation peak, and the multipath delay of the correlation peak position sequence is extracted. Merge module: Merge position difference is sequence length N zc The sequence of related peak positions that are integer multiples of the first path correlation peak position is used to obtain the root value sequence Q' of the SEQ sequence in the extraction module. Calculation module: Based on the amplitude and phase information of each relevant peak recorded in the extraction module, the phase sequence and amplitude sequence of each path in the merging module are obtained, the frequency offset and amplitude of each path are calculated, and the detection communication is completed.

[0048] In one embodiment of the present invention, a computer device is provided, comprising a processor and a memory. The memory stores a computer program, which includes program instructions. The processor executes the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing and control core of the terminal, suitable for implementing one or more instructions, specifically suitable for loading and executing one or more instructions in the computer storage medium to achieve a corresponding method flow or corresponding function. The processor described in this embodiment of the present invention can be used in the operation of a high-speed aircraft channel sounding communication method based on combined ZC sequences.

[0049] This invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device used to store programs and data. It is understood that the computer-readable storage medium here can include both the built-in storage medium in the computer device and extended storage media supported by the computer device. The computer-readable storage medium provides storage space that stores the operating system of the terminal. Furthermore, the storage space also stores one or more instructions suitable for loading and execution by a processor. These instructions can be one or more computer programs (including program code). It should be noted that the computer-readable storage medium here can be high-speed RAM or non-volatile memory, such as at least one disk storage device. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the corresponding steps of the high-speed aircraft channel sounding communication method based on combined ZC sequences in the above embodiments.

[0050] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0051] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0052] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0053] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0054] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and not to limit it. Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art should understand that modifications or equivalent substitutions can still be made to the specific implementation of the present invention. Any modifications or equivalent substitutions that do not depart from the spirit and scope of the present invention should be covered within the scope of protection of the claims of the present invention.

Claims

1. A high-speed aircraft channel sounding and communication method based on combined ZC sequences, characterized in that, Includes the following steps: S1, Select a sequence with length N zc E ZC sequences that satisfy the aperiodic autocorrelation property are mapped to the root value sequence Q. S2, Generate a ZC sequence of multiple root values ​​based on the root value sequence Q and name it the SEQ sequence; S3, the SEQ sequence is slide correlated with the ZC sequence in step S1, and the position of the correlation peak and the amplitude and phase information of the correlation peak are recorded after the slide correlation. The correlation peak position sequence is obtained based on the correlation peak position and the amplitude and phase information of the correlation peak, and the multipath delay of the correlation peak position sequence is extracted. S4, the merge position difference is the sequence length N. zc The sequence of related peak positions is an integer multiple of the sequence of related peak positions. Based on the root value information of the first path correlation peak positions, the root value sequence Q' of the SEQ sequence in step S3 is obtained. S5. Based on the amplitude and phase information of each relevant peak recorded in step S3, obtain the phase sequence and amplitude sequence of each path in step S4, calculate the frequency offset and amplitude of each path, and complete the detection communication.

2. The high-speed aircraft channel sounding communication method based on combined ZC sequences according to claim 1, characterized in that, The specific process of S1 is as follows: Select a sequence length of N zc The E ZC sequences that satisfy the aperiodic autocorrelation property are E ZC sequences whose main peak is much larger than the secondary peak of the aperiodic autocorrelation function, and the recorded root values ​​are q1, q2, ..., q. E ; The information bits required for communication are mapped to a root value sequence Q=[Q1,Q2,…,Q]. M ], where Q m =q1, q2, ..., q E m = 1, 2, ..., M, where M is the number of roots in the root sequence Q.

3. The high-speed aircraft channel sounding communication method based on combined ZC sequences according to claim 2, characterized in that, The specific process of S2 is as follows: Based on the root value sequence Q in step S1, generate a combination ZC sequence of multiple root values ​​and name it the SEQ sequence, SEQ=[ZC1(n),ZC2(n),…,ZC…]. M [(n)], where Q m =q1, q2, ..., q E m = 1, 2, ...; M, where M is the number of roots in the root sequence Q; Send the SEQ sequence to the receiving end.

4. The high-speed aircraft channel sounding communication method based on combined ZC sequences according to claim 3, characterized in that, The specific process of S3 is as follows: The receiving end receives the SEQ sequence; The E different implanted ZC sequences selected in step S1 are used as local sequences and slide correlated with the SEQ sequences respectively, and the positions and amplitude and phase information of all correlation peaks are recorded. The correlation peak positions obtained by sliding correlation of the local ZC sequence for each root value are merged in ascending order to form the correlation peak position sequence D. qe D qe =[D qe (1), D qe (2),…, D qe (T e )] indicates that the root value is q e The sequence of related peak positions; where e = 1, 2, ..., E; T e D represents the number of relevant peaks in the sequence. qe (t e ) represents the t-th element in the sequence. e The positions of the relevant peaks, t e =1,2,…,T e ;q e It is a variable, q E It is q e The maximum value, D qE The root value is q E The sequence of related peak positions; Extract the multipath delay of the relevant peak position sequence.

5. A high-speed aircraft channel sounding and communication method based on combined ZC sequences according to claim 4, characterized in that, The specific steps for extracting multipath delay based on the relevant peak position sequence are as follows: S31, note D qe (1) is D qe The position of the first relevant peak in D qe Extract all elements related to D qe (1) The position difference is the sequence length N zc Integer multiples of D qe (1) Forming positional sequence D qe 1 D qe 1 =[D qe 1 (1), D qe 1 (2), ..., D qe 1 (B)], where B represents D qe 1 The number of relevant peak positions in D qe 1 (b) indicates D qe 1 The position of the b-th related peak in the data, where b = 1, 2, ..., B; S32, find D qe 1 The positions of the two adjacent related peaks with the largest position difference are D qe 1 (b') and D qe 1 (b'+1), where b'=1, 2, ..., B; D qe 1 (b') and D qe 1 (b'+1) represents D qe 1 The positions of the two adjacent related peaks with the largest position difference; S33, Extract D qe Located in D qe 1 (b') and D qe 1 The sequence of related peak positions between (b'+1) belonging to the same path, D qe r =D qe 2 , ...D qe R , where D qe r The root value is q e In the relevant peak position sequence, it belongs to the r-th path and is located in D qe 1 (b') and D qe 1 For any sequence of related peak positions between (b'+1), r=2,…,R; S34, extract D qe 2 ,…, D qe R Position of the relevant peak at the end of the middle D qe r (end)=D qe 2 (end), ..., D qe R (end), D qe r (end) indicates D qe r The position of the relevant peak at the end of the middle is used to calculate D. qe r (end)-D qe 1 (b') can then be used to obtain the time delay difference between the r-th path and the first path; S35, calculate the path frequency offset using the phase sequence of each path correlation peak, using the following formula: ; Among them, F r Let W represent the frequency offset of the r-th path, and W be the bandwidth of the transmitted signal. Indicates to Perform a floor operation; Among them, P r =[P r (1), P r (2),…P r (m),P r [M] represents the phase sequence of the correlation peak of the r-th path, where r = 1, 2, ..., R.

6. The high-speed aircraft channel sounding communication method based on combined ZC sequences according to claim 5, characterized in that, The specific process of S4 is as follows: The D recorded in step S3 is merged q1 D q2 ,…,D qE The positional difference is the sequence length N. zc Multiples of D1, D2, ..., D2 are obtained. R D r =[D r (1), D r (2),…,D r [M] represents the sequence of correlation peak positions for the r-th path, r = 1, 2, ..., R; D r (m)=D qe (t e ) indicates that the position of the m-th correlation peak in this sequence originates from D. qe The tth e There are several related peak positions, m = 1, 2, ..., M; The root value information of each relevant peak position in the first path merging result D1 yields the root value sequence Q'=[Q1',Q2',…,Q] of the ZC sequence in the received signal. M ').

7. A high-speed aircraft channel sounding communication method based on combined ZC sequences according to claim 6, characterized in that, The specific process of S5 is as follows: Based on the amplitude and phase information of each related peak recorded in step S3, the phase sequence P of related peaks belonging to the same path is merged in step S4. r =P1, P2, ..., P R and amplitude sequence A r =A1, A2, ..., A R ; Among them, P r =[P r (1), P r (2),…P r (m),P r [M] represents the phase sequence of the correlation peak of the r-th path, r=1,2,…,R, P r (m) represents the phase of the m-th correlation peak in the sequence, m = 1, 2, ..., M; A r =[A r (1), A r (2), …A r (m), A r [M] represents the correlation peak amplitude sequence of the r-th path, r = 1, 2, ..., R; A r (m) represents the amplitude of the m-th related peak in the sequence, m=1,2,…,M; R represents that there are R different paths in total; The frequency offset and amplitude of each path are calculated using the phase sequence and amplitude sequence of the correlation peaks, thus completing the detection and communication.

8. A high-speed aircraft channel sounding communication system based on combined ZC sequences, employing the high-speed aircraft channel sounding communication method based on combined ZC sequences as described in any one of claims 1-7, characterized in that, include: Mapping module: Select sequence length N zc E ZC sequences that satisfy the aperiodic autocorrelation property are mapped to the root value sequence Q. Sending module: Generates a ZC sequence of multiple root values ​​based on the root value sequence Q and names it as the SEQ sequence; Extraction module: The SEQ sequence is correlated with the ZC sequence in step S1, and the position and amplitude and phase information of the correlation peak after the correlation are recorded. The correlation peak position sequence is obtained based on the correlation peak position and amplitude and phase information of the correlation peak, and the multipath delay of the correlation peak position sequence is extracted. Merge module: Merge position difference is sequence length N zc The sequence of related peak positions that are integer multiples of the first path correlation peak position is used to obtain the root value sequence Q' of the SEQ sequence in the extraction module. Calculation module: Based on the amplitude and phase information of each relevant peak recorded in the extraction module, the phase sequence and amplitude sequence of each path in the merging module are obtained, the frequency offset and amplitude of each path are calculated, and the detection communication is completed.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the method as described in any one of claims 1-7.

10. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements the steps of the method as described in any one of claims 1-7.

Images