CN108828596A - GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion - Google Patents
GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion Download PDFInfo
- Publication number
- CN108828596A CN108828596A CN201810667153.XA CN201810667153A CN108828596A CN 108828596 A CN108828596 A CN 108828596A CN 201810667153 A CN201810667153 A CN 201810667153A CN 108828596 A CN108828596 A CN 108828596A
- Authority
- CN
- China
- Prior art keywords
- channel
- individual
- signal
- indicate
- geo
- 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.)
- Pending
Links
Classifications
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
Abstract
The present invention discloses a kind of GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion, in order to solve how to rationally design multichannel configuration under the conditions of meeting signal-to-noise ratio change mark factor constraint to realize that optimal blurred signal compares problem;The present invention realizes that process includes:Firstly, determining GEO-SAR the select of satellite and system parameter;Secondly, the definition of channel allocation problem and modeling;Then, utilization feasibility criterion handles constraint condition, and solves channel in conjunction with difference progress algorithm and configure;Finally output optimal channel configuration;The channel configuration obtained by the method for the invention, can be realized the optimization of signal reconstruction performance, to instruct multichannel GEO star-machine Bistatic SAR practical application.
Description
Technical field
The invention belongs to Radar Technology field, in particular in the case of GEO star-machine multichannel SAR, channel configuration technology.
Background technique
Synthetic aperture radar (SAR) is a kind of round-the-clock, round-the-clock high-resolution imaging system, passes through the transmitting big time
The linear FM signal of bandwidth product, when reception, matched filtering obtained pulse compression signal, to obtain distance to high-resolution, benefit
The high-resolution of orientation is realized with synthetic aperture technique.Image quality is not influenced by weather condition (cloud layer, illumination) etc., is had
The characteristics of distant object is detected and is positioned.The typical application field of SAR includes disaster monitoring, resource exploration, geology
Mapping, military surveillance etc..
Geostationary orbit synthetic aperture radar (GEO SAR) has bigger mapping bandwidth and shorter revisiting period,
Make it possible to widely apply to disaster monitoring, tectonics imaging.GEO star-machine Bistatic SAR uses GEO satellite as irradiation source,
Airborne platform is as receiving station, since Receiver And Transmitter is spatially separating, so that the interference being subject to is greatly reduced.In addition, GEO
Star-machine Bistatic SAR has higher spatial resolution, so GEO star-machine Bistatic SAR is got growing concern for.
In GEO star-machine double-base SAR system, due to the introducing of airborne receiving station, it will lead to doppler bandwidth and become larger, when
When it is more than the transmitting pulse recurrence frequency of GEO-SAR, echo Doppler aliasing will cause, so as to cause occurring in imaging results
False target causes azimuth resolution reduction and image quality to deteriorate.In document " Digital beamforming on
receive:Techniques and optimization strategies for high-resolution wide-swath
SAR imaging, " in IEEE Trans.Aerosp.Electron.Syst., vol.45, no.2, pp.564-592,2009.
The reception information using multiple receivers is proposed to inhibit the method for azimuth ambiguity.But the reconstruction property of the method and multi-pass
Road configuration is related.Solve the conventional method of above-mentioned optimization problem, such as Newton method and steepest descent method.The iteration mistake of this method
Journey needs the derivative information of objective function.Due to channel response matrix be interchannel away from nonlinear function, derivative table
It is extremely difficult up to formula.In addition, this optimization method has biggish numerical instability, algorithm can be made to fall into if initial value setting is improper
Enter locally optimal solution and even result in algorithm and does not restrain.Therefore, numerical optimization is not suitable for solving above-mentioned constrained optimization problem.
Summary of the invention
In order to solve in GEO star-machine SAR, multichannel configures design problem;The present invention proposes a kind of based on feasible
Property criterion GEO star machine Bistatic SAR multi-pass channel-distribution method, while optimize multichannel blurred signal ratio and signal-to-noise ratio become mark because
Son.
The technical solution adopted by the present invention is:GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion,
Including:
S1, the constrained optimization problem for establishing multichannel configuration;The objective function of the constrained optimization problem is according to multichannel
Blurred signal ratio obtains, and the constraint condition of the constrained optimization problem is the maximum that signal-to-noise ratio becomes that the mark factor is less than or equal to setting
Signal-to-noise ratio becomes the mark factor;
S2, initialization differential evolution algorithm, including population size P, maximum evolution number Gmax, become mark factor F and intersection
Probability Cr;
S3, generate initial population, and solve in initial population the blurred signal ratio of each individual and signal-to-noise ratio become mark because
Son;The corresponding N kind channel configuration of each individual;
S4, when the number of iterations is less than maximum evolution number, work as using based on the selection operator of feasibility criterion to determine
The optimal solution of the corresponding constrained optimization problem of preceding individual;
Generate each corresponding composite vector of individual in current group according to mutation operator, then by composite vector with work as
Intersected to obtain trial vector between each channel configuration of preceding solution individual;
S5, step S4 is repeated, until the number of iterations is equal to maximum evolution number, obtains the blurred signal for meeting constraint condition
It is configured than corresponding channel.
Further, constrained optimization problem described in step S1 is:
min AASRb(x), x=(Δ x1,Δx2,...,ΔxN)T
Wherein, AASRb(x) the corresponding blurred signal ratio of channel configuration x is indicated;For maximum signal to noise ratio set by user
Become the mark factor;Φbf(x) indicate that configuration x corresponding signal-to-noise ratio in channel becomes the mark factor, Δ xjIndicate j-th of channel configuration, j=
1,...,N;N indicates receiving channel number;B is the biradical abbreviation of bistatic;Bf be beamforming network wave beam at
The abbreviation of type network;LaIndicate receiving antenna total length.
Further, the N is calculated according to the following formula:
Wherein, BaIndicate orientation bandwidth, PRF indicates transmitting pulse recurrence frequency.
Further, the blurred signal of certain individual described in step S3 is than calculating formula:
Wherein, f indicates Doppler frequency variable, psFor signal power, E [] indicates mean operator, eb(f) it indicates to correspond to
The remaining blurred signal of Doppler frequency.
Further, described in step S3 certain individual signal-to-noise ratio become mark factor calculating formula into:
Wherein, j indicates channel, j=1 ..., N;N indicates receiving channel number;Pj(f) each in representing matrix P (f)
Row, f indicate Doppler frequency variable.
Further, the solution procedure of the matrix P (f) includes the following steps:
A1, calculate according to the following formula j-th of channel apart from history:
Wherein, RR0Indicate receiving station's oblique distance, vsIndicate receiving station's flying speed, t indicates the orientation time;
A2, in order to solve the frequency spectrum in j-th of channel and the phase difference of single channel frequency spectrum, to transmitting station apart from history carry out four
Rank expansion carries out second order expension apart from history to step A1, obtains the time domain orientation phase of j-th of channel echo:
Wherein, RT0Oblique distance for transmitting station in synthetic aperture central instant, k1、k2、k3、k4Indicate satellite distance model
One arrives the expansion coefficient of quadravalence
A3, according to the time domain orientation phase of j-th obtained of channel echo of step A2, acquired by sequence inverting method
The doppler phase of j-th of channel echo:
Wherein, f indicates Doppler frequency variable;
A4, according to the doppler phase of j-th obtained of channel echo of step A3, calculate the channel of GEO star machine Bistatic SAR
Response matrix:
Wherein, Hb,j(f) respective function in j-th of channel, H are indicatedb,j(f)=exp { j [Φb,j(f)-Φb(f)]};
A5, restructuring matrix is obtained according to the channel response matrix of GEO star machine Bistatic SAR:
Further, step S4 is specially:
B1, for each of current group individual xi, random selection three is mutually different standby in current group
Gate road configurationWithComposite vector v is generated by following mutation operatori:
B2, composite vector is solved into x with currentiComponent between intersected to obtain trial vector ui;
Wherein, jrandIt is 1 to the integer being randomly generated between N, vj,iIndicate composite vector viIn component;xj,iIt indicates i-th
The component of individual;
If B3, i-th of individual and trial vector corresponding to i-th of individual are feasible solution, select wherein to obscure letter
It is number more individual as i-th of next iteration than more excellent corresponding feasible solution;
If one is feasible solution in trial vector corresponding to current individual and current individual, another is infeasible solutions,
Then select i-th individual of the feasible solution as next iteration;
If trial vector corresponding to current individual and current individual is infeasible solutions, wherein signal-to-noise ratio change mark is selected
Lesser i-th of the individual as next iteration of the factor.
Beneficial effects of the present invention:GEO star machine Bistatic SAR multichannel configuration side based on feasibility criterion of the invention
Method, it is first determined GEO star-machine Bistatic SAR multichannel selection parameter is defined and models to channel allocation problem;
Then utilization feasibility criterion handles constraint condition, and solves channel in conjunction with difference progress algorithm and configure, and this method is excellent simultaneously
Multichannel blurred signal ratio is changed and signal-to-noise ratio becomes the mark factor.Finally, simulation results show the validity of the method for proposition.
Detailed description of the invention
Fig. 1 is program flow chart provided in an embodiment of the present invention;
Fig. 2 is multichannel GEO star provided in an embodiment of the present invention-machine Bistatic SAR imaging arrangement schematic diagram;
Fig. 3 is that least confusion signal ratio changes schematic diagram in group in iterative process provided in an embodiment of the present invention;
Fig. 4 is multichannel configuration schematic diagram provided in an embodiment of the present invention;
Fig. 5 is single channel echo-wave imaging result before signal reconstruction provided in an embodiment of the present invention;
Fig. 6 is L provided in an embodiment of the present inventionaImaging results after multichannel reconstruct when=3m.
Specific embodiment
For convenient for those skilled in the art understand that technology contents of the invention, with reference to the accompanying drawing to the content of present invention into one
Step is illustrated.
In order to facilitate the description contents of the present invention, following term is explained first:
Term 1:Bistatic SAR
Bistatic SAR refers to the SAR system that system transmitting station and receiving station are placed in different platform, wherein at least one
Platform is motion platform, conceptually belongs to bistatic radar.
Term 2:GEO SAR
GEO SAR is geostationary orbit Synthetic Aperture Radar satellite, uses GEO satellite as irradiation source, operating in has
On the geostationary orbit at certain inclination angle, the cycle of operation is identical as earth rotation period.
Method of the invention has determined GEO star-machine Bistatic SAR multichannel selection parameter first, to channel allocation problem
It is defined and models;Then utilization feasibility criterion handles constraint condition, and solves channel in conjunction with difference progress algorithm
Configuration, this method optimizes multichannel blurred signal ratio simultaneously and signal-to-noise ratio becomes the mark factor.Finally, simulation results show proposing
The validity of method.Including following procedure:
It is as shown in Figure 1 the solution of the present invention flow chart, including:S1, the constrained optimization problem for establishing multichannel configuration;Institute
The objective function for stating constrained optimization problem is obtained according to multichannel blurred signal ratio, and the constraint condition of the constrained optimization problem is
Signal-to-noise ratio becomes the maximum signal to noise ratio change mark factor that the mark factor is less than or equal to setting;Specially:
It is illustrated in figure 2 multichannel GEO star-machine Bistatic SAR imaging arrangement, selects suitable GEO-SAR irradiation source, according to
The geographical location of required imaging and the wave cover characteristic of GEO-SAR determine its orbit parameter and system parameter, such as carrier frequency,
Signal bandwidth and transmitting pulse recurrence frequency.After obtaining the parameter of transmitting station, determine the flight parameter of airborne receiving station with
And beam angle, so as to which the parameters such as GEO star-machine Bistatic SAR doppler bandwidth are calculated according to biradical configuration.
Determine the parameter of airborne receiving station, i.e. receiving station's oblique distance RR0, receiving station flying speed vs, receiving antenna total length
La.According to the transmitting pulse recurrence frequency of the doppler bandwidth and GEO-SAR found out in the first step, receiving channel number is determined
N, such as following formula:
Wherein, Ba=Ta·ka,TaFor synthetic aperture time, kaFor orientation frequency modulation rate,
The constrained optimization problem of multichannel configuration is established according to (2):
Wherein, x is to include N number of channel position Δ xj, the decision vector of j=1 ..., N, channel position Δ xjNo more than
Overall antenna length degree.Become the mark factor for maximum signal to noise ratio set by user, corresponding to the channel configuration that constraint requirements obtain
Signal-to-noise ratio becomes the mark factor no more thanAASR is multichannel blurred signal ratio.
S2, initialization differential evolution algorithm, including population size P, maximum evolution number Gmax, become mark factor F and intersection
Probability Cr.When receiving antenna total length is respectively least confusion in 3m, 2.5m, 2.2m group in iterative process as shown in Figure 3
The variation schematic diagram of signal ratio, it is known that the variation of least confusion signal ratio tends to be steady substantially after the number of iterations is greater than 60 times
It is fixed.Fig. 3 abscissa Generation number indicates the number of iterations;Ordinate MinimunAASR in current
Generation indicates current iteration minimum AASR.
S3, generate initial population, and solve in initial population the blurred signal ratio of each individual and signal-to-noise ratio become mark because
Son;Channel configuration in the corresponding N of each individual;Specially:
It is random to generate P individual composition initial population X0(candidate channel configuration) xi, i=1,2 ..., P, xi=(Δ xi,1,
Δxi,2,…,Δxi,N).Each channel position Δ xi,jIt generates as follows:
Wherein, randij[0,1] is the equally distributed random number between 0 to 1, the jth corresponding to the configuration of i-th of channel
A channel position.Generating initial population X0Later, X is solved0In each alternative channel configure xiBlurred signal ratio AASRbWith
Signal-to-noise ratio becomes mark factor Φbf, and store its performance indicator.Such as following formula:
In definition, E [] indicate mean operator, molecule be current restructuring matrix input and output signal-to-noise ratio ratio,
Denominator is the ratio of the input and output signal-to-noise ratio when pulse recurrence frequency is uniform sampling.Wherein, PjIt (f) is the every of P (f)
A line, reconstruction factors of each same channel of behavior on different sub-band.It is specifically as follows to the solution procedure of P (f):
If j-th of receiving channel and the spacing of reference channel are Δ xj, instantaneous oblique distance of the GEO-SAR away from target point P is RT
(t), j-th of receiving channel and the instantaneous oblique distance of target point P are expressed as RRj(t), t indicates the orientation time, j-th channel away from
It can be expressed as from history:
For reference channel, Δ xj=0, second order expension is carried out to channel distance history is received, available j-th is logical
The time domain orientation phase of road echo:
Wherein, RT0It is oblique distance of the transmitting station in synthetic aperture central instant for hair.It can be in the hope of by sequence inverting method
The doppler phase of above-mentioned time domain echo:
Wherein, f indicates Doppler frequency variable.
In formula (5), if enabling Δ xjThe doppler phase of=0 available single channel GEO star machine Bistatic SAR, i.e.,Therefore, the Doppler frequency spectrum in j-th of channel and single channel Doppler frequency spectrum are divided by, are can be obtained
The band response function of GEO star machine Bistatic SAR:
Hb,j(f)=exp { j [Φb,j(f)-Φb(f)]} (8)
Then, H (f) is that channel response matrix can indicate as follows:
Restructuring matrix is the inverse matrix of H (f):
Signal-to-noise ratio can be found out through the above steps becomes the mark factor.
AASR is multichannel blurred signal ratio, is expressed as follows:
Wherein, psFor signal power,Sbs,k(f) it is k-th
The frequency spectrum of periodic extension.
S4, when the number of iterations is less than maximum evolution number, work as using based on the selection operator of feasibility criterion to determine
The optimal solution of the corresponding constrained optimization problem of preceding individual, as shown in figure 4, channel configuration result is converted to Channel spacing;
Generate each corresponding composite vector of individual in current group according to mutation operator, then by composite vector with work as
Intersected to obtain trial vector between each channel configuration of preceding solution individual;Specially:
Firstly, for each of current group individual xi, random selection three is mutually different in current group
Alternative channel configurationWithComposite vector v is generated by following mutation operatori:
F indicates to become the mark factor;
Then, by composite vector and current solution xiComponent between intersected to obtain trial vector ui。
Wherein, jrand∈ [1,2 ..., D] is 1 to the integer being randomly generated between D, ensure that uiIn at least one component
From composite vector viMiddle selection.As the number rand generated at random between 0 to 1ij[0,1] when being less than crossover probability, uj,iFrom
viMiddle value.Above-mentioned variation and crossover operator can effectively guide group to search for optimal solution in decision space, and guarantee to calculate
The convergence and population diversity of method.
Finally, each alternative channel configuration of calculating is corresponding to be obscured after generating trial vector to individual each in group
Signal ratio and signal-to-noise ratio become the mark factor.Then, using determining current individual x based on the selection operator of feasibility criterioniAnd examination
Test vector uiMiddle objective function and constraint condition more preferably enter the next generation.The selection method of feasibility criterion is as follows:
(1) work as xiWith ui(i.e. x when being feasible solutioniWith uiCorresponding signal-to-noise ratio becomes the mark factor and is both less than equal to),
Select AASRbPerformance more preferably individual x of the feasible solution as next iterationi。
(2) work as xiWith uiIn one be feasible solution, another is infeasible solutions, then selects feasible solution as next iteration
Individual xi。
(3) work as xiWith uiWhen being infeasible solutions, signal-to-noise ratio is selected to become the mark lesser individual as next iteration of the factor
xi。
Selection operator based on feasibility criterion can effectively instruct group close to feasible region of search, accelerate algorithm
Convergence rate.
S5, step S4 is repeated, until the number of iterations is equal to maximum evolution number, obtains the blurred signal for meeting constraint condition
It is configured than corresponding channel.From specific objective result can be seen that Fig. 5 be signal reconstruction before single channel echo-wave imaging as a result,
Fig. 6 is LaImaging results after multichannel reconstruct when=3.
Those of ordinary skill in the art will understand that the embodiments described herein, which is to help reader, understands this hair
Bright principle, it should be understood that protection scope of the present invention is not limited to such specific embodiments and embodiments.For ability
For the technical staff in domain, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made
Any modification, equivalent substitution, improvement and etc. should be included within scope of the presently claimed invention.
Claims (7)
1. the GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion, which is characterized in that including:
S1, the constrained optimization problem for establishing multichannel configuration;The objective function of the constrained optimization problem is fuzzy according to multichannel
Signal ratio obtains, and the constraint condition of the constrained optimization problem is the maximum noise that signal-to-noise ratio becomes that the mark factor is less than or equal to setting
The factor is marked than becoming;
S2, initialization differential evolution algorithm, including population size, maximum evolution number become the mark factor and crossover probability;
S3, initial population is generated, and solves the blurred signal ratio of each individual and signal-to-noise ratio in initial population and becomes the mark factor;Often
An individual corresponds to the configuration of N kind channel;
S4, when the number of iterations is less than maximum evolution number, determine using based on the selection operator of feasibility criterion when the one before
The optimal solution of the corresponding constrained optimization problem of body;
Each corresponding composite vector of individual in current group is generated according to mutation operator, then by composite vector and current solution
Intersected to obtain trial vector between each channel configuration of individual;
S5, step S4 is repeated, until the number of iterations is equal to maximum evolution number, the blurred signal for obtaining meeting constraint condition is compared
The channel configuration answered.
2. the GEO star machine Bistatic SAR multi-pass channel-distribution method according to claim 1 based on feasibility criterion, feature
It is, constrained optimization problem described in step S1 is:
min AASRb(x), x=(Δ x1,Δx2,...,ΔxN)T
Wherein, AASRb(x) the corresponding blurred signal ratio of channel configuration x is indicated;Become for maximum signal to noise ratio set by user and marks
The factor;Φbf(x) indicate that configuration x corresponding signal-to-noise ratio in channel becomes the mark factor, Δ xjIndicate the configuration of j-th channel, j=1 ...,
N;N indicates receiving channel number;LaIndicate receiving antenna total length.
3. the GEO star machine Bistatic SAR multi-pass channel-distribution method according to claim 2 based on feasibility criterion, feature
It is, the N is calculated according to the following formula:
Wherein, BaIndicate orientation bandwidth, PRF indicates transmitting pulse recurrence frequency.
4. the GEO star machine Bistatic SAR multi-pass channel-distribution method according to claim 2 based on feasibility criterion, feature
It is, the blurred signal of certain individual described in step S3 is than calculating formula:
Wherein, f indicates Doppler frequency variable, psFor signal power, E [] indicates mean operator, eb(f) indicate that how general correspondence is
Strangle the remaining blurred signal of frequency.
5. the GEO star machine Bistatic SAR multi-pass channel-distribution method according to claim 4 based on feasibility criterion, feature
Be, described in step S3 certain individual signal-to-noise ratio become mark factor calculating formula into:
Wherein, j indicates channel, j=1 ..., N;N indicates receiving channel number;Pj(f) every a line in representing matrix P (f), f
Indicate Doppler frequency variable.
6. the GEO star machine Bistatic SAR multi-pass channel-distribution method according to claim 5 based on feasibility criterion, feature
It is, the solution procedure of the matrix P (f) includes the following steps:
A1, calculate according to the following formula j-th of channel apart from history:
Wherein, RR0Indicate receiving station's oblique distance, vsIndicate receiving station's flying speed, t indicates the orientation time;
A2, in order to solve the frequency spectrum in j-th of channel and the phase difference of single channel frequency spectrum, to transmitting station apart from history carry out quadravalence exhibition
It opens, second order expension is carried out apart from history to step A1, obtains the time domain orientation phase of j-th of channel echo:
Wherein, RT0Oblique distance for transmitting station in synthetic aperture central instant, k1、k2、k3、k4Indicate that the one of satellite distance model arrives
The expansion coefficient of quadravalence
A3, according to the time domain orientation phase of j-th obtained of channel echo of step A2, jth is acquired by sequence inverting method
The doppler phase of a channel echo:
Wherein, f indicates Doppler frequency variable;
A4, according to the doppler phase of j-th obtained of channel echo of step A3, calculate the channel response of GEO star machine Bistatic SAR
Matrix:
Wherein, Hb,j(f) respective function in j-th of channel, H are indicatedb,j(f)=exp { j [Φb,j(f)-Φb(f)]};
A5, restructuring matrix is obtained according to the channel response matrix of GEO star machine Bistatic SAR:
7. the GEO star machine Bistatic SAR multi-pass channel-distribution method according to claim 6 based on feasibility criterion, feature
It is, step S4 is specially:
B1, for each of current group individual xi, three mutually different alternative channels are randomly choosed in current group
ConfigurationWithComposite vector v is generated by following mutation operatori:
Wherein, F indicates to become the mark factor;
B2, composite vector is solved into x with currentiComponent between intersected to obtain trial vector ui;
Wherein, jrandIt is 1 to the integer being randomly generated between N, vj,iIndicate composite vector viIn component;xj,iIndicate i-th of individual
Component;Cr indicates crossover probability;
If B3, i-th of individual and trial vector corresponding to i-th of individual are feasible solution, wherein blurred signal ratio is selected
I-th individual of the more excellent corresponding feasible solution as next iteration;
If one is feasible solution in trial vector corresponding to current individual and current individual, another is infeasible solutions, then selects
Select i-th individual of the feasible solution as next iteration;
If trial vector corresponding to current individual and current individual is infeasible solutions, the wherein signal-to-noise ratio change mark factor is selected
Lesser i-th of individual as next iteration.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810667153.XA CN108828596A (en) | 2018-06-26 | 2018-06-26 | GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810667153.XA CN108828596A (en) | 2018-06-26 | 2018-06-26 | GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108828596A true CN108828596A (en) | 2018-11-16 |
Family
ID=64137528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810667153.XA Pending CN108828596A (en) | 2018-06-26 | 2018-06-26 | GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108828596A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110018474A (en) * | 2019-01-25 | 2019-07-16 | 北京理工大学 | Three-D imaging method based on geostationary orbit synthetic aperture radar chromatographic technique |
CN110133749A (en) * | 2019-05-30 | 2019-08-16 | 中国地质科学院 | Gravity and magnetic data processing method and its system in a kind of exploration of geological resource |
CN113419241A (en) * | 2021-05-28 | 2021-09-21 | 中国科学院空天信息创新研究院 | Signal reconstruction method and equipment |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101609151A (en) * | 2009-07-17 | 2009-12-23 | 重庆大学 | A kind of moving target detection method that decomposes based on single-channel synthetic aperture radar (SAR) image sequence characteristic value |
CN101915911A (en) * | 2010-08-26 | 2010-12-15 | 哈尔滨工业大学 | Space arbitrary configuration distributed SAR moving target parameter estimation method based on cancellation accumulation space-time spectrum |
CN103630894A (en) * | 2013-06-04 | 2014-03-12 | 中国科学院电子学研究所 | Broadband multichannel coherent radar imaging system and control method thereof |
CN103885052A (en) * | 2014-03-25 | 2014-06-25 | 西安空间无线电技术研究所 | Polarization echo separation method of wide-range complete-polarization satellite-borne SAR |
CN106291625A (en) * | 2015-06-24 | 2017-01-04 | 法国国家太空研究中心 | There is the GNSS receiver of the ability for Solution operator carrier track fuzziness of improvement |
US9658341B2 (en) * | 2006-04-28 | 2017-05-23 | Telecommunication Systems, Inc. | GNSS long-code acquisition, ambiguity resolution, and signal validation |
CN107947846A (en) * | 2017-09-26 | 2018-04-20 | 江苏星基航讯科技发展有限公司 | Lead to all purpose aircraft mobile system and communication monitoring system of No.1 satellite based on day |
-
2018
- 2018-06-26 CN CN201810667153.XA patent/CN108828596A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9658341B2 (en) * | 2006-04-28 | 2017-05-23 | Telecommunication Systems, Inc. | GNSS long-code acquisition, ambiguity resolution, and signal validation |
CN101609151A (en) * | 2009-07-17 | 2009-12-23 | 重庆大学 | A kind of moving target detection method that decomposes based on single-channel synthetic aperture radar (SAR) image sequence characteristic value |
CN101915911A (en) * | 2010-08-26 | 2010-12-15 | 哈尔滨工业大学 | Space arbitrary configuration distributed SAR moving target parameter estimation method based on cancellation accumulation space-time spectrum |
CN103630894A (en) * | 2013-06-04 | 2014-03-12 | 中国科学院电子学研究所 | Broadband multichannel coherent radar imaging system and control method thereof |
CN103885052A (en) * | 2014-03-25 | 2014-06-25 | 西安空间无线电技术研究所 | Polarization echo separation method of wide-range complete-polarization satellite-borne SAR |
CN106291625A (en) * | 2015-06-24 | 2017-01-04 | 法国国家太空研究中心 | There is the GNSS receiver of the ability for Solution operator carrier track fuzziness of improvement |
CN107947846A (en) * | 2017-09-26 | 2018-04-20 | 江苏星基航讯科技发展有限公司 | Lead to all purpose aircraft mobile system and communication monitoring system of No.1 satellite based on day |
Non-Patent Citations (1)
Title |
---|
孙稚超: "基于GEO辐射源的星机SAR成像理论与方法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110018474A (en) * | 2019-01-25 | 2019-07-16 | 北京理工大学 | Three-D imaging method based on geostationary orbit synthetic aperture radar chromatographic technique |
CN110133749A (en) * | 2019-05-30 | 2019-08-16 | 中国地质科学院 | Gravity and magnetic data processing method and its system in a kind of exploration of geological resource |
CN113419241A (en) * | 2021-05-28 | 2021-09-21 | 中国科学院空天信息创新研究院 | Signal reconstruction method and equipment |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105182293B (en) | Based on relatively prime array MIMO radar DOA and DOD methods of estimation | |
CN108828596A (en) | GEO star machine Bistatic SAR multi-pass channel-distribution method based on feasibility criterion | |
CN103927551B (en) | Polarimetric SAR semi-supervised classification method based on superpixel correlation matrix | |
CN106772370B (en) | The sparse imaging method in Terahertz near field based on multiple-input multiple-output array | |
CN110187347A (en) | A kind of big breadth imaging method of the biradical synthetic aperture radar of geostationary orbit star machine | |
CN106872970B (en) | A kind of multiple target change data transfer rate tracking device and its method based on differential evolution | |
Addabbo et al. | HRR profile estimation using SLIM | |
CN110266616A (en) | A kind of channel estimation methods based on difference beam angle-measuring method | |
An et al. | Geosynchronous spaceborne–airborne multichannel bistatic SAR imaging using weighted fast factorized backprojection method | |
CN113162670A (en) | Digital multi-beam correction and synthesis method | |
CN111693993B (en) | Self-adaptive 1-bit data radar imaging method | |
CN115085774A (en) | Joint sensation fusion hybrid beam forming method based on Cramer-Lo boundary | |
Zhang et al. | Learning reflection beamforming codebooks for arbitrary RIS and non-stationary channels | |
Sun et al. | JPBA of ARN for target tracking in clutter | |
Chen et al. | Robust DoA estimation using denoising autoencoder and deep neural networks | |
CN107479055B (en) | A kind of self-focusing method and apparatus of high rail satellite-borne SAR | |
CN109946696B (en) | Radar staring correlated imaging method based on target prior information | |
CN115825944B (en) | Single-snapshot multi-target incoming wave direction estimation method based on external radiation source radar | |
McCallum et al. | Star scheduling mode—a new observing strategy for monitoring weak southern radio sources with the AuScope VLBI array | |
Quan et al. | Image encryption algorithm based on chaos system having Markov portion | |
Yuan et al. | Adaptive regularized sparse representation for weather radar echo super-resolution reconstruction | |
RU2431862C1 (en) | Polarisation independent direction finding method of multi-beam radio signals | |
Huang et al. | A novel unambiguous spectrum reconstruction algorithm for space-borne high resolution and wide swath SAR imaging based on iterative adaptive algorithm | |
Luwanga et al. | Automatic spread-F detection using deep learning | |
Huang et al. | Space-Based Electromagnetic Spectrum Sensing and Situation Awareness |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181116 |