CN108088358B - 一种基于多基线雷达轨道形变检测方法 - Google Patents
一种基于多基线雷达轨道形变检测方法 Download PDFInfo
- Publication number
- CN108088358B CN108088358B CN201711365760.2A CN201711365760A CN108088358B CN 108088358 B CN108088358 B CN 108088358B CN 201711365760 A CN201711365760 A CN 201711365760A CN 108088358 B CN108088358 B CN 108088358B
- Authority
- CN
- China
- Prior art keywords
- antenna
- scaling point
- signal
- swept
- linear
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
- B61K9/00—Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
- B61K9/08—Measuring installations for surveying permanent way
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种基于多基线雷达轨道形变检测方法,应用于雷达检测领域,本发明采用线性调频信号与扫频信号联合测距,将脉冲压缩技术与跳频技术相结合,通过对定标点的距离测量,实现了对轨道形变量的测量,并且采用线性调频信号与扫频信号联合测距的方式,很好地减小了测距误差影响。
Description
技术领域
本发明属于雷达检测领域,特别涉及一种轨道形变测量技术。
背景技术
随着中国经济的发展,轨道交通作为重要的运输手段之一,对我国国民经济发展具有重大价值。在列车行驶过程中,其作用在基面上的动荷载会引起地基变化,在长时间的过程中可能会导致轨道的沉降。如果高速铁路发生沉降,将会导致路基和线路形变等灾难性后果,并可能引发桥头跳车、路面沉陷等现象,轻则影响列车正常行驶,重则导致事故发生。由于轨道沉降对列车安全影响大,危害强,目前其已经成为轨道交通领域的研究热点。根据本发明人了解以及已发表的文献,例如:王雷霆:“论精密水准测量在地面沉降监测中的应用[J]”华北国土资源,2013(2):113-116和侯林山,王金龙,朱三妹,等:“利用差分GPS进行地面沉降监测的研究[J]”岩土力学,2006,27(5):811-815传统的沉降监测方法包括水准测量和GPS等。对于水准测量,尽管其测量误差小,但是测量需要大量的人力、物力和资源,易受天气等因素的干扰。GPS测量可对长时间发生沉降的区域进行测量,但是其精度差于水准测量。
在现代地形形变检测方法中,微波遥感技术是其中一种重要的测量手段。根据本发明人了解以及已发表的文献,例如,Amelung F,Galloway D L,Bell J W,et al.Sensingthe ups and downs of Las Vegas:“InSAR reveals structural control of landsubsidence and aquifer-system deformation[J]”Geology,1999,27(6):483-486。合成孔径雷达(SAR)测量具有全天候、全时段的特点,利用干涉SAR测量技术可获得地面微小形变变化,但该技术需要星载或者机载平台,并且数据处理复杂、计算量大。
发明内容
为解决上述技术问题,本发明提出一种基于多基线雷达轨道形变检测方法,利用布设于列车上的多基线雷达回波相位信息,通过求解定标点方程得到雷达的三维坐标信息,并获取轨道的形变量。
本发明采用的技术方案为:一种基于多基线雷达轨道形变检测方法,包括:
S1、接收回波信号,得到经各天线对定标点采样处理后的线性调频信号回波复向量以及扫频信号回波复向量;
S2、根据对应天线的线性调频信号回波复向量,经脉冲压缩处理,得到线性调频信号的定标点到该天线的延时;
S3、根据对应天线的扫频信号回波复向量,经快速傅里叶变换处理,得到扫频信号的定标点到该天线的模糊延时;
S4、根据对应天线的线性调频信号回波的相位,计算得到线性调频信号的定标点到该天线的相对距离延时;
S5、根据步骤S3经快速傅里叶变换处理后的扫频信号的相位,计算得到扫频信号的定标点到该天线的相对距离延时;
S6、重复步骤S2至步骤S5,分别得到线性调频信号的定标点到各天线的延时历史,以及扫频信号的定标点到各天线的模糊延时历史;
S7、根据线性调频信号的定标点到天线的延时历史与扫频信号的定标点到天线的模糊延时历史,计算得到定标点到天线的精确延时历史;
S8、对步骤S4得到的线性调频信号的定标点到天线的相对距离延时进行相位解缠,得到线性调频信号的定标点到天线的相对距离延时历史;对步骤S5得到的扫频信号的定标点到该天线的相对距离延时进行相位解缠,得到扫频信号的定标点到该天线的相对距离延时历史;
S9、根据定标点到天线的精确延时历史、线性调频信号的定标点到天线的相对距离延时历史以及扫频信号的定标点到天线的相对距离延时历史;得到定标点到对应天线校正后的精确距离延时历史;
S10、重复步骤S2至步骤S9,得到二次测量定标点到各天线校正后的精确距离延时历史;
S11、根据步骤S9得到的定标点到对应天线校正后的精确距离延时历史,以及步骤S10得到的二次测量定标点到各天线校正后的精确距离延时历史;计算局部坐标系形变量历史;
S12、将步骤S11得到的各局部坐标系形变量历史转换至东北天坐标系,得到轨道在东北天坐标系下的形变量。
进一步地,所述步骤S1之前还包括:S0、发射信号,具体为:根据已知的线性调频信号带宽,线性调频信号时宽和雷达载波频率,发射线性调频信号;
在线性调频信号发射完毕,延时线性调频信号和扫频信号时间差,发射扫频信号。
进一步地,所述天线包括:布设于测量列车原点处的基准天线;布设于测量列车原点且垂直于测量列车基准面,并到测量列车原点距离为天线基线长度b位置处的高度向天线;布设于测量列车原点且平行于测量列车基准面中x轴的方向,并到测量列车原点距离为天线基线长度b位置处的视向天线;布设于测量列车原点且平行于测量列车基准面中y轴的方向,并到测量列车原点距离为天线基线长度b位置处的运动向天线。
更进一步地,所述天线至少包括:基准天线、高度向天线、视向天线、运动向天线中的两种。
进一步地,步骤S3所述计算扫频信号的定标点到该天线的模糊延时,具体为:根据步骤得到的对应天线的扫频信号回波相连,进行k×N点的快速傅里叶变换,得到压缩后该天线扫频信号回波复向量,选择压缩后该天线扫频信号回波复向量的幅度最大值对应的位置,根据以下公式,计算得到定标点到该天线的模糊延时θ:
其中,i为压缩后该天线扫频信号回波复向量的幅度最大值对应的位置,k为自然数,N为基准天线的扫频信号回波复向量的维度,Δf为跳频间隔。
进一步地,步骤S4所述计算线性调频信号的定标点到该天线的相对距离延时,具体为:对该天线的线性调频信号回波复向量,采用脉冲压缩技术,得到压缩后该天线的线性调频信号回波复向量,寻找压缩后该天线的线性调频信号回波复向量的幅度最大的复数值,根据以下公式,计算得到线性调频信号的定标点到该天线的相对距离延时
其中,angle[·]为取复数的角度,fc为雷达载波频率,ρ为压缩后该天线的线性调频信号回波复向量的幅度最大的复数值。
进一步地,步骤S5所述计算扫频信号的定标点到该天线的相对距离延时,具体为:根据压缩后该天线的扫频信号回波复向量,寻找压缩后该天线的扫频信号回波复向量的幅度最大的复数值,根据以下公式,计算得到扫频信号的定标点到基准天线的相对距离延时φ;
φ=angle[μ]/(2πfc)
其中,angle[·]为取复数的角度,fc雷达载波频率,μ为压缩后该天线的扫频信号回波复向量的幅度最大的复数值。
进一步地,步骤S7所述计算定标点到天线的精确延时历史,具体为:根据线性调频信号的定标点到各天线的延时历史,以及扫频信号的定标点到各天线的模糊延时历史,结合下式,得到定标点到基准天线的精确延时历史o:
其中,floor[·]为向下取整,Δf为跳频间隔,θ为定标点到该天线的模糊延时。
进一步地,步骤S9所述计算定标点到对应天线校正后的精确距离延时历史,具体为:根据定标点到该天线的精确延时历史,解缠后的线性调频信号的定标点到基准天线的相对距离延时历史和解缠后的扫频信号的定标点到基准天线的相对距离延时历史,根据下式,计算得到定标点到基准天线校正后的精确距离延时历史:
其中,α为精确距离延时权系数,β为相对距离延时权系数,o(n)为定标点到基准天线的精确延时历史,ε(n)为解缠后的线性调频信号的定标点到基准天线的相对距离延时历史,η(n)为解缠后的扫频信号的定标点到高度向天线的相对距离延时历史,为定标点到基准天线校正后的精确距离延时历史。
本发明的有益效果:本发明的一种基于多基线雷达轨道形变检测方法,采用线性调频信号与扫频信号联合测距,将脉冲压缩技术与跳频技术相结合,通过对定标点的距离测量,实现了对轨道形变量的测量,并且采用线性调频信号与扫频信号联合测距的方式,很好地减小了测距误差影响。
附图说明
图1多基线雷达轨道形变检测系统结构示意图;
图2为发射机系统框图;
图3为接收机系统框图;
图4为多基线雷达天线示意图;
其中,1表示基准天线A,2表示高度向天线B,3表示视线向天线C,4表示运动向天线D,5表示测量列车基准面;
图5为本发明的方案流程图;
图6为线性调频信号与扫频信号时频图;
其中,1表示线性调频信号,2、3、4表示扫频信号各频率分量;
图7为本发明具体实施方式中采用形变量数据;
其中,1表示x轴真实形变量,2表示z轴真实形变量,3表示y轴真实形变量;
图8为本发明所提出方法得到的测量形变量数据;
其中,1表示x轴测量形变量,2表示z轴测量形变量,3表示y轴测量形变量。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
为了方便描述本发明的内容,首先作以下术语定义:
定义1、东北天坐标系
东北天坐标系也叫做站心坐标系,以站心(如GPS接收天线)为坐标系原点,z轴与椭圆法线重合,向上为正(天向),y轴与椭圆短半轴重合(北向),x轴与地球椭圆的长半轴重合(东向)所构成的直角坐标系。详细内容可参考“地理信息系统导论(第3版)”Kang-tsungChang,张康聪,and陈健飞.清华大学出版社,2009.
定义2、定标点
定标点是用来接收发射信号并将发射信号转发回天线的设备。定标点主要分为无源和有源两种。无源定标点,如定标球、角反射器等,具有成本低、便于维护等优点;有源定标点,如应答器,接收雷达信号并直接转发回雷达,能更好保持信号信噪比,并能提供自身位置等辅助信息,更利于后续数据处理。详细内容可参考“精密跟踪测量雷达技术[M]”王德纯,丁家会,程望东.电子工业出版社,2006.
定义3、测量列车
测量列车即用于铁路高速综合检测列车。测量列车基准面指测量列车上与轨道平面平行的雷达安装平面,本发明中定义测量列车速度方向为y轴,垂直于速度方向为x轴。详细内容参考“高速铁路工程测量规范[M]”中铁二院工程集团有限责任公司.中国铁道出版社,2010。
如图1所示为本发明实施例的多基线雷达轨道形变检测系统结构示意图,本发明的检测系统包括:发射机、接收机、天线以及定标点等,多基线雷达布设在测量列车上,随测量列车在轨道上运动,定标点位于轨道两侧。发射机系统框图如图2所示,依次包括:基带信号模块、载频调制模块、功率放大器模块以及发射天线;接收机系统框图如图3所示,依次包括:接收天线、下调制模块、前置放大器模块、数模转换器模块以及数据存储器。
本实施例的多基线雷达轨道形变检测系统初始化参数,包括:定标点在东北天坐标系下的位置Pw=(0,0,0)T;天线基线长度b=1m;雷达载波频率fc=10GHZ;跳频点数量M=1000;采样频率fs=2GHZ;线性调频信号带宽B=50MHZ;线性调频信号时宽,记作T1=10μs;扫频信号每个频率点的驻留时间,记作T2=10μs;跳频间隔Δf=50MHZ;采用公式k=50×106/10×10-6计算出雷达调频斜率k=5×1012(HZ)2;线性调频信号和扫频信号时间差T3=20μs。精确距离延时权系数α=0.5;相对距离延时权系数β=0.5;光传播速度c=3×108m/s。
如图4所示,本实施例布设4部天线。首先,在测量列车原点处布设一部天线,称为基准天线,记作A;然后,在经过测量列车原点且垂直于测量列车基准面,并到测量列车原点距离为天线基线长度b的位置布设一部天线,称为高度向天线,记作B;其次,在经过测量列车原点且平行于测量列车基准面中x轴的方向,并到测量列车原点距离为天线基线长度b的位置布设一部天线,称为视线向天线,记作C;最后,在经过测量列车原点且平行于测量列车基准面中y轴的方向,并到测量列车原点距离为天线基线长度b的位置布设一部天线,称为运动向天线,记作D。
工作原理为:发射机发射线性调频信号与扫频信号;接收机接收到的线性调频信号回波以及扫频信号回波;并对接收到的信号进行一系列的处理;具体处理过程如图5所示。线性调频信号与扫频信号如图6所示。
首先根据已知的线性调频信号带宽B=50MHZ,线性调频信号时宽T1=10μs和雷达载波频率fc=10GHZ,发射线性调频信号,记作F1;
其次,根据已知的线性调频信号和扫频信号时间差T3=20μs,跳频点数量M=1000,跳频间隔Δf=50MHZ,扫频信号每个频率点的驻留时间T2=10μs和雷达载波频率fc=10GHZ,在线性调频信号发射完毕,并延时线性调频信号和扫频信号时间差T3=20μs,发射扫频信号,记作F2。
本发明的技术方案为:一种基于多基线雷达轨道形变检测方法,包括:
S1、接收回波信号,得到经各天线对定标点采样处理后的线性调频信号回波复向量以及扫频信号回波复向量;具体的:
根据已知的线性调频信号时宽T1=10μs,采样频率fs=2GHZ和雷达载波频率fc=10GHZ,经过基准天线对定标点的反射回波进行采样,得到基准天线的线性调频信号回波复向量,记作s0;经过高度向天线对定标点的反射回波进行采样,得到高度向天线的线性调频信号回波复向量,记作s1;经过视线向天线对定标点的反射回波进行采样,得到视线向天线的线性调频信号回波复向量,记作s2;经过运动向天线对定标点的反射回波进行采样,得到运动向天线的线性调频信号回波复向量,记作s3。
根据已知的线性调频信号和扫频信号时间差T3=20μs,跳频点数量M=1000,跳频间隔Δf=50MHZ,扫频信号每个频率点的驻留时间T2=10μs和雷达载波频率fc=10GHZ,采样频率fs=2GHZ和雷达载波频率fc=10GHZ,经过基准天线对定标点的反射回波进行采样,得到基准天线的扫频信号回波复向量,记作w0;经过高度向天线对定标点的反射回波进行采样,得到高度向天线的扫频信号回波复向量,记作w1;经过视线向天线对定标点的反射回波进行采样,得到视线向天线的扫频信号回波复向量,记作w2;经过运动向天线对定标点的反射回波进行采样,得到运动向天线的扫频信号回波复向量,记作w3。
S2、根据对应天线的线性调频信号回波复向量,经脉冲压缩处理,得到线性调频信号的定标点到该天线的延时;具体的:
根据得到的基准天线的线性调频信号回波复向量s0,采用脉冲压缩技术,得到定标点到基准天线的延时,记作τ0;
根据得到的高度向天线的线性调频信号回波复向量s1,采用脉冲压缩技术,得到定标点到高度向天线的延时,记作τ1;
根据得到的视线向天线的线性调频信号回波复向量s2,采用脉冲压缩技术,得到定标点到视线向天线的延时,记作τ2;
根据得到的运动向天线的线性调频信号回波复向量s3,采用脉冲压缩技术,得到定标点到运动向天线的延时,记作τ3。
S3、根据对应天线的扫频信号回波复向量,经快速傅里叶变换处理,得到扫频信号的定标点到该天线的模糊延时;具体的:
对得到的基准天线的扫频信号回波向量w0,进行k×N点的快速傅里叶变换,其中,k=50,N=1000,得到压缩后基准天线扫频信号回波复向量,记作v0,选择压缩后基准天线扫频信号回波复向量v0的幅度最大值对应的位置,记作i0,利用公式得到定标点到基准天线的模糊延时,记作θ0;
对得到的高度向天线的扫频信号回波复向量w1,进行k×N点的快速傅里叶变换,得到压缩后高度向天线扫频信号回波复向量,记作v1,选择压缩后高度向天线扫频信号回波复向量v1的幅度最大值对应的位置,记作i1,利用公式得到定标点到高度向天线的模糊延时,记作θ1;
对得到的视线向天线的扫频信号回波复向量w2,进行k×N点的快速傅里叶变换,得到压缩后的视线向天线的扫频信号回波复向量,记作v2,选择压缩后视线向天线扫频信号回波复向量v2的幅度最大值对应的位置,记作i2,利用公式得到定标点到视线向天线的模糊延时,记作θ2;
对得到的运动向天线的扫频信号回波复向量w3,进行k×N点的快速傅里叶变换,得到压缩后运动向天线扫频信号回波复向量,记作v3,选择压缩后运动向天线扫频信号回波复向量v3的幅度最大值对应的位置,记作i3,利用公式得到定标点到运动向天线的模糊延时,记作θ3;
S4、根据对应天线的线性调频信号回波的相位,计算得到线性调频信号的定标点到该天线的相对距离延时;具体的:
根据得到的基准天线的线性调频信号回波复向量s0,采用脉冲压缩技术,得到压缩后基准天线的线性调频信号回波复向量,记作p0,寻找压缩后基准天线的线性调频信号回波复向量p0的幅度最大的复数值,记作ρ0,利用公式其中,angle[·]为取复数的角度,得到线性调频信号的定标点到基准天线的相对距离延时,记作
根据得到的高度向天线的线性调频信号回波复向量s1,采用脉冲压缩技术,得到压缩后高度向天线的线性调频信号回波复向量,记作p1,寻找压缩后高度向天线的线性调频信号回波复向量p1的幅度最大的复数值,记作ρ1,利用公式得到线性调频信号的定标点到高度向天线的相对距离延时,记作
根据得到的视线向天线的线性调频信号回波复向量s2,采用脉冲压缩技术,得到压缩后视线向天线的线性调频信号回波复向量,记作p2,寻找压缩后视线向天线的线性调频信号回波复向量p2的幅度最大的复数值,记作ρ2,利用公式得到线性调频信号的定标点到视线向天线的相对距离延时,记作
根据得到的运动向天线的线性调频信号回波复向量s3,采用脉冲压缩技术,得到压缩后运动向天线的线性调频信号回波复向量,记作p3,寻找压缩后运动向天线的线性调频信号回波复向量p3的幅度最大的复数值,记作ρ3,利用公式得到线性调频信号的定标点到运动向天线的相对距离延时,记作
S5、根据步骤S3经快速傅里叶变换处理后的扫频信号的相位,计算得到扫频信号的定标点到该天线的相对距离延时;具体的:
根据得到压缩后基准天线的扫频信号回波复向量v0,寻找压缩后基准天线的扫频信号回波复向量v0的幅度最大的复数值,记作μ0,利用公式φ=angle[μ]/(2π×10×109),得到扫频信号的定标点到基准天线的相对距离延时,记作φ0;
根据得到压缩后高度向天线的扫频信号回波复向量v1,寻找压缩后高度向天线的扫频信号回波复向量v1的幅度最大的复数值,记作μ1,利用公式φ=angle[μ]/(2π×10×109),得到扫频信号的定标点到高度向天线的相对距离延时,记作φ1;
根据得到压缩后视线向天线的扫频信号回波复向量v2,寻找压缩后视线向天线的扫频信号回波复向量v2的幅度最大的复数值,记作μ2,利用公式φ=angle[μ]/(2π×10×109),得到扫频信号的定标点到视线向天线的相对距离延时,记作φ2;
根据得到压缩后运动向天线的扫频信号回波复向量v3,寻找压缩后运动向天线的扫频信号回波复向量v3的幅度最大的复数值,记作μ3,利用公式φ=angle[μ]/(2π×10×109),得到扫频信号的定标点到运动向天线的相对距离延时,记作φ3;
S6、重复步骤S2至步骤S5,分别得到线性调频信号的定标点到各天线的延时历史,以及扫频信号的定标点到各天线的模糊延时历史;具体的:
根据测量列车运动规律,等距离间隔地发射线性调频信号及扫频信号,并接收回波信号、利用线性调频信号计算距离延时、利用扫频信号计算模糊距离延时、利用信号相位计算相对距离延时,得到定标点到基准天线的延时历史,记作τ0(n);定标点到高度向天线的延时历史,记作τ1(n);定标点到视线向天线的延时历史,记作τ2(n);定标点到运动向天线的延时历史,记作τ3(n);定标点到基准天线的模糊延时历史,记作θ0(n);定标点到高度向天线的模糊延时历史,记作θ1(n);定标点到视线向天线的模糊延时历史,记作θ2(n);定标点到运动向天线的模糊延时历史,记作θ3(n);线性调频信号的定标点到基准天线的相对距离延时历史,记作线性调频信号的定标点到高度向天线的相对距离延时历史,记作线性调频信号的定标点到视线向天线的相对距离延时历史,记作线性调频信号的定标点到运动向天线的相对距离延时历史,记作扫频信号的定标点到基准天线的相对距离延时历史,记作φ0(n);扫频信号的定标点到高度向天线的相对距离延时历史,记作φ1(n);扫频信号的定标点到视线向天线的相对距离延时历史,记作φ2(n);扫频信号的定标点到运动向天线的相对距离延时历史,记作φ3(n),其中n为发射信号序号,n=0,1,2,…。
S7、根据线性调频信号的定标点到天线的延时历史与扫频信号的定标点到天线的模糊延时历史,计算得到定标点到天线的精确延时历史;具体的:
利用定标点到基准天线的延时历史τ0(n)和定标点到基准天线的模糊延时历史θ0(n),根据公式其中floor[·]为向下取整,得到定标点到基准天线的精确延时历史,记作o0(n);利用定标点到高度向天线的延时历史τ1(n)和定标点到高度向天线的模糊延时历史θ1(n),根据公式得到定标点到高度向天线的精确延时历史,记作o1(n);利用定标点到视线向天线的延时历史τ2(n)和定标点到视线向天线的模糊延时历史θ2(n),根据公式得到定标点到视线向天线的精确延时历史,记作o2(n);利用定标点到运动向天线的延时历史τ3(n)和定标点到运动向天线的模糊延时历史θ3(n),根据公式得到定标点到运动向天线的精确延时历史,记作o3(n)。
S8、对步骤S4得到的线性调频信号的定标点到天线的相对距离延时进行相位解缠,得到线性调频信号的定标点到天线的相对距离延时历史;对步骤S5得到的扫频信号的定标点到该天线的相对距离延时进行相位解缠,得到扫频信号的定标点到该天线的相对距离延时历史;具体的:
对步骤S4得到的线性调频信号的定标点到基准天线的相对距离延时历史进行相位解缠,得到解缠后的线性调频信号的定标点到基准天线的相对距离延时历史,记作ε0(n);
对步骤S4得到的线性调频信号的定标点到高度向天线的相对距离延时历史进行相位解缠,得到解缠后的线性调频信号的定标点到高度向天线的相对距离延时历史,记作ε1(n);
对步骤S4得到的线性调频信号的定标点到视线向天线的相对距离延时历史进行相位解缠,得到解缠后的线性调频信号的定标点到视线向天线的相对距离延时历史,记作ε2(n);
对步骤S4得到的线性调频信号的定标点到运动向天线的相对距离延时历史进行相位解缠,得到解缠后的线性调频信号的定标点到运动向天线的相对距离延时历史,记作ε3(n);
对步骤S5得到的扫频信号的定标点到基准天线的相对距离延时历史φ0(n)进行相位解缠,得到解缠后的扫频信号的定标点到基准天线的相对距离延时历史,记作η0(n);
对步骤S5得到的扫频信号的定标点到高度向天线的相对距离延时历史φ1(n)进行相位解缠,得到解缠后的扫频信号的定标点到高度向天线的相对距离延时历史,记作η1(n);
对步骤S5得到的扫频信号的定标点到视线向天线的相对距离延时历史φ2(n)进行相位解缠,得到解缠后的扫频信号的定标点到视线向天线的相对距离延时历史,记作η2(n);
对步骤S5得到的扫频信号的定标点到运动向天线的相对距离延时历史φ3(n)进行相位解缠,得到解缠后的扫频信号的定标点到运动向天线的相对距离延时历史,记作η3(n);
S9、根据定标点到天线的精确延时历史、线性调频信号的定标点到天线的相对距离延时历史以及扫频信号的定标点到天线的相对距离延时历史;得到定标点到对应天线校正后的精确距离延时历史;具体的:
利用定标点到基准天线的精确延时历史o0(n),解缠后的线性调频信号的定标点到基准天线的相对距离延时历史ε0(n)和解缠后的扫频信号的定标点到基准天线的相对距离延时历史η0(n),根据公式得到定标点到基准天线校正后的精确距离延时历史,记作
利用定标点到高度向天线的精确延时历史o1(n),解缠后的线性调频信号的定标点到高度向天线的相对距离延时历史ε1(n)和解缠后的扫频信号的定标点到高度向天线的相对距离延时历史η1(n),根据公式得到定标点到高度向天线校正后的精确距离延时历史,记作
利用定标点到视线向天线的精确延时历史o2(n),解缠后的线性调频信号的定标点到视线向天线的相对距离延时历史ε2(n)和解缠后的扫频信号的定标点到视线向天线的相对距离延时历史η2(n),根据公式得到定标点到视线向天线校正后的精确距离延时历史,记作
利用定标点到运动向天线的精确延时历史o3(n),解缠后的线性调频信号的定标点到运动向天线的相对距离延时历史ε3(n)和解缠后的扫频信号的定标点到运动向天线的相对距离延时历史η3(n),根据公式得到定标点到运动向天线校正后的精确距离延时历史,记作
S10、重复步骤S2至步骤S9,得到二次测量定标点到各天线校正后的精确距离延时历史;具体的:重复步骤S2至步骤S9,进行二次测量,得到二次测量定标点到基准天线校正后的精确距离延时历史,记作二次测量定标点到高度向天线校正后的精确距离延时历史,记作二次测量定标点到视线向天线校正后的精确距离延时历史,记作二次测量定标点到运动向天线校正后的精确距离延时历史,记作
S11、根据步骤S9得到的定标点到对应天线校正后的精确距离延时历史,以及步骤S10得到的二次测量定标点到各天线校正后的精确距离延时历史;计算局部坐标系形变量历史;具体的:
利用已知的天线基线长度b=1米,光传播速度c=3×108米每秒,已知二次测量定标点到高度向天线校正后的精确距离延时历史二次测量定标点到基线天线校正后的精确距离延时历史定标点到高度向天线校正后的精确距离延时历史和定标点到基准天线校正后的精确距离延时历史根据公式得到形变量历史的x轴分量,记作u(n);
利用已知的天线基线长度b=1米,光传播速度c=3×108米每秒,已知二次观测定标点到视线向天线校正后的精确距离延时历史二次观测定标点到基线天线校正后的精确距离延时历史定标点到视线向天线校正后的精确距离延时历史和定标点到基准天线校正后的精确距离延时历史根据公式得到形变量历史的y轴分量,记作v(n);
利用已知的天线基线长度b=1米,光传播速度c=3×108米每秒,已知二次观测定标点到运动向天线校正后的精确距离延时历史二次观测定标点到基线天线校正后的精确距离延时历史定标点到运动向天线校正后的精确距离延时历史和定标点到基准天线校正后的精确距离延时历史根据公式得到形变量历史的z轴分量,记作w(n)。
形变量历史的x轴分量、形变量历史的y轴分量、形变量历史的z轴分量即构成了轨道在测量列车局部坐标系下的形变量。
S12、将步骤S11得到的各局部坐标系形变量历史转换至东北天坐标系,得到轨道在东北天坐标系下的形变量。具体的:
根据步骤12得到的形变量历史的x轴分量、形变量历史的y轴分量、形变量历史的z轴分量,并利用定义4的方法得到局部坐标系的三个坐标基(ξx、ξy、ξz)在东北天坐标系下的坐标,用公式采用公式(ENU)T=u(n)ξu(n)+v(n)ξv(n)+w(n)ξw(n),其中,ξu(n)、ξv(n)、ξw(n)为局部坐标系在东北天坐标系下的坐标,得到轨道在东北天坐标系下的形变量。
经过以上步骤,即可得到轨道在测量列车局部坐标系下的形变量和轨道在东北天坐标系下的形变量。
本实施例采用的形变量数据如图7所示,经过本发明方法得到的测量形变量数据如图8所示;从中可以看出,本发明所采用方法具有高精度测量能力,并保持了测量信号的完整性和信噪比。与干涉技术相比,本发明采用在地面布设有源或强定标点的方式,因此,信号不受点间串扰以及散射点散射系数的影响,能更好地保证测量信号的完整性和信噪比,以达到与干涉技术相似甚至更高的测量精度。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (9)
1.一种基于多基线雷达轨道形变检测方法,其特征在于,包括:
S1、接收回波信号,得到经各天线对定标点采样处理后的线性调频信号回波复向量以及扫频信号回波复向量;
S2、根据对应天线的线性调频信号回波复向量,经脉冲压缩处理,得到线性调频信号的定标点到该天线的延时;
S3、根据对应天线的扫频信号回波复向量,经快速傅里叶变换处理,得到扫频信号的定标点到该天线的模糊延时;
S4、根据对应天线的线性调频信号回波的相位,计算得到线性调频信号的定标点到该天线的相对距离延时;
S5、根据步骤S3经快速傅里叶变换处理后的扫频信号的相位,计算得到扫频信号的定标点到该天线的相对距离延时;
S6、重复步骤S2至步骤S5,分别得到线性调频信号的定标点到各天线的延时历史,以及扫频信号的定标点到各天线的模糊延时历史;
S7、根据线性调频信号的定标点到天线的延时历史与扫频信号的定标点到天线的模糊延时历史,计算得到定标点到天线的精确延时历史;
S8、对步骤S4得到的线性调频信号的定标点到天线的相对距离延时进行相位解缠,得到线性调频信号的定标点到天线的相对距离延时历史;对步骤S5得到的扫频信号的定标点到该天线的相对距离延时进行相位解缠,得到扫频信号的定标点到该天线的相对距离延时历史;
S9、根据定标点到天线的精确延时历史、线性调频信号的定标点到天线的相对距离延时历史以及扫频信号的定标点到天线的相对距离延时历史;得到定标点到对应天线校正后的精确距离延时历史;
S10、重复步骤S2至步骤S9,得到二次测量定标点到各天线校正后的精确距离延时历史;
S11、根据步骤S9得到的定标点到对应天线校正后的精确距离延时历史,以及步骤S10得到的二次测量定标点到各天线校正后的精确距离延时历史;计算局部坐标系形变量历史;
S12、将步骤S11得到的各局部坐标系形变量历史转换至东北天坐标系,得到轨道在东北天坐标系下的形变量。
2.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,所述步骤S1之前还包括:S0、发射信号,具体为:根据已知的线性调频信号带宽,线性调频信号时宽和雷达载波频率,发射线性调频信号;
在线性调频信号发射完毕,延时线性调频信号和扫频信号时间差,发射扫频信号。
3.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,所述天线包括:布设于测量列车原点处的基准天线;布设于测量列车原点且垂直于测量列车基准面,并到测量列车原点距离为天线基线长度b位置处的高度向天线;布设于测量列车原点且平行于测量列车基准面中x轴的方向,并到测量列车原点距离为天线基线长度b位置处的视向天线;布设于测量列车原点且平行于测量列车基准面中y轴的方向,并到测量列车原点距离为天线基线长度b位置处的运动向天线。
4.根据权利要求3所述的一种基于多基线雷达轨道形变检测方法,其特征在于,所述天线至少包括:基准天线、高度向天线、视向天线、运动向天线中的两种。
5.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,步骤S3所述得到扫频信号的定标点到该天线的模糊延时,具体为:根据步骤得到的对应天线的扫频信号回波相连,进行k×N点的快速傅里叶变换,得到压缩后该天线扫频信号回波复向量,选择压缩后该天线扫频信号回波复向量的幅度最大值对应的位置,根据以下公式,计算得到定标点到该天线的模糊延时θ:
其中,i为压缩后该天线扫频信号回波复向量的幅度最大值对应的位置,k为自然数,N为基准天线的扫频信号回波复向量的维度,Δf为跳频间隔。
6.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,步骤S4所述计算线性调频信号的定标点到该天线的相对距离延时,具体为:对该天线的线性调频信号回波复向量,采用脉冲压缩技术,得到压缩后该天线的线性调频信号回波复向量,寻找压缩后该天线的线性调频信号回波复向量的幅度最大的复数值,根据以下公式,计算得到线性调频信号的定标点到该天线的相对距离延时
其中,angle[·]为取复数的角度,fc为雷达载波频率,ρ为压缩后该天线的线性调频信号回波复向量的幅度最大的复数值。
7.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,步骤S5所述计算得到扫频信号的定标点到该天线的相对距离延时,具体为:根据压缩后该天线的扫频信号回波复向量,寻找压缩后该天线的扫频信号回波复向量的幅度最大的复数值,根据以下公式,计算得到扫频信号的定标点到基准天线的相对距离延时φ;
φ=angle[μ]/(2πfc)
其中,angle[·]为取复数的角度,fc雷达载波频率,μ为压缩后该天线的扫频信号回波复向量的幅度最大的复数值。
8.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,步骤S7所述计算得到定标点到天线的精确延时历史,具体为:根据线性调频信号的定标点到各天线的延时历史,以及扫频信号的定标点到各天线的模糊延时历史,结合下式,得到定标点到基准天线的精确延时历史o:
其中,floor[·]为向下取整,Δf为跳频间隔,θ为定标点到该天线的模糊延时。
9.根据权利要求1所述的一种基于多基线雷达轨道形变检测方法,其特征在于,步骤S9所述得到定标点到对应天线校正后的精确距离延时历史,具体为:根据定标点到该天线的精确延时历史,解缠后的线性调频信号的定标点到基准天线的相对距离延时历史和解缠后的扫频信号的定标点到基准天线的相对距离延时历史,根据下式,计算得到定标点到基准天线校正后的精确距离延时历史:
其中,α为精确距离延时权系数,β为相对距离延时权系数,o(n)为定标点到基准天线的精确延时历史,ε(n)为解缠后的线性调频信号的定标点到基准天线的相对距离延时历史,η(n)为解缠后的扫频信号的定标点到高度向天线的相对距离延时历史,为定标点到基准天线校正后的精确距离延时历史。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711365760.2A CN108088358B (zh) | 2017-12-18 | 2017-12-18 | 一种基于多基线雷达轨道形变检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711365760.2A CN108088358B (zh) | 2017-12-18 | 2017-12-18 | 一种基于多基线雷达轨道形变检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108088358A CN108088358A (zh) | 2018-05-29 |
CN108088358B true CN108088358B (zh) | 2019-08-20 |
Family
ID=62176994
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711365760.2A Active CN108088358B (zh) | 2017-12-18 | 2017-12-18 | 一种基于多基线雷达轨道形变检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108088358B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115447632A (zh) * | 2022-09-13 | 2022-12-09 | 中车工业研究院(青岛)有限公司 | 一种列车轨道检测系统、方法、装置以及介质 |
Citations (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09189762A (ja) * | 1996-01-08 | 1997-07-22 | Mitsubishi Electric Corp | レーダ装置を用いた地表変動観測方法並びにこの方法に用いる合成開口レーダ装置及びトランスポンダ |
US6914553B1 (en) * | 2004-11-09 | 2005-07-05 | Harris Corporation | Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the Faraday rotation, and associated methods |
JP2008232852A (ja) * | 2007-03-20 | 2008-10-02 | Yamaguchi Univ | 電磁波レーダによる鉄筋コンクリート構造物中の鉄筋径の非破壊計測装置及び非破壊計測方法 |
CN101706577A (zh) * | 2009-12-01 | 2010-05-12 | 中南大学 | InSAR监测高速公路路面沉降方法 |
EP2208016A1 (fr) * | 2007-11-09 | 2010-07-21 | Thales | Capteur optique pour mesurer la deformation au cours du temps d'une structure plane deformable |
CN102331575A (zh) * | 2011-06-23 | 2012-01-25 | 中国电子科技集团公司第三十八研究所 | 一种单航过全极化干涉合成孔径雷达 |
CN102608584A (zh) * | 2012-03-19 | 2012-07-25 | 中国测绘科学研究院 | 基于多项式反演模型的时间序列InSAR形变监测方法及装置 |
CN103048656A (zh) * | 2012-12-25 | 2013-04-17 | 河南东陆高科实业有限公司 | 基于连续波雷达的隧道围岩形变量测预警系统及方法 |
CN103091676A (zh) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | 矿区地表开采沉陷合成孔径雷达干涉测量的监测及解算方法 |
CN103308911A (zh) * | 2013-06-06 | 2013-09-18 | 重庆大学 | 基于距离欺骗技术的微位移测量方法及系统 |
CN103675790A (zh) * | 2013-12-23 | 2014-03-26 | 中国国土资源航空物探遥感中心 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
CN103792531A (zh) * | 2014-02-21 | 2014-05-14 | 重庆大学 | 基于天线阵的多目标微位移测量方法 |
CN104111456A (zh) * | 2014-07-23 | 2014-10-22 | 中国国土资源航空物探遥感中心 | 一种高速铁路沿线地表形变高分辨率InSAR监测方法 |
CN104122553A (zh) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | 一种集成多轨道、长条带CTInSAR的区域性地面沉降监测方法 |
CN104849712A (zh) * | 2015-04-22 | 2015-08-19 | 北京理工大学 | 一种基于多基地mimo-sar的三维形变监测系统 |
CN105043284A (zh) * | 2015-08-06 | 2015-11-11 | 中国电子科技集团公司第三十八研究所 | 一种大型雷达天线变形实时测量系统及方法 |
CN204925385U (zh) * | 2015-08-12 | 2015-12-30 | 桂林电子科技大学 | 一种山体滑坡的mimo雷达监测系统 |
CN106707281A (zh) * | 2017-01-05 | 2017-05-24 | 北京航空航天大学 | 一种基于多频数据处理的机载D‑InSAR形变检测方法 |
CN107389029A (zh) * | 2017-08-24 | 2017-11-24 | 北京市水文地质工程地质大队 | 一种基于多源监测技术融合的地面沉降集成监测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8587493B2 (en) * | 2010-09-23 | 2013-11-19 | North Carolina State University | Reversibly deformable and mechanically tunable fluidic antennas |
-
2017
- 2017-12-18 CN CN201711365760.2A patent/CN108088358B/zh active Active
Patent Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09189762A (ja) * | 1996-01-08 | 1997-07-22 | Mitsubishi Electric Corp | レーダ装置を用いた地表変動観測方法並びにこの方法に用いる合成開口レーダ装置及びトランスポンダ |
US6914553B1 (en) * | 2004-11-09 | 2005-07-05 | Harris Corporation | Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the Faraday rotation, and associated methods |
JP2008232852A (ja) * | 2007-03-20 | 2008-10-02 | Yamaguchi Univ | 電磁波レーダによる鉄筋コンクリート構造物中の鉄筋径の非破壊計測装置及び非破壊計測方法 |
JP4815607B2 (ja) * | 2007-03-20 | 2011-11-16 | 国立大学法人山口大学 | 電磁波レーダによる鉄筋コンクリート構造物中の鉄筋径の非破壊計測装置及び非破壊計測方法 |
EP2208016A1 (fr) * | 2007-11-09 | 2010-07-21 | Thales | Capteur optique pour mesurer la deformation au cours du temps d'une structure plane deformable |
CN101706577A (zh) * | 2009-12-01 | 2010-05-12 | 中南大学 | InSAR监测高速公路路面沉降方法 |
CN102331575A (zh) * | 2011-06-23 | 2012-01-25 | 中国电子科技集团公司第三十八研究所 | 一种单航过全极化干涉合成孔径雷达 |
CN102608584A (zh) * | 2012-03-19 | 2012-07-25 | 中国测绘科学研究院 | 基于多项式反演模型的时间序列InSAR形变监测方法及装置 |
CN103048656A (zh) * | 2012-12-25 | 2013-04-17 | 河南东陆高科实业有限公司 | 基于连续波雷达的隧道围岩形变量测预警系统及方法 |
CN103091676A (zh) * | 2013-01-22 | 2013-05-08 | 中国矿业大学 | 矿区地表开采沉陷合成孔径雷达干涉测量的监测及解算方法 |
CN103308911A (zh) * | 2013-06-06 | 2013-09-18 | 重庆大学 | 基于距离欺骗技术的微位移测量方法及系统 |
CN103675790A (zh) * | 2013-12-23 | 2014-03-26 | 中国国土资源航空物探遥感中心 | 一种基于高精度DEM提高InSAR技术监测地表形变精度的方法 |
CN103792531A (zh) * | 2014-02-21 | 2014-05-14 | 重庆大学 | 基于天线阵的多目标微位移测量方法 |
CN104111456A (zh) * | 2014-07-23 | 2014-10-22 | 中国国土资源航空物探遥感中心 | 一种高速铁路沿线地表形变高分辨率InSAR监测方法 |
CN104122553A (zh) * | 2014-07-23 | 2014-10-29 | 中国国土资源航空物探遥感中心 | 一种集成多轨道、长条带CTInSAR的区域性地面沉降监测方法 |
CN104849712A (zh) * | 2015-04-22 | 2015-08-19 | 北京理工大学 | 一种基于多基地mimo-sar的三维形变监测系统 |
CN105043284A (zh) * | 2015-08-06 | 2015-11-11 | 中国电子科技集团公司第三十八研究所 | 一种大型雷达天线变形实时测量系统及方法 |
CN204925385U (zh) * | 2015-08-12 | 2015-12-30 | 桂林电子科技大学 | 一种山体滑坡的mimo雷达监测系统 |
CN106707281A (zh) * | 2017-01-05 | 2017-05-24 | 北京航空航天大学 | 一种基于多频数据处理的机载D‑InSAR形变检测方法 |
CN107389029A (zh) * | 2017-08-24 | 2017-11-24 | 北京市水文地质工程地质大队 | 一种基于多源监测技术融合的地面沉降集成监测方法 |
Non-Patent Citations (5)
Title |
---|
D-InSAR技术在地震滑坡监测中的应用;任峰;《科技创新导报》;20170926(第19期);22-24 * |
GPS天线阵列变形监测系统设计;廖群 等;《工程勘察》;20060430(第4期);59-62 * |
基于BP的InSAR成像算法及多基线相位解缠算法研究;时代奇;《中国优秀硕士学位论文全文数据库信息科技辑》;20160315;I136-2364 * |
面向大型结构的调频连续波雷达位移测量系统VCO及回波的影响研究;陈丽;《中国优秀硕士学位论文全文数据库信息科技辑》;20160615;I136-356 * |
高精度InSAR成像及处理方法研究;邹光浩;《中国优秀硕士学位论文全文数据库信息科技辑》;20160315;I136-2283 * |
Also Published As
Publication number | Publication date |
---|---|
CN108088358A (zh) | 2018-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stewart et al. | HF radio measurements of surface currents | |
CN106353744B (zh) | 基于双基地fda-mimo雷达的多参数联合估计方法 | |
Gurgel et al. | Wellen Radar (WERA): A new ground-wave HF radar for ocean remote sensing | |
Long et al. | Mapping of North Atlantic winds by HF radar sea backscatter interpretation | |
CN106405556B (zh) | 车辆目标信息探测识别系统及其信号处理方法 | |
Stewart et al. | Radio measurements of oceanic winds at long ranges: An evaluation | |
CN101915911B (zh) | 基于相消积累空时谱的空间任意构型分布式sar动目标参数估计方法 | |
Ruegg et al. | Capabilities of dual-frequency millimeter wave SAR with monopulse processing for ground moving target indication | |
CN101666873A (zh) | 基于调变脉冲序列的高精度测距雷达的模糊处理方法 | |
CN106526583B (zh) | 一种基于天线方向图信息的地面运动目标定位方法 | |
CN102141611B (zh) | 斜视合成孔径雷达多普勒模糊数快速测定方法 | |
CN101363913B (zh) | 扇形波束圆锥扫描微波散射计信号处理方法 | |
Rajkumar et al. | Design and Development of DSP Interfaces and Algorithm for FMCW Radar Altimeter | |
CN108088358B (zh) | 一种基于多基线雷达轨道形变检测方法 | |
Arcone et al. | 1.4 GHz radar penetration and evidence of drainage structures in temperate ice: Black Rapids Glacier, Alaska, USA | |
JP2001501732A (ja) | 遅延補償型ドップラー・レーダー高度計 | |
Muller et al. | Phase center of L-band radar in polar snow and ice | |
Begum et al. | Analysis of rain cell size distribution for application in site diversity | |
RU2492504C1 (ru) | Способ определения нерадиальной проекции вектора скорости цели | |
Kovorotniy et al. | ATI SAR simulation shows signatures of complex objects | |
RU2686674C1 (ru) | Бесконтактный способ измерения пройденного пути | |
Li et al. | Ocean gravity wave phase velocity detection by HFSWR | |
Morey | Airborne sea ice thickness profiling using an impulse radar | |
RU2133480C1 (ru) | Радиолокационный способ определения параметров движения объекта | |
Tian et al. | Target localization based on doppler frequency and range sum measurement |
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 |