基于级数反演和数值计算的广义双基SAR距离徙动成像算法

李燕平; 张振华; 邢孟道; 保铮

doi:10.3724/SP.J.1146.2007.00810

基于级数反演和数值计算的广义双基SAR距离徙动成像算法

doi: 10.3724/SP.J.1146.2007.00810

基金项目:

教育部重点项目(20051009)和新世纪优秀人才支持计划项目(NCET -06-0861)资助课题

计量
- 文章访问数: 3840
- HTML全文浏览量: 76
- PDF下载量: 895
- 被引次数: 4
出版历程
- 收稿日期: 2007-05-28
- 修回日期: 2007-10-22
- 刊出日期: 2008-12-19

A Novel Range Migration Algorithm for General Bistatic SAR Imaging Based on Series Reversion and Numerical Computation

摘要

摘要: 收发平台速度大小和方向均不相同的双基SAR是一种广义的双基SAR，具有广阔的应用前景。该文首先分析指出，对于基线较长的广义双基SAR，很难简单地将其通过一定补偿而等效为单基SAR处理，接着提出了一种基于级数反演和数值计算的距离徙动算法，最后采用该方法对仿真数据进行处理，得到了良好的成像结果。理论分析和仿真结果表明，该方法可有效应用于长基线、大斜视角的广义双基SAR成像。
- 双基合成孔径雷达 /
- 距离徙动算法\雷达成像\级数反演\数值计算
Abstract: The bistatic configuration, when transmitter and receiver move along nonparallel tracks with unequal velocities, is called the general bistatic Synthetic Aperture Radar (SAR). This kind of bistatic SAR system has the potential for various applications. In this paper, it is shown that, for the general bistatic SAR with a long baseline, it is hard to convert the bistatic SAR data into a monostatic one for which the monostatic processor can be employed. Then, a Range Migration Algorithm (RMA) based on series reversion and numerical computation is presented. Finally, the accuracy of the proposed approach is verified with a simulation. The analysis and the simulation result show that this algorithm can be applied to the general bistatic case with a long baseline and a high squint angle.
- B

HTML全文

1. 引言

有限元分析(Finite Element Analysis, FEA)是将连续介质离散成有限的单元，使无限自由度的问题转化成有限自由度的问题，再使用计算机进行计算^[1]。在进行数据分析之前，需对待分析的物体做网格剖分，网格划分的质量直接影响到计算结果的好坏。一个高质量的网格剖分，要结合计算时间和工程精度要求，设定合理的网格数量和相对网格密度^[2,3]。

有限元离散误差是评估有限元模型仿真精度的主要依据。误差估计方法可分为“前误差估计”和“后误差估计”^[4]，关于后误差估计方法前人已经做了大量的研究，如文献[5]基于极小化总势能和极大化应变能等原理建立了基本框架对离散误差进行评估。文献[6]针对线性椭圆问题用局部误差指标定义了最优网格，并给出了误差的界限。文献[7]在误差估计方面做了大量富有成效的研究，并提出采用能量范数度量两应力场间的差值表示离散误差，文献[8]证明了这种方法的有效性。文献[9]提出有限元法通常在保证收敛性时是按能量收敛。

在有限元分析中网格尺寸大小的确定是一个较大的问题。本文通过在高压输电环境中建立不同复杂度的3维人体模型；对人体模型进行自适应剖分和不断细化网格单元大小的手动剖分的研究，以能量误差(EnergyError)为判据指导模型的建立和最佳剖分尺寸的设置。

2. 网格剖分

2.1 有限元网格自适应剖分

自适应网格剖分是一种场量计算方法，在生成的初始网格的基础上，通过自动估计网格剖分误差和网格细分得到合理的网格^[10,11]。它是通过循环实现网格加密，如图1所示。

图 1 自适应网络剖分流程

下载: 全尺寸图片幻灯片

2.2 有限元离散误差的基本理论

有限元法中可用范数进行误差度量，其中能量范数的误差度量法是指对模型中的每个单元进行能量误差估计^[12]。误差分为局部误差和全局误差^[13]。每个单元的误差称为局部误差，将所有的单元误差求和即得到全局误差。

对于线弹性力学问题，令 $\varOmega$ 为有界域， $\varGamma = {\varGamma _t} \cup {\varGamma _\mu }$ 为其边界条件，则基本方程可表示为

${{L\mu }} - {{f}} = {{{S}}^{\rm{T}}}{{DS\mu }} - {{f}} = 0$

(1)

其中，

$\;\sigma = {{DS\mu}} , {\text{在}}\varOmega{\text{区内}}$

(2)

$\left. \begin{aligned} & n \cdot \sigma = t,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\text{在边界}}{\varGamma _t}{\text{上}}\\ & \mu = \mu {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} ,{\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\text{在边界}}{\varGamma _\mu }{\text{上}} \end{aligned} \right\}$

(3)

其中，L为线性微分算子， ${{L}} = {{{S}}^{\rm{T}}}{{DS}}$ ; $\mu$ 为位移；f为体积力矢量；D为弹性常数矩阵；S为应变微分矩阵； $\sigma$ 为应力。

将式(1)离散后，有限元方程为

${{K}}\mu = {{R}}$

(4)

式中，K为总刚度矩阵；R为载荷列向量。

有限元单元i的能量误差为

${{{e}}_{{i}}} = \frac{1}{2}{\int\limits_\varOmega {\left( {{\sigma ^ * } - \hat \sigma } \right)} ^{\rm{T}}}{\left[ {{D}} \right]^{ - 1}}\left( {{\sigma ^ * } - \hat \sigma } \right){\rm{d}}{\varOmega _{{i}}}$

(5)

式中， ${\varOmega _i}$ 为单元i的有界域体积；[D]表示弹性矩阵； $\left( {{\sigma ^ * } - \hat \sigma } \right)$ 为有限元应力误差， $\hat \sigma$ 为应力近似解， ${\sigma ^ * }$ 为 $\hat \sigma$ 修匀得到的改进应力解。

将所有单元能量误差相加得单元总能量误差为

$e = \sum\limits_{i = 1}^{{N_r}} {{e_i}}$

(6)

式中， ${N_r}$ 为总的单元数。

在有限元软件中，采用全域能量误差作为有限元计算控制标准，计算公式为

$E = {\left( {\frac{e}{{U + e}}} \right)^{\frac{1}{2}}} \times 100\%$

(7)

式中， $U$ 为单元总应变能； $e$ 为单元总能量误差。

本文中，以式(7)计算得的能量误差作为精度评估依据，研究输电线路有限元模型在不同网格剖分方案下的能量误差，从而选择合适的网格剖分尺寸。

3. 仿真模型的建立

有限元分析的核心思想是结构的离散化，即有限元网格剖分。它是有限元分析中至关重要的一步，直接影响计算结构的精确度^[14]。由此可见，建立高压输电线路系统模型，并进行合理的网格剖分是高压输电线路有限元分析中最重要的一步。

3.1 输电线路有限元模型

本文根据国家电网公司配电网工程典型设计中的Z1-2单回直线水泥单杆杆头示意图(其中L为0.8 m)，如图2(a)所示。采用Ansoft Maxwell仿真软件建立相应的输电线路有限元模型，如图2(b)所示。在Ansoft Maxwell中进行瞬态电场研究。

图 2 三相输电线路布局图

下载: 全尺寸图片幻灯片

在无人情况下空间电场强度仿真结果分布图如图3所示，为了方便观察将立体图的场强最大值设置为30 kV/m，剖面图场强最大值设置为10 kV/m。从图中可以看出电场强度是从导线向四周不断减小的形式变化。

图 3 无人情况下空间中电位仿真结果分布图

下载: 全尺寸图片幻灯片

3.2 人体模型的建立

人体结构复杂，其内部含有大量可能受电场作用的带电粒子，这是因为人体内部有各种生物分子。在外加电场作用下，会使人体生物组织表现出传导性和介电性^[15]。本文并不着眼于人体某个部位的具体分析，因此建立3种不同复杂程度3维简化人体模型，如图4所示。

图 4 简化3D人体模型

下载: 全尺寸图片幻灯片

其中人体模型的尺寸如表1所示。

表 1 人体模型尺寸(m)

	人体模型I	人体模型II				人体模型III
	人体模型I	脚	身体	脖子	头	脚	身体	手臂	脖子	头
半径	0.28	0.19	0.24	0.10	0.17	0.05	0.12	0.03	0.04	0.20
高度	1.71	0.70	0.65	0.08	无	0.80	0.70	0.75	0.08	无

下载: 导出CSV

| 显示表格

3.3 实际案例的有限元模型仿真

根据实际环境和对象参数，在Ansoft Maxwell中建立输电线路有限元模型。考虑到空间中人体可能带来的影响，加入人体模型进行仿真。仿真结果如图5所示。为了方便观察将立体图的场强最大值设置为10² kV/m，剖面图的场强最大值设置为10 kV/m。将仿真结果与图3对比，可知放入人体模型后空间电场的衰减比无人时更快。电场强度分布规律均由输电导线向四周不断减小，在人体附近空气场强要比无人体处的场强度大，在靠近导线一侧的人体头部尤为突出，这表明人体的加入使其周围空气中的电场发生严重畸变。这是因为人体的相对介电常数非常大，可近似为一接地导体，且感应电荷随人体与正极板距离的增加而减少。因此，在靠近导线一侧的人体头部离正极板较近，会聚集有大量感应电荷，场强增大；随着人体各个部分离正极板距离的增加，场强也会随之减小。

图 5 有人情况下空间电场强度仿真结果分布图

下载: 全尺寸图片幻灯片

4. 网格剖分优化仿真

4.1 自适应网格剖分

对3种人体模型以身高为权重求半径的方差和进行自适应网格剖分仿真分析，得到不同人体模型的方差和能量误差，如表2所示。

表 2 不同人体模型的方差和自适应网格剖分能量误差

人体模型类型	人体模型I	人体模型II	人体模型III
能量误差(%)	21.0735	22.4244	26.8348
方差(m²)	0	0.036	0.050

下载: 导出CSV

| 显示表格

根据表2，可以拟合出能量误差随着人体模型方差的变化规律，如图6所示。

图 6 不同人体模型的能量误差的变化

下载: 全尺寸图片幻灯片

从图6中可知能量误差随着人体模型方差的增大而增大。因为人体模型方差越大表示建立的人体模型越复杂，在进行有限元网格剖分时模型越复杂产生的误差就越大。因此，可以根据人体模型能量误差来指导模型的选择。

4.2 手动网格剖分

对人体模型进行不同网格剖分设置，并仿真分析。从1.2节的有限元离散误差分析中可知，最大能量误差在一定程度上体现了计算结果的精度，因此本文选择最大的能量误差作为误差大小判断依据。本节仿真分析了3种人体模型进行手动网格剖分设置后的能量误差变化。如图7所示。从图中可知，电场中的能量误差随着计算部分网格剖分单元边长的增大而增大。

图 7 3种人体模型的不同网格单元大小的能量误差的变化

下载: 全尺寸图片幻灯片

图7拟合的对数函数分别为

$\left. \begin{aligned} & {y_{\rm{I}}} = 0.928\ln (x) + 22.093 \\ & {y_{\rm{II}}} = 0.0875\ln (x) + 21.937 \\ &{y_{\rm{III}}} = 0.3549\ln (x) + 24.625 \\ \end{aligned} \right\}$

(8)

分别对3个拟合的对数函数进行求导得

$\left. \begin{aligned} & y'_{\rm{I}} = \frac{{0.928}}{x}\\ & y'_{\rm{II}} = \frac{{0.0875}}{x}\\ & y'_{\rm{III}} = \frac{{0.3549}}{x} \end{aligned} \right\}$

(9)

从图7中可知能量误差越小，则模型的计算复杂度就越大，所占用的系统资源和分析时间越大。在实际中要结合计算时间、计算机内存和工程精度要求等因素，使得计算复杂度和能量误差达到一个均衡状态。因此，令 $y' = 1$ 为最佳网格剖分点，求得人体模型I、人体模型II和人体模型III的最佳剖分尺寸分别为0.928 m, 0.0875 m和0.3549 m。

对3种模型选择相同的网格剖分尺寸，如选择网格尺寸为0.2 m和0.5 m。它们的网格数差，如表3所示。从表中看出，全局网格数变化及人体网格数变化最小的是人体模型II，其次为人体模型III，最后为人体模型I。使得人体模型II的仿真结果能量误差变化最小，其次为人体模型III，最后为人体模型I。从而求得的x为人体模型II最小，其次为人体模型III，最后为人体模型I。

表 3 网格尺寸0.2 m和0.5 m的模型网格数变化

		人体模型I	人体模型II	人体模型III
全局网格数	网格尺寸0.2 m	281086	282154	285858
全局网格数	网格尺寸0.5 m	277204	280179	283806
网格数差	全局网格数差	3882	1975	2052
人体网格数	网格尺寸0.2 m	2646	2769	4987
人体网格数	网格尺寸0.5 m	566	1579	3615
网格数差	人体网格数差	2080	1190	1372

下载: 导出CSV

| 显示表格

5. 结束语

本文在分析了有限元网格自适应剖分及有限元离散误差的基础上，采用Ansoft Maxwell仿真软件，建立3种不同复杂程度的人体模型，并对3种人体模型进行自适应剖分和手动剖分的仿真，分析了能量误差的变化趋势，研究结果可用于对人体模型的选择提供参考和对不同模型的最佳网格剖分尺寸进行估计。本文的研究成果，对其它有限元网格剖分方案的研究也具有一定的参考价值。

参考文献(1)

[1] Krieger G and Moreira A. Spaceborne bi- and multistaticSAR: Potentials and challenges[J].IEE Proc., Radar SonarNavig.2006, 153(3):184-198 [2] 汤子跃, 张守融. 双站合成孔径雷达系统原理[M]. 北京: 科学出版社, 2003: 1-9. [3] Ender J, Walterscheid I, and Brenner A. Bistatic SARtranslationalinvariant processing and experimental results[J].IEE Proc.- Radar Sonar Navig.2006, 153(3):177-183 [4] Soumekh M. Bistatic synthetic aperture radar inverison withapplication in dynamic object imaging[J].IEEE Trans. onSignal Proccess.1991, 39(9):2044-2055 [5] 何峰, 梁甸农, 董臻. 适于大斜视角的星载双基SAR 波数域成像算法[J]. 电子学报, 2005, 33(6): 1011-1014.He Feng, Liang Dian-nong, and Dong Zhen. A wavenumberdomain algorithm for spaceborne bistatic SAR imaging withlarge squint angle[J]. Acta Electronica Sinica, 2005, 33(6):1011-1014. [6] 闫鸿慧, 王岩飞, 于海锋等. 一种基于距离补偿的分布式小卫星双基SAR 成像方法[J].电子与信息学报.2005, 27(5):771-774浏览 [7] Walterscheid I, Ender J, and Brenner A, et al.. Bistatic SARprocessing and experiments[J].IEEE Trans. on Geosc. andRemote Sens.2006, 44(10):2710-2717 [8] 朱振波, 汤子跃, 蒋兴舟. 机载双站合成孔径雷达ChirpScaling 成像算法[J].电子与信息学报.2006, 28(6):977-981浏览 [9] Barber B. Theory of digital imaging from orbital syntheticaperture radar[J].Int. J. Remote Sens.1985, 6(6):1009-1057 [10] Cheng Hu and Tao Zeng. Forward-looking bistatic SAR rangemigration algorithm[A]. Proceedings of 2006 CIEinternational conference on radar[C]. Shanghai: IEEE, 2006:724-727. [11] Loffeld O, Nies H, Peters V, and Knedlik S. Models anduseful relations for bistatic SAR processing[J].IEEE Trans.on Geosc. and Remote Sens.2004, 42(10):2031-2038 [12] Lam Neo Y, Wong F G, and Cumming I. A two-dimensionalspectrum for bistatic SAR processing using series reversion[J].IEEE Trans. on Geosc. and Remote Sens.2007, 4(1):93-96 [13] Pettersson M L. Extraction of moving ground targets by abistatic ultra-wideband SAR[J].IEE Proc.- Radar SonarNavi.2001, 148(1):35-40 [14] Morse P M and Feshbach H. Methods of theoreticalphysics[M]. New York: McGraw-Hill, 1953: 61-68.

施引文献

期刊类型引用(2)

1.	程姣，王利刚，高文斌，陶婷. 基于元分析模型的师生动机感染仿真研究. 现代电子技术. 2022(15): 109-112 . 百度学术
2.	闻小龙，杨鹏飞，储昭志，彭春荣，刘宇涛，吴双. 基于MEMS的距离自适应型非接触静电仪. 电子与信息学报. 2021(10): 3068-3074 . 本站查看

其他类型引用(2)

资源附件(0)

访问统计

计量

文章访问数: 3840
HTML全文浏览量: 76
PDF下载量: 895
被引次数: 4

1. 引言
2. 网格剖分
2.1 有限元网格自适应剖分
2.2 有限元离散误差的基本理论
3. 仿真模型的建立
3.1 输电线路有限元模型
3.2 人体模型的建立
3.3 实际案例的有限元模型仿真
4. 网格剖分优化仿真
4.1 自适应网格剖分
4.2 手动网格剖分
5. 结束语

1. 引言
2. 网格剖分
2.1 有限元网格自适应剖分
2.2 有限元离散误差的基本理论
3. 仿真模型的建立
3.1 输电线路有限元模型
3.2 人体模型的建立
3.3 实际案例的有限元模型仿真
4. 网格剖分优化仿真
4.1 自适应网格剖分
4.2 手动网格剖分
5. 结束语

参考文献(1)

施引文献

资源附件(0)

访问统计

姓名
邮箱
手机号码
标题
留言内容
验证码

留言板

基于级数反演和数值计算的广义双基SAR距离徙动成像算法

doi: 10.3724/SP.J.1146.2007.00810

计量

A Novel Range Migration Algorithm for General Bistatic SAR Imaging Based on Series Reversion and Numerical Computation

1. 引言

2. 网格剖分

2.1 有限元网格自适应剖分

2.2 有限元离散误差的基本理论

3. 仿真模型的建立

3.1 输电线路有限元模型

3.2 人体模型的建立

3.3 实际案例的有限元模型仿真

4. 网格剖分优化仿真

4.1 自适应网格剖分

4.2 手动网格剖分

5. 结束语

期刊类型引用(2)

其他类型引用(2)

计量

目录

1. 引言

2. 网格剖分

2.1 有限元网格自适应剖分

2.2 有限元离散误差的基本理论

3. 仿真模型的建立

3.1 输电线路有限元模型

3.2 人体模型的建立

3.3 实际案例的有限元模型仿真

4. 网格剖分优化仿真

4.1 自适应网格剖分

4.2 手动网格剖分

5. 结束语

留言板

基于级数反演和数值计算的广义双基SAR距离徙动成像算法

doi: 10.3724/SP.J.1146.2007.00810

计量

出版历程

A Novel Range Migration Algorithm for General Bistatic SAR Imaging Based on Series Reversion and Numerical Computation

1. 引言

2. 网格剖分

2.1 有限元网格自适应剖分

2.2 有限元离散误差的基本理论

3. 仿真模型的建立

3.1 输电线路有限元模型

3.2 人体模型的建立

3.3 实际案例的有限元模型仿真

4. 网格剖分优化仿真

4.1 自适应网格剖分

4.2 手动网格剖分

5. 结束语

期刊类型引用(2)

其他类型引用(2)

计量

出版历程

目录

1. 引言

2. 网格剖分

2.1 有限元网格自适应剖分

2.2 有限元离散误差的基本理论

3. 仿真模型的建立

3.1 输电线路有限元模型

3.2 人体模型的建立

3.3 实际案例的有限元模型仿真

4. 网格剖分优化仿真

4.1 自适应网格剖分

4.2 手动网格剖分

5. 结束语