CN110031838A - A kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform - Google Patents
A kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform Download PDFInfo
- Publication number
- CN110031838A CN110031838A CN201910387567.1A CN201910387567A CN110031838A CN 110031838 A CN110031838 A CN 110031838A CN 201910387567 A CN201910387567 A CN 201910387567A CN 110031838 A CN110031838 A CN 110031838A
- Authority
- CN
- China
- Prior art keywords
- wall
- observation position
- vector
- front surface
- signal
- 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.)
- Granted
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/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
- G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform, reflection coefficient sequence and transmitting signal are obtained simultaneously from wall echo time domain sampled data first with sparse blind deconvolution method, to accurately estimate the delay inequality of wall front surface and rear surface back wave, it is then based on the time delay estimation value and delay inequality theoretical value construction objective function of multiple observation positions, the thickness and relative dielectric constant of wall are accurately estimated by minimizing objective function realization, finally single base transceiver is being combined accurately to be estimated with the conductivity for setting lower wall rear surface and front surface reflection wave amplitude comparison wall.The estimation method is especially suitable for the unknown estimation with reference in the case of transmitted waveform to wall parameter, it can be avoided the problem of wall parameter estimation procedure needs to refer to transmitted waveform prior information, the efficiency of wall parameter Estimation is not only increased, and improves the accuracy of wall parameter Estimation.
Description
Technical field
The present invention relates to through-wall radar technical field more particularly to a kind of through-wall radar walls without necessarily referring to transmitted waveform
Body method for parameter estimation, this method is mainly used in the fields such as city law enforcement, disaster assistance and military operation, especially suitable for ginseng
Examine the estimation under transmitted waveform unknown situation to wall parameter.
Background technique
Through-wall radar has the characteristics that distance resolution is high, penetration capacity is strong, it can penetrate cob wall, wood using electromagnetic wave
The nonmetallic wall such as wall, brick wall and concrete is detected and is imaged to human body after wall or constructure inner structure, regardless of
On military or civilian, through-wall radar is all widely used.When electromagnetic wave is by air incidence wall, in wall
Spread speed can reduce, while also occur that phenomena such as refraction, reflection and signal are decayed, so electromagnetic wave penetrates the ability of wall
The accurate estimation of wall parameter not only can be improved into three thickness, dielectric constant and conductivity relating to parameters with wall
Image quality amount, moreover it is possible to guarantee the accuracy of target position positioning.
In recent years, it is contemplated that the importance of wall parameter Estimation has many scholars to be dedicated to the research of wall parameter Estimation.
Currently, wall method for parameter estimation is mostly that the information comprising wall parameter is extracted from wall front and rear surfaces echo, utilizes phase
It closes function formula and calculates determining wall parameter, measurement wall front and rear surfaces reflection echo delay inequality is most commonly seen under normal conditions,
Existing delay time estimation method is generallyd use based on Fast Fourier Transform (FFT) and subspace super-resolution method, however passes through these
The premise that method carrys out estimation time delay is that important affair first passes through metal plate acquisition with reference to transmitted waveform, and which results in the more of actual measurement
Secondary property, it is both complex and time-consuming.
Therefore, how through-wall radar wall parameter accurately to be estimated under the premise of not considering transmitted waveform, is become
Urgent problem to be solved.
Summary of the invention
The technical problem to be solved by the present invention is in view of the above shortcomings of the prior art, provide one kind without necessarily referring to transmitting
Waveform through-wall radar wall method for parameter estimation, realization accurately estimate wall parameter in through-wall radar detection process.
To achieve the above object, the present invention adopts the following technical scheme: a kind of thunder through walls without necessarily referring to transmitted waveform
Up to wall method for parameter estimation, comprising the following steps:
Step 1: the distance of through-wall radar transmitting antenna and receiving antenna to wall front surface is r, and transmitting antenna is kept
Motionless, receiving antenna is moved M times along horizontal line direction by fixed step size, M observation position is obtained, in each observation bit
It sets using Δ τ as time interval and reception signal is sampled, m (m=0,1 ..., M-1) a observation position records N number of uniform
The measurement data T of samplingm(n Δ τ), n=0,1 ..., N-1, Δ τ are sampling time interval；
Step 2: dual-mode antenna being placed in free space, the corresponding each and identical dual-mode antenna spacing of step 1, every
A observation position samples reception signal using Δ τ as time interval, a observation position note of m (m=0,1 ..., M-1)
Record the antenna direct wave measurement data b of N number of uniform samplingm(n Δ τ), n=0,1 ..., N-1；
Step 3: the antenna direct-path signal in measurement data obtained using background cancel method removal step 1 obtains m
The wall echometric measurement data y of (m=0,1 ..., M-1) a N number of uniform sampling of observation positionm(n Δ τ)=Tm(nΔτ)-bm(n
Δ τ), the convolution signal model of n=0,1 ..., N-1, m-th of observation position wall echometric measurement data are expressed as:
ym(n Δ τ)=am(nΔτ)*sm(nΔτ)+vm(nΔτ) (1)
In formula, am(n Δ τ) is the corresponding reference transmitted signal of m-th of observation position, sm(n Δ τ) is m-th of observation position
Corresponding reflection coefficient signals, vm(n Δ τ) is the measurement noise of m-th of observation position, and * indicates convolution algorithm；
Convolution signal model described in formula (1) is expressed as the form of matrix-vector, as shown in formula (2):
ym=Amsm+vm (2)
In formula, ym=[ym(0), ym(Δ τ) ..., ym((N-1)Δτ)]TFor m-th of observation position corresponding N × 1 Wei Qiang
Body echometric measurement data vector, Am=[am0, am1..., am(Q-1)] it is N × Q dimension transmitting signal convolution matrix, sm=[sm(0), sm(Δ
τ) ..., sm((Q-1)Δτ)]TReflection coefficient vector, v are tieed up for Q × 1m=[vm(0), vm(Δ τ) ..., vm((N-1)Δτ)]TFor
The dimension measurement noise vector of N × 1；
Step 4: in m (m=0,1 ..., M-1) a observation position, wall echo being adopted using sparse blind deconvolution method
Sample data are deconvoluted, and obtain reflection coefficient sequence and transmitting signal, while obtaining wall front surface and rear surface back wave
Time delay estimation valueSpecific steps are as follows:
(1) centered on corresponding to the moment at wall echometric measurement data amplitude maximum choose L (L=N/10) a data as
Reference transmitted signal is expressed as the dimensional vector of L × 1And initialize transmitting letter
Number vectorThe number of iterations h=0；
(2) emission signal vector is utilizedIt generates N × Q and ties up convolution matrixWherein Q
=N-L+1, is embodied as:
(3) it is deconvoluted, is obtained sparse to wall echometric measurement data using the sparse algorithm for reconstructing of orthogonal matching pursuit
Reflection coefficient vectorWherein, Q=N-L+1, specific steps are as follows:
1. initializing residual errorSupported collectionFor empty set, the number of iterations k=0；
2. calculating residual errorWith convolution matrixColumn vector inner product in maximum value corresponding to indexed set, i.e.,Wherein, related coefficientQ=0,1 ..., Q-1；
3. updating supported collectionIt calculates
4. updating residual error
5. ifThen stop iteration, exports sparse vectorNon-zero entry
Plain position is by supported collectionIt is determined, corresponding nonzero element value isOtherwise, the number of iterations k adds 1, return step
②；
(4) the reflection coefficient vector obtained using step (3)It generates N × L and ties up convolution matrixIt is embodied as:
(5) emission signal vector is updated
(6) whenWhen stop iteration, otherwise, the number of iterations h adds 1 and returns
It returns step (2)；
(7) output reflection coefficient vectorAnd emission signal vectorObtain wall front surface and rear surface back wave
Round trip transmission time delay difference estimated value, be denoted as
Step 5: the theoretical delay inequality Δ t of wall front surface and rear surface back wave is calculated by geometrical modelm(d,εr),
In m-th of observation position, the theoretical delay inequality of wall front surface and rear surface back wave is expressed as follows:
Wherein, 2LmIt is the distance of m-th observation position transmitting antenna and receiving antenna, c is that electromagnetic wave is propagated in a vacuum
Speed, d are thickness of wall body, εrFor the relative dielectric constant of wall, xmIndicate the position of the corresponding refraction point P of m-th of observation position
It sets, is expressed as
Step 6: construction objective function f (d, εr), obtain thickness of wall body d and relative dielectric constant εrEstimated value；
Utilize the time delay estimation value for the M observation position that step 4 obtainsThe M observation position obtained with step 5
Delay inequality theoretical value Δ tm(d,εr) construction objective function it is as follows:
Minimum value by solving objective function shown in formula (7) obtains thickness of wall body d and relative dielectric constant εrEstimation
Value；
Step 7: using the conductivityσ of the solving result estimation wall of step 6, the specific method is as follows:
Transmitting-receiving is set antenna together to be placed at wall front surface r, obtains the back wave of wall front surface and rear surface
Amplitude R1And R2, therefore, the Amplitude Ratio of the back wave of wall rear surface and front surface is
Solution formula (8) obtains wall loss attenuation rate expression formula
The thickness of wall body d and relative dielectric constant ε that step 6 is estimatedrIt brings formula (9) into, acquires wall loss attenuation rate α,
Lower wall is lost for electromagnetic wave, the conductivityσ of wall is calculated using following formula
Wherein, free space wave impedance η0=120 π.
The invention firstly uses sparse blind deconvolution methods to obtain reflection system simultaneously from wall echo time domain sampled data
Number Sequence and transmitting signal, so that the accurately delay inequality of estimation wall front surface and rear surface back wave, is then based on multiple sights
The time delay estimation value and delay inequality theoretical value that location is set construct objective function, are realized by minimizing objective function to wall
Thickness and relative dielectric constant are accurately estimated, finally set lower wall rear surface together in the single base transceiver of combination and front surface is anti-
Ejected wave Amplitude Ratio accurately estimates the conductivity of wall.
The beneficial effects of the present invention are: in wall parameter estimation procedure using sparse blind deconvolution method from wall return
Reflection coefficient sequence and reference transmitted signal are obtained simultaneously in wave time domain sampled data, then obtains wall under multiple measurement positions
The time delay estimation value of front surface and rear surface back wave, to realize the accurate estimation to wall parameter.It is provided by the invention
Wall method for parameter estimation avoids wall parameter estimation procedure and needs to refer to transmitted waveform prior information, not only increases wall
The efficiency of parameter Estimation, and improve the accuracy of wall parameter estimation result.Through-wall radar wall ginseng provided by the invention
Number estimation method, especially suitable for the unknown estimation with reference in the case of transmitted waveform to wall parameter.
Detailed description of the invention
Fig. 1 is a kind of the total of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform provided by the invention
Body flow chart:
Fig. 2 is a kind of the dilute of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform provided by the invention
Dredge the flow chart of blind deconvolution algorithm；
Fig. 3 is to measure schematic diagram of a scenario under different dual-mode antenna spacing provided by the invention；
Fig. 4 is that same set of single base transceiver provided by the invention measures schematic diagram of a scenario under antenna；
Fig. 5 is using a kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform provided by the invention
The wall of different-thickness carries out the relative error result figure of wall parameter Estimation when being 20dB to signal-to-noise ratio；
Fig. 6 is using a kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform provided by the invention
The relative error result figure of wall parameter Estimation is carried out under the conditions of different signal-to-noise ratio to the wall with a thickness of 0.20m.
Specific embodiment
With reference to the accompanying drawings and examples, specific embodiments of the present invention will be described in further detail.Implement below
Example is not intended to limit the scope of the invention for illustrating the present invention.
Embodiment
As depicted in figs. 1 and 2, a kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform, the party
Method is realized by following step:
Step 1: the distance of through-wall radar transmitting antenna and receiving antenna to wall front surface is r, and transmitting antenna is kept
Motionless, receiving antenna is moved M times along horizontal line direction by fixed step size, M observation position is obtained, in each observation bit
It sets using Δ τ as time interval and reception signal is sampled, m (m=0,1 ..., M-1) a observation position records N number of uniform
The measurement data T of samplingm(n Δ τ), n=0,1 ..., N-1, Δ τ are sampling time interval；
Step 2: dual-mode antenna being placed in free space, the corresponding each and identical dual-mode antenna spacing of step 1, every
A observation position samples reception signal using Δ τ as time interval, a observation position note of m (m=0,1 ..., M-1)
Record the antenna direct wave measurement data b of N number of uniform samplingm(n Δ τ), n=0,1 ..., N-1；
Step 3: the antenna direct-path signal in measurement data obtained using background cancel method removal step 1 obtains m
The wall echometric measurement data y of (m=0,1 ..., M-1) a N number of uniform sampling of observation positionm(n Δ τ)=Tm(nΔτ)-bm(n
Δ τ), the convolution signal model of n=0,1 ..., N-1, m-th of observation position wall echometric measurement data are expressed as:
ym(n Δ τ)=am(nΔτ)*sm(nΔτ)+vm(nΔτ) (1)
In formula, am(n Δ τ) is the corresponding reference transmitted signal of m-th of observation position, sm(n Δ τ) is m-th of observation position
Corresponding reflection coefficient signals, vm(n Δ τ) is the measurement noise of m-th of observation position, and * indicates convolution algorithm；
Convolution signal model described in formula (1) is expressed as the form of matrix-vector, as shown in formula (2):
ym=Amsm+vm (2)
In formula, ym=[ym(0), ym(Δ τ) ..., ym((N-1)Δτ)]TFor m-th of observation position corresponding N × 1 Wei Qiang
Body echometric measurement data vector, Am=[am0, am1..., am(Q-1)] it is N × Q dimension transmitting signal convolution matrix, sm=[sm(0), sm
(Δ τ) ..., sm((Q-1)Δτ)]TReflection coefficient vector, v are tieed up for Q × 1m=[vm(0), vm(Δ τ) ..., vm((N-1)Δτ)]TFor
The dimension measurement noise vector of N × 1；
Step 4: in m (m=0,1 ..., M-1) a observation position, wall echo being adopted using sparse blind deconvolution method
Sample data are deconvoluted, and obtain reflection coefficient sequence and transmitting signal, while obtaining wall front surface and rear surface back wave
Time delay estimation valueSpecific steps are as follows:
(1) centered on corresponding to the moment at wall echometric measurement data amplitude maximum choose L (L=N/10) a data as
Reference transmitted signal is expressed as the dimensional vector of L × 1And initialize transmitting letter
Number vectorThe number of iterations h=0；
(2) emission signal vector is utilizedIt generates N × Q and ties up convolution matrixWherein Q
=N-L+1, is embodied as:
(3) it is deconvoluted, is obtained sparse to wall echometric measurement data using the sparse algorithm for reconstructing of orthogonal matching pursuit
Reflection coefficient vectorWherein, Q=N-L+1, specific steps are as follows:
1. initializing residual errorSupported collectionFor empty set, the number of iterations k=0；
2. calculating residual errorWith convolution matrixColumn vector inner product in maximum value corresponding to indexed set, i.e.,Wherein, related coefficientQ=0,1 ..., Q-1；
3. updating supported collectionIt calculates
4. updating residual error
5. ifThen stop iteration, exports sparse vectorNon-zero entry
Plain position is by supported collectionIt is determined, corresponding nonzero element value isOtherwise, the number of iterations k adds 1, and return step is 2.；
(4) the reflection coefficient vector obtained using step (3)It generates N × L and ties up convolution matrixIt is embodied as:
(5) emission signal vector is updated
(6) whenWhen stop iteration, otherwise, the number of iterations h adds 1 and returns
It returns step (2)；
(7) output reflection coefficient vectorAnd emission signal vectorObtain wall front surface and rear surface back wave
Round trip transmission time delay difference estimated value, be denoted as
Step 5: the theoretical delay inequality Δ t of wall front surface and rear surface back wave is calculated by geometrical modelm(d,εr),
In m-th of observation position, the theoretical delay inequality of wall front surface and rear surface back wave is expressed as follows:
Wherein, 2LmIt is the distance of m-th observation position transmitting antenna and receiving antenna, c is that electromagnetic wave is propagated in a vacuum
Speed, d are thickness of wall body, εrFor the relative dielectric constant of wall, xmIndicate the position of the corresponding refraction point P of m-th of observation position
It sets, is expressed as
Step 6: construction objective function f (d, εr), obtain thickness of wall body d and relative dielectric constant εrEstimated value；
Utilize the time delay estimation value for the M observation position that step 4 obtainsThe M observation position obtained with step 5
Delay inequality theoretical value Δ tm(d,εr) construction objective function it is as follows:
Minimum value by solving objective function shown in formula (7) obtains thickness of wall body d and relative dielectric constant εrEstimation
Value；
Step 7: using the conductivityσ of the solving result estimation wall of step 6, the specific method is as follows:
Transmitting-receiving is set antenna together to be placed at wall front surface r, obtains the back wave of wall front surface and rear surface
Amplitude R1And R2, therefore, the Amplitude Ratio of the back wave of wall rear surface and front surface is
Solution formula (8) obtains wall loss attenuation rate expression formula
The thickness of wall body d and relative dielectric constant ε that step 6 is estimatedrIt brings formula (9) into, acquires wall loss attenuation rate α,
Lower wall is lost for electromagnetic wave, the conductivityσ of wall is calculated using following formula
Wherein, free space wave impedance η0=120 π.
In the present embodiment, using simulation model respectively to a thickness of 0.10m, 0.15m, 0.20m, 0.25m, 0.30m,
0.35m, 6 groups of walls that relative dielectric constant is 6 and conductivity is 0.01S/m carry out parametric inversion.As shown in figure 3, transmitting day
Line and receiving antenna are placed in parallel at a distance of 0.3m at wall 0.5m, and transmitting antenna is motionless, and receiving antenna presses step-length 0.1m
It is moved 9 times along orientation, is corresponding with 10 observation positions, be 1GHz in each observation position driving source centre frequency, when sampling
Between between be divided into 0.0047ns.When thickness of wall body and relative dielectric constant are estimated, 10 observation positions and 4000 sampled points are chosen
Measurement data be used for sparse blind deconvolution, in sparse blind deconvolution solution procedure, most with wall echometric measurement data amplitude
400 data are chosen centered on the sampled point of general goal to be used to initialize transmitting signal.As shown in figure 4, transmitting-receiving is placed with antenna is set
At wall 0.5m, wall front surface and rear surface reflex amplitude R are being measured1And R2, utilize thickness of wall body and opposite dielectric
The estimated value of constant calculates the conductivity of wall.
In the present embodiment, through-wall radar wall parameter estimation result is as shown in Table 1 and Table 2, is recorded in table 1 when signal-to-noise ratio is
20dB, thickness of wall body, relative dielectric constant when thickness of wall body is respectively 0.10m, 0.15m, 0.20m, 0.25m, 0.30m, 0.35m
The data result estimated with conductivity, the relative error of corresponding estimated value are as shown in Figure 5；It is recorded in table 2 when thickness of wall body is
0.20m, thickness of wall body, relative dielectric constant and conductance when signal-to-noise ratio is respectively 10dB, 15dB, 20dB, 25dB, 30dB, 35dB
The data result of rate estimation, the relative error of corresponding estimated value are as shown in Figure 6, it can be seen that the present invention can be unknown with reference to hair
Wall parameter is accurately estimated in the case of ejected wave shape.
Table 1
Table 2
Using through-wall radar wall method for parameter estimation of the invention, avoids wall parameter estimation procedure and need to refer to hair
Ejected wave shape prior information, not only increases the efficiency of wall parameter Estimation, and improves the accuracy of wall parameter Estimation.
Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention., rather than its limitations；To the greatest extent
Pipe present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: its according to
So be possible to modify the technical solutions described in the foregoing embodiments, or to some or all of the technical features into
Row equivalent replacement；And these are modified or replaceed, it does not separate the essence of the corresponding technical solution, and the claims in the present invention are limited
Fixed range.
Claims (1)
1. a kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform, which is characterized in that including walking as follows
It is rapid:
Step 1: the distance of through-wall radar transmitting antenna and receiving antenna to wall front surface is r, and transmitting antenna remains stationary,
Receiving antenna is moved M times along horizontal line direction by fixed step size, and M observation position is obtained, each observation position with
Δ τ is sampled as time interval to signal is received, and m (m=0,1 ..., M-1) a observation position records N number of uniform sampling
Measurement data Tm(n Δ τ), n=0,1 ..., N-1, Δ τ are sampling time interval；
Step 2: dual-mode antenna being placed in free space, the corresponding each and identical dual-mode antenna spacing of step 1, in each sight
Location is set using Δ τ as time interval and is sampled to reception signal, and m (m=0,1 ..., M-1) a observation position records N number of
The antenna direct wave measurement data b of uniform samplingm(n Δ τ), n=0,1 ..., N-1；
Step 3: the antenna direct-path signal in measurement data obtained using background cancel method removal step 1 obtains m (m=
0,1 ..., M-1) a N number of uniform sampling of observation position wall echometric measurement data ym(n Δ τ)=Tm(nΔτ)-bm(nΔτ),
The convolution signal model of n=0,1 ..., N-1, m-th of observation position wall echometric measurement data are expressed as:
ym(n Δ τ)=am(nΔτ)*sm(nΔτ)+vm(nΔτ) (1)
In formula, am(n Δ τ) is the corresponding reference transmitted signal of m-th of observation position, sm(n Δ τ) is that m-th of observation position is corresponding
Reflection coefficient signals, vm(n Δ τ) is the measurement noise of m-th of observation position, and * indicates convolution algorithm；
Convolution signal model described in formula (1) is expressed as the form of matrix-vector, as shown in formula (2):
ym=Amsm+vm (2)
In formula, ym=[ym(0),ym(Δτ),…,ym((N-1)Δτ)]TWall is tieed up for m-th of observation position corresponding N × 1 to return
Wave measurement data vector, Am=[am0,am1,…,am(Q-1)] it is N × Q dimension transmitting signal convolution matrix, sm=[sm(0),sm(Δ
τ),…,sm((Q-1)Δτ)]TReflection coefficient vector, v are tieed up for Q × 1m=[vm(0),vm(Δτ),…,vm((N-1)Δτ)]TFor
The dimension measurement noise vector of N × 1；
Step 4: in m (m=0,1 ..., M-1) a observation position, using sparse blind deconvolution method to wall echo samples number
According to deconvoluting, obtain reflection coefficient sequence and transmitting signal, at the same obtain wall front surface and rear surface back wave when
Prolong poor estimated valueSpecific steps are as follows:
(1) L (L=N/10) a data are chosen centered on corresponding to the moment at wall echometric measurement data amplitude maximum as reference
Emit signal, is expressed as the dimensional vector of L × 1And initialize transmitting signal to
AmountThe number of iterations h=0；
(2) emission signal vector is utilizedIt generates N × Q and ties up convolution matrixWherein, Q=N-L
+ 1, it is embodied as:
(3) it is deconvoluted using the sparse algorithm for reconstructing of orthogonal matching pursuit to wall echometric measurement data, obtains sparse reflection
Coefficient vectorWherein, Q=N-L+1, specific steps are as follows:
1. initializing residual errorSupported collectionFor empty set, the number of iterations k=0；
2. calculating residual errorWith convolution matrixColumn vector inner product in maximum value corresponding to indexed set, i.e.,Wherein, related coefficient
3. updating supported collectionIt calculates
4. updating residual error
5. ifThen stop iteration, exports sparse vector Nonzero element position
By supported collectionIt is determined, corresponding nonzero element value isOtherwise the number of iterations k adds 1, and return step is 2.；
(4) the reflection coefficient vector obtained using step (3)It generates N × L and ties up convolution matrixIt is embodied as:
(5) emission signal vector is updated
(6) whenWhen stop iteration, otherwise, the number of iterations h adds 1 and returns to step
Suddenly (2)；
(7) output reflection coefficient vectorAnd emission signal vectorObtain pair of wall front surface and rear surface back wave
Journey transmission time delay difference estimated value, is denoted as
Step 5: the theoretical delay inequality Δ t of wall front surface and rear surface back wave is calculated by geometrical modelm(d,εr), in m
The theoretical delay inequality of a observation position, wall front surface and rear surface back wave is expressed as follows:
Wherein, 2LmIt is the distance of m-th observation position transmitting antenna and receiving antenna, c is electromagnetic wave spread speed in a vacuum,
D is thickness of wall body, εrFor the relative dielectric constant of wall, xmIt indicates the position of the corresponding refraction point P of m-th of observation position, indicates
For
Step 6: construction objective function f (d, εr), obtain thickness of wall body d and relative dielectric constant εrEstimated value；
Utilize the time delay estimation value for the M observation position that step 4 obtainsThe time delay of the M observation position obtained with step 5
Poor theoretical value Δ tm(d,εr) construction objective function it is as follows:
Minimum value by solving objective function shown in formula (7) obtains thickness of wall body d and relative dielectric constant εrEstimated value；
Step 7: using the conductivityσ of the solving result estimation wall of step 6, the specific method is as follows:
Transmitting-receiving is set antenna together to be placed at wall front surface r, obtains the reflex amplitude R of wall front surface and rear surface1
And R2, therefore, the Amplitude Ratio of the back wave of wall rear surface and front surface is
Solution formula (8) obtains wall loss attenuation rate expression formula
The thickness of wall body d and relative dielectric constant ε that step 6 is estimatedrIt brings formula (9) into, wall loss attenuation rate α is acquired, for electricity
Lower wall is lost in magnetic wave, and the conductivityσ of wall is calculated using following formula
Wherein, free space wave impedance η0=120 π.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910387567.1A CN110031838B (en) | 2019-05-10 | 2019-05-10 | Through-wall radar wall parameter estimation method without reference to transmitted waveform |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910387567.1A CN110031838B (en) | 2019-05-10 | 2019-05-10 | Through-wall radar wall parameter estimation method without reference to transmitted waveform |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110031838A true CN110031838A (en) | 2019-07-19 |
| CN110031838B CN110031838B (en) | 2021-03-12 |
Family
ID=67241706
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910387567.1A Active CN110031838B (en) | 2019-05-10 | 2019-05-10 | Through-wall radar wall parameter estimation method without reference to transmitted waveform |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110031838B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110568439A (en) * | 2019-09-25 | 2019-12-13 | 中南大学 | Impulse type through-wall radar antenna ringing suppression method based on deconvolution |
| CN111398687A (en) * | 2020-03-06 | 2020-07-10 | 浙江省交通运输科学研究院 | Test method for estimating dielectric constant of asphalt pavement |
| CN112147236A (en) * | 2020-09-21 | 2020-12-29 | 大连理工大学 | Ultrasonic signal resolution improving method based on sparse blind deconvolution |
| CN114152943A (en) * | 2021-12-15 | 2022-03-08 | 电子科技大学 | Two-stage wall parameter estimation method based on ultra-wideband through-wall radar |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001021651A (en) * | 1999-07-06 | 2001-01-26 | Hitachi Engineering & Services Co Ltd | Method and device for inspecting inner wall of tunnel |
| CN105911544A (en) * | 2016-05-09 | 2016-08-31 | 西安理工大学 | Spread spectrum through-the-wall radar imaging method based on compression perception technology |
| CN106504209A (en) * | 2016-10-27 | 2017-03-15 | 西安电子科技大学 | The iteration of passive millimeter wave radar image weights blind deconvolution method again |
| CN106772270A (en) * | 2017-01-16 | 2017-05-31 | 哈尔滨工业大学 | The method of sampling and reconstructing method of a kind of radar echo signal |
| CN107589413A (en) * | 2017-08-21 | 2018-01-16 | 南京理工大学 | A kind of ground high resolution radar target multicycle detection method based on compressed sensing |
| CN108562897A (en) * | 2018-01-26 | 2018-09-21 | 桂林电子科技大学 | A kind of sparse imaging method of structure and device of MIMO through-wall radars |
| CN109696672A (en) * | 2018-12-29 | 2019-04-30 | 南京汇君半导体科技有限公司 | A kind of high-resolution through-wall radar imaging method based on space structure relevance |
-
2019
- 2019-05-10 CN CN201910387567.1A patent/CN110031838B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001021651A (en) * | 1999-07-06 | 2001-01-26 | Hitachi Engineering & Services Co Ltd | Method and device for inspecting inner wall of tunnel |
| CN105911544A (en) * | 2016-05-09 | 2016-08-31 | 西安理工大学 | Spread spectrum through-the-wall radar imaging method based on compression perception technology |
| CN106504209A (en) * | 2016-10-27 | 2017-03-15 | 西安电子科技大学 | The iteration of passive millimeter wave radar image weights blind deconvolution method again |
| CN106772270A (en) * | 2017-01-16 | 2017-05-31 | 哈尔滨工业大学 | The method of sampling and reconstructing method of a kind of radar echo signal |
| CN107589413A (en) * | 2017-08-21 | 2018-01-16 | 南京理工大学 | A kind of ground high resolution radar target multicycle detection method based on compressed sensing |
| CN108562897A (en) * | 2018-01-26 | 2018-09-21 | 桂林电子科技大学 | A kind of sparse imaging method of structure and device of MIMO through-wall radars |
| CN109696672A (en) * | 2018-12-29 | 2019-04-30 | 南京汇君半导体科技有限公司 | A kind of high-resolution through-wall radar imaging method based on space structure relevance |
Non-Patent Citations (3)
| Title |
|---|
| CHRISTOPHER THAJUDEEN,ET AL: "A Hybrid Bistatic–Monostatic Radar Technique for Calibration-Free Estimation of Lossy Wall Parameters", 《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》 * |
| 孙延鹏等: "基于 OFDM 的探地雷达压缩感知成像方法", 《沈阳航空航天大学学报》 * |
| 朱江: "超宽带穿墙雷达未知墙体参数估计研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110568439A (en) * | 2019-09-25 | 2019-12-13 | 中南大学 | Impulse type through-wall radar antenna ringing suppression method based on deconvolution |
| CN110568439B (en) * | 2019-09-25 | 2023-08-22 | 中南大学 | Deconvolution-based impulse type through-wall radar antenna ringing suppression method |
| CN111398687A (en) * | 2020-03-06 | 2020-07-10 | 浙江省交通运输科学研究院 | Test method for estimating dielectric constant of asphalt pavement |
| CN112147236A (en) * | 2020-09-21 | 2020-12-29 | 大连理工大学 | Ultrasonic signal resolution improving method based on sparse blind deconvolution |
| CN114152943A (en) * | 2021-12-15 | 2022-03-08 | 电子科技大学 | Two-stage wall parameter estimation method based on ultra-wideband through-wall radar |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110031838B (en) | 2021-03-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110031838A (en) | A kind of through-wall radar wall method for parameter estimation without necessarily referring to transmitted waveform | |
| CN109799499B (en) | Wall parameter estimation method of through-wall radar | |
| CN103713288B (en) | Sparse Bayesian reconstruct linear array SAR formation method is minimized based on iteration | |
| CN107976660B (en) | Missile-borne multi-channel radar ultra-low-altitude target analysis and multi-path echo modeling method | |
| CN106872969B (en) | Radar target angle estimation method based on MTD pulse accumulation and sliding processing | |
| CN109001297B (en) | Method for measuring acoustic reflection coefficient of large-sample underwater acoustic material based on single-vector hydrophone | |
| CN109375206B (en) | Moving target speed measurement method based on speed search | |
| Hamdollahzadeh et al. | Moving target localization in bistatic forward scatter radars: Performance study and efficient estimators | |
| CN109597021A (en) | A kind of Wave arrival direction estimating method and device | |
| CN107300694A (en) | A kind of unknown wall method for parameter estimation based on Electromgnetically-transparent coefficient | |
| CN110531309A (en) | There are coherent signal angle-of- arrival estimation methods when amplitude phase error based on atom norm | |
| Cui et al. | Ice radar investigation at Dome A, East Antarctica: Ice thickness and subglacial topography | |
| CN105548959A (en) | Multi-sensor multi-target positioning method based on sparse reconstruction | |
| CN106093927B (en) | Target based on the big pulse width signal of radar tests the speed distance measuring method | |
| Crowgey et al. | Comparison of UWB short-pulse and stepped-frequency radar systems for imaging through barriers | |
| CN109752633B (en) | Method and system for positioning partial discharge signal of transformer substation | |
| CN109782249A (en) | A kind of two target associated time delays algorithm for estimating | |
| CN115932824A (en) | FMCW radar ranging method and system based on multiple antennas | |
| CN111443333B (en) | Multi-azimuth matching signal generation method based on spectrum synthesis | |
| CN111965617B (en) | GPU-based time division MIMO radar signal processing method | |
| CN110231590B (en) | Array target angle measurement method based on DFT (discrete Fourier transform) | |
| CN103308912B (en) | Coherent tracking-before-detection method based on compressive sensing | |
| Chen et al. | An improved BP algorithm for high-resolution MIMO imaging radar | |
| CN113126029A (en) | Multi-sensor pulse sound source positioning method suitable for deep sea reliable acoustic path environment | |
| CN117148273B (en) | Underwater sound source positioning method for matching probability measure |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |