CN103197285B - 一种用于合成孔径雷达成像的导航数据拟合方法 - Google Patents
一种用于合成孔径雷达成像的导航数据拟合方法 Download PDFInfo
- Publication number
- CN103197285B CN103197285B CN201310093771.5A CN201310093771A CN103197285B CN 103197285 B CN103197285 B CN 103197285B CN 201310093771 A CN201310093771 A CN 201310093771A CN 103197285 B CN103197285 B CN 103197285B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- mfrac
- fitting
- lgd
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000003384 imaging method Methods 0.000 title claims abstract description 15
- 238000012545 processing Methods 0.000 claims abstract description 18
- 239000011159 matrix material Substances 0.000 claims description 75
- 230000001133 acceleration Effects 0.000 claims description 20
- 238000010561 standard procedure Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 5
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 2
- 238000007781 pre-processing Methods 0.000 abstract description 10
- 230000001629 suppression Effects 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000007547 defect Effects 0.000 abstract description 2
- 239000000243 solution Substances 0.000 description 40
- 238000005516 engineering process Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000012086 standard solution Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出了一种用于合成孔径雷达成像的导航数据拟合方法,它是通过采用分段高阶拟合方法很好地保持了导航数据预处理结果的光滑性,降低了合成孔径雷达处理过程中的相位噪声,具有很强的噪声抑制能力;在进行速度数据拟合后,通过对输出的速度数据进行高阶积分,得到高精度的位置数据,克服了传统导航数据预处理方法中,通过对速度数据进行简单累加方法得到位置数据,而导致的数值计算误差。本发明可用于各种机载合成孔径雷达导航数据预处理。
Description
技术领域
本发明属于雷达技术领域,它特别涉及了合成孔径雷达成像技术领域。
背景技术
导航数据是处理非匀速直线运动合成孔径雷达数据的重要信息,对提高合成孔径雷达成像处理质量和精度具有重要意义。根据本发明人了解以及已发表的文献,例如:JianchengFangandXiaolinGong,“PredictiveIteratedKalmanFilterforINS/GPSIntegrationandItsApplicationtoSARMotionCompensation”,IEEETransactionsonInstrumentationandMeasurement,vol.59,no.4,April2010,导航数据已经成为合成孔径雷达非直线运动误差补偿的重要方法之一。
但是,在某些情况下,GPS、IMU等导航***获得的导航数据由于其数据输出格式差异,例如,某些IMU***提供瞬时速度而非瞬时位置数据,导航***测量误差,导航数据更新率与合成孔径雷达脉冲重复周期不匹配等问题,无法直接用于合成孔径雷达成像处理,而需要对导航数据进行预处理。目前的预处理方法主要采用线性插值技术。但是,由于线性插值方法不具有光滑特性,而合成孔径雷达实际运动轨迹,由于惯性的作用一般为光滑曲线,使得线性插值结果与合成孔径雷达实际天线相位中心轨迹之间存在处理误差。另外,由于线性插值方法不具有噪声抑制能力,该方法的处理结果可能会随着测量误差的增加而变差,影响高精度合成孔径雷达成像处理***的稳定性。
发明内容
为了克服传统导航数据预处理方法存在的处理误差以及不具备噪声抑制能力等问题,本发明提出了一种用于合成孔径雷达成像的导航数据拟合方法。该方法采用分段高阶拟合方法,在很好的保持了导航数据预处理结果的光滑性的同时还具有很好的噪声抑制能力。
为了方便描述本发明的内容,首先作以下术语定义:
定义1、导航数据
导航数据是指由GPS、INS、IMU等导航设备提供的,记录合成孔径雷达运动规律的数据。对于某些导航设备,其提供的导航数据为合成孔径雷达瞬时三维位置数据;而某些导航设备提供的导航数据为合成孔径雷达瞬时三维速度数据。详细内容可参考文献:王巍主编.光纤陀螺惯性***.北京:中国宇航出版社,2010。
定义2、勒让德多项式函数
勒让德多项式函数是在区间[-1,1]之间相互正交的一组多项式,记作:Lgd(t;k),其中,t为勒让德多项式函数的变量,k为勒让德多项式函数的阶数,其定义和数学表达式详见文献:《数学手册》编写组,数学手册,高等教育出版社,1979年。
另外,本发明还涉及勒让德多项式函数的导函数,记作DLgd(t;k),其定义为:
以及勒让德多项式函数的原函数,记作ILgd(t;k),其定义为:
定义3、求解解空间的正交基矩阵的方法
解空间是指满足线性方程组Ax=0的所有解构成的空间。解空间的正交基矩阵是指张成解空间的坐标基构成的矩阵。解空间的正交基矩阵可以利用matlab提供的null直接计算得到。详细内容可参考文献:张德丰,Matlab数值计算方法,机械工业出版社,2011年。
定义4、标准的求解线性方程组特解的方法
线性方程组的特解是指满足线性方程组Ax=b的任意一个解向量。标准的求解线性方程组特解的方法是指高斯消元法等线性代数中采用的计算线性方程特解的方法,详细内容可参考文献:王坤,周岩主编,线性代数,机械工业出版社,2012年。
定义5、标准的多项式最小二乘拟合方法
标准的多项式最小二乘拟合方法是指求解以多项式函数作为拟合基函数的最小二乘拟合方法。标准的多项式最小二乘拟合问题的求解可采用最优化方法实现,详细内容可参考文献:张德丰,Matlab数值计算方法,机械工业出版社,2011年。由于多项式函数具有不同的表示形式和定义区间,在拟合前需要确定拟合基函数。
由于本发明中前一拟合段的速度、加速度等作为当前拟合段的附加最优化约束条件,因此标准的多项式最小二乘拟合方法无法直接用于本发明的分段多项式拟合。但在分段迭代的第一步,由于缺少前一拟合段的速度、加速度等信息,可以采用标准的多项式最小二乘拟合方法进行导航数据拟合。
本发明提供了一种用于合成孔径雷达成像的导航数据拟合方法,它包括以下几个步骤:
步骤1、初始化导航数据处理参数和数据
为了对导航数据进行预处理,需要提供如下初始化参数和待处理数据,包括:导航数据更新率,记作FIMU;合成孔径雷达脉冲重复频率,记作PRF;导航速度数据,记作v(n),其中,n为导航数据更新序号,n=0,1,2,...,N,N为需要处理的导航数据数目;拟合阶数,记作Q,Q≥3,拟合段长度,记作L,L>Q;导航轨迹初始位置,记作s0。
步骤2、构造拟合矩阵和边界向量
利用公式
构造拟合矩阵,Φ为(L+1)×(Q+1)的矩阵,其中,表示k阶勒让德多项式函数在处的值,k为勒让德多项式函数的阶数,k=0,1,...,Q,i为拟合矩阵的行序号,i=0,1,...,L。
利用公式 构造边界矩阵,Mag为2×(Q+1)的矩阵,其中,DLgd(-1;k')为k'阶勒让德多项式函数的导函数在-1处的值,k'=0,1,...,Q。
步骤3、计算解空间的正交基矩阵
根据步骤2得到的边界矩阵Mag,采用求解解空间的正交基矩阵的传统方法,计算边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵,记作Ψ,Ψ为(Q+1)×(Q-1)的矩阵。
步骤4、初始段导航数据拟合
根据步骤1得到的导航速度数据v(n),选择导航速度数据中导航数据更新序号n=0,1,...,L的L+1个数据,以0阶到Q阶勒让德函数构成的函数簇作为拟合基函数,利用传统标准的多项式最小二乘拟合方法,对导航速度数据的前L+1个数据进行拟合,得到初始段导航数据的多项式拟合系数向量,记作C[1], 其中,为初始段导航数据的多项式拟合系数向量的第q个分量,q=0,1,...,Q,上标T表示矩阵的转置。
步骤5、计算重采样初始段位置数据
利用公式 其中,表示对的四舍五入取整,计算拟合段数据对应的合成孔径雷达脉冲重复周期数目,记作LSAR。根据步骤4得到的初始段导航数据的多项式拟合系数向量C[1],利用公式得到初始段与合成孔径雷达脉冲重复周期对应位置数据,记作s[1](n'),其中,ILgd(n';q')为q'阶勒让德函数的原函数,q'=0,1,...,Q,n'为与合成孔径雷达脉冲重复周期对应的归一化时刻序列,
其中,j表示与合成孔径雷达脉冲重复周期对应的归一化时刻序列序号,j=0,1,...,LSAR。
步骤6、选择需要拟合的导航数据
令当前拟合段序号l=2,根据步骤1得到的导航速度数据v(n),采用标准的文件读取方法,读取导航速度数据中导航数据更新序号n=(l-1)L,((l-1)L+1),...,lL的L+1个数据,得到当前段需要拟合的导航数据,记作v[l](n),并利用公式 将当前段需要拟合的导航数据排列为向量形式,得到当前段需要拟合的导航数据向量,记作v[l],v[l]为(L+1)×1的向量。
步骤7、计算边界方程特解
如果当前拟合段序号l=2,则利用步骤4得到的初始段导航数据的多项式拟合系数向量C[1],利用公式计算初始段的末速度,记作v[1](1),利用公式计算初始段的末加速度,记作a[1](1)。利用步骤2得到的边界矩阵Mag以及初始段的末速度v[1](1)和初始段的末加速度a[1](1),构造当前段的边界线性方程组
其中,c[l]为线性方程组的未知列向量;
采用传统标准的求解线性方程组特解的方法,求解当前段的边界线性方程组,得到当前段的边界线性方程组的特解,记作为(Q+1)×1的向量。
如果当前拟合段序号l>2,则利用前次迭代得到的导航数据的多项式拟合系数向量C[l-1],利用公式计算前次迭代段的末速度,记作v[l-1](1),利用公式计算前次迭代段的末加速度,记作a[l-1](1)。利用步骤2得到的边界矩阵Mag以及前次迭代段的末速度v[l-1](1)和前次迭代段的末加速度a[l-1](1),构造当前段的边界线性方程组
采用传统标准的求解线性方程组特解的方法,求解当前段的边界线性方程组,得到当前段的边界线性方程组的特解,记作为(Q+1)×1的向量。
步骤8、计算最优化系数
根据步骤2得到的拟合矩阵Φ,步骤3得到的边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ,步骤6得到的当前段需要拟合的导航数据向量v[l],步骤7得到的当前段的边界线性方程组的特解利用公式 计算当前段的最优化系数,记作 为(Q-1)×1的向量,其中,上标-1表示矩阵的逆。
步骤9、计算多项式拟合系数
根据步骤3得到的边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ,步骤7得到的当前段的边界线性方程组的特解步骤8得到的当前段的最优化系数利用公式计算当前段的导航数据的多项式拟合系数,记作C[l],C[l]为(Q+1)×1的向量, 其中,为第l段导航数据的多项式拟合系数向量的第q个分量,q=0,1,...,Q。步骤10、计算拟合/重采样后的位置数据
根据步骤5得到的与合成孔径雷达脉冲重复周期对应的归一化时刻序列n',第9步得到的当前段的导航数据的多项式拟合系数C[l],利用公式得到当前拟合段与合成孔径雷达脉冲重复周期对应位置数据,记作s[l](n'),其中,s[l-1](1)表示前一个拟合段得到的末端时刻的位置数据。
步骤11、拟合段迭代
令当前拟合段序号l分别为:3,4,...,...,对于当前拟合段序号l分别为:3,4,...,...,时,分别重复步骤6~步骤10,处理完所有导航数据,其中,表示对向下取整,得到满足合成孔径雷达成像处理需要的导航数据。
需要说明的是,为了叙述简洁,本发明只考虑了三维导航数据沿某个坐标轴的分量,选择三维导航数据沿其他坐标轴的分量,采用与本发明相同的方法即可实现对三维导航数据的预处理。另外,本发明步骤2中采用勒让德多项式函数构造拟合矩阵,将勒让德多项式函数改为其它函数簇,如三角函数簇等,不影响本发明的有效性。
本发明的创新点在于:
1、将分段插值的运算量小、适宜处理大批量数据的优点与拟合技术具有噪声抑制能力的优点相结合,提出了分段拟合的导航数据处理方法,克服了分段线性插值方法的固有缺陷。
本发明的优点在于:
1、分段高阶拟合方法很好地保持了导航数据预处理结果的光滑性,降低了合成孔径雷达处理过程中的相位噪声;
2、拟合方法的采用使得本发明方法具有很好的噪声抑制能力;
3、在进行速度数据拟合后,通过对输出的速度数据进行高阶积分,得到高精度的位置数据,克服了传统导航数据预处理方法中,通过对速度数据进行简单累加方法得到位置数据,而导致的数值计算误差。
本发明可用于各种机载合成孔径雷达导航数据预处理。
附图说明
图1本发明具体实施方式中采用的导航速度数据
图2本发明所提出方法得到的拟合后的导航速度数据
可以看出,拟合后的导航数据保持了实际合成孔径雷达运动过程中的光滑特征,并有效降低了导航速度数据中的噪声。
图3本发明方法流程图
具体实施方式
为了验证本发明所提出方法的有效性,申请人对该方法进行了仿真实验分析,具体实施步骤如下:
步骤1、初始化导航数据处理参数和数据
为了对导航数据进行预处理,需要提供如下初始化参数和待处理数据,包括:导航数据更新率FIMU=40Hz;合成孔径雷达脉冲重复频率PRF=1000Hz;导航速度数据v(n),其中,n为导航数据更新序号,n=0,1,2,...,1000;拟合阶数Q=9,拟合段长度L=100;导航轨迹初始位置s0。
步骤2、构造拟合矩阵和边界向量
利用公式
构造拟合矩阵,Φ为(L+1)×(Q+1)的矩阵,其中,表示k阶勒让德多项式函数在处的值,k为勒让德多项式函数的阶数,k=0,1,...,Q,i为拟合矩阵的行序号,i=0,1,...,L。
利用公式 构造边界矩阵,Mag为2×(Q+1)的矩阵,其中,DLgd(-1;k')为k'阶勒让德多项式函数的导函数在-1处的值,k'=0,1,...,Q。
步骤3、计算解空间的正交基矩阵
根据步骤2得到的边界矩阵Mag,采用求解解空间的正交基矩阵的现有方法,计算边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ。
步骤4、初始段导航数据拟合
根据步骤1得到的导航速度数据v(n),选择导航速度数据中导航数据更新序号n=0,1,...,100的101个数据,以0阶到9阶勒让德函数构成的函数簇作为拟合基函数,利用标准的多项式最小二乘拟合方法,对导航速度数据的前101个数据进行拟合,得到初始段导航数据的多项式拟合系数向量C[1], 其中,为初始段导航数据的多项式拟合系数向量的第q个分量,q=0,1,...,9,上标T表示矩阵的转置。
步骤5、计算重采样初始段位置数据
利用公式 其中,表示对的四舍五入取整,计算拟合段数据对应的合成孔径雷达脉冲重复周期数目LSAR=2500。根据步骤4得到的初始段导航数据的多项式拟合系数向量C[1],利用公式得到初始段与合成孔径雷达脉冲重复周期对应位置数据s[1](n'),其中,ILgd(n';q')为q'阶勒让德函数的原函数,q'=0,1,...,Q,n'为与合成孔径雷达脉冲重复周期对应的归一化时刻序列,
其中,j表示与合成孔径雷达脉冲重复周期对应的归一化时刻序列序号,j=0,1,...,2500。
步骤6、选择需要拟合的导航数据
令当前拟合段序号l=2,根据步骤1得到的导航速度数据v(n),,采用标准的文件读取方法,读取导航速度数据中导航数据更新序号n=100(l-1),(100(l-1)+1),...,100l的101个数据,得到当前段需要拟合的导航数据
v[l](n),并利用公式 将当前段需要拟合的导航数据排列为向量形式,得到当前段需要拟合的导航数据向量v[l]。
步骤7、计算边界方程特解
如果当前拟合段序号l=2,则利用步骤4得到的初始段导航数据的多项式拟合系数向量C[1],利用公式计算初始段的末速度,记作v[1](1),利用公式计算初始段的末加速度,记作a[1](1)。利用步骤2得到的边界矩阵Mag以及初始段的末速度v[1](1)和初始段的末加速度a[1](1),构造当前段的边界线性方程组
其中,c[l]为线性方程组的未知列向量。
采用标准的求解线性方程组特解的方法,求解当前段的边界线性方程组,得到当前段的边界线性方程组的特解,记作为10×1的向量。
如果当前拟合段序号l>2,则利用前次迭代得到的导航数据的多项式拟合系数向量C[l-1],利用公式计算前次迭代段的末速度,记作v[l-1](1),利用公式计算前次迭代段的末加速度,记作a[l-1](1)。利用步骤2得到的边界矩阵Mag以及前次迭代段的末速度v[l-1](1)和前次迭代段的末加速度a[l-1](1),构造当前段的边界线性方程组
其中,c[l]为线性方程组的未知列向量。
采用标准的求解线性方程组特解的方法,求解当前段的边界线性方程组,得到当前段的边界线性方程组的特解,记作为10×1的向量。
步骤8、计算最优化系数
根据步骤2得到的拟合矩阵Φ,步骤3得到的边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ,步骤6得到的当前段需要拟合的导航数据向量v[l],步骤7得到的当前段的边界线性方程组的特解利用公式 计算当前段的最优化系数其中,上标-1表示矩阵的逆。
步骤9、计算多项式拟合系数
根据步骤3得到的边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ,步骤7得到的当前段的边界线性方程组的特解步骤8得到的当前段的最优化系数利用公式计算当前段的多项式拟合系数C[l]。
步骤10、计算拟合/重采样后的位置数据
根据步骤5得到的与合成孔径雷达脉冲重复周期对应的归一化时刻序列n',第9步得到的当前段的导航数据的多项式拟合系数C[l],利用公式得到当前拟合段与合成孔径雷达脉冲重复周期对应位置数据s[l](n'),其中,s[l-1](1)表示前一个拟合段得到的末端时刻的位置数据。
步骤11、拟合段迭代
令当前拟合段序号l=3,4,...,100,重复步骤6-10,直到处理完所有导航数据。
经过以上步骤,即可得到满足合成孔径雷达成像处理需要的导航数据。
从中可以看出,本发明所采用方法具有很好的噪声抑制能力,并保持了实际合成孔径雷达运动过程中的光滑特征。
Claims (1)
1.一种用于合成孔径雷达成像的导航数据拟合方法,它包括以下几个步骤:
步骤1、初始化导航数据处理参数和数据
为了对导航数据进行预处理,需要提供如下初始化参数和待处理数据,包括:导航数据更新率,记作FIMU;合成孔径雷达脉冲重复频率,记作PRF;导航速度数据,记作v(n),其中,n为导航数据更新序号,n=0,1,2,...,N,N为需要处理的导航数据数目;拟合阶数,记作Q,Q≥3,拟合段长度,记作L,L>Q;导航轨迹初始位置,记作s0;
步骤2、构造拟合矩阵和边界向量
利用公式
构造拟合矩阵,Φ为(L+1)×(Q+1)的矩阵,其中,表示k阶勒让德多项式函数在处的值,k为勒让德多项式函数的阶数,k=0,1,...,Q,i为拟合矩阵的行序号,i=0,1,...,L;
利用公式 构造边界矩阵,Mag为2×(Q+1)的矩阵,其中,DLgd(-1;k')为k'阶勒让德多项式函数的导函数在-1处的值,k'=0,1,...,Q;
步骤3、计算解空间的正交基矩阵
根据步骤2得到的边界矩阵Mag,采用求解解空间的正交基矩阵的传统方法,计算边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵,记作Ψ,Ψ为(Q+1)×(Q-1)的矩阵;
步骤4、初始段导航数据拟合
根据步骤1得到的导航速度数据v(n),选择导航速度数据中导航数据更新序号n=0,1,...,L的L+1个数据,以0阶到Q阶勒让德函数构成的函数簇作为拟合基函数,利用传统标准的多项式最小二乘拟合方法,对导航速度数据的前L+1个数据进行拟合,得到初始段导航数据的多项式拟合系数向量,记作C[1], 其中,为初始段导航数据的多项式拟合系数向量的第q个分量,q=0,1,...,Q,上标T表示矩阵的转置;
步骤5、计算重采样初始段位置数据
利用公式 其中,表示对的四舍五入取整,计算拟合段数据对应的合成孔径雷达脉冲重复周期数目,记作LSAR;根据步骤4得到的初始段导航数据的多项式拟合系数向量C[1],利用公式得到初始段与合成孔径雷达脉冲重复周期对应位置数据,记作s[1](n'),其中,ILgd(n';q')为q'阶勒让德函数的原函数,q'=0,1,...,Q,n'为与合成孔径雷达脉冲重复周期对应的归一化时刻序列,
其中,j表示与合成孔径雷达脉冲重复周期对应的归一化时刻序列序号,j=0,1,...,LSAR;
步骤6、选择需要拟合的导航数据
令当前拟合段序号l=2,根据步骤1得到的导航速度数据v(n),采用标准的文件读取方法,读取导航速度数据中导航数据更新序号n=(l-1)L,((l-1)L+1),...,lL的L+1个数据,得到当前段需要拟合的导航数据,记作v[l](n),并利用公式 将当前段需要拟合的导航数据排列为向量形式,得到当前段需要拟合的导航数据向量,记作v[l],v[l]为(L+1)×1的向量;
步骤7、计算边界方程特解
如果当前拟合段序号l=2,则利用步骤4得到的初始段导航数据的多项式拟合系数向量C[1],利用公式计算初始段的末速度,记作v[1](1),利用公式计算初始段的末加速度,记作a[1](1);利用步骤2得到的边界矩阵Mag以及初始段的末速度v[1](1)和初始段的末加速度a[1](1),构造当前段的边界线性方程组
其中,c[l]为线性方程组的未知列向量;
采用传统标准的求解线性方程组特解的方法,求解当前段的边界线性方程组,得到当前段的边界线性方程组的特解,记作为(Q+1)×1的向量;
如果当前拟合段序号l>2,则利用前次迭代得到的导航数据的多项式拟合系数向量C[l-1],利用公式计算前次迭代段的末速度,记作v[l-1](1),利用公式计算前次迭代段的末加速度,记作a[l-1](1);利用步骤2得到的边界矩阵Mag以及前次迭代段的末速度v[l-1](1)和前次迭代段的末加速度a[l-1](1),构造当前段的边界线性方程组
采用传统标准的求解线性方程组特解的方法,求解当前段的边界线性方程组,得到当前段的边界线性方程组的特解,记作为(Q+1)×1的向量;
步骤8、计算最优化系数
根据步骤2得到的拟合矩阵Φ,步骤3得到的边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ,步骤6得到的当前段需要拟合的导航数据向量v[l],步骤7得到的当前段的边界线性方程组的特解利用公式 计算当前段的最优化系数,记作为(Q-1)×1的向量,其中,上标-1表示矩阵的逆;
步骤9、计算多项式拟合系数
根据步骤3得到的边界矩阵Mag为系数矩阵的线性方程组的解空间的正交基矩阵Ψ,步骤7得到的当前段的边界线性方程组的特解步骤8得到的当前段的最优化系数利用公式计算当前段的导航数据的多项式拟合系数,记作C[l],C[l]为(Q+1)×1的向量, 其中,为第l段导航数据的多项式拟合系数向量的第q个分量,q=0,1,...,Q;
步骤10、计算拟合/重采样后的位置数据
根据步骤5得到的与合成孔径雷达脉冲重复周期对应的归一化时刻序列n',第9步得到的当前段的导航数据的多项式拟合系数C[l],利用公式得到当前拟合段与合成孔径雷达脉冲重复周期对应位置数据,记作s[l](n'),其中,s[l-1](1)表示前一个拟合段得到的末端时刻的位置数据;步骤11、拟合段迭代
令当前拟合段序号l分别为:3,4,...,...,对于当前拟合段序号l分别为:3,4,...,...,时,分别重复步骤6~步骤10,处理完所有导航数据,其中,表示对向下取整,得到满足合成孔径雷达成像处理需要的导航数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310093771.5A CN103197285B (zh) | 2013-03-22 | 2013-03-22 | 一种用于合成孔径雷达成像的导航数据拟合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310093771.5A CN103197285B (zh) | 2013-03-22 | 2013-03-22 | 一种用于合成孔径雷达成像的导航数据拟合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103197285A CN103197285A (zh) | 2013-07-10 |
CN103197285B true CN103197285B (zh) | 2014-08-13 |
Family
ID=48719962
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310093771.5A Expired - Fee Related CN103197285B (zh) | 2013-03-22 | 2013-03-22 | 一种用于合成孔径雷达成像的导航数据拟合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103197285B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103744057B (zh) * | 2013-12-24 | 2016-06-01 | 河海大学 | 基于输出相关自适应卡尔曼滤波的弹道轨迹形成方法 |
CN104459660B (zh) * | 2014-06-19 | 2017-02-15 | 西安电子科技大学 | 基于数据拟合的阵元多普勒域空时二维自适应处理方法 |
CN113567939B (zh) * | 2021-05-17 | 2024-05-10 | 陕西航天技术应用研究院有限公司 | 一种合成孔径雷达***的预失真补偿方法、模块及收发*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6212132B1 (en) * | 1998-08-04 | 2001-04-03 | Japan Radio Co., Ltd. | Three-dimensional radar apparatus and method for displaying three-dimensional radar image |
CN101329400A (zh) * | 2008-07-30 | 2008-12-24 | 电子科技大学 | 一种基于拟合优度检验的雷达目标恒虚警检测方法 |
CN101776755A (zh) * | 2010-02-05 | 2010-07-14 | 西安电子科技大学 | 高速平台超高分辨率sar数据预处理方法 |
CN102193086A (zh) * | 2010-02-09 | 2011-09-21 | 中国科学院电子学研究所 | 对合成孔径雷达***一体化数据处理装置和方法 |
CN102279386A (zh) * | 2011-05-12 | 2011-12-14 | 西安电子科技大学 | 基于fpga的sar成像信号处理数据转置方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101177050B1 (ko) * | 2005-05-09 | 2012-08-24 | 엘타 시스템즈 리미티드 | 감소된 탐색 시간을 갖는 페이즈드 어레이 레이더 안테나및 이의 사용방법 |
-
2013
- 2013-03-22 CN CN201310093771.5A patent/CN103197285B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6212132B1 (en) * | 1998-08-04 | 2001-04-03 | Japan Radio Co., Ltd. | Three-dimensional radar apparatus and method for displaying three-dimensional radar image |
CN101329400A (zh) * | 2008-07-30 | 2008-12-24 | 电子科技大学 | 一种基于拟合优度检验的雷达目标恒虚警检测方法 |
CN101776755A (zh) * | 2010-02-05 | 2010-07-14 | 西安电子科技大学 | 高速平台超高分辨率sar数据预处理方法 |
CN102193086A (zh) * | 2010-02-09 | 2011-09-21 | 中国科学院电子学研究所 | 对合成孔径雷达***一体化数据处理装置和方法 |
CN102279386A (zh) * | 2011-05-12 | 2011-12-14 | 西安电子科技大学 | 基于fpga的sar成像信号处理数据转置方法 |
Non-Patent Citations (8)
Title |
---|
利用非线性拟合抑制前视探地雷达杂波;樊勇等;《现代雷达》;20081231;第30卷(第12期);43-48 * |
基于压缩传感的线阵SAR三维成像方法研究;韦顺军等;《宇航学报》;20111130;第32卷(第11期);2403-2409 * |
提高结构光三维形貌测量精度的分段拟合滤波方法;许中奇等;《天津理工大学学报》;20101231;第26卷(第6期);45-50 * |
樊勇等.利用非线性拟合抑制前视探地雷达杂波.《现代雷达》.2008,第30卷(第12期),43-48. |
许中奇等.提高结构光三维形貌测量精度的分段拟合滤波方法.《天津理工大学学报》.2010,第26卷(第6期),45-50. |
非移变模式记载双基地SAR成像算法及实验结果;龚镇强等;《电子与信息学报》;20090131;第31卷(第1期);173-176 * |
韦顺军等.基于压缩传感的线阵SAR三维成像方法研究.《宇航学报》.2011,第32卷(第11期),2403-2409. |
龚镇强等.非移变模式记载双基地SAR成像算法及实验结果.《电子与信息学报》.2009,第31卷(第1期),173-176. |
Also Published As
Publication number | Publication date |
---|---|
CN103197285A (zh) | 2013-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106803271A (zh) | 一种视觉导航无人机的摄像机标定方法及装置 | |
CN106767780B (zh) | 基于Chebyshev多项式插值逼近的扩展椭球集员滤波方法 | |
CN114018274B (zh) | 车辆定位方法、装置及电子设备 | |
CN103344260B (zh) | 基于rbckf的捷联惯导***大方位失准角初始对准方法 | |
CN107765244B (zh) | 基于机载双天线InSAR基线测量方法和装置 | |
CN105136145A (zh) | 一种基于卡尔曼滤波的四旋翼无人机姿态数据融合的方法 | |
CN103699798B (zh) | 一种实现地震波场数值模拟方法 | |
CN102121847A (zh) | 一种瞬态声场重建方法 | |
CN103197285B (zh) | 一种用于合成孔径雷达成像的导航数据拟合方法 | |
CN104597488B (zh) | 非等边长网格波动方程有限差分模板优化设计方法 | |
CN103115625A (zh) | 一种浮体横纵荡及升沉运动的测量方法及*** | |
CN108508463B (zh) | 基于Fourier-Hermite正交多项式扩展椭球集员滤波方法 | |
CN113218391A (zh) | 一种基于ewt算法的姿态解算方法 | |
CN102183761B (zh) | 星载干涉合成孔径雷达数字高程模型重建方法 | |
CN107741223B (zh) | 磁罗盘校准方法、装置及无人机 | |
CN111578931B (zh) | 基于在线滚动时域估计的高动态飞行器自主姿态估计方法 | |
CN110672127B (zh) | 阵列式mems磁传感器实时标定方法 | |
Dai et al. | Geomagnetic field aided inertial navigation using the SITAN algorithm | |
CN106709150B (zh) | 一种电流分布和近场电磁分布三维可视化方法 | |
CN110722798B (zh) | 一种基于组合平面切片的轨迹规划方法及产品 | |
CN104614751B (zh) | 基于约束信息的目标定位方法 | |
CN112987054B (zh) | 标定sins/dvl组合导航***误差的方法和装置 | |
CN115270579A (zh) | 二阶声波方程有限差分数值模拟参数选取方法 | |
CN111435543B (zh) | 电子地图中道路的绘制方法、装置及绘制设备、存储介质 | |
CN102800100A (zh) | 基于距离势场和自适应气球力的图像分割方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140813 |