CN105182308A - 一种机载gnss海洋反射信号的生成方法 - Google Patents
一种机载gnss海洋反射信号的生成方法 Download PDFInfo
- Publication number
- CN105182308A CN105182308A CN201510532780.9A CN201510532780A CN105182308A CN 105182308 A CN105182308 A CN 105182308A CN 201510532780 A CN201510532780 A CN 201510532780A CN 105182308 A CN105182308 A CN 105182308A
- Authority
- CN
- China
- Prior art keywords
- signal
- time delay
- reflection
- correlation
- reflected
- 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
- 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
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种机载GNSS海洋反射信号的生成方法,GNSS接收机接收到的海面反射信号是许多条具有不同时延反射信号的叠加。因此,海面反射信号的时延相关曲线就是由许多个不同时延反射信号的相关曲线叠加而成的。根据这一原理,可假设有n条等时延间隔的反射信号,运用近似拟合的方法,求得其各自的幅值,进而得到其相对于直射信号的功率。再求出n条反射信号相对于直射信号的时延,反射信号的多普勒频率可直接用直射信号的多普勒频率。这样就可以通过时延、多普勒频率、功率进行GNSS海洋反射信号的生成。该方法对复杂的海面反射信号用近似拟合的方法进行简化,使得反射信号的模拟切实可行,为GNSS反射信号模拟器的研制奠定了基础。
Description
技术领域
本发明涉及一种机载GNSS海洋反射信号的生成方法,属于GNSS-R技术领域。
背景技术
GNSS-R技术是一种新兴的遥感技术,它是在GNSS***得到广泛应用的基础上发展起来,主要利用GNSS***提供的信号源,通过其反射信号来探测海面风场、海面高度、有效波高、海水盐度、土壤湿度等地球物理参数,并进行相应的反演。与其他海洋遥感技术相比,GNSS-R具有高时间和空间分辨率、低成本、信号源丰富、高机动性等优点,与其他探测手段优势互补,可以增加海洋探测手段的多样性,弥补局部探测手段不足的状况。因此,GNSS-R技术在海洋环境监测中具有重要的应用前景。
目前,国内外在反射信号建模和应用及专用接收机研制方面取得了大量的成果,并开展了大量的试验,包括机载、气球、星载试验以及岸基试验,其中以机载试验居多。由于GNSS-R的试验通常需要对测试场景有苛刻的要求,需要耗费大量的人力物力才能进行。而且,由于受到气候的影响,海情海面的不同,使得试验可重复性非常差。因此,需研制一套可用于实验室测试的GNSS-R海洋反射信号模拟***。
发明内容
本发明的目的是为了解决上述问题,提出了一种机载GNSS海洋反射信号的生成方法,为GNSS反射信号模拟器的研制提供保障。
本发明的一种机载GNSS海洋反射信号的生成方法,包括以下几个步骤:
(1)通过ZV模型得到GNSS海洋反射信号的相关功率曲线,以直射信号作归一化处理,再开方得到归一化的反射信号相关曲线。
(2)对归一化的反射信号相关曲线进行近似计算,得到n个等间隔时延的反射信号幅值的数学解。
(3)对n个反射信号幅值的数学解进行拟合,求得双指数函数的待定系数。
(4)根据拟合函数和直射信号的功率,就可求得n条不同时延反射信号的功率。
(5)计算镜点反射信号相对于直射信号的时延,则可知道n条反射信号相对于直射信号的时延。忽略反射信号多普勒频率的影响,则可设n条反射信号的多普勒频率与直射信号的多普勒频率相同。这样,就得到了生成n条反射信号所需的时延、多普勒频率、功率3个参数。
本发明的优点在于:
采用近似拟合的方法对GNSS信号经过海面反射的复杂现象进行了简化,模拟适当数量的反射信号,既能反映出海洋风场的特性,又易于工程实现。
附图说明
图1是GNSS海洋反射信号模拟的流程图。
图2是归一化的反射信号相关曲线图。
图3是等时延反射信号自相关曲线组成的示意图。
图4是GNSS-R的几何关系图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是一种机载GNSS海洋反射信号的生成方法,流程如图1所示,包括以下几个步骤:
步骤一、通过ZV模型(Zavorotny和Voronovich建立的海面风场模型)得到GNSS海洋反射信号的相关功率曲线,以直射信号作归一化处理,再开方得到归一化的反射信号相关曲线;
具体为:
(1)设置接收机和卫星的基本参数。
在设置参数前,选择本地散射坐标系,即设镜面反射点为坐标系的原点,导航卫星、接收机和镜面反射点在一个平面内。关于接收机的参数是:接收机高度、接收机速度矢量,关于导航卫星的参数是:卫星高度、卫星速度矢量、卫星高度角、卫星信号功率。
(2)设置海洋相关参数。
GNSS导航信号经过海面反射后,反射信号承载了海面特性的特性信息,而海面特性与海洋参数是相关联的。关于海洋的主要参数有:海面风速、风向、海温、盐度、介电常数。(3)计算ZV模型时延一维相关功率曲线。
海面反射信号的时延—多普勒二维相关功率形式如下:
其中,Ti是相干积分时间;D为接收机天线的增益;Λ是GPS伪随机码的自相关函数;S为多普勒滤波函数;σ0为粗糙海面的标准化散射截面;为镜像反射点到平均海面上某一点的向量;R0,R分别为GPS卫星和接收机到端点的距离。τ为时延,fc为镜点处的多普勒频率,fD为不同反射点处的多普勒频率。
海面反射信号相关功率模型为时间延迟τ和多普勒频率fD两变量的函数,当fD固定时,为时延一维相关功率曲线。
(4)对时延相关功率曲线归一化。
为了表达结果的相对值,相关功率可以相对于直射信号进行归一化,可表示为
其中:Rd为GPS卫星到接收机的距离,Ti为是相干积分时间。PcorrR为反射信号的相关功率,PcorrD为直射信号的相关功率。
对归一化的时延相关功率曲线开方,得到反射信号相关曲线,如图2所示。
步骤二、对反射信号相关曲线进行幅值近似计算。
接收机接收到的反射信号是多条具有不同时延反射信号的叠加。忽略多普勒频率的影响,那么反射信号相关曲线就是本地伪码与不同时延信号的相关值。假设接收到的反射信号由时延间隔1/4码片的16条反射信号叠加而成。假设每个码片被采样成16个点,两条反射信号之间的时延差4个采样点。
对归一化的时延相关功率曲线开方,得到反射信号相关曲线。众所周知,CA码自相关曲线是个三角形。那么,16个自相关曲线乘以其相对应的幅值叠加后就是反射信号的相关曲线。图3是反射信号自相关曲线组成的示意图。
图3中线条a被假设为镜面反射点的信号,在第3个采样点处可假设只有镜点信号,通过反射信号相关曲线第3个采样点的相关值与CA码自相关值的比可得到镜点信号的幅值。反射信号相关曲线第7个采样点的值由两条信号的相关值叠加,则第7个采样点的相关值减去第1路信号在该采样点时的相关值,再除以CA码自相关值可得到第2路的幅值。以此类推,可得到16路反射信号的幅值。
步骤三、对反射信号幅值的结果进行拟合。
近似计算得到的结果中可能有负值,因此需要其进行拟合。拟合函数采用双指数的数学模型,其表达式如下:
A=a×exp(-b·τ)+c×exp(-d·τ)(3)
其中,A表示信号的幅值,τ表示时延。a、b、c和d为待定系数。τ=0为镜点反射信号的时延。
步骤四、用拟合函数计算反射信号功率。
由于相关功率曲线的归一化因子是直射信号,因此根据拟合公式可得到直射信号幅值为1时的不同时延反射信号幅值。然后根据幅值与功率的关系可以计算得到不同时延反射信号的功率。
假如在步骤一中卫星信号功率设为-158.5dBW(即1.4×10-16W),也就是直射信号的功率为1.4×10-16W,由信号功率和幅值的关系可得直射信号的幅值为1.6×10-8V。假如τ=0时,由拟合函数(式3)计算A=0.25,则反射信号的幅值为1.6×10-8V*0.25=4.2×10-9V,再计算反射信号的功率为9×10-18W,也即-170dBW。
步骤五、计算镜点反射信号相对于直射信号的时延。
图4反映了GNSS-R反射事件的几何关系,图4中T代表导航卫星,R代表接收机,S为镜面反射点,hR为接收机平台相对基准水平面的高度,θ为反射信号对应卫星的仰角。
利用GNSS-R反射事件的几何关系可以获取镜点反射信号相对于直射信号的时延值为:
ρE=2hRsinθ(4)
然后根据镜点反射信号相对于直射信号的时延,就可以得到各条反射信号相对于直射信号的时延。另外,由于忽略多普勒频率的影响,可假设各条反射信号的多普勒频率与直射信号的多普勒频率相同。这样,就可以得到16条反射信号的时延、多普勒频率、功率3个参数,从而实现机载GNSS海洋反射信号的生成。
Claims (1)
1.一种机载GNSS海洋反射信号的生成方法,包括以下几个步骤:
步骤一、通过ZV模型,得到GNSS海洋反射信号的相关功率曲线,对直射信号作归一化处理,最终得到归一化的反射信号相关曲线;
具体为:
(1)设置接收机和卫星的基本参数;
设置接收机参数:接收机高度、接收机速度矢量;
设置卫星参数:卫星高度、卫星速度矢量、卫星高度角、卫星信号功率;
(2)设置海洋相关参数;
海洋相关参数:海面风速、风向、海温、盐度、介电常数;
(3)计算ZV模型时延一维相关功率曲线;
海面反射信号相关功率模型:
其中,Ti是相干积分时间;D为接收机天线的增益;Λ是GPS伪随机码的自相关函数;S为多普勒滤波函数;σ0为粗糙海面的标准化散射截面;为镜像反射点到平均海面上某一点的向量;R0,R分别为GPS卫星和接收机到端点的距离;τ为时延,fc为镜点处的多普勒频率,fD为不同反射点处的多普勒频率;
根据海面反射信号相关功率模型,给定fD,得到时延一维相关功率曲线;
(4)对时延相关功率曲线进行归一化;
对时延相关功率曲线进行归一化:
其中:Rd为GPS卫星到接收机的距离,Ti为是相干积分时间;PcorrR为反射信号的相关功率,PcorrD为直射信号的相关功率;
对归一化的时延相关功率曲线进行开方,得到反射信号相关曲线;
步骤二、对反射信号相关曲线进行幅值近似计算;
设反射信号相关曲线为本地伪码与不同时延信号的相关值,假设接收到的反射信号由时延间隔1/4码片的16条反射信号叠加而成,每个码片被采样成16个点,两条反射信号之间的时延差4个采样点;
设置镜面反射点的信号,假设在第3个采样点处只有镜点信号,通过反射信号相关曲线第3个采样点的相关值与CA码自相关值的比,得到镜点信号的幅值;反射信号相关曲线第7个采样点的值由两条信号的相关值叠加,则第7个采样点的相关值减去第1路信号在该采样点时的相关值,再除以CA码自相关值,得到第2条的幅值,以此类推,得到16条反射信号的幅值;
步骤三、对反射信号幅值的结果进行拟合;
对反射信号幅值的结果进行拟合,拟合函数为:
A=a×exp(-b·τ)+c×exp(-d·τ)(3)
其中,A表示信号的幅值,τ表示时延;a、b、c和d为待定系数;τ=0为镜点反射信号的时延;
步骤四、用拟合函数计算反射信号功率;
根据拟合函数,得到直射信号幅值为1时的不同时延反射信号幅值,根据幅值与功率的关系,计算得到不同时延反射信号的功率;
步骤五、计算镜点反射信号相对于直射信号的时延;
利用GNSS-R反射事件的几何关系,获取镜点反射信号相对于直射信号的时延值为:
ρE=2hRsinθ(4)
其中,hR为接收机平台相对基准水平面的高度,θ为反射信号对应卫星的仰角;
根据镜点反射信号相对于直射信号的时延,得到各条反射信号相对于直射信号的时延,假设各条反射信号的多普勒频率与直射信号的多普勒频率相同,最终,得到16条反射信号的时延、多普勒频率、功率,完成机载GNSS海洋反射信号的生成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510532780.9A CN105182308B (zh) | 2015-08-27 | 2015-08-27 | 一种机载gnss海洋反射信号的生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510532780.9A CN105182308B (zh) | 2015-08-27 | 2015-08-27 | 一种机载gnss海洋反射信号的生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105182308A true CN105182308A (zh) | 2015-12-23 |
CN105182308B CN105182308B (zh) | 2018-01-12 |
Family
ID=54904508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510532780.9A Expired - Fee Related CN105182308B (zh) | 2015-08-27 | 2015-08-27 | 一种机载gnss海洋反射信号的生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105182308B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107037460A (zh) * | 2017-05-27 | 2017-08-11 | 金华航大北斗应用技术有限公司 | 一种用于海面测高的北斗meo卫星机载反射信号建模方法 |
CN107037461A (zh) * | 2017-05-27 | 2017-08-11 | 金华航大北斗应用技术有限公司 | 一种可用于海面测风的北斗meo卫星机载反射信号模拟器 |
CN107272033A (zh) * | 2017-05-27 | 2017-10-20 | 金华航大北斗应用技术有限公司 | 一种北斗meo卫星机载反射信号的随机时延建模方法 |
CN110440685A (zh) * | 2019-07-16 | 2019-11-12 | 金华航大北斗应用技术有限公司 | 一种gnss-r岸基海面溢油遥感探测方法 |
CN112817017A (zh) * | 2021-02-06 | 2021-05-18 | 北京航空航天大学 | 一种考虑涌浪影响的gnss反射信号建模方法 |
CN112817018A (zh) * | 2021-02-06 | 2021-05-18 | 北京航空航天大学 | 一种考虑洋流影响的gnss反射信号建模方法 |
CN113050134A (zh) * | 2021-03-19 | 2021-06-29 | 中国人民解放军92859部队 | 基于卫星导航信息的海面风场反演观测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000022454A1 (en) * | 1998-10-15 | 2000-04-20 | Northrop Grumman Corporation | Moving target simulator |
JP2000147112A (ja) * | 1998-11-12 | 2000-05-26 | Sekisui Chem Co Ltd | 相関式探査装置 |
US20030164793A1 (en) * | 2002-03-01 | 2003-09-04 | Jurgen Ludewig | Device for generating a transit time delay of a pulsed radar signal and method for operation thereof |
CN104678386A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种利用gnss海面反射信号相关功率探测目标的方法 |
-
2015
- 2015-08-27 CN CN201510532780.9A patent/CN105182308B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2000022454A1 (en) * | 1998-10-15 | 2000-04-20 | Northrop Grumman Corporation | Moving target simulator |
JP2000147112A (ja) * | 1998-11-12 | 2000-05-26 | Sekisui Chem Co Ltd | 相関式探査装置 |
US20030164793A1 (en) * | 2002-03-01 | 2003-09-04 | Jurgen Ludewig | Device for generating a transit time delay of a pulsed radar signal and method for operation thereof |
CN104678386A (zh) * | 2015-02-11 | 2015-06-03 | 北京航空航天大学 | 一种利用gnss海面反射信号相关功率探测目标的方法 |
Non-Patent Citations (1)
Title |
---|
白伟华: "GNSS_R海洋遥感技术研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107037460A (zh) * | 2017-05-27 | 2017-08-11 | 金华航大北斗应用技术有限公司 | 一种用于海面测高的北斗meo卫星机载反射信号建模方法 |
CN107037461A (zh) * | 2017-05-27 | 2017-08-11 | 金华航大北斗应用技术有限公司 | 一种可用于海面测风的北斗meo卫星机载反射信号模拟器 |
CN107272033A (zh) * | 2017-05-27 | 2017-10-20 | 金华航大北斗应用技术有限公司 | 一种北斗meo卫星机载反射信号的随机时延建模方法 |
CN107272033B (zh) * | 2017-05-27 | 2019-07-02 | 金华航大北斗应用技术有限公司 | 一种北斗meo卫星机载反射信号的随机时延建模方法 |
CN107037460B (zh) * | 2017-05-27 | 2019-07-02 | 金华航大北斗应用技术有限公司 | 一种用于海面测高的北斗meo卫星机载反射信号建模方法 |
CN110440685A (zh) * | 2019-07-16 | 2019-11-12 | 金华航大北斗应用技术有限公司 | 一种gnss-r岸基海面溢油遥感探测方法 |
CN112817017A (zh) * | 2021-02-06 | 2021-05-18 | 北京航空航天大学 | 一种考虑涌浪影响的gnss反射信号建模方法 |
CN112817018A (zh) * | 2021-02-06 | 2021-05-18 | 北京航空航天大学 | 一种考虑洋流影响的gnss反射信号建模方法 |
CN113050134A (zh) * | 2021-03-19 | 2021-06-29 | 中国人民解放军92859部队 | 基于卫星导航信息的海面风场反演观测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105182308B (zh) | 2018-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105182308A (zh) | 一种机载gnss海洋反射信号的生成方法 | |
CN103345759B (zh) | 一种海底大型复杂沙波地貌的精确探测方法 | |
CN102749622B (zh) | 基于多波束测深的声速剖面及海底地形的联合反演方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN104062660A (zh) | 一种基于时域离散InSAR干涉对的矿区地表时序形变监测方法 | |
CN103023586B (zh) | 一种天波超视距雷达电离层信道仿真方法 | |
CN103093057A (zh) | 一种船舶导航雷达信号仿真方法 | |
CN104698457A (zh) | 一种迭代曲面预测InSAR成像及高度估计方法 | |
CN102914773B (zh) | 一种多航过圆周sar三维成像方法 | |
CN103576137B (zh) | 一种基于成像策略的多传感器多目标定位方法 | |
CN103293521B (zh) | 一种利用x波段雷达探测近海海域水深的方法 | |
CN102841385A (zh) | 一种基于多重分形克里金法的局部地磁图构建方法 | |
CN105021843A (zh) | 600kHZ宽带声学多普勒流速剖面仪及其实现方法 | |
CN103344653B (zh) | 一种基于双gps接收机的土壤湿度实时测量***及方法 | |
CN102749067A (zh) | 自适应多级计算的精确反演海面高程遥测装置 | |
CN111950438B (zh) | 基于深度学习的天宫二号成像高度计有效波高反演方法 | |
CN104730518A (zh) | 一种基于高斯拟合的雷达多普勒谱估计海面流场的方法 | |
CN104318593A (zh) | 一种雷达海杂波的仿真方法和*** | |
CN105204022A (zh) | 海面风场的反演方法及装置 | |
Park et al. | New approach to sea surface wind retrieval from GNSS-R measurements | |
Cario et al. | Analysis of error sources in underwater localization systems | |
CN103778633B (zh) | 确定数字高程模型单元网格遮挡的方法及装置 | |
CN103105603B (zh) | 一种x波段测波雷达海流反演预处理方法 | |
CN109632258A (zh) | 一种基于矢量传感器的收发分离的海洋内波声学检测方法 | |
CN103487793A (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 | ||
GR01 | Patent grant | ||
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: 20180112 Termination date: 20190827 |