CN107271981B - 一种目标电磁散射数据模型的建立方法 - Google Patents
一种目标电磁散射数据模型的建立方法 Download PDFInfo
- Publication number
- CN107271981B CN107271981B CN201710637865.2A CN201710637865A CN107271981B CN 107271981 B CN107271981 B CN 107271981B CN 201710637865 A CN201710637865 A CN 201710637865A CN 107271981 B CN107271981 B CN 107271981B
- Authority
- CN
- China
- Prior art keywords
- wavelet
- target
- scale
- coefficient
- electromagnetic scattering
- 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013499 data model Methods 0.000 title claims abstract description 9
- 238000004458 analytical method Methods 0.000 claims abstract description 24
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000013598 vector Substances 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 238000004088 simulation Methods 0.000 abstract description 4
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
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
- 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种目标电磁散射数据模型的建立方法,包括:根据目标电磁散射的实测数据进行多重分形理论分析,初步建立目标电磁散射多重分形建模,对数据进行小波分解和重建,比较重建数据与实测数据的分形结果并对小波分析的维数进行调整,以精确建立目标电磁散射模型。本发明能利用目标电磁散射的实测数据建立散射模型,利用仿真方法反演出目标电磁散射的数据。
Description
技术领域
本发明涉及一种目标电磁散射数据模型的建立方法。
背景技术
当物体被电磁波照射时,能量朝各个方向散射,散射能量的空间分布取决于物体的形状、大小、结构、入射波的频率、极化等,产生电磁散射的物体通常称为目标或散射体。在复杂多变的电磁散射环境中,由于被测目标位置和姿态不断地变化,需要时刻对目标的RCS特性、速度及加速度变化、角变化、目标振动所引起的回波幅度起伏等进行模拟,这就需要对目标瞬时电磁散射特性进行实时模拟,以满足雷达信号处理机系统的研制与开发需要。
工程上可以通过实测和建模仿真获取目标电磁散射的相关数据,但目标的电磁散射特性会影响雷达的回波信号,实际情况下目标的散射点数目有限且散射强度是视角的函数,采用经典的随机模型不足以反映目标散射的物理本质和复杂性,且实测结果费用高,受实际条件的限制较多,很难得到完备的目标电磁散射数据。
发明内容
本发明的目的在于克服现有技术中的不足,提供一种目标电磁散射数据模型的建立方法,解决现有技术中经典随机模型不足以反应目标散射的物力本质和复杂性、实测结果费用高、目标电磁散射数据精度不高的技术问题。
为解决上述技术问题,本发明所采用的技术方案是:一种目标电磁散射数据模型的建立方法,包括如下步骤:
(1):根据目标电磁散射的实测数据初步建立目标电磁散射的多重分形模型,确定小波分解级数、计算分形参数;
(2)对实测数据进行小波分解,计算每个尺度下的小波系数和尺度系数;
(3)根据小波系数和尺度系数进行小波重建;
(4)比较重建数据与实测数据的分形结果,并对小波分析的分形维数进行调整,以精确建立目标电磁散射模型。
计算分形参数的具体方法如下:
定义目标电磁散射系数的测量值为
其中:x(M,N)表示目标电磁散射系数在(M,N)处的测量值,M,N则分别表示目标电磁散射系数测量值的行数和列数;
设空间内存在一正方柱,所述正方柱的底边宽度为ε,正方柱的高度N(ε)由其范围内目标电磁散射幅度的最高点和最低点之差决定,具体定义为:
N(ε)=max|x(m,n)-x(m,,n,)| (2)
其中,点(m,n)和(m,,n,)都在底边为ε的正方形内,改变ε的大小逐步移动正方柱遍及所有目标散射点;
用分形维数来描述局部与整体的自相似程度,其分形维数定义为:
D=ln(N(ε))/ln(1/ε) (3)。
计算每个尺度下的小波系数和尺度系数的方法如下:
对目标的电磁散射系数的测量值x进行行变换,得到一维序列y:
y=[x(1,1) x(1,2) … x(1,N) x(2,1) x(2,2) … x(2,N) … x(M,1) x(M,2)… x(M,N)] (4)
小波函数φ(t)表示为:
利用dj[k]=<x,φj,k(t)>和式(9),小波系数dj[k]表示为:
式(10)和(11)称为小波变换的分解方程,尺度j的小波展开系数cj[k]和dj[k]由尺度j+1的展开系数cj+1[k]经过滤波器组h0[-n]和h1[-n]后再抽取而得到;同理,通过cj[k]按此过程进一步得到展开系数cj-1[k]和dj-1[k],如此重复,最后由cj+1[k]经过滤波器组得到一系列不同尺度的小波展开系数。
根据小波系数和尺度系数进行小波重建的具体方法如下:
小波逆变换的重建方程可以表示为:
进行分形维数调整的具体方法如下:
规定常数J0为满足上述多重分形模型反演目标电磁散射的数据平方误差积分的最大值,则有:
J(ε)≤J0 (18)
当满足式(18)时,说明该建模过程成立,若(18)不成立,则改变ε的大小,重复上述过程,直到式(18)成立为止。
所述目标电磁散射实测数据来源于用逆合成孔径雷达分别对目标成像所获得目标各散射点的散射系数,所述散射系数是由实部和虚部组成的矢量。
与现有技术相比,本发明所达到的有益效果是:
能利用目标电磁散射的实测数据建立散射模型,利用仿真方法反演出目标电磁散射的数据,利用分形技术的优点准确地建立目标电磁散射模型。
附图说明
图1是目标的电磁散射系数强度的示意图;
图2是本发明的流程图。
具体实施方式
本发明是一种利用多重分形理论和小波分析对目标电磁散射特性进行建模的方法,包括:根据目标电磁散射的实测数据进行多重分形理论分析,初步建立目标电磁散射多重分形建模,对数据进行小波分解和重建,比较重建数据与实测数据的分形结果并对小波分析的维数进行调整,以精确建立目标电磁散射模型,具体为:计算正方柱底边的正方形为ε时正方柱的高度N(ε)和分形维数D,其中D=ln(N(ε))/ln(1/ε),改变ε的大小逐步移动正方柱,遍及所有目标散射点;分析小波变换的分解方程和得到小波展开系数cj[k]和dj[k];利用小波逆变换的重建方程对实测目标散射数据进行数据重建;然后对重建数据进行多重分形分析和相应分形参数的计算;通过分析实测数据的多重分形参数和重建数据的多重分形参数的平方误差积分准则逐步调整ε值的大小以达到最精确的。本发明能很好地利用目标电磁散射的实测数据建立散射模型,利用仿真方法反演出目标电磁散射的数据。
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
目标散射实测数据来源于用逆合成孔径雷达分别对目标成像所获得目标各散射点的散射系数,散射系数是由实部和虚部组成的矢量。图1是目标的散射系数强度的示意图,从图中可以看出实际的复杂分布式目标由较少的强散射点和较多的弱散射点以及大量的散射强度很小的散射点组成。首先分析实测数据的多重分形特征,然后根据分形特征参数建立合理、有效的目标散射多重分形模型,利用目标散射多重分形模型反演出目标的散射系数。
对于如图1所示的目标电磁散射数据,定义目标电磁散射系数的测量值为
其中:x(M,N)表示目标电磁散射系数在(M,N)处的测量值,M,N则分别表示目标电磁散射系数测量值的行数和列数;
步骤一:设空间内存在一正方柱,该正方柱的底边宽度为ε,正方柱的高度N(ε)由其范围内目标电磁散射幅度的最高点和最低点之差决定,具体定义为:
N(ε)=max|x(m,n)-x(m,,n,)| (2)
其中,点(m,n)和(m,,n,)都在底边为ε的正方形内,改变ε的大小逐步移动正方柱遍及所有目标散射点计算。
可以用分形维数来描述局部与整体的自相似程度,其分形维数定义为:
D=ln(N(ε))/ln(1/ε) (3)
分形在区间上满足相似性且计算出的分维数具有一定的稳定性。工程上用ln(N(ε))和ln(1/ε)的拟合曲线的线性度来反映分形的自相似性程度,所以作ln(N(ε))和ln(1/ε)的拟合曲线,则分形维数D的物理意义就是拟合曲线的线性部分的斜率。
目标在雷达的不同观测通道上的电磁散射数据具有自相似性,为利用实测数据重建各观测通道上的目标电磁散射数据模型提供了理论基础,具体流程如图2所示,产生雷达各观测通道上的目标电磁散射数据。
步骤二:对目标的电磁散射系数模值x进行行变换,得到1维序列y:
y=[x(1,1) x(1,2) … x(1,n) x(2,1) x(2,2) … x(2,N) … x(M,1) x(M,2)… x(M,N)](1×MN) (4)
式(5)中,尺度系数为cj[k],小波系数为dj[k],j为尺度下标,j越大,分辨率越高;k为小波变换点数;t为时间参数;
小波函数φ(t)可以表示为:
利用和式(8),系数cj[k]可表示为:
利用dj[k]=<x,φj,k(t)>和式(9),系数dj[k]可表示为:
式(10)和(11)称为小波变换的分解方程,尺度j的小波展开系数cj[k]和dj[k]都可以由尺度j+1的展开系数cj+1[k]经过滤波器组h0[-l]和h1[-l]后再抽取而得到。同理,可以通过cj[k]按此过程进一步得到展开系数cj-1[k]和dj-1[k],如此重复,最后可以由cj+1[k]经过滤波器组得到一系列不同尺度的小波展开系数,小波变换的展开级数为N。
小波重建算法为分解算法的逆过程,即由展开系数cj[k]和dj[k]重建展开系数cj+1[k],从而实现小波逆变换。
步骤三:小波逆变换的重建方程可以表示为:
步骤五:将实测数据的多重分形参数D和重建数据的多重分形参数进行平方误差积分准则分析,设分析的结果J(ε)为:
分析可知,J(ε)的值越小,说明重建目标电磁散射数据模型的结果越精准,不完全分集MIMO雷达目标检测效果也就越好。规定常数J0为满足上述多重分形模型反演目标电磁散射的数据平方误差积分的最大值,则有:
J(ε)≤J0 (18)
当满足式(18)时,说明该建模过程成立,若(18)不成立,则返回到第一步,改变ε的大小,重复上述过程,直到式(18)成立为止。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (2)
1.一种目标电磁散射数据模型的建立方法,其特征在于,包括如下步骤:
(1):根据目标电磁散射的实测数据初步建立目标电磁散射的多重分形模型,确定小波分解级数、计算分形参数;
(2)对实测数据进行小波分解,计算每个尺度下的小波系数和尺度系数;
(3)根据小波系数和尺度系数进行小波重建;
(4)比较重建数据与实测数据的分形结果,并对小波分析的分形维数进行调整,以精确建立目标电磁散射模型;
计算分形参数的具体方法如下:
定义目标电磁散射系数的测量值为
其中:x(M,N)表示目标电磁散射系数在(M,N)处的测量值,M,N则分别表示目标电磁散射系数测量值的行数和列数;
设空间内存在一正方柱,所述正方柱的底边宽度为ε,正方柱的高度N(ε)由其范围内目标电磁散射幅度的最高点和最低点之差决定,具体定义为:
N(ε)=max|x(m,n)-x(m’,n’)| (2)
其中,点(m,n)和(m’,n’)都在底边为ε的正方形内,改变ε的大小逐步移动正方柱遍及所有目标散射点;
用分形维数来描述局部与整体的自相似程度,其分形维数定义为:
D=ln(N(ε))/ln(1/ε) (3);
计算每个尺度下的小波系数和尺度系数的方法如下:
对目标的电磁散射系数的测量值x进行行变换,得到一维序列y:
y=[x(1,1)x(1,2)…x(1,N)x(2,1)x(2,2)…x(2,N)…x(M,1)x(M,2)…x(M,N)](4)
小波函数φ(t)表示为:
利用dj[k]=<x,φj,k(t)>和式(9),小波系数dj[k]表示为:
式(10)和(11)称为小波变换的分解方程,尺度j的小波展开系数cj[k]和dj[k]由尺度j+1的展开系数cj+1[k]经过滤波器组h0[-n]和h1[-n]后再抽取而得到;同理,通过cj[k]按此过程进一步得到展开系数cj-1[k]和dj-1[k],如此重复,最后由cj+1[k]经过滤波器组得到一系列不同尺度的小波展开系数;
根据小波系数和尺度系数进行小波重建的具体方法如下:
小波逆变换的重建方程可以表示为:
重建后的信号为
进行分形维数调整的具体方法如下:
规定常数J0为满足上述多重分形模型反演目标电磁散射的数据平方误差积分的最大值,则有:
J(ε)≤J0 (18)
当满足式(18)时,说明该建模过程成立,若(18)不成立,则改变ε的大小,重复上述过程,直到式(18)成立为止。
2.根据权利要求1所述的目标电磁散射数据模型的建立方法,其特征在于,所述目标电磁散射实测数据来源于用逆合成孔径雷达分别对目标成像所获得目标各散射点的散射系数,所述散射系数是由实部和虚部组成的矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710637865.2A CN107271981B (zh) | 2017-07-31 | 2017-07-31 | 一种目标电磁散射数据模型的建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710637865.2A CN107271981B (zh) | 2017-07-31 | 2017-07-31 | 一种目标电磁散射数据模型的建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107271981A CN107271981A (zh) | 2017-10-20 |
CN107271981B true CN107271981B (zh) | 2020-02-21 |
Family
ID=60075921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710637865.2A Active CN107271981B (zh) | 2017-07-31 | 2017-07-31 | 一种目标电磁散射数据模型的建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107271981B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109471084B (zh) * | 2018-11-23 | 2022-10-25 | 西北工业大学 | 一种电磁散射机理分解方法 |
CN111859704B (zh) * | 2020-07-30 | 2024-06-18 | 上海无线电设备研究所 | 一种分布式多视角下非刚体目标电磁散射建模方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6847925B2 (en) * | 2000-06-21 | 2005-01-25 | Hrl Laboratories, Llc | Method and apparatus for modeling three-dimensional electromagnetic scattering from arbitrarily shaped three-dimensional objects |
CN102999908A (zh) * | 2012-11-19 | 2013-03-27 | 西安电子科技大学 | 基于改进视觉注意模型的sar机场分割方法 |
CN104573289A (zh) * | 2015-02-06 | 2015-04-29 | 西安电子科技大学 | 一维导体粗糙海面与二维导体目标的电磁散射仿真方法 |
CN106650048A (zh) * | 2016-12-05 | 2017-05-10 | 中国舰船研究设计中心 | 一种基于斜率分布的船海互耦散射预测方法 |
CN106872978A (zh) * | 2017-01-16 | 2017-06-20 | 西北工业大学 | 一种复杂场景的电磁建模仿真方法 |
-
2017
- 2017-07-31 CN CN201710637865.2A patent/CN107271981B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107271981A (zh) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Franceschetti et al. | On ocean SAR raw signal simulation | |
Zhang et al. | Resolution enhancement for large-scale real beam mapping based on adaptive low-rank approximation | |
Zhang et al. | A TV forward-looking super-resolution imaging method based on TSVD strategy for scanning radar | |
CN111965632B (zh) | 一种基于黎曼流形降维的雷达目标检测方法 | |
Nouguier et al. | Nonlinear ocean wave reconstruction algorithms based on simulated spatiotemporal data acquired by a flash LIDAR camera | |
CN101915943B (zh) | 均匀背景介质的介电常数和隐蔽目标参数的联合反演方法 | |
Knospe et al. | Covariance estimation for dInSAR surface deformation measurements in the presence of anisotropic atmospheric noise | |
CN112906300B (zh) | 基于双通道卷积神经网络的极化sar土壤湿度反演方法 | |
CN109683134B (zh) | 一种面向旋转声源的高分辨率定位方法 | |
CN107390216A (zh) | 基于波数域相干因子的高速超分辨率驻点扫描成像方法 | |
CN108008389A (zh) | 一种基于gpu的快速频域后向投影三维成像方法 | |
CN107271981B (zh) | 一种目标电磁散射数据模型的建立方法 | |
CN104698462A (zh) | 基于变分的合成孔径雷达海面风场融合方法 | |
CN102944872A (zh) | 雷达散射截面近场到近场的变换方法 | |
CN112444773A (zh) | 基于空域融合的压缩感知二维doa估计方法 | |
CN110232342A (zh) | 基于卷积神经网络的海况等级判断方法及装置 | |
CN114428249B (zh) | 非高斯海面sar成像方法 | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
CN113608218A (zh) | 一种基于后向投影原理的频域干涉相位稀疏重构方法 | |
CN117436329A (zh) | 一种基于机器学习的合成孔径雷达图像风向反演方法 | |
Johnson et al. | Surface shape reconstruction from phaseless scattered acoustic data using a random forest algorithm | |
Refice et al. | On the use of anisotropic covariance models in estimating atmospheric DInSAR contributions | |
Wang et al. | Interference image registration combined by enhanced scale-invariant feature transform characteristics and correlation coefficient | |
Wan et al. | Azimuth cutoff compensation method for SAR wave observation based on multiview wave spectrum data fusion | |
CN104167019B (zh) | 一种sar图像精确度计算方法 |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231123 Address after: Room 6-231, Tengfei Building, No. 88 Jiangmiao Road, Research and Innovation Park, Nanjing Area, China (Jiangsu) Pilot Free Trade Zone, Nanjing, Jiangsu Province, 210000 (Information Application) Patentee after: Jiangsu Haozhifeng High tech Co.,Ltd. Address before: 211100 No. 8 West Buddha Road, Jiangning District, Jiangsu, Nanjing Patentee before: HOHAI University |
|
TR01 | Transfer of patent right |