CN109215072A - 基于tanh-sinh积分法的箔条云RCS获取方法 - Google Patents
基于tanh-sinh积分法的箔条云RCS获取方法 Download PDFInfo
- Publication number
- CN109215072A CN109215072A CN201810818113.0A CN201810818113A CN109215072A CN 109215072 A CN109215072 A CN 109215072A CN 201810818113 A CN201810818113 A CN 201810818113A CN 109215072 A CN109215072 A CN 109215072A
- Authority
- CN
- China
- Prior art keywords
- chaff
- formula
- point
- cloud
- coordinate system
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了基于tanh‑sinh积分法的箔条云RCS获取方法。传统的箔条云RCS算法的计算效率低，当计算大规模箔条云时，其计算量将无法满足实际需求。本发明如下：一、获取箔条云及特征箔条的信息，并建立全局坐标系、局部坐标系。二、获取特征箔条的双站散射系数。三、计算被积函数σ⊥to⊥、σ⊥to//、σ//to//的奇异点。四、剖分σ⊥to⊥、σ⊥to//、σ//to//的待积分区域。五、分别积分各小区域的雷达散射截面积。六、通过叠加的方式计算特征箔条的雷达散射截面积。七、通过乘上箔条数量计算箔条云的雷达散射截面积。本发明在计算精度接近蒙特卡罗法的同时，大大增强了计算效率。
Description
技术领域
本发明属于计算机辅助分析与设计以及软件设计技术领域，具体涉及一种基于tanh-sinh积分法的箔条云RCS获取方法。
背景技术
箔条云通常是由若干条细金属线所组成，当目标感知到危险时，它会喷射出若干条细金属线，形成一个云状结构。有一种典型的箔条形式是偶极子，其工作原理是一定长度的偶极子能够在雷达载波频率上产生共振，从而隐藏待打击目标，有效保护真实目标。因此箔条及箔条云雷达散射截面积(RCS)的计算一直以来都是科学界的热点问题。目前为止，传统的算法是利用蒙特卡罗法进行计算，该方法的优点是计算正确性高，缺点是计算效率低，当计算大规模箔条云时，其计算量将无法满足实际需求。因此，需要一种新的RCS获取方法来提升箔条云雷达散射截面积的获取速率。
发明内容
本发明的目的在于提供一种基于tanh-sinh积分法的箔条云RCS获取方法。
本发明的具体步骤如下：
步骤一、雷达对箔条云进行检测，得到箔条云内箔条的数量x，任取箔条云中的一条箔条作为特征箔条；并获取特征箔条的长度L、入射电场的垂直极化部分水平极化部分反射电场的垂直极化部分水平极化部分
建立全局坐标系和局部坐标系。全局坐标系的坐标原点为点T；点T为任意一个不在特征箔条上的点。局部坐标系的坐标原点为特征箔条的几何中心点D。全局坐标系的三根坐标轴分别为X轴、Y轴、Z轴。局部坐标系的三根坐标轴分别为X'轴、Y'轴、Z'轴。点D作为源点D，取空间中的一个与点T、点D不重合的点作为场点R。局部坐标系的X'-Y'平面与点T、点D、点R所成平面重合。局部坐标系的X'轴平分∠TDR。
步骤二、计算矩阵如式(1)所示，矩阵如式(2)所示、矩阵如式(3)所示：
式(1)、(2)、(3)中，为矩阵转换后的入射电场垂直极化部分、水平极化部分；为矩阵转换后的接收电场垂直极化部分、水平极化部分；矩阵[A]为箔条散射矩阵，表达式为矩阵[B]为入射端的极化转换矩阵；矩阵[C]为接收端的极化转换矩阵。d11的表达式如式(4a)所示，d12的表达式如式(4b)所示，d21的表达式如式(4c)所示，d22的表达式如式(4d)所示
d11＝BA0(β/2)cos2θd (4a)
式(4a)、(4b)、(4c)、(4d)中，θd为箔条中心轴线与局部坐标系内Z'轴的偏置角度；为箔条与局部坐标系内X'轴的偏置角度；β为线TD与线DR之间的夹角；B的表达式如式(5)所示；A0(β/2)的表达式如式(6)所示；
B＝[-jη/(2λzradr2)]exp(-j2πr2/λ) (5)
A0(β/2)＝A(π/2,β/2)A(π/2,-β/2) (6)
式(5)、(6)中，j为虚数单位；η为空气的特征阻抗；λ为对箔条云进行检测的雷达的工作频率；zrad为箔条的特性阻抗；r2为源点D与场点R之间的距离。
计算特征箔条雷达散射截面积的垂直同极化值σ⊥to⊥的表达式如式(7a)所示，平行同极化值σ//to//的表达式如式(7b)所示，交叉极化值σ⊥to//的表达式如式(7c)所示。
式(7a)、(7b)、(7c)中，为特征箔条的双站散射系数，其表达式如式(8)所示。
步骤三、根据式(1)、式(2)、式(3)及式(8)计算被积函数σ⊥to⊥、σ⊥to//、σ//to//的奇异点，分别是(π/2,β/2)、(π/2,π-β/2)、(π/2,π+β/2)、(π/2,2π-β/2)。
步骤四、将σ⊥to⊥、σ⊥to//、σ//to//的待积分区域分别剖分为十六个小区域，使得每个奇异点均成为其中四个小区域的顶点，且任意一个小区域上均有且仅有一个顶点为奇异点。
步骤五、对步骤四所得σ⊥to⊥、σ⊥to//、σ//to//的所有小区域分别进行积分，分别得到所有小区域的雷达散射截面积。
步骤六、将垂直同极化值σ⊥to⊥对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的垂直同极化值σ⊥to⊥。将平行同极化值σ//to//对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的平行同极化值σ//to//。将交叉极化值σ⊥to//对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的交叉极化值σ⊥to//。
步骤七、计算箔条云雷达散射截面积的垂直同极化值σ⊥to⊥总＝x·σ⊥to⊥；计算箔条云雷达散射截面积的交叉极化值σ⊥to//总＝x·σ⊥to//；计算箔条云雷达散射截面积的平行同极化值σ//to//总＝x·σ//to//。
进一步地，步骤四中，对σ⊥to⊥、σ⊥to//、σ//to//的所有小区域分别进行积分的方法采用利用tanh-sinh积分算法。
本发明具有的有益效果是：
1、本发明在计算大规模箔条云的雷达散射截面积时，通过提高单根箔条的双站散射系数的计算效率，可以在同等时间下，计算更大规模的箔条云。
2、本发明在计算同等规模箔条云，在保证同样精度的前提下，能够大幅减少计算时间，实现大规模箔条云雷达散射截面积的快速计算，继而满足实际应用场景。
附图说明
图1为本发明建立的特征箔条双站散射的全局坐标系；
图2为本发明建立的特征箔条双站散射的局部坐标系；
图3为本发明剖分待积区域的示意图；
图4为本发明与蒙特卡洛法对第一种箔条云进行雷达散射截面积计算的结果比较图；
图5为本发明与蒙特卡洛法对第二种箔条云进行雷达散射截面积计算的结果比较图；
图6为蒙特卡洛法计算不同规模箔条云雷达散射截面积所需时间随相对误差的变化曲线图；
图7为本发明计算不同规模箔条云雷达散射截面积所需时间随相对误差的变化曲线图；
图8为本发明与蒙特卡洛法计算箔条云雷达散射截面积所需时间随箔条云规模变化的对比曲线图；
图9为本发明与蒙特卡洛法计算箔条云雷达散射截面积所需时间随相对误差的对比曲线图。
具体实施方式
以下结合附图对本发明作进一步说明。
一种基于tanh-sinh积分法的箔条云RCS获取方法具体步骤如下：
步骤一、如图1和2所示，雷达对箔条云进行检测，得到箔条云内箔条的数量x，任取箔条云中的一条箔条作为特征箔条；并获取特征箔条的长度L、入射电场的垂直极化部分水平极化部分反射电场的垂直极化部分水平极化部分
建立特征箔条的双站散射的全局坐标系和局部坐标系。全局坐标系的坐标原点为点T；点T为任意一个不在特征箔条上的点。局部坐标系的坐标原点为特征箔条的几何中心点D。全局坐标系的三根坐标轴分别为X轴、Y轴、Z轴。局部坐标系的三根坐标轴分别为X'轴、Y'轴、Z'轴。点D作为源点D，取空间中的一个与点T、点D不重合的点作为场点R。局部坐标系的X'-Y'平面(X'轴、Y'轴所成平面)与点T、点D、点R所成平面重合。局部坐标系的X'轴平分∠TDR。∠TDR即线TD与线TR所成的角。
步骤二、计算矩阵如式(1)所示，矩阵如式(2)所示、矩阵如式(3)所示：
式(1)、(2)、(3)中，为矩阵转换后的入射电场垂直极化部分、水平极化部分；为矩阵转换后的接收电场垂直极化部分、水平极化部分；矩阵[A]为箔条散射矩阵，表达式为矩阵[B]为入射端的极化转换矩阵；矩阵[C]为接收端的极化转换矩阵。d11的表达式如式(4a)所示，d12的表达式如式(4b)所示，d21的表达式如式(4c)所示，d22的表达式如式(4d)所示
d11＝BA0(β/2)cos2θd (4a)
式(4a)、(4b)、(4c)、(4d)中，θd为箔条中心轴线与局部坐标系内Z'轴的偏置角度(夹角)；为箔条与局部坐标系内X'轴的偏置角度；β为线TD(点T与源点D之间的连线)与线DR(源点D与场点R之间的连线)之间的夹角；B的表达式如式(5)所示；A0(β/2)的表达式如式(6)所示；
B＝[-jη/(2λzradr2)]exp(-j2πr2/λ) (5)
A0(β/2)＝A(π/2,β/2)A(π/2,-β/2) (6)
式(5)、(6)中，j为虚数单位；η为空气的特征阻抗，其值为120πΩ；λ为对箔条云进行检测的雷达的工作频率；zrad为箔条的特性阻抗；r2为源点D与场点R之间的距离；exp(-j2πr2/λ)为自然常数e的(-j2πr2/λ)次方。
计算特征箔条雷达散射截面积的垂直同极化值σ⊥to⊥的表达式如式(7a)所示，平行同极化值σ//to//的表达式如式(7b)所示，交叉极化值σ⊥to//的表达式如式(7c)所示。
式(7a)、(7b)、(7c)中，为特征箔条的双站散射系数，其表达式如式(8)所示。
步骤三、根据式(1)、式(2)、式(3)及式(8)计算被积函数σ⊥to⊥、σ⊥to//、σ//to//的奇异点(奇异点获取为常规技术)，分别是(π/2,β/2)、(π/2,π-β/2)、(π/2,π+β/2)、(π/2,2π-β/2)。
步骤四、如图3所示，根据tanh-sinh积分算法，将σ⊥to⊥、σ⊥to//、σ//to//的待积分区域分别剖分为十六个小区域，使得每个奇异点均成为其中四个小区域的顶点，且任意一个小区域上均有且仅有一个顶点为奇异点。
步骤五、对步骤四所得σ⊥to⊥、σ⊥to//、σ//to//的所有小区域分别利用tanh-sinh积分算法进行积分，分别得到所有小区域的雷达散射截面积。
步骤六、将垂直同极化值σ⊥to⊥对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的垂直同极化值σ⊥to⊥。将平行同极化值σ//to//对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的平行同极化值σ//to//。将交叉极化值σ⊥to//对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的交叉极化值σ⊥to//。
步骤七、计算箔条云雷达散射截面积的垂直同极化值σ⊥to⊥总＝x·σ⊥to⊥；计算箔条云雷达散射截面积的交叉极化值σ⊥to//总＝x·σ⊥to//；计算箔条云雷达散射截面积的平行同极化值σ//to//总＝x·σ//to//。由于箔条云是离散的，即箔条云里的每一根箔条之间相互没有影响。因此，将单根箔条的雷达散射截面积乘以箔条云的数目即可得到整个箔条云的雷达散射截面积。
为了验证本发明的正确性和有效性，申请人分别使用蒙特卡洛法和本发明计算两种不同尺寸箔条所形成箔条云的雷达散射截面积。第一种箔条云的计算结果图如图4所示。第二种箔条云的计算结果图如图5所示。图4及5中，正方形点、圆点、正置三角形分别对应本发明计算出的雷达散射截面积垂直同极化值、交叉极化值、平行同极化值；倒置的三角形点、朝向左侧的三角形点朝向右侧的三角形点分别对应蒙特卡洛法明计算出的雷达散射截面积垂直同极化值、交叉极化值、平行同极化。结果表明，本发明所提算法得到的数据和蒙特卡洛法一致，其正确性得以保证。
为了显示本发明的高速性，申请人分别使用蒙特卡洛法和本发明计算不同规模箔条云的雷达散射截面积。用蒙特卡洛法计算规模为20根、2000根、20000根的箔条云雷达散射截面积所消耗的时间折线图如图6所示。
用本发明计算规模为20根、2000根、20000根的箔条云雷达散射截面积所消耗的时间折线图如图7所示。对比图6和图7可以发现，本发明计算不同规模箔条云的雷达散射截面积所需的时间均低于蒙特卡洛法。
本发明与蒙特卡洛法计算箔条云雷达散射截面积所需时间随箔条云规模变化的对比曲线图如图8所示，可以看出，本发明计算不同规模箔条云雷达散射截面积所需时间均低于蒙特卡洛法
本发明与蒙特卡洛法计算箔条云雷达散射截面积所需时间随相对误差的对比曲线图如图9所示，可以看出，本发明在不同相对误差要求下，计算雷达散射截面积所需时间均低于蒙特卡洛法。
综上所示，本发明的计算效率在同等精度条件下，远高于传统算法。
Claims (2)
1.基于tanh-sinh积分法的箔条云RCS获取方法，其特征在于：步骤一、雷达对箔条云进行检测，得到箔条云内箔条的数量x，任取箔条云中的一条箔条作为特征箔条；并获取特征箔条的长度L、入射电场的垂直极化部分水平极化部分反射电场的垂直极化部分水平极化部分
建立全局坐标系和局部坐标系；全局坐标系的坐标原点为点T；点T为任意一个不在特征箔条上的点；局部坐标系的坐标原点为特征箔条的几何中心点D；全局坐标系的三根坐标轴分别为X轴、Y轴、Z轴；局部坐标系的三根坐标轴分别为X'轴、Y'轴、Z'轴；点D作为源点D，取空间中的一个与点T、点D不重合的点作为场点R；局部坐标系的X'-Y'平面与点T、点D、点R所成平面重合；局部坐标系的X'轴平分∠TDR；
步骤二、计算矩阵如式(1)所示，矩阵如式(2)所示、矩阵如式(3)所示：
式(1)、(2)、(3)中，为矩阵转换后的入射电场垂直极化部分、水平极化部分；为矩阵转换后的接收电场垂直极化部分、水平极化部分；矩阵[A]为箔条散射矩阵，表达式为矩阵[B]为入射端的极化转换矩阵；矩阵[C]为接收端的极化转换矩阵；d11的表达式如式(4a)所示，d12的表达式如式(4b)所示，d21的表达式如式(4c)所示，d22的表达式如式(4d)所示
d11＝BA0(β/2)cos2θd (4a)
式(4a)、(4b)、(4c)、(4d)中，θd为箔条中心轴线与局部坐标系内Z'轴的偏置角度；为箔条与局部坐标系内X'轴的偏置角度；β为线TD与线DR之间的夹角；B的表达式如式(5)所示；A0(β/2)的表达式如式(6)所示；
B＝[-jη/(2λzradr2)]exp(-j2πr2/λ) (5)
A0(β/2)＝A(π/2,β/2)A(π/2,-β/2) (6)
式(5)、(6)中，j为虚数单位；η为空气的特征阻抗；λ为对箔条云进行检测的雷达的工作频率；zrad为箔条的特性阻抗；r2为源点D与场点R之间的距离；
计算特征箔条雷达散射截面积的垂直同极化值σ⊥to⊥的表达式如式(7a)所示，平行同极化值σ//to//的表达式如式(7b)所示，交叉极化值σ⊥to//的表达式如式(7c)所示；
式(7a)、(7b)、(7c)中，为特征箔条的双站散射系数，其表达式如式(8)所示；
步骤三、根据式(1)、式(2)、式(3)及式(8)计算被积函数σ⊥to⊥、σ⊥to//、σ//to//的奇异点，分别是(π/2,β/2)、(π/2,π-β/2)、(π/2,π+β/2)、(π/2,2π-β/2)；
步骤四、将σ⊥to⊥、σ⊥to//、σ//to//的待积分区域分别剖分为十六个小区域，使得每个奇异点均成为其中四个小区域的顶点，且任意一个小区域上均有且仅有一个顶点为奇异点；
步骤五、对步骤四所得σ⊥to⊥、σ⊥to//、σ//to//的所有小区域分别进行积分，分别得到所有小区域的雷达散射截面积；
步骤六、将垂直同极化值σ⊥to⊥对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的垂直同极化值σ⊥to⊥；将平行同极化值σ//to//对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的平行同极化值σ//to//；将交叉极化值σ⊥to//对应的十六个小区域的雷达散射截面积相累加，得到特征箔条的交叉极化值σ⊥to//；
步骤七、计算箔条云雷达散射截面积的垂直同极化值σ⊥to⊥总＝x·σ⊥to⊥；计算箔条云雷达散射截面积的交叉极化值σ⊥to//总＝x·σ⊥to//；计算箔条云雷达散射截面积的平行同极化值σ//to//总＝x·σ//to//。
2.根据权利要求1所述的基于tanh-sinh积分法的箔条云RCS获取方法，其特征在于：步骤四中，对σ⊥to⊥、σ⊥to//、σ//to//的所有小区域分别进行积分的方法采用利用tanh-sinh积分算法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810818113.0A CN109215072B (zh) | 2018-07-24 | 2018-07-24 | 一种箔条云rcs获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810818113.0A CN109215072B (zh) | 2018-07-24 | 2018-07-24 | 一种箔条云rcs获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109215072A true CN109215072A (zh) | 2019-01-15 |
| CN109215072B CN109215072B (zh) | 2021-07-20 |
Family
ID=64990130
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810818113.0A Active CN109215072B (zh) | 2018-07-24 | 2018-07-24 | 一种箔条云rcs获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109215072B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110502782A (zh) * | 2019-07-04 | 2019-11-26 | 西安电子科技大学 | 基于分区矢量输运的大体量箔条云电磁散射测定方法 |
| CN111199097A (zh) * | 2019-12-25 | 2020-05-26 | 西安电子科技大学 | 瞬态矢量辐射传输理论箔条云散射处理方法 |
| CN112364435A (zh) * | 2020-11-19 | 2021-02-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 多次抛撒箔条云气动融合轨道再入实时跟踪统计模拟方法 |
| CN113759322A (zh) * | 2021-09-08 | 2021-12-07 | 中国人民解放军海军航空大学 | 一种用于对抗雷达的箔条云rcs预测方法及系统 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4633255A (en) * | 1983-05-16 | 1986-12-30 | The United States Of America As Represented By The Secretary Of The Navy | Method for sea surface high frequency radar cross-section estimation using Dopler spectral properties |
| US6809682B1 (en) * | 2003-01-16 | 2004-10-26 | The United States Of America As Represented By The Secretary Of The Army | Method and device for the detection and track of targets in high clutter |
| CN102508220A (zh) * | 2011-10-24 | 2012-06-20 | 西瑞克斯(北京)通信设备有限公司 | 均匀双各向同性媒质物体的雷达散射截面获取方法 |
| CN102902841A (zh) * | 2012-08-28 | 2013-01-30 | 同济大学 | 一种计算机载复杂结构雷达横截面积的方法 |
| US20130038486A1 (en) * | 2010-04-05 | 2013-02-14 | Raytheon Company | Generating radar cross-section signatures |
| CN103530469A (zh) * | 2013-10-23 | 2014-01-22 | 武汉大学 | 一种角反射器雷达截面积的计算方法 |
| CN103675772A (zh) * | 2013-11-26 | 2014-03-26 | 北京宇航系统工程研究所 | 一种多功能sar复杂电磁环境模拟装置 |
| CN103870654A (zh) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | 基于并行矩量法与物理光学混合的电磁散射仿真方法 |
| CN106355641A (zh) * | 2016-08-26 | 2017-01-25 | 西安电子科技大学 | 三维粗糙目标雷达散射截面的仿真方法 |
| CN108614251A (zh) * | 2016-12-11 | 2018-10-02 | 南京理工大学 | 多涂覆目标的雷达散射截面可视化计算方法 |
-
2018
- 2018-07-24 CN CN201810818113.0A patent/CN109215072B/zh active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4633255A (en) * | 1983-05-16 | 1986-12-30 | The United States Of America As Represented By The Secretary Of The Navy | Method for sea surface high frequency radar cross-section estimation using Dopler spectral properties |
| US6809682B1 (en) * | 2003-01-16 | 2004-10-26 | The United States Of America As Represented By The Secretary Of The Army | Method and device for the detection and track of targets in high clutter |
| US20130038486A1 (en) * | 2010-04-05 | 2013-02-14 | Raytheon Company | Generating radar cross-section signatures |
| CN102508220A (zh) * | 2011-10-24 | 2012-06-20 | 西瑞克斯(北京)通信设备有限公司 | 均匀双各向同性媒质物体的雷达散射截面获取方法 |
| CN102902841A (zh) * | 2012-08-28 | 2013-01-30 | 同济大学 | 一种计算机载复杂结构雷达横截面积的方法 |
| CN103530469A (zh) * | 2013-10-23 | 2014-01-22 | 武汉大学 | 一种角反射器雷达截面积的计算方法 |
| CN103675772A (zh) * | 2013-11-26 | 2014-03-26 | 北京宇航系统工程研究所 | 一种多功能sar复杂电磁环境模拟装置 |
| CN103870654A (zh) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | 基于并行矩量法与物理光学混合的电磁散射仿真方法 |
| CN106355641A (zh) * | 2016-08-26 | 2017-01-25 | 西安电子科技大学 | 三维粗糙目标雷达散射截面的仿真方法 |
| CN108614251A (zh) * | 2016-12-11 | 2018-10-02 | 南京理工大学 | 多涂覆目标的雷达散射截面可视化计算方法 |
Non-Patent Citations (5)
| Title |
|---|
| CHUAN YIN等: ""An improved model for computing the scattering amplitude coefficient of the chaff cloud"", 《2014 IEEE WORKSHOP ON ADVANCED RESEARCH AND TECHNOLOGY IN INDUSTRY APPLICATIONS》 * |
| JINGSHUANG LU等: ""Research on meliorative average RCS of single chaff"", 《WORLD AUTOMATION CONGRESS 2012》 * |
| TANG, B.等: ""The correlation characteristics of channel matrix of chaff-supported MIMO system"", 《IEEE ANT. & WIRELESS PROPAGAT》 * |
| 曲长文 等: ""单根任意长度箔条双站散射截面积研究"", 《海军航空工程学院学报》 * |
| 沈允春 等: ""箔条的雷达截面积 "", 《系统工程与电子技术》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110502782A (zh) * | 2019-07-04 | 2019-11-26 | 西安电子科技大学 | 基于分区矢量输运的大体量箔条云电磁散射测定方法 |
| CN110502782B (zh) * | 2019-07-04 | 2022-04-26 | 西安电子科技大学 | 基于分区矢量输运的大体量箔条云电磁散射测定方法 |
| CN111199097A (zh) * | 2019-12-25 | 2020-05-26 | 西安电子科技大学 | 瞬态矢量辐射传输理论箔条云散射处理方法 |
| CN111199097B (zh) * | 2019-12-25 | 2022-09-13 | 西安电子科技大学 | 瞬态矢量辐射传输理论箔条云散射处理方法 |
| CN112364435A (zh) * | 2020-11-19 | 2021-02-12 | 中国空气动力研究与发展中心超高速空气动力研究所 | 多次抛撒箔条云气动融合轨道再入实时跟踪统计模拟方法 |
| CN113759322A (zh) * | 2021-09-08 | 2021-12-07 | 中国人民解放军海军航空大学 | 一种用于对抗雷达的箔条云rcs预测方法及系统 |
| CN113759322B (zh) * | 2021-09-08 | 2023-08-11 | 中国人民解放军海军航空大学 | 一种用于对抗雷达的箔条云rcs预测方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109215072B (zh) | 2021-07-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109215072A (zh) | 基于tanh-sinh积分法的箔条云RCS获取方法 | |
| CN102411136B (zh) | 一种扩展基线解模糊的相位干涉仪测向方法 | |
| CN101692257B (zh) | 一种复杂曲面的配准方法 | |
| Fan et al. | An improved backward SBR-PO/PTD hybrid method for the backward scattering prediction of an electrically large target | |
| CN103530469B (zh) | 一种角反射器雷达截面积的计算方法 | |
| CN102419430B (zh) | 一种基于平行基线的圆阵相位干涉仪二维测向方法 | |
| CN104036078B (zh) | 一种基于安装高度的阵列天线辐射和散射性能综合设计方法 | |
| CN104122527A (zh) | 一种基于查表法的圆阵相位干涉仪宽带瞬时测向方法 | |
| CN108663653A (zh) | 基于l形电磁矢量传感器阵列的波达方向估计方法 | |
| CN104363650B (zh) | 一种野外条件下无线传感器网络定位优化方法 | |
| CN104200074A (zh) | 快速获取目标电磁散射特性的多层复波束方法 | |
| CN104239678B (zh) | 一种实现干涉仪测向定位的方法和装置 | |
| CN104933213A (zh) | 基于空间映射的大规模相控天线阵列宽角扫描优化方法 | |
| CN103177193A (zh) | 薄介质涂敷的金属旋转对称目标电磁散射快速计算方法 | |
| US9864042B2 (en) | Optimizing storage and usage of angle-of-arrival heatmaps | |
| CN103487791A (zh) | 一种基于散射中心矩阵的rcs转换方法 | |
| CN109633527A (zh) | 基于低秩及几何约束的嵌入式平面麦克风阵列声源测向方法 | |
| CN201576308U (zh) | 一种复杂曲面的配准系统 | |
| CN110907925B (zh) | 一种高频地波雷达双站模型下的权重定位方法 | |
| CN105184033B (zh) | 基于阶数步进金属目标的宽频带电磁特性快速预估方法 | |
| CN106997044B (zh) | 一种雷达模拟中天线方向性加权因子计算方法和装置 | |
| CN104778293B (zh) | 非均匀介质目标电磁散射的体积分Nystrom分析方法 | |
| CN110083904A (zh) | 基于gpu加速的量子雷达散射截面计算方法 | |
| CN109813968A (zh) | 一种阵列天线方向图的测量方法、设备、系统以及计算机可读存储介质 | |
| CN105242264B (zh) | 一种平面阵列多目标角度高分辨实现方法 |
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: 20230113 Address after: 312000 No. 16, Meixi Road, Xinchang County, Shaoxing City, Zhejiang Province (Building 1 and 2) Patentee after: ZHEJIANG JINMA AUTOMATION TECHNOLOGY Co.,Ltd. Address before: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang Patentee before: HANGZHOU DIANZI University |
|
| TR01 | Transfer of patent right |