CN104698444B - 确定阵元位置误差对阵列天线三维成像质量影响的方法 - Google Patents
确定阵元位置误差对阵列天线三维成像质量影响的方法 Download PDFInfo
- Publication number
- CN104698444B CN104698444B CN201310716836.7A CN201310716836A CN104698444B CN 104698444 B CN104698444 B CN 104698444B CN 201310716836 A CN201310716836 A CN 201310716836A CN 104698444 B CN104698444 B CN 104698444B
- Authority
- CN
- China
- Prior art keywords
- peak
- sigma
- side lobe
- spread function
- point target
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000001914 filtration Methods 0.000 claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000009825 accumulation Methods 0.000 claims description 10
- 230000001427 coherent effect Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 230000001186 cumulative effect Effects 0.000 claims description 2
- 230000001902 propagating effect Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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/40—Means for monitoring or calibrating
-
- 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/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种确定阵元位置误差对阵列天线三维成像质量影响的方法。该方法包括:在三维空间中分布阵列天线和点目标P,接收输入的阵元位置测量误差(Δxnm,Δynm,Δznm);基于该阵列天线及该阵元位置测量误差(Δxnm,Δynm,Δznm)建立阵列天线信号模型,生成所述点目标P的回波信号srnm(t);对所述点目标P的回波信号srnm(t)进行匹配滤波处理得到sonm(t);根据匹配滤波后的回波信号sonm(t)计算点目标P在跨航向、航迹向和高程向上的点扩展函数、峰值旁瓣比和积分旁瓣比。本发明能够总体上定量分析阵元位置误差对三维图像成像质量的影响,从而为测量设备及运动补偿算法的选择提供参考依据。
Description
技术领域
本发明涉及雷达技术领域,尤其涉及一种确定阵元位置误差对阵列天线成像质量影响的方法。
背景技术
合成孔径雷达在遥感领域得到了广泛的应用。在实际应用中,大多采用机载合成孔径雷达进行遥感成像。但是由于飞机在飞行的过程中受大气扰动等外部环境影响,不能实现直线匀速运动,这会导致合成孔径雷达的回波相位噪声和误差增大,相应地会引起合成孔径雷达成像指标的恶化。虽然目前已经存在许多考虑这些误差并相应地进行运动补偿的技术方案,但是器件本身存在的阵元位置误差也会导致图像质量恶化。现有技术中一般都考虑器件在一维或二维成像时阵元位置误差对成像质量的影响。
然而,在实现本发明的过程中,申请人发现现有技术还无法确定阵元位置误差对阵列天线三维成像质量影响。
发明内容
(一)要解决的技术问题
鉴于上述技术问题,本发明提供了一种确定阵元位置误差对阵列天线三维成像质量影响的方法。
(二)技术方案
根据本发明的一个方面,提供了一种确定阵元位置误差对阵列天线三维成像质量影响的方法。该方法包括:步骤A,在三维空间中分布阵列天线和点目标P,接收输入的阵元位置测量误差(Δxnm,Δynm,Δznm);步骤B,基于该阵列天线及该阵元位置测量误差(Δxnm,Δynm,Δznm)建立阵列天线信号模型,生成点目标P的回波信号srnm(t);步骤C,对点目标P的回波信号srnm(t)进行匹配滤波处理得到sonm(t);步骤D,根据匹配滤波后的回波信号sonm(t)计算点目标P在跨航向上的点扩展函数Sx(xi)、峰值旁瓣比PSLRx和积分旁瓣比ISLRx;步骤E,根据匹配滤波后的回波信号sonm(t)计算点目标P在航迹向上的点扩展函数Sy(yk)、峰值旁瓣比PSLRy和积分旁瓣比ISLRy;以及步骤F,根据匹配滤波后的回波信号sonm(t)计算点目标P在高程向上的点扩展函数Sz(zl)、峰值旁瓣比PSLRz和积分旁瓣比ISLRz。其中,其中,点目标P在跨航向上的点扩展函数Sx(xi)及峰值旁瓣比PSLRx和积分旁瓣比ISLRx,在航迹向上的点扩展函数Sy(yk)及峰值旁瓣比PSLRy和积分旁瓣比ISLRy,在高程向上的点扩展函数Sz(zl)及峰值旁瓣比PSLRz和积分旁瓣比ISLRz即为阵元位置测量误差(Δxnm,Δynm,Δznm)对阵列天线成像质量的影响。
(三)有益效果
从上述技术方案可以看出,本发明确定阵元位置误差对阵列天线成像质量影响的方法能够总体上定量分析阵元位置误差对三维图像成像质量的影响,从而为测量设备及运动补偿算法的选择提供参考依据,以便得到符合要求的三维图像。
附图说明
图1是本发明提供的阵元位置误差对阵列天线成像质量影响的分析流程图;
图2是阵列天线雷达成像阵元位置误差几何示意图;
图3是点目标P在跨航向的点扩展函数求解几何示意图;
图4是点目标P在航迹向的点扩展函数求解几何示意图;
图5是点目标P在高程向的点扩展函数求解几何示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。需要说明的是,在附图或说明书描述中,相似或相同的部分都使用相同的图号。附图中未绘示或描述的实现方式,为所属技术领域中普通技术人员所知的形式。另外,虽然本文可提供包含特定值的参数的示范,但应了解,参数无需确切等于相应的值,而是可在可接受的误差容限或设计约束内近似于相应的值。
在本发明的一个示例性实施例中,提供了一种确定阵元位置误差对阵列天线成像质量影响的方法。图1是本发明提供的阵元位置误差对阵列天线成像质量影响的分析流程图。请参照图1,本实施例包括:
步骤A,在三维空间中分布阵列天线和点目标P,接收输入的阵元位置测量误差(Δxnm,Δynm,Δznm);
图2是阵列天线雷达成像阵元位置误差几何示意图。如图2所示,x轴代表跨航向,y轴代表航迹向,z轴代表高程向;阵列天线包含N个沿跨航向分布、间距为dv的阵元An,该N个阵元An以速度V沿y轴正方向运动并在y轴上有M个采样位置,采样间距为du。在第m个航迹向采样位置处,每个阵元位置的测量值为(xnm,ynm,znm),输入的测量误差在跨航向、航迹向和高程向的分量分别为Δxnm、Δynm和Δznm,则每个阵元的实际位置为(xnm+Δxnm,ynm+Δynm,znm+Δznm),m=1,2,…,M,n=1,2,…,N;观测场景中包含一个点目标P,假定其坐标为(xp,yp,zp)。
步骤B,基于该阵列天线及该阵元位置测量误差(Δxnm,Δynm,Δznm)建立阵列天线信号模型,生成点目标P的回波信号srnm(t);
回波信号srnm(t)的表达式如下:
其中,σ为点目标P的散射系数;t为信号传播方向的时间变量;Rnm为天线阵元到点目标P的距离;T为脉冲宽度;fc为信号中心频率;c为电磁波在空气中传播的速度;K为调制频率,等于信号带宽B除以脉冲宽度T;rect(·)为矩形窗函数:当t∈[-T/2,T/2],当t为其他值,
步骤C,对点目标P的回波信号srnm(t)进行匹配滤波处理得到sonm(t);
该步骤C进一步包括:
子步骤C1,将回波信号srnm(t)与本振信号sl(t)做下混频(即与sl(t)的共轭sl *(t)相乘),得到基带回波信号srbnm(t);
sl(t)=exp(2jπfct) (3)
sl *(t)=exp(-2jπfct) (4)
子步骤C2,对基带回波信号srbnm(t)进行傅里叶变换,得到srbnm(f);
子步骤C3,将srbnm(f)乘以参考信号sMF(f)进行匹配滤波,得到sonm(f),其中:
子步骤C4,对sonm(f)再进行傅里叶逆变换,得到匹配滤波后的信号sonm(t),其中:
sonm(t)=σ sin c(t-2Rnm/c)exp(-4jπfcRnm/c) (9)
步骤D,根据匹配滤波后的回波信号sonm(t)计算点目标P在跨航向上的点扩展函数Sx(xi)、峰值旁瓣比PSLRx和积分旁瓣比ISLRx;
该步骤D进一步包括如下步骤:
子步骤D1,计算点扩展函数Sx(xi)的位置坐标(xi=xp+i·dx,yp,zp),其中,dx=dv·Nx/2,i∈[-Nx·N,Nx·N-1],i∈Z,ceil(·)表示对小数向上取整,H=max(znm),为阵元在跨航向的波束宽度,如图3所示;
子步骤D2,沿x轴对匹配滤波后的回波信号sonm(t)进行相干累积,计算点扩展函数Sx(xi)在x轴上的分量Sxm(xi);
该子步骤D2具体为:根据Sx(xi)位置(xi,yp,zp)到天线阵元的距离计算这些位置在sonm(t)中对应的数值,并乘以补偿相位再沿着x轴累积求和,得到Sxm(xi)如下:
其中,
子步骤D3,沿y轴方向对点扩展函数Sx(xi)在x轴上的分量Sxm(xi)进行累积求和,得到点目标P在跨航向上的点扩展函数Sx(xi):
子步骤D4,由点目标P在跨航向上的点扩展函数Sx(xi)计算点目标P在跨航向上的峰值旁瓣比PSLRx及积分旁瓣比ISLRx;
该子步骤D4具体包括:
子分步骤D4a,对点目标P在跨航向上的点扩展函数Sx(xi)取模,得|Sx(xi)|,找出其最大值Px=|Sx(xxpeak)|即主瓣峰值,xxpeak即主瓣峰值位置;在xxpeak左侧找出离xxpeak最近的|Sx(xi)|的一个谷底(极小值),对应位置记为xxleft;在xxpeak右侧找出离xxpeak最近的|Sx(xi)|的一个谷底(极小值),对应位置记为xxright;在xi(i∈[-Nx·N,xleft]∪[xright,Nx·N-1])范围内找出|Sx(xi)|的最大值Lx=|Sx(Lxlobe)|即副瓣峰值最大值,xxlobe即副瓣峰值最大值位置。则峰值旁瓣比
子分步骤D4b,分别在xi(i∈[xleft,xright])和xi(i∈[-Nx·N,Nx·N-1])范围内对|Sx(xi)|2进行求和,得到主瓣功率和总功率则积分旁瓣比
步骤E,根据匹配滤波后的回波信号sonm(t)计算点目标P在航迹向上的点扩展函数Sy(yk)、峰值旁瓣比PSLRy和积分旁瓣比ISLRy;
该步骤E包括如下步骤:
子步骤E1:计算点扩展函数Sy(yk)的位置坐标(xp,yk=yp+k·dy,zp),其中dy=du/2My,k∈[-My·M,My·M-1],k∈Z,My>0,My∈Z,如图4所示,其中Z代表整数;
子步骤E2:沿x轴对匹配滤波后的回波信号sonm(t)进行相干累积,计算点扩展函数Sy(yk)在x轴上的分量Sym(yk):根据Sy(yk)位置(xp,yk,zp)到天线阵元的距离计算这些位置在sonm(t)中对应的数值,并乘以补偿相位再沿着x轴累积求和,得到Sym(yk)如下:
其中,
子步骤E3:沿y轴方向对Sym(yk)进行累积求和,得到点目标P在航迹向上的点扩展函数Sy(yk):
子步骤E4:计算点目标P在航迹向上的峰值旁瓣比PSLRy及积分旁瓣比ISLRy;
该子步骤E4具体包括:
子分步骤E4a,对Sy(yk)取模,得|Sy(yk)|,找出其最大值Py=|Sy(yypeak)|即主瓣峰值,yypeak即主瓣峰值位置。在yypeak左侧找出离yypeak最近的|Sy(yk)|的一个谷底(极小值),对应位置记为yyleft,在yyleft右侧找出离yypeak最近的|Sy(yk)|的一个谷底(极小值),对应位置记为yyright。在yk(k∈[-My·M,yleft]∪[yright,My·M-1])范围内找出|Sy(yk)|的最大值Ly=|Sy(yylobe)|即副瓣峰值最大值,yylobe即副瓣峰值最大值位置。则峰值旁瓣比
子分步骤E4b,分别在yk(k∈[yleft,yright])和yk(k∈[-My·M,My·M-1])范围内对|Sy(yk)|2进行求和,得到主瓣功率和总功率则积分旁瓣比
步骤F,根据匹配滤波后的回波信号sonm(t)计算点目标P在高程向上的点扩展函数Sz(zl)、峰值旁瓣比PSLRz和积分旁瓣比ISLRz;
该步骤F具体包括如下子步骤:
子步骤F1:计算Sz(zl)的位置坐标,x轴坐标为xp,y轴坐标为yp,z轴坐标为zl=zp+l·dz,其中,dz=(H-zp)/2Mz,H=max(znm),Mz>0,Mz∈Z,l∈[-Mz,Mz-1],l∈Z,如图5所示;
子步骤F2:沿x轴对sonm(t)进行相干累积,计算Sz(zl)在x轴上的分量Szm(zl):根据Sz(zl)位置(xp,yp,zl)到天线阵元的距离计算这些位置在sonm(t)中对应的数值,并乘以补偿相位再沿着x轴累积求和,得到Szm(zl)如下:
其中,
子步骤F3:沿y轴方向对Szm(zl)进行累积求和,得到点目标P在高程向上的点扩展函数Sz(zl):
子步骤F4:计算点目标P在高程向上的峰值旁瓣比PSLRz及积分旁瓣比ISLRz;
该子步骤F4具体包括:
子分步骤F4a,对Sz(zl)取模得|Sz(zl)|,找出其最大值Pz=|Sz(zzpeak)|即主瓣峰值,zzpeak即主瓣峰值位置。在zzpeak左侧找出离zzpeak最近的|Sz(zl)|的一个谷底(极小值),对应位置记为zzleft,在zzleft右侧找出离zzpeak最近的|Sz(zl)|的一个谷底(极小值),对应位置记为zzright。在zl(l∈[-Mz,zleft]∪[zright,Mz-1])范围内找出|Sz(zl)|的最大值Lz=|Sz(zzlobe)|即副瓣峰值最大值,zzlobe即副瓣峰值最大值位置。则峰值旁瓣比
子分步骤F4b,分别在zl(l∈[zleft,zright])和zl(l∈[-Mz,Mz-1])范围内对|Sz(zl)|2进行求和,得到主瓣功率和总功率则积分旁瓣比
其中,点目标P在跨航向上的点扩展函数Sx(xi)及峰值旁瓣比PSLRx和积分旁瓣比ISLRx,在航迹向上的点扩展函数Sy(yk)及峰值旁瓣比PSLRy和积分旁瓣比ISLRy,在高程向上的点扩展函数Sz(zl)及峰值旁瓣比PSLRz和积分旁瓣比ISLRz即为阵元位置测量误差(Δxnm,Δynm,Δznm)对阵列天线成像质量的影响。
本发明能够总体上定量分析阵元位置误差对三维图像成像质量的影响,从而为测量设备及运动补偿算法的选择提供参考依据,以便得到符合要求的三维图像。
至此,已经结合附图对本实施例进行了详细描述。依据以上描述,本领域技术人员应当对本发明确定阵元位置误差对阵列天线成像质量影响的方法有了清楚的认识。
此外,上述对各元件和方法的定义并不仅限于实施方式中提到的各种具体结构、形状或方式,本领域的普通技术人员可对其进行简单地熟知地替换,例如:
(1)天线回波信号还可以用其他非脉冲LFM信号形式表示;
(2)主瓣右侧零点位置可以用主瓣右侧3dB幅度处来代替。
综上所述,本发明提供一种可利用点目标成像的点扩展函数及其峰值旁瓣比和积分旁瓣比来分析阵元位置误差对阵列天线成像质量影响的方式。即使测量设备在测量阵列天线阵元位置时有未知的测量误差,本发明也可以分析该误差对成像质量的影响,从而为测量设备及后续运动补偿算法的选择提供参考依据,可应用于飞行导航、灾害监测、战场指挥、城市规划等领域。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种确定阵元位置误差对阵列天线三维成像质量影响的方法,其特征在于,包括:
步骤A,在三维空间中分布阵列天线和点目标P,接收输入的阵元位置测量误差(Δxnm,Δynm,Δznm);
步骤B,基于该阵列天线及该阵元位置测量误差(Δxnm,Δynm,Δznm)建立阵列天线信号模型,生成所述点目标P的回波信号srnm(t);
步骤C,对所述点目标P的回波信号srnm(t)进行匹配滤波处理得到sonm(t);
步骤D,根据匹配滤波后的回波信号sonm(t)计算所述点目标P在跨航向上的点扩展函数Sx(xi)、峰值旁瓣比PSLRx和积分旁瓣比ISLRx;
步骤E,根据匹配滤波后的回波信号sonm(t)计算点目标P在航迹向上的点扩展函数Sy(yk)、峰值旁瓣比PSLRy和积分旁瓣比ISLRy;以及
步骤F,根据匹配滤波后的回波信号sonm(t)计算点目标P在高程向上的点扩展函数Sz(zl)、峰值旁瓣比PSLRz和积分旁瓣比ISLRz;
其中,所述点目标P在跨航向上的点扩展函数Sx(xi)及峰值旁瓣比PSLRx和积分旁瓣比ISLRx,在航迹向上的点扩展函数Sy(yk)、峰值旁瓣比PSLRy和积分旁瓣比ISLRy,在高程向上的点扩展函数Sz(zl)及峰值旁瓣比PSLRz和积分旁瓣比ISLRz即为阵元位置测量误差(Δxnm,Δynm,Δznm)对阵列天线成像质量的影响;其中,
所述步骤A中,所述阵列天线包含N个沿跨航向分布、间距为dv的阵元An,该N个阵元An以速度V沿y轴正方向运动并在y轴上有M个采样位置,采样间距为du,测量获得阵列天线各个阵元的位置(xnm,ynm,znm),并假定点目标P的坐标为(xp,yp,zp);
所述步骤B中,所述点目标P的回波信号srnm(t)为:
其中,σ为点目标P的散射系数;t为信号传播方向的时间变量;Rnm为天线阵元到点目标P的距离;T为脉冲宽度;fc为信号中心频率;c为电磁波在空气中传播的速度;K为调制频率,等于信号带宽B除以脉冲宽度T;rect(·)为矩形窗函数:当t∈[-T/2,T/2],当t为其他值,
所述步骤C进一步包括:
子步骤C1,将回波信号srnm(t)与本振信号sl(t)做下混频,得到基带回波信号srbnm(t);
sl(t)=exp(2jπfct)
sl *(t)=exp(-2jπfct)
子步骤C2,对基带回波信号srbnm(t)进行傅里叶变换,得到srbnm(f);
子步骤C3,将srbnm(f)乘以参考信号sMF(f)进行匹配滤波,得到sonm(f),其中:
子步骤C4,对sonm(f)再进行傅里叶逆变换,得到匹配滤波后的信号sonm(t),其中:
sonm(t)=σsinc(t-2Rnm/c)exp(-4jπfcRnm/c);
所述步骤D具体包括:
子步骤D1,计算点扩展函数Sx(xi)的位置坐标(xi=xp+i·dx,yp,zp),其中,dx=dv·Nx/2,i∈[-Nx·N,Nx·N-1],i∈Z,ceil(·)表示对小数向上取整,H=max(znm),为阵元在跨航向的波束宽度;
子步骤D2,沿x轴对匹配滤波后的回波信号sonm(t)进行相干累积,计算点扩展函数Sx(xi)在x轴上的分量Sxm(xi),其具体为:根据Sx(xi)位置(xi,yp,zp)到天线阵元的距离Rxinm,计算这些位置在sonm(t)中对应的数值,并乘以补偿相位再沿着x轴累积求和,得到Sxm(xi)如下:
其中,
子步骤D3,沿y轴方向对点扩展函数Sx(xi)在x轴上的分量Sxm(xi)进行累积求和,得到点目标P在跨航向上的点扩展函数Sx(xi):
子步骤D4,由点目标P在跨航向上的点扩展函数Sx(xi)计算点目标P在跨航向上的峰值旁瓣比PSLRx及积分旁瓣比ISLRx。
2.根据权利要求1所述的方法,其特征在于,所述子步骤D3包括:
子分步骤D4a,对点目标P在跨航向上的点扩展函数Sx(xi)取模,得|Sx(xi)|,找出其最大值Px=|Sx(xxpeak)|即主瓣峰值,xxpeak即主瓣峰值位置;在xxpeak左侧找出离xxpeak最近的|Sx(xi)|的一个谷底,对应位置记为xxleft;在xxpeak右侧找出离xxpeak最近的|Sx(xi)|的一个谷底,对应位置记为xxright;在xi(i∈[-Nx·N,xleft]∪[xright,Nx·N-1])范围内找出|Sx(xi)|的最大值Lx=|Sx(xxlobe)|即副瓣峰值最大值,xxlobe即副瓣峰值最大值位置,则峰值旁瓣比
子分步骤D4b,分别在xi(i∈[xleft,xright])和xi(i∈[-Nx·N,Nx·N-1])范围内对|Sx(xi)|2进行求和,得到主瓣功率和总功率则积分旁瓣比
3.根据权利要求1所述的方法,其特征在于,所述步骤E具体包括:
子步骤E1:计算点扩展函数Sy(yk)的位置坐标(xp,yk=yp+k·dy,zp),其中,dy=du/2My,k∈[-My·M,My·M-1],k∈Z,My>0;
子步骤E2:沿x轴对匹配滤波后的回波信号sonm(t)进行相干累积,计算点扩展函数Sy(yk)在x轴上的分量Sym(yk);根据Sy(yk)位置(xp,yk,zp)到天线阵元的距离计算这些位置在sonm(t)中对应的数值,并乘以补偿相位再沿着x轴累积求和,得到Sym(yk)如下:
其中,
子步骤E3:沿y轴方向对Sym(yk)进行累积求和,得到点目标P在航迹向上的点扩展函数Sy(yk):
子步骤E4:计算点目标P在航迹向上的峰值旁瓣比PSLRy及积分旁瓣比ISLRy。
4.根据权利要求3所述的方法,其特征在于,所述子步骤E4具体包括:
子分步骤E4a,对Sy(yk)取模,得|Sy(yk)|,找出其最大值Py=|Sy(yypeak)|即主瓣峰值,yypeak即主瓣峰值位置,在yypeak左侧找出离yypeak最近的|Sy(yk)|的一个谷底,对应位置记为yyleft,在yyleft右侧找出离yypeak最近的|Sy(yk)|的一个谷底,对应位置记为yyright;在yk(k∈[-My·M,yleft]∪[yright,My·M-1])范围内找出|Sy(yk)|的最大值Ly=|Sy(yylobe)|即副瓣峰值最大值,yylobe即副瓣峰值最大值位置,则峰值旁瓣比
子分步骤E4b,分别在yk(k∈[yleft,yright])和yk(k∈[-My·M,My·M-1])范围内对|Sy(yk)|2进行求和,得到主瓣功率和总功率则积分旁瓣比
5.根据权利要求1所述的方法,其特征在于,所述步骤F具体包括:
子步骤F1:计算Sz(zl)的位置坐标,x轴坐标为xp,y轴坐标为yp,z轴坐标为zl=zp+l·dz,其中,dz=(H-zp)/2Mz,H=max(znm),Mz>0,Mz∈Z,l∈[-Mz,Mz-1],l∈Z;
子步骤F2:沿x轴对sonm(t)进行相干累积,计算Sz(zl)在x轴上的分量Szm(zl):根据Sz(zl)位置(xp,yp,zl)到天线阵元的距离计算这些位置在sonm(t)中对应的数值,并乘以补偿相位再沿着x轴累积求和,得到Szm(zl)如下:
其中,
子步骤F3:沿y轴方向对Szm(zl)进行累积求和,得到点目标P在高程向上的点扩展函数Sz(zl):
子步骤F4:计算点目标P在高程向上的峰值旁瓣比PSLRz及积分旁瓣比ISLRz。
6.根据权利要求5所述的方法,其特征在于,所述子步骤F4具体包括:
子分步骤F4a,对Sz(zl)取模得|Sz(zl)|,找出其最大值Pz=|Sz(zzpeak)|即主瓣峰值,zzpeak即主瓣峰值位置,在zzpeak左侧找出离zzpeak最近的|Sz(zl)|的一个谷底,对应位置记为zzleft,在zzleft右侧找出离zzpeak最近的|Sz(zl)|的一个谷底,对应位置记为zzright,在zl(l∈[-Mz,zleft]∪[zright,Mz-1])范围内找出|Sz(zl)|的最大值Lz=|Sz(zzlobe)|即副瓣峰值最大值,zzlobe即副瓣峰值最大值位置,则峰值旁瓣比
子分步骤F4b,分别在zl(l∈[zleft,zright])和zl(l∈[-Mz,Mz-1])范围内对|Sz(zl)|2进行求和,得到主瓣功率和总功率则积分旁瓣比
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310716836.7A CN104698444B (zh) | 2013-12-23 | 2013-12-23 | 确定阵元位置误差对阵列天线三维成像质量影响的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310716836.7A CN104698444B (zh) | 2013-12-23 | 2013-12-23 | 确定阵元位置误差对阵列天线三维成像质量影响的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104698444A CN104698444A (zh) | 2015-06-10 |
CN104698444B true CN104698444B (zh) | 2017-04-26 |
Family
ID=53345767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310716836.7A Active CN104698444B (zh) | 2013-12-23 | 2013-12-23 | 确定阵元位置误差对阵列天线三维成像质量影响的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104698444B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108132466B (zh) * | 2017-12-22 | 2021-11-09 | 内蒙古工业大学 | 一种机载阵列天线下视三维成像方法和系统 |
CN109671153B (zh) * | 2018-12-05 | 2022-12-30 | 安徽慧软科技有限公司 | 一种复杂器官轮廓线的三维重建方法 |
CN116086586B (zh) * | 2023-04-11 | 2023-06-20 | 中船重工海目测试技术(海南)有限公司 | 一种基于粒子群优化阵列处理的船舶辐射噪声测量方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1664611A (zh) * | 2005-03-24 | 2005-09-07 | 武汉大学 | 一种基于非直线天线阵列的无源通道校正方法 |
GB2489072A (en) * | 2011-02-09 | 2012-09-19 | Kongsberg Maritime As | A transmit/receive array arranged for the suppression of sidelobes |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6954173B2 (en) * | 2003-07-02 | 2005-10-11 | Raytheon Company | Techniques for measurement of deformation of electronically scanned antenna array structures |
-
2013
- 2013-12-23 CN CN201310716836.7A patent/CN104698444B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1664611A (zh) * | 2005-03-24 | 2005-09-07 | 武汉大学 | 一种基于非直线天线阵列的无源通道校正方法 |
GB2489072A (en) * | 2011-02-09 | 2012-09-19 | Kongsberg Maritime As | A transmit/receive array arranged for the suppression of sidelobes |
Non-Patent Citations (2)
Title |
---|
MIMO雷达成像算法研究;王怀军;《中国博士学位论文全文数据库 信息科技辑》;20110815(第8期);第37-38页 * |
基于空间谱域分析的雷达成像技术研究;韩兴斌;《中国博士学位论文全文数据库 信息科技辑》;20090715(第7期);第92-93、101-105页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104698444A (zh) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102645651B (zh) | 一种sar层析超分辨成像方法 | |
CN103487802B (zh) | 扫描雷达角超分辨成像方法 | |
CN106908787A (zh) | 一种实波束扫描雷达前视角超分辨率成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN103207380B (zh) | 基于二维频域稀疏约束的宽带目标测向方法 | |
EP3477331A1 (en) | Below-noise after transmit (bat) chirp radar | |
CN104536000A (zh) | 一种实波束扫描雷达角超分辨方法 | |
CN103777186B (zh) | 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 | |
CN103885058B (zh) | 一种利用稀疏统计特性的机载雷达前视超分辨成像方法 | |
CN104122549B (zh) | 基于反卷积的雷达角超分辨成像方法 | |
CN108226891B (zh) | 一种扫描雷达回波计算方法 | |
CN104698444B (zh) | 确定阵元位置误差对阵列天线三维成像质量影响的方法 | |
Huang et al. | Detection and fast motion parameter estimation for target with range walk effect based on new axis rotation moving target detection | |
Martín et al. | Aircraft tracking by means of the Acoustical Doppler Effect | |
EP3418770B1 (en) | Synthetic-aperture-radar signal processing device | |
Chen et al. | Compensation for high-frequency vibration of SAR imaging in the terahertz band based on linear chirplet transform and empirical mode decomposition | |
CN114814838A (zh) | 一种基于正交基斜距展开的合成孔径无源定位方法及装置 | |
Zhang et al. | An airborne/missile-borne array radar forward-looking imaging algorithm based on super-resolution method | |
CN106646418A (zh) | 一种基于改进omp的机载雷达杂波空时谱快速估计方法 | |
CN112859018A (zh) | 一种基于图像几何校正的视频sar成像方法 | |
CN113093187A (zh) | 一种无道路信息辅助的csar地面动目标跟踪方法 | |
Sjanic et al. | Simultaneous navigation and synthetic aperture radar focusing | |
Tang et al. | Forward-Looking Super-Resolution Imaging Based on Echo Denoising and Noise Weighting at Low SNR | |
CN111190173A (zh) | 一种基于预测值量测转换的相控阵雷达目标跟踪方法 | |
Wang et al. | Fast Iterative Wiener Filter-Based ISAR Imaging and Cross-Range Scaling With Periodically Gapped CPI |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |