三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究

何姿 陈如山

何姿, 陈如山. 三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究[J]. 雷达学报, 2019, 8(3): 318–325. doi:  10.12000/JR19032
引用本文: 何姿, 陈如山. 三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究[J]. 雷达学报, 2019, 8(3): 318–325. doi:  10.12000/JR19032
HE Zi and CHEN Rushan. Study on efficient high frequency method of electromagnetic scattering from ships above sea surface[J]. Journal of Radars, 2019, 8(3): 318–325. doi:  10.12000/JR19032
Citation: HE Zi and CHEN Rushan. Study on efficient high frequency method of electromagnetic scattering from ships above sea surface[J]. Journal of Radars, 2019, 8(3): 318–325. doi:  10.12000/JR19032

三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究

doi: 10.12000/JR19032
基金项目: 国家自然科学基金(61701232, 61431006),江苏省自然科学基金(BKs20170854),中国科协青年人才托举工程(2017QNRC001),中国博士后基金(2017M620861, 2018T110127)
详细信息
    作者简介:

    何 姿(1988–),女,博士,讲师。2016年6月在南京理工大学获得博士学位,2015年9月到2016年9月在美国伊利诺伊大学香槟分校(UIUC)交流访问。主要研究方向为计算电磁学、特征提取、目标识别,近年来在学术期刊、会议上发表论文60余篇,其中被SCI收录期刊论文30余篇。博士论文《目标电磁散射特性的抛物线方程方法研究》荣获2016年度中国电子学会优秀博士学位论文优秀奖;2017年度成功入选中国科协信息科技学会联合体青年人才托举工程;2018年度荣获International Applied Computational Electromagnetics Society (ACES) Symposium国际会议组颁发的“Young Scientist Award Scheme”

    陈如山(1965–),男,博士、教授、博士生导师,江苏省政协委员,教育部“长江学者”特聘教授,国家杰出青年科学基金获得者,享受国务院政府特殊津贴,国务院学位委员会学科评议组“电子科学与技术组”成员,“新世纪百千万人才工程”国家级人选,“十二五”国防基础研究专家组成员,教育部“长江学者和创新团队发展计划”创新团队带头人,国家重点学科“电磁场与微波技术”学科带头人,江苏省“青蓝工程”科技创新团队带头人,国际电子学报副主编,现为南京理工大学教授,博士生导师,电子工程与光电技术学院副院长。在电磁脉冲与瞬态电磁场、电磁兼容与雷达隐身技术、计算电磁学等方面发表论文160多篇,科学引文索引(SCI)收录论文140篇,(EI)收录论文120篇。曾获2000年度工程类科学引文索引(SCI)收录论文数为全国第一,南京市自然科学优秀学术论文三等奖,江苏省科技进步三等奖,国家教委科技进步三等奖,国家教委科技进步二等奖,江苏省国防工办科技进步一等奖,兵器工业集团公司科技进步二等奖,军队科技进步奖二等奖2项、三等奖1项。2015年江苏省研究生培养模式改革成果一等奖。目前已获得国防发明专利14项,实用新型发明专利8项

    通讯作者:

    陈如山 eerschen@njust.edu.cn

  • 中图分类号: TN011

Study on Efficient High Frequency Method of Electromagnetic Scattering from Ships above Sea Surface

Funds: The National Natural Science Foundation of China (61701232, 61431006), Jiangsu Province Natural Science Foundation (BKs20170854), The Young Elite Scientists Sponsorship Program (2017QNRC001), China Postdoctoral Science Foundation (2017M620861, 2018T110127)
More Information
    Corresponding author: CHEN Rushan, eerschen@njust.edu.cn
  • 摘要: 3维随机粗糙海面与其上方复杂目标复合电磁(EM)散射特性的建模与分析在微波遥感、目标识别、雷达成像、导弹制导等领域中有着重要的研究价值。该文主要研究了基于高频算法的随机粗糙海面及舰船的复合电磁散射特性,开发了PO-IPO混合方法,为3维随机粗糙海面与复杂目标一体化高效求解提供了新思路。文中分别使用了物理光学方法(PO)、迭代物理光学方法(IPO)、PO-PO以及PO-IPO混合方法对海面及舰船进行了建模与仿真,其中,引入锥形波来代替平面波作为发射源,锥形波可以更好地抑制粗糙面在边缘位置被突然截断而形成的电磁反射和边缘绕射等效应。从数值仿真结果中可以看出,PO-IPO混合方法可将复杂物体本身面元间以及粗糙海面与物体间的耦合作用考虑在内,因此PO-IPO可以作为一种有效的途径来快速获取随机粗糙海面及舰船的复合电磁散射特性。
  • 图  1  多路径模型基本原理

    Figure  1.  Theory of multipath model

    图  2  PO-IPO算法流程图

    Figure  2.  Flow chart for PO-IPO

    图  3  PM谱粗糙面模型

    Figure  3.  Model of PM rough surface

    图  4  入射角度为${\theta _{\rm i}}$=30°粗糙海面模型的双站RCS结果对比

    Figure  4.  Bistatic RCS results for the PM rough surface when ${\theta _{\rm i}}$=30°

    图  5  入射角度为${\theta _{\rm i}}$=60°粗糙海面模型的双站RCS结果对比

    Figure  5.  Bistatic RCS results for the PM rough surface when ${\theta _{\rm i}}$=60°

    图  6  舰船模型示意图

    Figure  6.  Model of ship

    图  7  舰船模型的双站RCS结果对比

    Figure  7.  Bistatic RCS results for the ship model

    图  8  粗糙海面与缩比舰船组合模型示意图

    Figure  8.  Model of ship over rough sea surface

    图  9  粗糙海面与缩比舰船组合模型的双站RCS结果对比

    Figure  9.  Bistatic RCS results for the ship model over rough sea surface

    图  10  粗糙海面与舰船组合模型示意图

    Figure  10.  Model of full-size ship over rough sea surface

    图  11  粗糙海面与舰船组合模型的双站RCS结果对比

    Figure  11.  Bistatic RCS results for the full-size ship model over rough sea surface

    表  1  PO与JMCFIE计算资源对比

    Table  1.   Comparison of computational resources between PO and JMCFIE

    数值方法时间(s)内存(MB)
    SBR27176
    PO23174
    JMCFIE19565912
    下载: 导出CSV

    表  2  SBR, IPO与MoM-MLFMA计算资源对比

    Table  2.   Comparison of computational resources among SBR, IPO and MoM-MLFMA

    数值方法时间(s)内存(MB)
    SBR1558
    IPO259374
    MoM-MLFMA1125701178
    下载: 导出CSV

    表  3  PO-IPO与JMCFIE计算资源对比

    Table  3.   Comparison of computational resources between PO-IPO and JMCFIE

    数值方法时间(s)内存(MB)
    PO-IPO471.842.5
    JMCFIE8994.021645.0
    下载: 导出CSV
  • [1] THORSOS E I. The validity of the Kirchhoff approximation for rough surface scattering using a Gaussian roughness spectrum[J]. The Journal of the Acoustical Society of America, 1988, 83(1): 78–92. doi:  10.1121/1.396188
    [2] LEADER J. The relationship between the Kirchhoff approach and small perturbation analysis in rough surface scattering theory[J]. IEEE Transactions on Antennas and Propagation, 1971, 19(6): 786–788. doi:  10.1109/TAP.1971.1140044
    [3] BASS F G and FUKS I M. Wave Scattering from Statistically Rough Surfaces[M]. Oxford: Pergamon Press, 1979.
    [4] FUKS I M. Wave diffraction by a rough boundary of an arbitrary plane-layered medium[J]. IEEE Transactions on Antennas and Propagation, 2001, 49(4): 630–639. doi:  10.1109/8.923325
    [5] 魏鹏博, 江旺强, 张晓梨, 等. 基于小斜率近似方法的沥青路面电磁散射[C]. 第十二届全国电波传播学术讨论会论文集, 青岛, 2013.WEI Pengbo, JIANG Wangqiang, ZHANG Xiaoli, et al. Electromagnetic L scattering of asphalt pavement based on small slope approximation method[C]. National Symposium on Radio Communication, Qingdao, 2013.
    [6] JOHNSON J T, SHIN R T, KONG Jinau, et al. A numerical study of the composite surface model for ocean backscattering[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(1): 72–83. doi:  10.1109/36.655319
    [7] CHEN M F and FUNG A K. A numerical study of the regions of validity of the Kirchhoff and small-perturbation rough surface scattering models[J]. Radio Science, 1988, 23(2): 163–170. doi:  10.1029/RS023i002p00163
    [8] THORSOS E I and JACKSON D R. The validity of the perturbation approximation for rough surface scattering using a Gaussian roughness spectrum[J]. The Journal of the Acoustical Society of America, 1989, 86(1): 261–277. doi:  10.1121/1.398342
    [9] JOHNSON J T. A numerical study of scattering from an object above a rough surface[J]. IEEE Transactions on Antennas and Propagation, 2002, 50(10): 1361–1367. doi:  10.1109/TAP.2002.802152
    [10] VORONOVICH A. Small-slope approximation for electromagnetic wave scattering at a rough interface of two dielectric half-spaces[J]. Waves in Random Media, 1994, 4(3): 337–367. doi:  10.1088/0959-7174/4/3/008
    [11] NIETO-VESPERINAS M. Depolarization of electromagnetic waves scattered from slightly rough random surfaces: A study by means of the extinction theorem[J]. Journal of the Optical Society of America, 1982, 72(5): 539–547. doi:  10.1364/JOSA.72.000539
    [12] BAHAR E and LEE B S. Radar scatter cross sections for two-dimensional random rough surfaces - full wave solutions and comparisons with experiments[J]. Waves in Random Media, 1996, 6(1): 1–23. doi:  10.1080/13616679609409792
    [13] WINEBRENNER D P and ISHIMARU A. Application of the phase-perturbation technique to randomly rough surfaces[J]. Journal of the Optical Society of America A, 1985, 2(12): 2285–2294. doi:  10.1364/JOSAA.2.002285
    [14] LIANG Y, ZENG X H, GUO L X, et al. Review on the pre-study about specific target and rough soil surface composite electromagnetic scattering basing on an efficient numerical strategy[C]. 2016 Progress in Electromagnetic Research Symposium (PIERS), Shanghai, China, 2016: 3198.
    [15] 叶红霞. 随机粗糙面与目标复合电磁散射的数值计算方法[D]. [博士论文], 复旦大学, 2007.YE Hongxia. Numerical methods for the composite electromagnetic scattering of random rough surface and targets[D]. [Ph.D. dissertation], Fudan University, 2007.
    [16] 朱国强, 孙劲, 郑立志, 等. 平板目标与随机粗糙面对电磁波的复合散射[J]. 武汉大学学报(自然科学版), 2000, 46(1): 99–103.ZHU Guoqiang, SUN Jin, ZHENG Lizhi, et al. Composite scattering of electromagnetic waves by the plate target and random rough surface[J]. Journal of Wuhan University (Natural Science Edition), 2000, 46(1): 99–103.
    [17] 杨伟. 三维复杂粗糙海面电磁散射建模研究与特性分析[D]. [博士论文], 电子科技大学, 2012.YANG Wei. Charecteristics analysis and modeling research on electromagnetic scattering from 3D complicated rough sea surface[D]. [Ph.D. dissertation], University of Electronic Science and Technology of China, 2012.
    [18] 王强, 郭立新, 刘忠玉. 粗糙海面与上方多目标瞬态电磁散射混合算法[J]. 西安电子科技大学学报, 2018, 45(5): 57–63. doi:  10.3969/j.issn.1001-2400.2018.05.010WANG Qiang, GUO Lixin, and LIU Zhongyu. Hybrid method for transient electromagnetic scattering from multiple targets above the rough sea surface[J]. Journal of Xidian University, 2018, 45(5): 57–63. doi:  10.3969/j.issn.1001-2400.2018.05.010
    [19] 叶红霞, 金亚秋. 三维介质粗糙面上导体目标散射的解析-数值混合算法[J]. 电波科学学报, 2008, 23(6): 1144–1153, 1187. doi:  10.3969/j.issn.1005-0388.2008.06.025YE Hongxia and JIN Yaqiu. A hybrid analytical-numerical algorithm for computation of scattering from a 3D PEC target above a dielectric rough surface[J]. Chinese Journal of Radio Science, 2008, 23(6): 1144–1153, 1187. doi:  10.3969/j.issn.1005-0388.2008.06.025
    [20] 张民, 郭立新, 聂丁, 等. 海面目标雷达散射特性与电磁成像[M]. 北京: 科学出版社, 2015.ZHANG Min, GUO Lixin, NIE Ding, et al. Radar Scattering Characteristics and Electromagnetic Imaging of Targets on the Sea Surface[M]. Beijing: Science Press, 2015.
    [21] 金亚秋, 刘鹏, 叶红霞. 随机粗糙面与目标复合散射数值模拟理论与方法[M]. 北京: 科学出版社, 2008.JIN Yaqiu, LIU Peng, and YE Hongxia. Theory and Method of Numerical Simulation of Composite Scattering from the Object and Randomly Rough Surface[M]. Beijing: Science Press, 2008.
    [22] 郭立新, 王蕊, 吴振森. 随机粗糙面散射的基本理论与方法[M]. 北京: 科学出版社, 2010.GUO Lixin, WANG Rui, and WU Zhensen. Basics Theory and Method for the Scattering from the Random Rough Surface[M]. Beijing: Science Press, 2010.
    [23] 康士峰, 王显德. 粗糙面与目标电磁散射统计特性分析[J]. 微波学报, 2004, 20(3): 43–46. doi:  10.3969/j.issn.1005-6122.2004.03.010KANG Shifeng and WANG Xiande. Analysis of statistical characteristics for EM scattering from rough surface and target[J]. Journal of Microwaves, 2004, 20(3): 43–46. doi:  10.3969/j.issn.1005-6122.2004.03.010
    [24] ZHU G Q, CAO Q F, and YANG H L. Analysis of double bounce scattering between target and rough surface[J]. Chinese Journal of Radio Science, 1996, 11(1): 80–84.
    [25] 邹高翔, 童创明, 高飞, 等. 复杂陆地粗糙面及其上方坦克目标复合散射研究[J]. 电波科学学报, 2017, 32(3): 261–272. doi:  10.13443/j.cjors.2017042801ZOU Gaoxiang, TONG Chuangming, GAO Fei, et al. Composite electromagnetic scattering from tank target above complicated ground rough surface[J]. Chinese Journal of Radio Science, 2017, 32(3): 261–272. doi:  10.13443/j.cjors.2017042801
    [26] 潘小敏, 盛新庆, 孔雷. 电特大舰船电磁环境的并行多层快速多极子研究[J]. 装备环境工程, 2008, 5(1): 48–52, 60. doi:  10.3969/j.issn.1672-9242.2008.01.012PAN Xiaomin, SHENG Xinqing, and KONG Lei. Analysis on electromagnetic environment of an extremely large ship by parallel MLFMA[J]. Equipment Environmental Engineering, 2008, 5(1): 48–52, 60. doi:  10.3969/j.issn.1672-9242.2008.01.012
    [27] 刁桂杰, 许小剑, 倪虹, 等. 舰船目标宽带单脉冲雷达三维成像建模与仿真[J]. 系统仿真学报, 2018, 30(7): 2515–2524. doi:  10.16182/j.issn1004731/x.joss.201807011DIAO Guijie, XU Xiaojian, NI Hong, et al. Three-dimensional imaging modeling and simulation of ships based on wide-band monopulse radar[J]. Journal of System Simulation, 2018, 30(7): 2515–2524. doi:  10.16182/j.issn1004731/x.joss.201807011
    [28] YE H X and JIN Y Q. Fast iterative approach to electromagnetic scattering from the target above a rough surface[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(1): 108–115.
    [29] DING D Z, LI G M, AN Y Y, et al. Application of hierarchical two-level spectral preconditioning method for electromagnetic scattering from the rough surface[J]. International Journal of Antennas and Propagation, 2014, 2014: 1–10. doi:  10.1155/2014/752418
    [30] AN Yuyuan, WANG Daoxiang, and CHEN Rushan. A fast numerical algorithm for calculating electromagnetic scattering from an object above a rough surface[J]. Electromagnetics, 2013, 33(1): 10–22. doi:  10.1080/02726343.2013.751006
  • [1] 张新勋, 周生华, 刘宏伟.  目标极化散射特性对极化分集雷达检测性能的影响 . 雷达学报, 2019, 8(4): 510-518. doi: 10.12000/JR18112
    [2] 郭忠义, 汪彦哲, 郑群, 尹超逸, 杨阳, 宫玉彬.  涡旋电磁波天线技术研究进展 . 雷达学报, 2019, 8(5): 631-655. doi: 10.12000/JR19091
    [3] 闫华, 陈勇, 李胜, 胡利平, 李焕敏, 殷红成.  基于弹跳射线法的海面舰船目标三维散射中心快速建模方法 . 雷达学报, 2019, 8(1): 107-117. doi: 10.12000/JR18078
    [4] 陈珲, 徐亮, 张言明, 周小阳, 崔铁军.  超电大复杂目标太赫兹散射特性建模微波方法延拓研究 . 雷达学报, 2018, 7(1): 108-118. doi: 10.12000/JR17097
    [5] 高敬坤, 邓彬, 秦玉亮, 王宏强, 黎湘.  THz全尺寸凸体粗糙目标雷达回波散射建模与成像仿真 . 雷达学报, 2018, 7(1): 97-107. doi: 10.12000/JR17086
    [6] 赵华, 郭立新.  分形粗糙表面涂覆目标太赫兹散射特性 . 雷达学报, 2018, 7(1): 91-96. doi: 10.12000/JR17091
    [7] 牟媛, 吴振森, 赵豪, 武光玲.  粗糙金属和介质目标的太赫兹散射特性分析 . 雷达学报, 2018, 7(1): 83-90. doi: 10.12000/JR17094
    [8] 陈刚, 党红杏, 谭小敏, 陈珲, 崔铁军.  太赫兹低频段随机粗糙金属板散射特性研究 . 雷达学报, 2018, 7(1): 75-82. doi: 10.12000/JR17093
    [9] 文贡坚, 朱国强, 殷红成, 邢孟道, 杨虎, 马聪慧, 闫华, 丁柏圆, 钟金荣.  基于三维电磁散射参数化模型的SAR目标识别方法 . 雷达学报, 2017, 6(2): 115-135. doi: 10.12000/JR17034
    [10] 白杨, 吴洋, 殷红成, 阙肖峰.  无人机极化散射特性室内测量研究 . 雷达学报, 2016, 5(6): 647-657. doi: 10.12000/JR16032
    [11] 艾小锋, 曾勇虎, 高磊, 王肖洋, 汪连栋.  飞机目标全极化双基地散射特性研究 . 雷达学报, 2016, 5(6): 639-646. doi: 10.12000/JR16070
    [12] 韩壮志, 宋春吉, 侯建强.  伪码族复合连续波信号的多分辨率特性分析 . 雷达学报, 2016, 5(3): 278-283. doi: 10.12000/JR15100
    [13] 路东伟, 阙肖峰, 綦鑫, 聂在平.  米波频段飞机目标全极化散射仿真与特性分析 . 雷达学报, 2016, 5(2): 182-189. doi: 10.12000/JR16030
    [14] 王佳宁, 许小剑.  二维线性与非线性海面的宽带散射特性仿真及分析 . 雷达学报, 2015, 4(3): 343-350. doi: 10.12000/JR15053
    [15] 邢相薇, 计科峰, 康利鸿, 詹明.  HRWS SAR图像舰船目标监视技术研究综述 . 雷达学报, 2015, 4(1): 107-121. doi: 10.12000/JR14144
    [16] 余凡, 袁捷.  一种改进的随机粗糙地表的双尺度微波散射模型研究(英文) . 雷达学报, 2015, 4(5): 560-570. doi: 10.12000/JR15067
    [17] 魏仪文, 郭立新, 殷红成.  海洋内波对海面电磁散射特性的影响分析 . 雷达学报, 2015, 4(3): 326-333. doi: 10.12000/JR15060
    [18] 孙红梅, 陈广东, 张弓.  基于测控电磁波极化特征信息的飞行器姿态估计 . 雷达学报, 2013, 2(4): 466-475. doi: 10.3724/SP.J.1300.2013.13087
    [19] 张问一, 胡东辉, 丁赤飚.  基于FABEMD 和Goldstein 滤波器的SAR 舰船尾迹图像增强方法 . 雷达学报, 2012, 1(4): 426-435. doi: 10.3724/SP.J.1300.2012.20059
    [20] 李宁, 汪玲, 张弓.  多基ISAR 舰船侧视及俯视高分辨率成像方法 . 雷达学报, 2012, 1(2): 163-170. doi: 10.3724/SP.J.1300.2012.20021
  • 加载中
图(11) / 表 (3)
计量
  • 文章访问数:  1287
  • HTML全文浏览量:  600
  • PDF下载量:  138
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-02-28
  • 修回日期:  2019-04-04
  • 网络出版日期:  2019-05-09
  • 刊出日期:  2019-06-01

三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究

doi: 10.12000/JR19032
    基金项目:  国家自然科学基金(61701232, 61431006),江苏省自然科学基金(BKs20170854),中国科协青年人才托举工程(2017QNRC001),中国博士后基金(2017M620861, 2018T110127)
    作者简介:

    何 姿(1988–),女,博士,讲师。2016年6月在南京理工大学获得博士学位,2015年9月到2016年9月在美国伊利诺伊大学香槟分校(UIUC)交流访问。主要研究方向为计算电磁学、特征提取、目标识别,近年来在学术期刊、会议上发表论文60余篇,其中被SCI收录期刊论文30余篇。博士论文《目标电磁散射特性的抛物线方程方法研究》荣获2016年度中国电子学会优秀博士学位论文优秀奖;2017年度成功入选中国科协信息科技学会联合体青年人才托举工程;2018年度荣获International Applied Computational Electromagnetics Society (ACES) Symposium国际会议组颁发的“Young Scientist Award Scheme”

    陈如山(1965–),男,博士、教授、博士生导师,江苏省政协委员,教育部“长江学者”特聘教授,国家杰出青年科学基金获得者,享受国务院政府特殊津贴,国务院学位委员会学科评议组“电子科学与技术组”成员,“新世纪百千万人才工程”国家级人选,“十二五”国防基础研究专家组成员,教育部“长江学者和创新团队发展计划”创新团队带头人,国家重点学科“电磁场与微波技术”学科带头人,江苏省“青蓝工程”科技创新团队带头人,国际电子学报副主编,现为南京理工大学教授,博士生导师,电子工程与光电技术学院副院长。在电磁脉冲与瞬态电磁场、电磁兼容与雷达隐身技术、计算电磁学等方面发表论文160多篇,科学引文索引(SCI)收录论文140篇,(EI)收录论文120篇。曾获2000年度工程类科学引文索引(SCI)收录论文数为全国第一,南京市自然科学优秀学术论文三等奖,江苏省科技进步三等奖,国家教委科技进步三等奖,国家教委科技进步二等奖,江苏省国防工办科技进步一等奖,兵器工业集团公司科技进步二等奖,军队科技进步奖二等奖2项、三等奖1项。2015年江苏省研究生培养模式改革成果一等奖。目前已获得国防发明专利14项,实用新型发明专利8项

    通讯作者: 陈如山 eerschen@njust.edu.cn
  • 中图分类号: TN011

摘要: 3维随机粗糙海面与其上方复杂目标复合电磁(EM)散射特性的建模与分析在微波遥感、目标识别、雷达成像、导弹制导等领域中有着重要的研究价值。该文主要研究了基于高频算法的随机粗糙海面及舰船的复合电磁散射特性,开发了PO-IPO混合方法,为3维随机粗糙海面与复杂目标一体化高效求解提供了新思路。文中分别使用了物理光学方法(PO)、迭代物理光学方法(IPO)、PO-PO以及PO-IPO混合方法对海面及舰船进行了建模与仿真,其中,引入锥形波来代替平面波作为发射源,锥形波可以更好地抑制粗糙面在边缘位置被突然截断而形成的电磁反射和边缘绕射等效应。从数值仿真结果中可以看出,PO-IPO混合方法可将复杂物体本身面元间以及粗糙海面与物体间的耦合作用考虑在内,因此PO-IPO可以作为一种有效的途径来快速获取随机粗糙海面及舰船的复合电磁散射特性。

English Abstract

何姿, 陈如山. 三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究[J]. 雷达学报, 2019, 8(3): 318–325. doi:  10.12000/JR19032
引用本文: 何姿, 陈如山. 三维随机粗糙海面与舰船的复合电磁特性的高频方法分析研究[J]. 雷达学报, 2019, 8(3): 318–325. doi:  10.12000/JR19032
HE Zi and CHEN Rushan. Study on efficient high frequency method of electromagnetic scattering from ships above sea surface[J]. Journal of Radars, 2019, 8(3): 318–325. doi:  10.12000/JR19032
Citation: HE Zi and CHEN Rushan. Study on efficient high frequency method of electromagnetic scattering from ships above sea surface[J]. Journal of Radars, 2019, 8(3): 318–325. doi:  10.12000/JR19032
    • 实现粗糙面上方目标的电磁散射特性的快速分析在低空飞行器识别与监测、微波遥感以及雷达制导等领域有重要的研究意义。为了简化模型的计算,通常将背景(例如海洋、草地及沙漠等)和散射目标(例如飞机、舰船及导弹等)分为两个独立的部分进行计算,因此散射目标和背景之间的影响不能够被考虑在内。

      对于2维随机粗糙面模型的建立,除了可以使用入射波波长来衡量随机粗糙面的粗糙程度,还可以使用其它统计量来描述,例如均方根高度、相关函数和表面相关长度、功率谱密度、均方根斜率等等。由于2维随机粗糙面的电尺寸大,很难使用全波分析方法进行仿真与计算,因此,很多学者都是在一定的假定条件下得出近似解。基尔霍夫近似方法(KA)[1,2]是计算粗糙表面电磁散射问题的最为常见的一种近似方法,它认为粗糙面上任意一点的场强都可以由这一点上的切平面反射波来确定,因此,粗糙表面上任意一点的总场可近似为入射场与这点处相切的无穷大平面的反射场之和。微扰法(SPM)[3,4]中可将散射场近似等效为由远离边界而传播的不同振幅的平面波相加而成,这种方法只适用于小尺度的表面高度波动情况比入射波波长小和均方根斜率远远小于1的情形。小斜率近似方法(SSA)[5,6]相对于上面两种方法,使用范围更加广泛,并且对于大角度入射的时候,计算精度也相对比较高。除此之外,常见的近似方法还有双尺度法(Two Scale Model, TSM)[7]、消光定理(Extinction Theorem, ET)[8]、全波法(Full Wave Algorithm, FWA)[9]以及相位微扰法(Phase Perturbation Technique, PPT)[10]等。以上几种方法虽然计算速度快,占用内存小,但是没有考虑目标各个单元之间的互耦以及目标的边缘绕射和电磁波的多径传播。因此,开发有效的数值仿真方法分析这类问题是非常必要和急迫的。全波数值仿真方法主要有积分方程方法[11]和微分方程方法[12,13],它们可以精确地求解麦克斯韦方程,并求得空间的电磁场分布,但是,随着目标的电尺寸的增加,全波数值仿真方法所需要的计算资源也会剧烈地增长,因此不能够有效地分析大型的随机粗糙面。今年,很多学者开发出了一些快速算法来加速随机粗糙海平面与其上方物体电磁散射的快速建模分析[1420],例如广义前后向迭代法(GFBM)[14],广义前后向算法与谱加速方法(GFBM/SAA)[15],二次弹射法[16],基于KA分析的快速多级子方法(KA-MLFMA)[17]等方法,进一步降低了计算资源的消耗。

      同期,国内的许多学者分别采用多种数值方法,系统地开展了海面等粗糙面与目标复合的电磁散射特性研究。金亚秋等人采用双向解析射线追踪(BART)算法、FBM和共轭梯度法(CG)相结合的快速互耦迭代混合算法(FBM-CG)、FBM和扩展传播内层展开(EPILE)混合算法、KA和MoM相结合的混合算法(KA-MoM)等技术计算了舰船目标与粗糙海面的复合散射[21]。郭立新等人研究了物理光学近似,以及时频域数值混合算法分析了目标与粗糙面复合散射问题[22]。康士峰等人应用矩量法研究了粗糙面与其上方目标的散射回波特征[23]。朱国强等人基于物理光学法的双次弹射法、多层UV分解法以及利用谱域积分和高频近似混合算法分析了粗糙面上目标的电磁波复合散射[24]。童创明等人提出一种快速精确求解复杂陆地粗糙面及其上方坦克目标复合散射的混合方法[25]。盛新庆等人提出基尔霍夫近似结合多层快速多极子技术研究1维介质粗糙面与上方目标的后向散射[26]。许小剑等人采用混合高频方法研究了海面与目标形成的角反射器结构,建立了时变海面舰船目标动态雷达特征信号的数值仿真模型[27]。张民等人运用时域表面积分方程方法研究了锥形波入射下粗糙海面和3维介质目标的瞬时散射特征,并基于结合小面元模型的多路径复合散射理论、IPO-MEC等算法广泛开展了海上目标电磁环境建模及SAR成像仿真工作[28]。陈如山等人基于矩量法,采用特征基函数(CBF)结合矩阵分解算法(MDA),以及谱预条件技术实现粗糙面目标电磁散射特性快速分析[29,30]。市场上主流的电磁仿真商用软件有FEKO, HFSS, CST以及FDTD Solutions等,但是这些软件的核心代码被国外垄断,价格偏高,可拓展性差,因此,仍需要开发更为有效的数值计算方法用于分析海上目标电磁散射问题。

      本文从高频近似方法出发,分别介绍了物理光学方法(PO)、迭代物理光学方法(IPO)以及PO-IPO混合方法,实现了3维随机粗糙海平面及其上方物体一体化建模与电磁散射特性的分析。可以看出,PO方法的计算速度最快,IPO方法相对于PO方法,当分析结构较为复杂的物体时,计算精度得到显著的提高,但是计算效率会有所降低。因此,本文采用PO-IPO的混合方法,使用IPO分析粗糙海面上方的物体,考虑到粗糙面上方物体的互耦作用,使用PO分析海面以及海面与物体间的作用,从而能够有效地对3维随机粗糙海面与舰船的复合目标的电磁散射特性进行分析。

    • 当物体接受到电磁波的照射时,物体表面就会感应出电磁流,从而可以获得远区的散射场。物理光学方法不考虑单元之间的互耦,因此不需要求解矩阵方程,因此计算效率高,在工程中有很强的实用性。由Stratton-Chu公式,可得到散射场公式如下

      $$ \begin{align} {{\text{E}}_{\rm s}}\left( {\text{r}} \right) =& \int_{\rm s}\{ - {\rm{j}}\omega {\mu _0}[ {\text{J}}\left( {{\text{r}}'} \right)g\left( {{\text{r}},{\text{r}}'} \right) \\ {\rm{}}& \left.+ \frac{1}{{{k^2}}}\nabla ' \cdot {\text{J}}\left( {{\text{r}}'} \right)\nabla g\left( {{\text{r}},{\text{r}}'} \right) \right] \\ {\rm{}}& + {\text{M}}\left( {{\text{r}}'} \right) \times g\left( {{\text{r}},{\text{r}}'} \right) \} {\rm{d}}{\rm s}' \end{align} $$ (1)

      其中,${\text{r}}$表示场点,${\text{r}}'$表示源点。$g\left( {{\text{r}},{\text{r}}'} \right) = {{{{\rm e}^{ - {\rm j}{\rm k}R}}}}\bigr/ $$\left({{4{\text{π}} R}}\right)$, $R = \left| {{\text{r}} - {\text{r}}'} \right|$

      根据远场近似条件,远区的散射场公式可简化为

      $$ \begin{align} {{\text{E}}_{\rm s}}\left( {\text{r}} \right) =& \frac{{{\rm jk}}}{{4{\text{π}} R}}{{\rm e}^{ - {\rm jk}R}}\int_{\rm s} \left[ {{{\hat {\rm k}}}_{\rm s}} \times {\text{M}}\left( {{\text{r}}'} \right)\right. \\ {\rm{}}& \left.+ \eta {{{\hat {\rm k}}}_{\rm s}} \times \left( {{{{\hat {\rm k}}}_{\rm s}} \times {\text{J}}\left( {{\text{r}}'} \right)} \right) \right]{{\rm e}^{{\rm jk} \cdot {{{\hat {\rm k}}}_{\rm s}}\cdot{\text{r}}'}} {\rm d}{\rm s}' \end{align} $$ (2)

      其中,$\rm k$代表入射波的波数,$\eta = 120{\text{π}} $自由空间中的本征阻抗。

      假设入射电场为${{\text{E}}_{\rm i}} = \hat a{E_0}{{\rm e}^{ - {\rm jk }\cdot {{{\hat {\rm k}}}_{\rm i}} \cdot {\text{r}}'}}$,建立局部正交坐标系$\left( {{\hat \tau} ,{\hat \delta} ,{{{\hat {\rm k}}}_{\rm i}}} \right)$,其中${\hat \tau} = {{\hat {\rm k}}_{\rm i}} \times {\hat n}\Bigr/\left| {{{{\hat {\rm k}}}_{\rm i}} \times {\hat n}} \right|$, ${\hat \delta} = $${{\hat {\rm k}}_{\rm i}} \times {\hat \tau} $。基于基尔霍夫近似,粗糙面的感应电磁流可表示为如下形式

      $$ \begin{align} {\text{J}}\left( {{\text{r}}'} \right) = \hat n \times {\text{H}} = & - \frac{1}{\eta }\left[ \left( {1 + {R_{\rm VV}}} \right)\left( {\hat a \cdot \hat \delta } \right)\left( {\hat n \times \hat \tau } \right) \right.\\ {\rm{}}& \left.+ \left( {1 - {R_{\rm HH}}} \right)\left( {{{{\hat {\rm k}}}_{\rm i}} \cdot \hat n} \right)\left( {\hat a \cdot \hat \tau } \right)\hat \tau \right]\\ {\rm{}} & \cdot{E_0}{{\rm e}^{ -{\rm jk} \cdot {{{\hat {\rm k}}}_{\rm i}} \cdot {\text{r}}'}} \end{align} $$ (3)
      $$ \begin{align} {\text{M}}\left( {{\text{r}}'} \right) = {\text{E}} \times \hat n = & - \Bigr[ \left( {1 + {R_{\rm HH}}} \right)\left( {\hat a \cdot \hat \tau } \right)\left( {\hat n \times \hat \tau } \right) \\ {\rm{}}& \left.- \left( {1 - {R_{\rm VV}}} \right)\left( {{{{\hat {\rm k}}}_{\rm i}} \cdot \hat n} \right)\left( {\hat a \cdot \hat \delta } \right)\hat \tau \right] \hspace{10pt}\\ {\rm{}}& \cdot{E_0}{{\rm e}^{ - {\rm{jk}} \cdot {{{\hat {\rm k}}}_{\rm i}} \cdot {\text{r}}'}} \end{align} $$ (4)

      其中,${R_{\rm VV}}$${R_{\rm HH}}$分别表示VV极化与HH极化的Fresnel反射系数。当物体为理想导体时,${R_{\rm VV}} = $$ - {R_{\rm HH}} = 1$,因此有${\text{J}} = {\hat n} \times {\text{H}} = 2{\hat n} \times {{\text{H}}\!_{\rm i}}$, ${\text{M}} = 0$

    • 物理光学方法不考虑各个面元之间的耦合作用,可以将各个面元视为单独的辐射单元。因此,物理光学方法可以较为有效地分析电大简单的模型,但是,对于含有强耦合结构的散射体,物理光学方法就不能够满足计算精度的要求。1995年,Feamndo Obelleiro-Basteiro, Jose Luis Rodriguez和Robert J.Burkholder等人首次提出了利用迭代物理光学法分析含腔物体的电磁散射特性,迭代物理光学法通过不断更新表面电磁流,从而更好地模拟了各个面元间的互耦作用,相对于物理光学方法获得了较高的计算精度。

      对于理想导体目标,可得某一个面元上的表面电流如下

      $$ \begin{align} {\text{J}} =& 2{\hat n} \times {\text{H}}_{\rm i}^{\rm total} = 2{\hat n} \times \left( {{\text{H}}\!_{\rm i} + {{\text{H}}^{\rm s}}} \right) \\ =& 2{\hat n} \times {\text{H}}\!_{\rm i} + 2{\hat n} \times P.V.\int_{\rm s} \!\!{\left({\text{J}} \times \nabla g\left( {{\text{r}},{\text{r}}'} \right)\right){\rm d}{\rm s}'} \quad\ \ \end{align} $$ (5)

      其中,$P.V.\displaystyle\int_{} {} $代表去除了自身面元在内的主值积分,通过迭代的方式来不断更新表面电流,式(5)可改写为

      $$ \begin{align} {{\text{J}}_{k + 1}} = &2{\hat n} \times {\text{H}}_{\rm i}^{\rm total} = 2{\hat n} \times {\text{H}}\!_{\rm i}+ 2{\hat n} \\ {\rm{}}&\times P.V.\int_{\rm s} \left({{\text{J}}_k} \times \nabla g\left( {{\text{r}},{\text{r}}'} \right)\right){\rm d}{\rm s}' \end{align} $$ (6)

      当式(6)中的电流达到稳定后,再由式(2)计算可得远区散射场。总之,迭代物理光学方法(IPO)可以将各个面元间的相互作用考虑在内,而弹跳射线方法(SBR)仅可以将射线路径上的面元的相互作用考虑在内,因此,对于结构相对较为复杂的目标,迭代物理光学方法的计算精度会优于弹跳射线方法,但是,所需的计算资源以及计算时间也会有一定程度的增加。因此,本文考虑使用迭代物理光学方法分析海面上的复杂目标,并结合物理光学方法来分析目标与海面的相互作用。

    • 对于粗糙面与目标的复合模型,当物体模型相对复杂的时候,则需要考虑粗糙面与目标之间的相互作用。如图1所示,${{\hat k}_{\rm i}}$表示入射波的方向,${\text{J}}_1^{\rm s}$,${\text{M}}_1^{\rm s}$表示粗糙面上由入射波的照射而引发的1阶电流和磁流,${\text{J}}_1^{\rm t}$表示目标物体上入射波的照射而引发的1阶电流。下标1为电流的阶数,上标s表示粗糙面,t表示目标。图1(a)表示了入射波直接作用于粗糙面和目标所产生的感应电磁流,不考虑粗糙面与目标之间的耦合作用。对于金属目标,1阶感应电流可由${\text{J}} = 2{\hat n} \times {{\text{H}}_{\rm i}}$计算得到。图1(b)代表粗糙面与目标之间的耦合作用,首先,入射波照射到目标表面产生1阶电流${\text{J}}_1^{\rm t}$,目标表面的1阶电流作为新的源产生了1阶散射电场${\text{E}}_1^{\rm t}$,从而作为新的入射场作用到粗糙面上,在粗糙面上引发了2阶感应电磁流${\text{J}}_2^{\rm s}$,${\text{M}}_2^{\rm s}$,通过式(3)和式(4)可以计算出粗糙面的2阶感应电磁流。同理,入射波照射到粗糙面上产生1阶电磁流${\text{J}}_1^{\rm s}$,${\text{M}}_1^{\rm s}$,然后产生1阶散射磁场${\text{H}}_1^{\rm s}$,从而可以在物体表面生成2阶感应电流${\text{J}}_2^{\rm t}$。最终,粗糙面上的2阶感应电磁流和物体表面上的2阶感应电流将会在空间产生2阶散射电磁场${\text{E}}_2^{\rm t}$${\text{H}}_2^{\rm s}$。依次类推,可以得到粗糙面与目标上的各阶感应电磁流,最终可累加得到总感应电磁流

      图  1  多路径模型基本原理

      Figure 1.  Theory of multipath model

      $$ {\text{J}}_{\rm s {\tiny{-}} total}^n = {\text{J}}_1^{\rm s} + {\text{J}}_2^{\rm s} +·\!·\!· + {\text{J}}_n^{\rm s} \hspace{15pt} $$ (7)
      $$ {\text{M}}_{\rm s {\tiny{-}} total}^n = {\text{M}}_1^{\rm s} + {\text{M}}_2^{\rm s} + ·\!·\!· + {\text{M}}_n^{\rm s} $$ (8)
      $$ {\text{J}}_{\rm t {\tiny{-}} total}^n = {\text{J}}_1^{\rm t} + {\text{J}}_2^{\rm t} + ·\!·\!· + {\text{J}}_n^{\rm t} \hspace{20pt} $$ (9)

      对于分析复杂目标与粗糙面复合模型时,由于物体结构复杂,物体各个面元间的相互耦合作用不可忽略,PO不考虑面元间的相互作用,SBR仅可以将射线路径上面元的相互作用考虑在内,因此使用PO, SBR等方法分析结构较为复杂的目标时会有较大的误差。对于PO-IPO混合方法,粗糙海面上的各个阶的电磁流由PO算法计算得到,假设复杂目标为电大金属目标,则其各个阶的电流由IPO算法得到,粗糙海面与复杂目标之间的耦合作用则使用图1中的处理方式。虽然PO-IPO混合方法的精度相比于PO, SBR的精度要高,但是由于IPO方法需要考虑目标各个面元间的耦合作用,因此PO-IPO混合方法需要的计算资源相对于PO, SBR要多。如图2给出了PO-IPO算法的流程框图。

      图  2  PO-IPO算法流程图

      Figure 2.  Flow chart for PO-IPO

    • 假设一个高斯谱生成的3维粗糙面,大小为${L_x} = {L_y} = 0.6\;{\rm{m}}$,相对介电常数为$ {\varepsilon _{\rm r}} =(5.4,{\rm j}0.04)$,粗糙面轮廓如图3所示。粗糙面的相关长度为${l_x} = {l_y} = 0.06\;{\rm{m}}$,均方根高度为$h = 0.1\;{\rm{m}}$。锥形波入射,频率为3 GHz。入射角度分别为$ {\theta _{\rm i}}{\rm{ = }}{30^ {\circ} },$${60^ {\circ} }$, ${\varphi _{\rm i}}{\rm{ = }}{0^ {\circ} }$,散射角为$ - {90^ {\circ} } \le {\theta _{\rm s}} \le {90^ {\circ} }$, ${\varphi _{\rm s}}{\rm{ = }}{0^ {\circ} }$。分别使用了物理光学方法(PO)、弹跳射线方法(SBR)以及全波数值仿真方法(JMCFIE)对粗糙海面进行了数值仿真,图4是PO, SBR和JMCFIE的入射角度为${\theta _{\rm i}}$=30°, VV极化双站RCS的结果对比,图5是PO, SBR和JMCFIE的入射角度为${\theta _{\rm i}}$=60°, VV极化双站RCS的结果对比,表1是两种方法消耗资源对比。可以看出,PO能够在保证一定精度的同时,大大节省了计算资源,当浪情较低时,SBR和PO的计算精度基本吻合。因此,PO可作为一种有效的方法来分析粗糙海平面的电磁散射问题。

      图  3  PM谱粗糙面模型

      Figure 3.  Model of PM rough surface

      图  4  入射角度为${\theta _{\rm i}}$=30°粗糙海面模型的双站RCS结果对比

      Figure 4.  Bistatic RCS results for the PM rough surface when ${\theta _{\rm i}}$=30°

      图  5  入射角度为${\theta _{\rm i}}$=60°粗糙海面模型的双站RCS结果对比

      Figure 5.  Bistatic RCS results for the PM rough surface when ${\theta _{\rm i}}$=60°

      表 1  PO与JMCFIE计算资源对比

      Table 1.  Comparison of computational resources between PO and JMCFIE

      数值方法时间(s)内存(MB)
      SBR27176
      PO23174
      JMCFIE19565912
    • 考虑一个全尺寸舰船模型的电磁散射特性,如图6所示,舰船长为80 m,宽为12 m,高为14 m。入射波为均匀平面波,频率为500 MHz,入射角为${\theta _{\rm i}} = {30^ {\circ} }$, ${\varphi _{\rm i}} = {180^ {\circ} }$,极化方式为垂直极化,观察角度为$ - {90^ {\circ} } \le {\theta _{\rm s}} \le {90^ {\circ} }$, ${\varphi _{\rm s}}{\rm{ = }}{0^ {\circ} }$。分别使用了弹跳射线法(SBR)、迭代物理光学方法(IPO)以及全波数值仿真方法(MoM-MLFMA)对舰船进行了数值仿真,如图7所示,给出了3种方法的双站RCS的对比图,表2对3种方法进行了计算资源的对比,可以看出对于这种较为复杂的模型,IPO算法的计算精度比SBR方法要高,但是所需的计算资源要比SBR方法多。因此,分析复杂目标的电磁散射特性时,IPO方法是一个很好的选择。

      图  6  舰船模型示意图

      Figure 6.  Model of ship

      表 2  SBR, IPO与MoM-MLFMA计算资源对比

      Table 2.  Comparison of computational resources among SBR, IPO and MoM-MLFMA

      数值方法时间(s)内存(MB)
      SBR1558
      IPO259374
      MoM-MLFMA1125701178

      图  7  舰船模型的双站RCS结果对比

      Figure 7.  Bistatic RCS results for the ship model

    • 首先,仿真了一个舰船与粗糙海面的缩比模型,分别使用了PO-IPO混合方法以及全波数值仿真方法(JMCFIE)对舰船进行了数值仿真,假设一个高斯谱生成的大小为$51\;{\rm{m}} \times 9\;{\rm{m}}$的海面,均方根高度为$0.6\;{\rm{m}}$,相关长度为$1.2\;{\rm{m}}$,相对介电常数为${\varepsilon _{\rm r}} = (5.4,{\rm j}0.04)$,上方目标为一长26 m,宽4 m,高4.6 m的理想导体舰船,组合模型示意图如图8所示,入射角为${\theta _{\rm i}}{\rm{ = }}{30^ {\circ} }$, ${\varphi _{\rm i}}{\rm{ = 18}}{0^ {\circ} }$,入射波频率为0.5 GHz,散射角为$ - {90^ {\circ} } \le {\theta _{\rm s}} \le {90^ {\circ} }$, ${\varphi _{\rm s}}{\rm{ = }}{0^ {\circ} }$,极化方式为VV极化。如图9所示,给出了两种方法的双站RCS的对比图,表3对两种方法进行了计算资源的对比,可以看出对于这种较为复杂的组合模型,PO-IPO算法的计算精度仍然比较高。

      图  8  粗糙海面与缩比舰船组合模型示意图

      Figure 8.  Model of ship over rough sea surface

      表 3  PO-IPO与JMCFIE计算资源对比

      Table 3.  Comparison of computational resources between PO-IPO and JMCFIE

      数值方法时间(s)内存(MB)
      PO-IPO471.842.5
      JMCFIE8994.021645.0

      图  9  粗糙海面与缩比舰船组合模型的双站RCS结果对比

      Figure 9.  Bistatic RCS results for the ship model over rough sea surface

      最后,分析一个全尺寸的带有精细结构的舰船与海面的复合模型,精确的数值仿真方法已经不能够有效地分析这类问题。假设粗糙海面的大小为250 m×40 m,海况设为3级,相对介电常数为${\varepsilon _{\rm r}} = (32.3,{\rm j}30.2)$,上方有一舰船,舰船长165 m,宽21 m,高59 m,模型如图10所示。入射波为锥形波,频率为30 GHz,入射角为$ - {75^ {\circ} } \le {\theta _{\rm i}} \le {75^ {\circ} }$, ${\varphi _{\rm i}}{\rm{ = }}{0^ {\circ} }$,垂直极化。由3.1小节和3.2小节的仿真数值结果分析可得,对于舰船模型可使用IPO方法对其进行仿真计算,对于粗糙海面模型可以使用PO对其进行建模计算。如图11所示,使用了PO-IPO混合方法对舰船与粗糙海平面模型进行了单站雷达散射截面积的数值仿真,模型的总未知量个数为320亿,总共使用了100个核进行并行计算,需要内存约为11 GB,每个角度上的计算时间为77 s。

      图  10  粗糙海面与舰船组合模型示意图

      Figure 10.  Model of full-size ship over rough sea surface

      图  11  粗糙海面与舰船组合模型的双站RCS结果对比

      Figure 11.  Bistatic RCS results for the full-size ship model over rough sea surface

    • 对于粗糙海面与简单目标的复合模型的电磁散射问题,可以使用PO, SBR等高频方法对其实现高效分析,但是随着目标结构变得复杂,传统高频方法的精度就会受到限制,因此,引入迭代物理光学方法对复杂目标进行计算,从而提出了PO-IPO混合算法。使用PO方法计算粗糙海面的表面电磁流,使用IPO方法计算复杂目标的表面电流,并通过惠更斯原理不断迭代更新两者的表面电磁流,直至达到稳定状态,便可得到远区的散射场。PO-IPO混合方法既考虑了复杂物体本身面元间的互作用,同时考虑了粗糙海面与物体间的耦合,因此可以保证计算效率的同时获得较高的计算精度。

参考文献 (30)

目录

    /

    返回文章
    返回