CN105099976B - A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal - Google Patents
A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal Download PDFInfo
- Publication number
- CN105099976B CN105099976B CN201510450571.XA CN201510450571A CN105099976B CN 105099976 B CN105099976 B CN 105099976B CN 201510450571 A CN201510450571 A CN 201510450571A CN 105099976 B CN105099976 B CN 105099976B
- Authority
- CN
- China
- Prior art keywords
- signal
- radar
- communication
- integration
- modulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/10—Frequency-modulated carrier systems, i.e. using frequency-shift keying
- H04L27/103—Chirp modulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
Abstract
A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal, the frequency curve shape of the integration signal is asymmetrical triangle, it is made of two sections of opposite linear FM signal (LFM) signals of tune frequency symbol, wherein first portion's signal is that frequency modulation rate is k1For LFM signals, the communication information is not modulated, for carrying out the synchronization of signal of communication or capture, second portion signal is that frequency modulation rate is k2LFM signals (wherein k1×k2<0), with the form modulation communication code element information of (1,1) code on second portion signal.Based on the parameterized model of asymmetric triangle frequency modulated(FM) radar communicating integral signal, carry out the parameter of optimization design integration signal to characterize the ambiguity function of radar data reduction and characterize the object function of communication performance.Institute's extracting method of the present invention ensure that the probability for the successfully capture that communicates using shortest synchronizing signal, and ensure that on the premise of radar data reduction, ensure the rate of information throughput under certain signal-to-noise ratio.
Description
Technical field
The present invention relates to a kind of parameter optimization methods of asymmetric triangle frequency modulated(FM) radar communicating integral signal, can be applicable in
In the radar-communication integration Waveform Design shared based on signal, belong to radar communication technology interdisciplinary field, more particularly to
Radar-communication integration technology.
Background technology
It is resisted with the development and the information-based needs fought, the form of war of technology from single weapon platform to system pair
Anti-rotation becomes, the comprehensive integrated trend for being implemented as following military electronic systems development of multifunction electronic system.Radar communication
Signal integration has great significance for promoting the comprehensive performance of electronic system.Exploring the design method of integrated waveform is
The primary study content of radar communication signal integrated technique.During radar-communication integration is realized, in order to realize thunder
Up to function and communication function, the mode of generally use time sharing transmissions or multi-beam space division.But both modes do not possess adjustable
Section property, can not can be regarded as real integrated technique, can only be called highly integrated.Only realize Radar-Communication letter signal
The integration of waveform is only real integration.
At present, domestic and international Duo Jia research institutions have carried out the integrated correlative study of radar communication signal, such as the cloth of Sweden
Lai Jin Polytechnics, Xian Electronics Science and Technology University, University of Electronic Science and Technology, early warning institute of air force and the Chinese Academy of Space Technology
Deng the achievement in research of existing literature is mainly the content of the following aspects:1) integration signal of carrier wave is become based on Step Frequency,
The integration signal is mainly used in radar network and realizes communication function, is specially provided with a radar and is used to implement synchronous thunder
It reaches;2) the radar-communication integration signal based on random Step Frequency between Orthogonal Frequency Division Multiplexing (OFDM) arteries and veins, the integration signal lead to
Data message is loaded into radar signal by overfrequency agile, data is transmitted using random step frequency, so as to realize simultaneously
Radar detection and data communication function;3) integration signal based on multicarrier direct sequence signal, the integration signal use OFDM
With the mutual-complementing code multicarrier direct expansion (MC-DSSS) that spread spectrum is combined;4) it is based on Ao Poman sequences (Oppermann
Sequences, OS) weighted pulse burst integration signal, which is weighted train of pulse with OS, optimizes thunder
The ambiguity function reached;5) signal integral method is realized using based on linear FM signal (LFM) spread spectrum.Preceding three classes one
It is mainly separated by the progress such as orthogonal on time-sharing multiplex or frequency domain to change radar signal and the signal of communication of signal, without abundant
Utilize time resource and band resource.Then two class integration signals do not account for communication synchronization mainly using spread spectrum
Trapped problems.
The content of the invention
The present invention solves the problems, such as:A kind of asymmetric two sections of frequency modulation rates of triangle frequency modulated(FM) radar communicating integral signal sampling
The opposite LFM signals composition of symbol, first portion's LFM signals are used for the synchronization of signal of communication, the modulation of second portion LFM signal subsections
Communication data.For the design of this integration signal waveform, a kind of asymmetric triangle frequency modulated(FM) radar communicating integral letter is proposed
Number parameter optimization method.Based on signal model in Parametric designing procedure using shortest synchronizing signal ensure that communication into
The probability of work(capture, and on the premise of having taken into account radar data reduction simultaneously, ensure that the information under certain error rate condition passes
Signal-to-noise ratio under defeated rate or certain information rate can be applied in radar communication signal integral system.
The present invention technical solution be:A kind of parameter optimization of asymmetric triangle frequency modulated(FM) radar communicating integral signal
Method comprises the following steps:
(1) parameterized model of asymmetric triangle frequency modulated(FM) radar communicating integral signal, i.e. integration signal s are establishedin
(t), it is expressed as:
In formula:B and T is respectively signal bandwidth and time width, s1(t) and s2(t) it is respectively
Positive frequency modulation part and the signal expression of negative frequency modulation part, s1(t) as synchronization acquistion sequence, s2(t) it is logical for average segmentation modulation
Believe bit information, c (t) is in signal s2(t) the communication bit information modulated on;k1And k2Respectively s1(t) and s2(t) frequency modulation
Rate, k1>0, k2<0, and haveT1And T2Respectively s1(t) and s2(t) time span of signal, and have T=T1+
T2;N is in signal s2(t) he number for the communication information modulated on;C (n) is n-th of communication bit information, c (n) for -1 or
1;Δ T is the time width shared by each modulation communication information code element;
(2) unoptimizable parameter in step (1) is initialized, i.e., the slant range resolution according to needed for radar determines radar
Signal bandwidth or Setting signal bandwidth;Meanwhile according to radar transmission power, determine signal time width T;Meanwhile set thresholding
Signal to Noise Ratio (SNR)th;
(3) the signal s in step (1)in(t) and s1(t) it is related to do slip, i.e. communication synchronization sequence capturing, does not examine
Consider modulation communication bit information c (n), obtain the result y (τ) of communication capture, formula is as follows:
In formula, Fresnel integral C (x) and S (x) are as follows:
(4) communication in step (3) captures result y (τ) and 1 He of constraints that communication correctly captures is implemented as described below
Constraints 2, to s1(t) the time span T of signal1It optimizes, obtains the minimum T for meeting communication contact conditions1, it is denoted as
T1_min;
Constraints 1:That is the maximum side petal and step of the y (τ) of step (3)
(3) ratio of the relevant peak point of slip of y (τ), less than η1;
Constraints 2:
η in formula1For the thresholding of setting, N0For the noise power spectral density of setting, B is signal bandwidth, N0B is noise power;
(5) choose and meet T1_min≤T1<T2The T of condition1, according to formula T2=T-T1, T is calculated2, according to formula k1=
B/T1And k2=-B/T2, s is calculated1(t) and s2(t) the frequency modulation rate k of signal1、k2;
(6) according to following constraints 3, to the time width Δ T shared by each modulation communication information code element in step (1)
It optimizes, obtains the Δ T, i.e. Δ T for meeting communication information demodulation condition minimummin;
Constraints 3:
N=T in formula2/ Δ T, N represents the communication code element information number of modulation in formula;SNRrealFor the actual letter of signal of communication
It makes an uproar ratio;
(7) selection meets Δ Tmin≤ΔT≤T2Δ T, and calculate N=T corresponding with Δ T2/ΔT;
(8) if the he number N of current modulation communication information>6 enter step (9), and N≤6 enter step (10);Setting
P is counting variable, initializes p=0,0≤p<P;P is the maximum of p;
(9) P=100 is made, randomly generates the modulation communication information sequence c that one group of code element number is N, formula is as follows:
C=2randint (1, N, [0,1]) -1
Randint (1, N, [0,1]) represents 1 row N row in formula, i.e. length is N, is worth the random sequence for 0 or 1;P at this time
After value increases by 1, p is assigned, is entered step (11);
(10) P=2 is madeN, the communication modulation information c_p that P groups code element number is N is generated, from the communication tune that P groups code element number is N
In information c_p processed, one group of c=c_p is chosen as communication modulation information;After the value of p increases by 1 at this time, p is assigned;Into (12);
C_p is whole combinations of N-bit data, and each bit data is 1 or -1;
(11) by the T in step (5)1With the Δ T in step (7) and communication modulation information c (n) generations in step (9)
Enter in the integration signal model in step (1), modulated the integration signal s of the communication informationin_p(t);According to constraint
Condition 4 judges following ambiguity functionBidimensional auto-correlation function, including apart from auto-correlation function and speed from phase
Close function, if meet radar performance, if being unsatisfactory for constraints 4, return to step (7) if meeting constraints 4, enters
Step (13);
The ambiguity function is by formula
It provides, f in formuladFor Doppler frequency shift;
Constraints 4:
It provides, wherein η2For the thresholding of setting, fd1And fd2The respectively Doppler frequency tolerance of radar;
(12) by the T in step (5)1With the Δ T in step (7) and the communication modulation information c_p (n) in step (10)
It substitutes into the integration signal model in step (1), has been modulated the integration signal s of the communication informationin_p(t), according to step
Suddenly the constraints 4 of (11), the bidimensional auto-correlation function of the ambiguity function of judgment step (11), including apart from auto-correlation function and
Velocity autocorrelation function, if meet radar performance, if being unsatisfactory for constraints 4, return to step (7), if meeting constraint item
Part 4 enters step (14);
(13) if p meets following condition:p<P, then return to step (9), until p=P terminates, the integration letter optimized
Number, which realizes signal waveform integration on the premise of radar ambiguity function and communication signal-to-noise ratio is met;
(14) if p meets following condition:p<P, then return to step (10), until p=P terminates, the integration optimized
Signal, the integration signal realize signal waveform integration on the premise of radar ambiguity function and communication signal-to-noise ratio is met.
According to the radar slant-range resolution ratio ρ of settingr, determine the bandwidth B of system, formula is as follows:
In formula, c=3 × 108M/s is the light velocity;ρrFor slant range resolution;
Slant range resolution ρrValue be 0.15 meter to 50 meters
According to mean power, peak power and pulse recurrence frequency PRF that the radar emission of setting exports, system is determined
Time width T, formula are as follows:
In formula, PtFor the peak power of radar transmitter output, PavFor the mean power of radar transmitter output, PRF is
Pulse recurrence frequency.
The mean power of the radar transmitter output set is 100 watt to 1000 watts.
For the peak power of the radar transmitter output set as 1 kilowatt to 10 kilowatts, the size of peak power set point value is big
In the value of mean power.
The pulse recurrence frequency PRF set is 1000Hz to 9000Hz.
The present invention compared with prior art the advantages of be:
(1) present invention is asymmetrical cam design by frequency curve shape, realizes the synchronizing sequence for communication
Tuning rate between signal and the signal for having modulated communication data has certain orthogonality, can realize the synchronization acquistion of communication
With data modulation, demodulation.And integration signal takes full advantage of system time and frequency resource.
(2) present invention proposes the object function based on radar ambiguity function and communication performance come optimization design integration signal
Parameter method, this method ensure that the probability for the successfully capture that communicates using shortest synchronizing signal, and ensure that radar
On the premise of detection performance, ensure the rate of information throughput under certain error rate condition or the noise under certain information rate
Than.
(3) integration signal that the parameterization design method carried by the present invention designs, while taken into account radar performance
And communication performance, it is effective and feasible in the design of radar-communication integration signal, suitable for radar-communication integration system
Practical engineering application.
(4) on the basis of the present invention is for the FM signal being used widely in field of radar, letter is obscured based on radar
Number and the parameterized model of communication capture and data transmission performance propose a kind of asymmetric triangle frequency modulated(FM) radar communicating integral letter
Number parameter optimization method, just (or negative) of the integration signal adjusts the FM signal of frequency-portions for synchronizing sequence, bear (or
The FM signal of frequency-portions just) is adjusted to be used for modulation communication information.Therefore, the present invention takes full advantage of time resource and frequency band money
Source can meet the design needs that integration signal takes into account radar detection and communication data transfer function simultaneously.
(5) slant range resolution ρ of the present inventionrValue for more than 0.15, by the hardware water of the signal generator on low orbit satellite
Flat limitation, bandwidth B are less than 1GHz, and preferred scheme is ρrValue for 0.15 to 50 meters, suitable for generally should for radar detection
With resolution ratio can reach the performance requirement of current radar detection.
(6) mean power for the radar transmitter output that the present invention is set is 100 watt to 1000 watts.The radar emission of setting
The peak power of machine output is 1 kilowatt to 10 kilowatts, and the size of peak power set point value is more than the value of mean power, preferred side
Case is 1 kilowatt to 10 kilowatts of the peak power in the radar-communication integration system of low rail (400 to 800 kilometers), and radar receives
The signal-to-noise ratio of machine can reach more than 10dB, realize the capability improving of radar detection target, and detection target is accurate.
(7) for the pulse recurrence frequency PRF that the present invention is set as 1000Hz to 9000Hz, PRF generates meeting in more than 9000Hz
Lead to the problem of range ambiguity, preferred scheme is PRF in 1000Hz to 9000Hz, can realize the purpose for detecting target.
Description of the drawings
Fig. 1 is the process flow block diagram of the present invention;
Fig. 2 is the frequency curve schematic diagram of non-symmetric triangular frequency modulated(FM) radar communicating integral signal;
Fig. 3 is the correct acquisition probability of capture of communicating under the conditions of thresholding signal-to-noise ratio;
Fig. 4 is the peak sidelobe ratio of communication acquisition procedure with the change curve of synchronizing sequence length;
Fig. 5 be the ambiguity function of integration signal along reference axis cutting drawing (communication information c_5={ -1,1,1 } of modulation),
Fig. 5 (a) is cutting drawing of the ambiguity function along y-coordinate axis, i.e., apart from auto-correlation function, Fig. 5 (b) is ambiguity function along x coordinate axis
Cutting drawing, i.e. velocity autocorrelation function;
Fig. 6 is the ambiguity function of integration signal along fd=200Hz axis cutting drawing (modulation communication information c_5=-
1,1,1});
Fig. 7 is to be based on the integrated radar-communication integration system block diagram of signal waveform.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
The present invention basic ideas be:A kind of parameter optimization side of asymmetric triangle frequency modulated(FM) radar communicating integral signal
Method, the frequency curve shape of the integration signal is asymmetrical triangle, is believed by two sections of opposite linear frequency modulations of tune frequency symbol
Number (LFM) signal composition, wherein first portion's signal is that frequency modulation rate is k1For LFM signals, the communication information is not modulated, for carrying out
The synchronization or capture of signal of communication, second portion signal are that frequency modulation rate is k2LFM signals (wherein k1×k2<0), at second
With the form modulation communication code element information of (- 1,1) code on sub-signal.Based on asymmetric triangle frequency modulated(FM) radar communicating integral signal
Parameterized model, with characterize the ambiguity function of radar data reduction and characterize communication performance object function carry out optimization design one
The parameter of body signal.Institute's extracting method of the present invention ensure that using shortest synchronizing signal communicate successfully capture probability, and
It ensure that on the premise of radar data reduction, ensure the rate of information throughput under certain signal-to-noise ratio.
As shown in fig. 7, based on the integrated radar-communication integration system of signal waveform, it is including sending system and receiving
System, transmission system include:Integrated waveform generating module, radar transmitter, duplexer and antenna;
Communications encoder, receive communication data, and carry out coding send to encryption processing module be encrypted after obtain it is pending
Send signal of communication;
Radar RF pulse generators generate bandwidth B and the radar radio-frequency pulse of time width T (i.e. system bandwidth B and time width T), send
Radar signal to be sent is obtained after to radar modulator;
Radar signal to be sent is divided into N sections, signal of communication to be sent is believed respectively with radar to be sent after segmentation
Number carry out dot product, obtain integration signal, sent after radar transmitter is modulated up-conversion by duplexer to radio frequency day
Line sends out to space, and duplexer can realize the reception of antenna and the switching of emission function.Due to the letter designed by the present invention
Number frequency curve is an asymmetrical triangle, therefore is referred to as a kind of asymmetric triangle frequency modulated(FM) radar communicating integral signal wave
Shape.
Reception system includes Radar Receiver System and communication receiving system;
Radar Receiver System, including radar receiver, Radar Signal Processing module;
Radar receiver carries out down-converted to received signal, is filtered and send to Radar Signal Processing mould
Block, Radar Signal Processing module carry out radar performance processing analysis, i.e. Range resolution and speed is differentiated, and realization is tested the speed and ranging.
Communication receiving system, including communication control processor antenna, communication control processor, signal of communication processing module, stay of two nights module;
Communication control processor receives signal by communication control processor antenna, send to signal of communication processing module, at signal of communication
Reason module is demodulated received signal, decrypts, and solution coding is sent to stay of two nights module, the symbol letter after finally being demodulated
Breath.
It is as shown in Figure 7 based on the integrated radar-communication integration system block diagram of signal waveform.
As shown in Figure 1, a kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal, including following
Step:
(1) parameterized model of asymmetric triangle frequency modulated(FM) radar communicating integral signal, i.e. integration signal s are establishedin
(t), it is expressed as:
In formula:B and T is respectively signal bandwidth and time width, s1(t) and s2(t) it is respectively
Positive frequency modulation part and the signal expression of negative frequency modulation part, s1(t) as synchronization acquistion sequence, s2(t) it is logical for average segmentation modulation
Believe bit information, c (t) is in signal s2(t) the communication bit information modulated on;k1And k2Respectively s1(t) and s2(t) frequency modulation
Rate, k1>0, k2<0, and haveT1And T2Respectively s1(t) and s2(t) time span of signal, and have T=T1+
T2;N is in signal s2(t) he number for the communication information modulated on;C (n) is n-th of communication bit information, c (n) for -1 or
1;Δ T is the time width shared by each modulation communication information code element;
It is asymmetrical cam design by frequency curve shape, realizes the synchronization sequence signals for communication and modulation
Tuning rate between the signal of communication data has certain orthogonality, can realize the synchronization acquistion of communication and data tune
System, demodulation.And integration signal takes full advantage of system time and frequency resource.
(2) as shown in Fig. 2, being initialized to unoptimizable parameter in step (1), i.e., the oblique distance according to needed for radar is differentiated
Rate determines the signal bandwidth of radar or Setting signal bandwidth;Meanwhile according to radar transmission power, determine signal time width T;Together
When, set thresholding Signal to Noise Ratio (SNR)th;According to the radar slant-range resolution ratio ρ of settingr, determine the bandwidth B of system, formula is as follows:
In formula, c=3 × 108M/s is the light velocity;ρrFor slant range resolution.
Slant range resolution ρrValue for more than 0.15, limited by the level of hardware of the signal generator on low orbit satellite,
Bandwidth B is less than 1GHz, and preferred scheme is ρrValue for 0.15 to 50 meters, suitable for the commonly used of radar detection, resolution ratio
The performance requirement of current radar detection can be reached.
According to mean power, peak power and pulse recurrence frequency PRF that the radar transmitter of setting exports, system is determined
Time width T, formula is as follows:
In formula, PtFor the peak power of radar transmitter output, PavFor the mean power of radar transmitter output, PRF is
Pulse recurrence frequency.The mean power of the radar transmitter output set is 100 watt to 1000 watts.The radar transmitter of setting is defeated
The peak power gone out is 1 kilowatt to 10 kilowatts, and the size of peak power set point value is more than the value of mean power, and preferred scheme is
In the radar-communication integration system of low rail (400 to 800 kilometers), 1 kilowatt to 10 kilowatts of peak power, radar receiver
Signal-to-noise ratio can reach more than 10dB, realize the capability improving of radar detection target, and detection target is accurate.The pulse of setting repeats
Frequency PRF is 1000Hz to 9000Hz, and PRF can cause target to lead to the problem of range ambiguity, preferred scheme in more than 9000Hz
It is PRF 1000 to 9000, can realizes the accuracy of the function of detection target, i.e., accurately detect target.
(3) the signal s in step (1)in(t) and s1(t) it is related to do slip, i.e. communication synchronization sequence capturing, does not examine
Consider modulation communication bit information c (n), obtain the result y (τ) of communication capture, formula is as follows:
In formula, Fresnel integral C (x) and S (x) are as follows:
(4) communication in step (3) captures result y (τ) and 1 He of constraints that communication correctly captures is implemented as described below
Constraints 2, to s1(t) the time span T of signal1It optimizes, obtains the minimum T for meeting communication contact conditions1, it is denoted as
T1_min;
Constraints 1:That is the maximum side petal and step of the y (τ) of step (3)
(3) ratio of the relevant peak point of slip of y (τ), less than η1;
Constraints 2:
η in formula1For the thresholding of setting, N0For the noise power spectral density of setting, B is signal bandwidth, N0B is noise power;
Noise power N0B is made an uproar to obtain by test communication control processor bottom.
Step (3) and step (4) can be succeeded the constraint that captures parameter T by synchronizing sequence1Magnitude range.
(5) choose and meet T1_min≤T1<T2The T of condition1, according to formula T2=T-T1, T is calculated2, according to formula k1=
B/T1And k2=-B/T2, s is calculated1(t) and s2(t) the frequency modulation rate k of signal1、k2;
(6) according to following constraints 3, to the time width Δ T shared by each modulation communication information code element in step (1)
It optimizes, obtains the Δ T, i.e. Δ T for meeting communication information demodulation condition minimummin;
Constraints 3:
N=T in formula2/ Δ T, N represents the communication code element information number of modulation in formula;SNRrealFor the actual letter of signal of communication
It makes an uproar ratio;
Step (5) by modulated communication code element information linear FM signal can correctly demodulate in the case of signal-to-noise ratio
SNRrealConstraints under obtain T2Value range.
(7) selection meets Δ Tmin≤ΔT≤T2Δ T, and calculate N=T corresponding with Δ T2/ΔT;
(8) if the he number N of current modulation communication information>6 enter step (9), and N≤6 enter step (10);Setting
P is counting variable, initializes p=0,0≤p<P;P is the maximum of p;
(9) P=100 is made, randomly generates the modulation communication information sequence c that one group of code element number is N, formula is as follows:
C=2randint (1, N, [0,1]) -1
Randint (1, N, [0,1]) represents 1 row N row in formula, i.e. length is N, is worth the random sequence for 0 or 1;P at this time
After value increases by 1, p is assigned, is entered step (11);
(10) P=2 is madeN, the communication modulation information c_p that P groups code element number is N is generated, from the communication tune that P groups code element number is N
In information c_p processed, one group of c=c_p is chosen as communication modulation information;After the value of p increases by 1 at this time, p is assigned;Into (12);
C_p is whole combinations of N-bit data, and each bit data is 1 or -1;
(11) by the T in step (5)1With the Δ T in step (7) and communication modulation information c (n) generations in step (9)
Enter in the integration signal model in step (1), modulated the integration signal s of the communication informationin_p(t);According to constraint
Condition 4 judges following ambiguity functionBidimensional auto-correlation function, including apart from auto-correlation function and speed from phase
Close function, if meet radar performance, if being unsatisfactory for constraints 4, return to step (7) if meeting constraints 4, enters
Step (13);
The ambiguity function is by formula
It provides, wherein fdFor Doppler frequency shift;Described judges whether to meet radar performance by formula
Constraints 4:
It provides, wherein η2For the thresholding of setting, fd1And fd2The respectively Doppler frequency tolerance of radar;
(12) by the T in step (5)1With the Δ T in step (7) and the communication modulation information c_p (n) in step (10)
It substitutes into the integration signal model in step (1), has been modulated the integration signal s of the communication informationin_p(t), according to step
Suddenly the constraints 4 of (11), the bidimensional auto-correlation function of the ambiguity function of judgment step (11), including apart from auto-correlation function and
Velocity autocorrelation function, if meet radar performance, if being unsatisfactory for constraints 4, return to step (7) otherwise enters step
(14);
(13) if p meets following condition:p<P, then return to step (9), until p=P terminates, the integration letter optimized
Number, which realizes signal waveform integration on the premise of radar ambiguity function and communication signal-to-noise ratio is met;
(14) if p meets following condition:p<P, then return to step (10), until p=P terminates, the integration optimized
Signal, the integration signal realize signal waveform integration on the premise of radar ambiguity function and communication signal-to-noise ratio is met.
Step (7) meets the parameter of radar system performance to step (14) by the constraints optimization of radar performance
T1, T2, Δ T.Finally met communication system performance and the integration signal waveform of radar system performance.
Embodiment:
A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal, specific implementation step are as follows:
(1) parameterized model of asymmetric triangle frequency modulated(FM) radar communicating integral signal is established;
Integration signal can be expressed as:
Wherein:
The B and T is respectively signal bandwidth and time width;The s1(t) and s2(t) it is respectively positive frequency modulation part and negative frequency modulation
Part, wherein s1(t) it is synchronization acquistion sequence, modulate communications information is in s2(t) on signal;The c (t) is in signal s2(t) on
The communication information of modulation;The k1And k2Respectively s1(t) and s2(t) frequency modulation rate (wherein k1>0, k2<0, and have
);The T1And T2Respectively s1(t) and s2(t) signal length, and have T=T1+T2;The N is in signal s2(t) modulated on
The communication information he number;The c (n) is n-th of communication information, is worth for -1 or 1;The Δ T is each code element institute
The time width accounted for;
(2) unoptimizable parameter in step (1) is initialized, signal bandwidth B=70MHz;According to thunder signal time width T=
500us;Assuming that Signal to Noise Ratio (SNR)=- 25dB of practical communication environment, gives thresholding Signal to Noise Ratio (SNR)th=-25dB;
(3) signal in step (1), for simplicity, without considering communication modulation information c (n), according to communication capture
Process obtains the result y of communication capture, and the y is by formula
It provides, wherein, Fresnel integral C (x) and S (x)
(4) constraints that the communication capture result y in step (3) and realization communication correctly capture, to parameter T1
It optimizes, obtains the minimum T for meeting communication contact conditions1, it is denoted as T1_min, the constraints is by formula
Constraints 1:
Constraints 2:
It provides, the thresholding η being provided with1=-13dB, N0For noise power spectral density, B is signal bandwidth, N0B is noise
Power is made an uproar acquisition in practice by testing the bottom of communication control processor;
Fig. 3 gives SNRthUnder the conditions of=- 25dB, meet the acquisition probability of constraints 2, from the figure 3, it may be seen that believing in thresholding
It makes an uproar than a timing, the correct acquisition probability of synchronizing sequence is with T in communication1The increase of/T and increase.Work as T1During/T >=0.25, lead to
The synchronization acquistion probability of letter is 100%.Fig. 4 gives synchronizing sequence acquisition procedure, peak sidelobe ratio (maximum secondary lobe and capture peak
The ratio of value) with synchronizing sequence length change curve, as seen from the figure, work as T1During/T >=0.03, peak sidelobe ratio is less than-
13dB.In thresholding η1During=- 13dB, meet the T of constraints 11≥0.03T.Therefore T1_min=0.25T.
(5) choose and meet T1_min≤T1<T2The T of condition1=0.25T=125us, according to formula T2=T-T1T is calculated2
=375us, according to formula k1=B/T1And k2=-B/T2K is calculated1=0.56 × 1012、k2=0.187 × 1012;
(6) constraints demodulated according to the communication information is realized, to the time width Δ T shared by each modulation communication information code element
(or he number N to modulating the communication information) is optimized, and obtains the Δ for meeting communication information demodulation condition minimum
T, i.e. Δ Tmin, the constraints is by formula
Constraints 3:
By optimization, meet the Δ T that the communication information demodulates condition minimummin=T1=125us;
(7) selection meets Δ Tmin≤ΔT≤T2Δ T=125us, and calculate corresponding N=3.
(8) because the he number N of current modulation communication information is unsatisfactory for following condition:N>6, into (10);Initialize p
=0;
(9) P=100 is made, randomly generates the modulation communication information sequence that one group of code element number is N, it is described to randomly generate sequence
By formula
C=2randint (1, N, [0,1]) -1
It provides, wherein randint (1, N, [0,1]) is the subfunction in Matlab, represents that it is that N values are 0 or 1 to generate length
Random sequence;P=p+1, into (11);
(10) P=2 is madeN=8, generate all combinations for the communication modulation information that 8 groups of code element numbers are 3:C_p (n)=- 1,
1 }, wherein p=1,2 ..., 8, i.e. c_1={ 1,1,1 }, c_2={ 1,1, -1 }, c_3={ 1, -1,1 }, c_4={ 1, -1, -1 },
C_5={ -1,1,1 }, c_6={ -1,1, -1 }, c_7={ -1, -1,1 } and c_8={ -1, -1, -1 } choose one group of c=c_p
(n) it is used as communication modulation information;P=p+1;Into (12);
(11) by the T in step (5)1With the Δ T in step (7) and in step (9) communication modulation information c (n=1,
2 ..., N) it substitutes into the integration signal model in step (1), modulated the integration signal s of the communication informationin_p(t),
Whether the bidimensional auto-correlation function (including apart from auto-correlation function and velocity autocorrelation function) of analysis ambiguity function meets radar
Performance if being unsatisfactory for radar performance, into (7), otherwise enters (13);
The ambiguity function is by formula
It provides, wherein fdFor Doppler frequency shift;Described judges whether to meet radar performance by formula
Constraints 4:
It provides, wherein η2For the thresholding of setting, fd1And fd2The Doppler frequency tolerance that respectively radar considers;
(12) by the T in step (5)1Communication modulation information c=c_p with the Δ T in step (7) and in step (10)
(n) (n=1,2 ..., N;P=1,2 ..., P) it substitutes into the integration signal model in step (1), modulated communication letter
The integration signal s of breathin_p(t), analyze ambiguity function bidimensional auto-correlation function (including apart from auto-correlation function and speed from
Correlation function) whether meet radar performance, if being unsatisfactory for radar performance, into (7), otherwise into (14), described is fuzzy
Function is by formula
It provides, wherein fdFor Doppler frequency shift;Described judges whether to meet radar performance by formula
Constraints 4:
It provides, the thresholding η being provided with2=-13dB, fd1=-200Hz and fd2=200Hz is respectively the more of radar consideration
General Le frequency tolerance;
(13) if p meets following condition:p<P, then return to step (9), otherwise terminate.
(14) if p meets following condition:p<P, then return to step (10), otherwise terminate.
Fig. 5 (a) and (b) give the integration signal s of modulation communication information c_5={ -1,1,1 }in_p(t) (p=5) phase
The distance by radar and velocity autocorrelation function answered.From figure 5 it can be seen that the integration signal apart from auto-correlation function zero
Point nearby has higher secondary lobe (peak sidelobe ratio is about -13.4715dB), and the secondary lobe away from zero point is relatively low but very steady, and fast
Spending auto-correlation function has extremely low secondary lobe, is approximately the impulse in origin.Fig. 6 gives Doppler frequency shift fdDuring=200Hz
The distance dimension ambiguity function of integration signal, can obtain peak sidelobe ratio as -13.4109dB from Fig. 6.Radar considers more
General Le frequency tolerance is fd∈ (- 200Hz, 200Hz), thresholding η2During=- 13dB, in constraints 4 corresponding 8 groups it is different
The MSR of modulation communication information be respectively -13.2870dB, -13.2624dB, -13.2625dB, -13.2854dB, -
13.3122dB, -13.2892dB, -13.2873dB and -13.3121dB, wherein maximum are -13.2624dB, are satisfied by constraining
Condition 4.
It is of the invention mainly to be verified that all steps, conclusion are all verified just on Matlab7.12.0 using emulation data
Really.
By the specific embodiment of the invention as can be seen that a kind of asymmetric triangle frequency modulated(FM) radar provided by the present invention leads to
Believe the parameter optimization method of integration signal, the probability for the successfully capture that communicates is ensure that using shortest synchronizing signal, and protecting
On the premise of having demonstrate,proved radar data reduction, ensure under the rate of information throughput or the certain information rate under certain error rate condition
Signal-to-noise ratio, which has taken into account radar performance and communication performance simultaneously, in setting for radar-communication integration signal
It is effective and feasible in meter, suitable for the practical engineering application of radar-communication integration system.
The content not being described in detail in description of the invention belongs to the known technology of those skilled in the art.
Claims (5)
1. a kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal, it is characterised in that:For being based on
The integrated radar-communication integration system of signal waveform includes including sending system and reception system, transmission system:Integration
Waveform generating module, radar transmitter, duplexer and antenna;
Communications encoder receives communication data, and carry out coding send to encryption processing module be encrypted after obtain it is to be sent logical
Believe signal;
Radar RF pulse generators generate bandwidth B and the radar radio-frequency pulse of time width T, are obtained after sending to radar modulator to be sent
Radar signal;
It is as follows that the integrated step of signal waveform is realized to sent radar signal:
(1) parameterized model of asymmetric triangle frequency modulated(FM) radar communicating integral signal, integration signal s are establishedin(t), represent
For:
In formula:B and T is respectively signal bandwidth and time width, is differentiated according to the radar slant-range of setting
Rate ρr, determine the bandwidth B of system, formula is as follows:C=3 × 108M/s is the light velocity;ρrFor slant range resolution;According to setting
Mean power, peak power and the pulse recurrence frequency PRF of fixed radar emission output, determine the time width T of system, formula is such as
Under:PtFor the peak power of radar transmitter output, PavFor radar transmitter output mean power,
PRF is pulse recurrence frequency;
s1(t) and s2(t) be respectively positive frequency modulation part and negative frequency modulation part signal expression, s1(t) it is used as synchronization acquistion sequence
Row, s2(t) to be averagely segmented modulation communication bit information, c (t) is in signal s2(t) the communication bit information modulated on;k1With
k2Respectively s1(t) and s2(t) frequency modulation rate, k1>0, k2<0, and haveT1And T2Respectively s1(t) and s2(t)
The time span of signal, and have T=T1+T2;N is in signal s2(t) he number for the communication information modulated on;C (n) is n-th
A communication bit information, c (n) are -1 or 1;Δ T is the time width shared by each modulation communication information code element;
(2) unoptimizable parameter in step (1) is initialized, the slant range resolution according to needed for radar determines the signal of radar
Bandwidth or Setting signal bandwidth;Meanwhile according to radar transmission power, determine signal time width T;Meanwhile set thresholding signal-to-noise ratio
SNRth;
(3) the signal s in step (1)in(t) and s1(t) it is related to do slip, communication synchronization sequence capturing, without considering modulation
Communication bit information c (n) obtains the result y (τ) of communication capture, and formula is as follows:
In formula, Fresnel integral C (x) and S (x) are as follows:
(4) communication in step (3) captures result y (τ) and constraints 1 and constraint that communication correctly captures is implemented as described below
Condition 2, to s1(t) the time span T of signal1It optimizes, obtains the minimum T for meeting communication contact conditions1, it is denoted as T1_min;
Constraints 1:The maximum side petal of the y (τ) of step (3) is with the y's (τ) of step (3)
The ratio of relevant peak point is slided, less than η1;
Constraints 2:
η in formula1For the thresholding of setting, N0For the noise power spectral density of setting, B is signal bandwidth, N0B is noise power;
(5) choose and meet T1_min≤T1<T2The T of condition1, according to formula T2=T-T1, T is calculated2, according to formula k1=B/T1
And k2=-B/T2, s is calculated1(t) and s2(t) the frequency modulation rate k of signal1、k2;
(6) according to following constraints 3, the time width Δ T shared by each modulation communication information code element in step (1) is carried out
Optimization obtains the Δ T, Δ T for meeting communication information demodulation condition minimummin;
Constraints 3:
N=T in formula2/ Δ T, N represents the communication he number of modulation in formula;SNRrealFor the actual signal-to-noise ratio of signal of communication;
(7) selection meets Δ Tmin≤ΔT≤T2Δ T, and calculate N=T corresponding with Δ T2/ΔT;
(8) if the he number N of current modulation communication information>6 enter step (9), and N≤6 enter step (10);Set p as
Counting variable initializes p=0,0≤p<P;P is the maximum of p;
(9) P=100 is made, randomly generates the modulation communication information sequence c that one group of code element number is N, formula is as follows:
C=2randint (1, N, [0,1]) -1
Randint (1, N, [0,1]) represents 1 row N row in formula, and length N is worth the random sequence for 0 or 1;The value of p increases at this time
After 1, p is assigned, is entered step (11);
(10) P=2 is madeN, the communication modulation information c_p that P groups code element number is N is generated, from the communication modulation information that P groups code element number is N
In c_p, one group of c=c_p is chosen as communication modulation information;After the value of p increases by 1 at this time, p is assigned;Into (12);
C_p is whole combinations of N-bit data, and each bit data is 1 or -1;
(11) by the T in step (5)1Step is substituted into the Δ T in step (7) and the communication modulation information c (n) in step (9)
(1) in the integration signal model in, the integration signal s of the communication information has been modulatedin_p(t);According to constraints 4,
Judge following ambiguity function χsin(τ,fd) bidimensional auto-correlation function, including apart from auto-correlation function and velocity autocorrelation function,
Whether radar performance is met, if being unsatisfactory for constraints 4, return to step (7) if meeting constraints 4, enters step
(13);
The ambiguity function is by formula
It provides, f in formuladFor Doppler frequency shift;
Constraints 4:
It provides, wherein η2For the thresholding of setting, fd1And fd2The respectively Doppler frequency tolerance of radar;
(12) by the T in step (5)1It is substituted into the Δ T in step (7) and the communication modulation information c_p (n) in step (10)
In integration signal model in step (1), the integration signal s of the communication information has been modulatedin_p(t), according to step
(11) constraints 4, the bidimensional auto-correlation function of the ambiguity function of judgment step (11), including apart from auto-correlation function and speed
Spend auto-correlation function, if meet radar performance, if being unsatisfactory for constraints 4, return to step (7), if meeting constraints
4, it enters step (14);
(13) if p meets following condition:p<P, then return to step (9), until p=P terminates, the integration signal optimized,
The integration signal realizes signal waveform integration on the premise of radar ambiguity function and communication signal-to-noise ratio is met;
(14) if p meets following condition:p<P, then return to step (10), until p=P terminates, the integration signal optimized,
The integration signal realizes signal waveform integration, obtains on the premise of radar ambiguity function and communication signal-to-noise ratio is met
Finally meet communication system performance and the integration signal waveform of radar system performance.
2. a kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal according to claim 1,
It is characterized in that:Slant range resolution ρrValue be 0.15 meter to 50 meters.
3. a kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal according to claim 1,
It is characterized in that:The mean power of the radar transmitter output set is 100 watt to 1000 watts.
4. a kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal according to claim 1,
It is characterized in that:The peak power of the radar transmitter output set is 1 kilowatt to 10 kilowatts, the size of peak power set point value
More than the value of mean power.
5. a kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal according to claim 1,
It is characterized in that:The pulse recurrence frequency PRF set is 1000Hz to 9000Hz.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510450571.XA CN105099976B (en) | 2015-07-28 | 2015-07-28 | A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510450571.XA CN105099976B (en) | 2015-07-28 | 2015-07-28 | A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105099976A CN105099976A (en) | 2015-11-25 |
CN105099976B true CN105099976B (en) | 2018-06-01 |
Family
ID=54579545
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510450571.XA Active CN105099976B (en) | 2015-07-28 | 2015-07-28 | A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105099976B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106990392B (en) * | 2017-02-15 | 2019-06-25 | 中国人民解放军国防科学技术大学 | A kind of extraterrestrial target fine motion information acquisition method based on random stepped frequency signal |
CN106911605A (en) * | 2017-03-01 | 2017-06-30 | 哈尔滨工业大学 | A kind of radar-communication integration Waveform generation method based on FMCW |
CN109085543B (en) * | 2018-08-24 | 2020-08-18 | 斯凯瑞利(北京)科技有限公司 | Linear frequency modulation signal generator |
CN109507661B (en) * | 2018-09-28 | 2023-04-07 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Radar and communication integrated signal processing method |
CN113238198B (en) * | 2021-03-31 | 2023-12-15 | 南京航空航天大学 | Radar communication integrated system radiation time control method based on time division multiplexing |
CN113687344B (en) * | 2021-07-20 | 2023-08-11 | 西安空间无线电技术研究所 | Triangular wave modulation linear frequency modulation continuous wave radar speed measuring method |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103019744A (en) * | 2012-12-31 | 2013-04-03 | 清华大学 | Computing middleware-based radar signal processing module library construction method and application thereof |
CN104569927A (en) * | 2015-01-27 | 2015-04-29 | 中国人民解放军国防科学技术大学 | Radar pulse transmitter and generation method of orthogonal frequency division multiplexing (OFDM) pulse signal |
CN104580052A (en) * | 2013-10-20 | 2015-04-29 | 杨瑞娟 | Novel application of multi-carrier direct spread signal in radar communication integration |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100939276B1 (en) * | 2008-04-22 | 2010-01-29 | 인하대학교 산학협력단 | UWB distance measurement system and driving method thereof |
ES2351868T3 (en) * | 2008-05-23 | 2011-02-11 | Thales Nederland B.V. | PROCEDURE TO EVALUATE THE POSITION AND SPEED OF AN OBJECTIVE WITH A RADAR THAT ISSUING A OFDM WAVE FORM. |
-
2015
- 2015-07-28 CN CN201510450571.XA patent/CN105099976B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103019744A (en) * | 2012-12-31 | 2013-04-03 | 清华大学 | Computing middleware-based radar signal processing module library construction method and application thereof |
CN104580052A (en) * | 2013-10-20 | 2015-04-29 | 杨瑞娟 | Novel application of multi-carrier direct spread signal in radar communication integration |
CN104569927A (en) * | 2015-01-27 | 2015-04-29 | 中国人民解放军国防科学技术大学 | Radar pulse transmitter and generation method of orthogonal frequency division multiplexing (OFDM) pulse signal |
Also Published As
Publication number | Publication date |
---|---|
CN105099976A (en) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105099976B (en) | A kind of parameter optimization method of asymmetric triangle frequency modulated(FM) radar communicating integral signal | |
CN108768446B (en) | The signal waveform design method of low probability of intercept radar communication integrated system | |
CN105162742B (en) | A kind of asymmetric triangle frequency modulated(FM) radar communicating integral signal waveform determines method | |
CN106161310B (en) | A kind of multicarrier difference chaotic shift keying modulation demodulation method and modem | |
CN109507661A (en) | Radar and communicating integral signal processing method | |
CN107026810A (en) | The PN synchronization method of burst directly-enlarging system and its DS waveform that happens suddenly | |
CN114124238B (en) | OTFS communication radar integrated waveform design method based on time division system | |
CN103905085B (en) | One is burst hybrid spread spectrum underwater sound concealed communication method | |
CN105122080A (en) | Ranging and positioning system | |
CN106911605A (en) | A kind of radar-communication integration Waveform generation method based on FMCW | |
JP2006174411A (en) | Transmission device, reception device, and communication system | |
CN101729089A (en) | Transmitter and receiver of communication system and synchronization method thereof | |
CN109655817B (en) | Relay attack defense method and system based on pulse flight time ranging | |
CN107171994B (en) | Radio fuze signal identification and reconstruction system and method | |
CN109412644A (en) | A kind of doppler frequency estimation method of direct expansion msk signal | |
CN112710992B (en) | Radio frequency stealth radar communication integrated system and method based on four-dimensional antenna array | |
CN107728110A (en) | One kind is based on the warbled shared Design of Signal method of LFM signal phases/tune | |
CN104215954A (en) | Networking method for stepping variable-frequency radar communication integration | |
CN104062641B (en) | The new method that a kind of OFDM radar range profile processes | |
CN104954042A (en) | MSK (minimum shift keying)/WALSH joint demodulating and de-spreading method on basis of FPGA (field programmable gate arrays) | |
CN102769476A (en) | PSWF (prolate spheroidal wave function) pulse-modulated signal synchronous capture method | |
CN1599273B (en) | Synchronous device and method of power line carrier communication transceiver | |
CN102780666B (en) | A kind of anti-timing jitter method of PSWF pulse-modulated signal | |
Wang et al. | Chirp noise waveform aided fast acquisition approach for large Doppler shifted TT&C system | |
CN101873293B (en) | Non-sinusoidal time domain orthogonal demodulation signal synchronization method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |