CN109725333A - 一种场景自适应的卫星信号接收及处理方法 - Google Patents

一种场景自适应的卫星信号接收及处理方法 Download PDF

Info

Publication number
CN109725333A
CN109725333A CN201811543803.6A CN201811543803A CN109725333A CN 109725333 A CN109725333 A CN 109725333A CN 201811543803 A CN201811543803 A CN 201811543803A CN 109725333 A CN109725333 A CN 109725333A
Authority
CN
China
Prior art keywords
satellite
signal
scene
measurement
adaptive
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
Application number
CN201811543803.6A
Other languages
English (en)
Other versions
CN109725333B (zh
Inventor
李久超
陈东
张千
柴源
陈明章
边炳秀
孙治国
高新宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Academy of Space Technology CAST
Original Assignee
China Academy of Space Technology CAST
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Academy of Space Technology CAST filed Critical China Academy of Space Technology CAST
Priority to CN201811543803.6A priority Critical patent/CN109725333B/zh
Publication of CN109725333A publication Critical patent/CN109725333A/zh
Application granted granted Critical
Publication of CN109725333B publication Critical patent/CN109725333B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radio Relay Systems (AREA)

Abstract

一种场景自适应的卫星信号接收及处理方法,首先产生测试用例,测量环境判定,并自适应调整增益,然后测量基础参数,完成天线自适应调整,获取直射分量、目标测量场景卫星信号,最后提取卫星信号特征,进而完成信道模型拟合。本发明与现有技术相比,面向S波段卫星,在不同经纬度和地表地貌情况下,测量从卫星轨位到地表的卫星信号传输信道参数,解决快速场景变化下的卫星信道参数测量问题,具有很好的使用价值。

Description

一种场景自适应的卫星信号接收及处理方法
技术领域
本发明涉及卫星通信测量领域,特别是一种场景自适应的卫星信号接收及处理方法。
背景技术
卫星信道特征是卫星通信体制设计的基础。而卫星移动信道衰落主要包括路径损耗、阴影衰落和多径衰落。理论上,路径衰落主要与传输距离和频率相关,阴影衰落服从Lognormal分布,多径衰落服从Rayleigh分布。美国、日本、加拿大和澳大利亚等国家针对测试点的经纬度、地形地貌以及电磁状态等,通过信号的特定参数测量,表征信号接收的实际变化情况。当前,卫星信道模型主要包括C.Loo、Corraza和Lutz模型,其测量方法是利用从气球、高塔或飞艇中放置信号发生装置,发射模拟卫星的信号,通过模拟的信道接收,测定卫星信号包络并进行数值统计,拟合后得到卫星信号特征参数。现有卫星信号接收及处理方法依赖于模拟的空中平台,其不能正确反映信号从卫星轨位到达地表的实际信道特征变化,同时空中平台的成本高难以支持长时间和大区域测量。
现有卫星信号接收及处理方法主要是以基站作为信号发射源发射测试信号,以测试仪表作为信号接收终端接收衰落信号以及分析特征。其设置接收节点获取待测量的码本组中码本数量,对信道进行测量,并向发射节点上报接收到的码本进行比对。此测量方法可扩展至多个接收节点,信道测量装置向至少两个用户设备发送天线域训练序列,信道测量装置对接收的至少两个用户设备的信道响应测量序列进行联合处理,得到下行信道状态信息。测试信号中可设定无线帧,其携带信道测量请求方标识信息、信道测量反馈方标识信息以及信道测量反馈指示信息,指示信道测量反馈方根据所述信道测量反馈指示信息进行信道测量。针对卫星条件下,主要考虑降雨、大气吸收以及大气闪烁等因素对卫星信道建模的影响,采用主成分分析和模糊聚类分析方法建立了卫星信道多状态模型。
综上所述,现有方法主要存在的缺点是:(1)缺乏S波段卫星的信道特征测量:卫星是以极化波束的形式辐射,现有测量方法面向地面移动通信系统采用全向辐射方式,其电磁波的传输模式上存在较大的差异性。同时,针对卫星的信道测量是以L波段和Ka波段为主,不同波段在传输损耗和大气衰减中均具有较大的差异性,难以直接应用;(2)缺乏灵活的快速场景变化的系统适应:现有卫星信道测量的装置集成度低,无法灵活装配和快速调试,而卫星测量环境往往随着测量设备载体运动而变化,需要快速调节设备状态,准确测定目标场景的信道参数。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种场景自适应的卫星信号接收及处理方法,通过增益自适应、天线追踪自适应、拟合算法自适应的方式实现对卫星信号的接收及处理。该方法面向S波段卫星,在不同经纬度和地表地貌情况下,通过场景自适应方法得到卫星信道参数,解决快速场景变化下的卫星信道参数测量问题。
本发明的技术解决方案是:一种场景自适应的卫星信号接收及处理方法,包括如下步骤:
(1)产生测试用例;
(2)测量环境判定,然后自适应调整增益;
(3)测量基础参数,然后自适应调整天线;
(4)获取直射分量;
(5)获取目标测量场景卫星信号;
(6)提取卫星信号特征;
(7)拟合信道模型。
所述的产生测试用例的方法为:
生成长码高斯伪随机序列,通过BPSK方式进行调制,经上变频向卫星发射,接收采集卫星信号的I、Q两路数据,通过IQ数据,测量得到卫星信号的包络、相位、多径效应、相干带宽、时延功率谱密度。
所述的BPSK调制后,经上变频向卫星发射的信号需要使得信号到达地表时功率不小于120dBm。
所述的测量环境判定,然后自适应调整增益的方法为:
接收采集卫星信号的I、Q两路数据时,对测试场景进行场景图片采集,然后对场景图片进行分析,若测试场景为空旷地,将信号接收增益调为50dB,若测试场景为阴影遮蔽地区,将信号接收增益调为60dB,若测试场景为城市密集区,将信号接收增益调为55dB。
所述的测量基础参数,然后自适应调整天线的方法为:
计算测量位置的理论卫星仰角E、方位角A、极化角P为
然后得到天线当前状态的实测卫星仰角E′、方位角A′、极化角P′为
其中
x=cosA·cosE
y=sinA·cosE
z=sinE
ω为自动追踪天线的俯仰角、δ为自动追踪天线的横滚角、为自动追踪天线的航向角、为卫星轨位的经度、为测量场景所在的经度、β为测量场景所在的纬度;
计算出天线误差信号ΔE=E-E′、ΔA=A-A′、ΔP=P-P′,采用自动控制算法,通过比例放大、微分、积分、比例放大来改变天线的方向使其对准卫星。
所述的获取直射分量的方法为:
将地面测试设备移动至测试场景周围的开阔地,对卫星信号进行滤波处理,得到信号的样本根据n=max(i)得到信号在当前场景下的LOS分量,i为正整数。
所述的获取目标测量场景卫星信号的方法为:
将地面测试设备移动至目标测量场景,接收卫星信号,并进行降噪预处理,然后采用二阶巴特沃斯滤波器,对频带外的噪声进行滤波,使得幅频响应为其中,m为滤波器阶数,ωc为截止频率,ω为通带频率。
所述的提取卫星信号特征的方法为:
将滤波后的I、Q两路数据通过的运算得到信号的包络,再利用进行归一化处理,计算出r′概率分布p(r′),然后由 计算r′残余概率分布CCDF,得到卫星信号的一阶统计量,然后生成与发射的长码高斯伪随机序列相同的序列,与接收到的下变频信号进行滑动互相关,得到Rpq(τ)=IFFT(FFT(p)×conj[FFT(q)]),进而得到卫星信号的二阶统计量;其中,Rpq(τ)为复数互相关值,p和q分别为接收到的下变频信号和数据处理单元生成的长码高斯伪随机序列。
所述的拟合信道模型的方法为:
步骤1:基于非线性最小二乘法的初值遍历自适应算法对卫星信号的一阶参量进行拟合,根据场景确定算法参数迭代初值范围以及搜索窗,b0∈[bmin,bmax],d0∈[dmin,dmax],μ0∈[μmin,μmax],若测试场景为空旷地,b0∈[0.1,0.3],d0∈[0.1,0.3],μ0∈[0.1,0.3],若测试场景为阴影遮蔽环境,b0∈[0.1,0.3],d0∈[0.6,0.8],μ0∈[-0.6,-0.4],若测试场景为城市密集区b0∈[0.5,0.7],d0∈[0.1,0.3],μ0∈[-0.6,-0.4],然后确定迭代初值b0=bmin,d0=dmin,μ0=μmin,给定拟合误差门限e,最大迭代次数it,梯度最小门限gmin,步长最小门限hmin,初始化阻尼系数λ0
步骤2:计算函数值为
其中,yi为由卫星信号得到的第i个实测数据,ri为根据地i个I、Q两路数据计算得到信号包络,f(b0,d0,μ0,ri)是由初值得到的模型函数值,若f2(b,d,μ)≤e,则退出迭代;若f2(b,d,μ)>e,则转入步骤3;
步骤3:计算拟合函数的雅可比矩阵
步骤4:根据Jk、fk计算得到梯度gk=Jk Tfk,若||gk||≤gmin,则退出迭代;若||gk||>gmin,则转入步骤5,其中,k为正整数,且初值为1;
步骤5:根据公式(Jk TJkkI)h=-gk,I为3×3单位矩阵,求解迭代步长h,若||h||≤hmin,则退出迭代;若||h||>hmin,则执行步骤6;其中,λk为第k次迭代时的阻尼系数,黄色部分为Jk
步骤6:根据得到新的参数若k<it,则返回步骤2,再次执行;否则执行步骤7;
步骤7:若b0∈[bmin,bmax],d0∈[dmin,dmax],μ0∈[μmin,μmax],则本初值的迭代最优解为fn 2(b,d,μ)min,令b0=b0+0.01,d0=d0+0.01,μ0=μ0+0.01,返回步骤2,进行下一组初值迭代;若b0>bmax,d0>dmax,μ0>μmax,则执行步骤8;fn 2是不同迭代初值下第n次拟合的最优解;
步骤8:将fn 2(b,d,μ)min中的最小值作为迭代最优解,并得到对应的b,d,μ,带入信道模型公式f(b,d,μ,r),得到符合当前测试场景的信道模型。
一种场景自适应的卫星信号接收及处理系统,包括信号发射模块、主控模块、环境监视模块、信号接收模块,信号接收模块包括自动追踪天线、伺服控制单元;其中:
信号发射模块生成长码高斯伪随机序列,调制后经上变频向卫星发射,信号接收模块中的自动追踪天线采集卫星信号的I、Q两路数据,测量得到接收此信号卫星信号的包络、相位、多径效应、相干带宽、时延功率谱密度,主控模块控制环境监视模块对测试场景进行场景图片采集,然后对采集的场景图片进行分析,确定测试场景,进而调整信号接收增益,信号接收模块中的自动追踪天线对测试所处地理位置进行测量与记录,得到测量位置的理论卫星仰角、方位角、极化角,信号接收模块中的伺服控制单元控制自动追踪天线在静止或移动过程中均能实现自动调整对星,主控模块获取直射分量、目标测量场景卫星信号,提取卫星信号特征,进而完成信道模型拟合。
本发明与现有技术相比的优点在于:
(1)本发明面向S波段卫星,在不同经纬度和地表地貌情况下,测量从卫星轨位到地表的卫星信号传输信道参数,解决快速场景变化下的卫星信道参数测量问题;
(2)本发明实现了卫星信号的自适应处理,提高了卫星信号处理的效率和精度。
附图说明
图1为本发明自适应的卫星信号接收及处理方法的实现流程图;
图2为实现卫星天线自适应追踪的流程图;
图3为自适应初值遍历拟合方法流程图。
具体实施方式
本发明克服现有技术的不足,提供了一种场景自适应的卫星信号接收及处理方法,通过增益自适应、天线追踪自适应、拟合算法自适应的方式实现对卫星信号的接收及处理。该方法面向S波段卫星,在不同经纬度和地表地貌情况下,通过场景自适应方法得到卫星信道参数,解决快速场景变化下的卫星信道参数测量问题。
一种自适应的卫星信号接收处理方法,实现步骤如下:
(1)信号发射模块产生测试用例,针对S波段卫星信道的测量,本方法设计单频正弦波和长码高斯伪随机序列等两种测量用例,可以对卫星信号的包络、相位、多径效应等特征进行测量;
(2)环境采集及增益自适应调节:环境监视单元将测试场景图片采集上传到主控单元,测量场景包括开阔地、阴影遮蔽地区、城市密集区等三类,根据不同的场景,信号接收单元调节不同的信号接收增益;
(3)基础参数的测量及天线追踪自适应调整:导航模块对测试所处地理位置进行测量与记录,并计算测量位置的理论卫星仰角(E)、方位角(A)、极化角(P),自适应天线根据卫星仰角、方位角、极化角,计算误差信号,采用自动控制算法,经过例放大、积分、微分、比例放大,通过伺服控制单元使天线在静止或移动过程中均能实现自适应调整对星;
(4)卫星信号特征提取:将测量装置移动到目标测量场景周围的空旷地,以获取卫星信号的直射分量,然后将测量装置移动到目标测量场景,对卫星信号进行测量。经滤波降噪处理和归一化处理后计算卫星信号的概率分布(PDF)参与概率分布(CCDF),得到卫星信号的一阶统计量。若测试信号为长码高斯伪随机序列,则可通过与本地产生的相同序列进行滑动互相关得到卫星信号的时延功率谱密度(PDP)。
(5)利用自适应的初值遍历拟合算法进行信道模型拟合:利用基于非线性最小二乘法的初值遍历自适应算法对卫星信号的一阶参量进行拟合,得到卫星信道模型参数。根据测试场景设置合适的参数初值搜索窗T,对参数进行迭代得到算法的局部最优解,对该搜索窗内的初值进行遍历,选取最小的拟合误差函数所对应的参数作为最佳拟合参数。
下面以静止状态下测量长码高斯伪随机序列为例,结合卫星信道测量流程,对本发明的实施方式进行进一步说明
(1)测试用例产生
信号发射模块生成长码高斯伪随机序列,通过BPSK方式进行调制,经上变频向卫星发射,接收采集卫星信号的I、Q两路数据,通过IQ数据,测量得到接收此信号卫星信号的包络、相位、多径效应、相干带宽、时延功率谱密度等特征,信号发射功率要求信号到达地表时功率不小于120dBm。
(2)测量环境判定及增益自适应调整
当接收采集卫星信号的I、Q两路数据时,主控单元控制环境监视单元(地面)对测试场景进行场景图片采集,然后主控单元对环境监视单元采集的场景图片进行分析,确定测试场景。若测试场景为空旷地,信号接收单元将信号接收增益调为50dB;若测试场景为阴影遮蔽地区,信号接收单元将信号接收增益调为60dB;若测试场景为城市密集区,信号接收单元将信号接收增益调为55dB。
(3)基础参数的测量及天线自适应调整
信号接收单元中的卫星自动追踪天线的导航模块对测试所处地理位置进行测量与记录,并计算测量位置的理论卫星仰角(E)、方位角(A)、极化角(P)。
其计算公式如下:
自动追踪天线根据卫星仰角、方位角、极化角,通过伺服控制单元使天线在静止或移动过程中均能实现自动调整对星。天线当前状态的实测卫星仰角(E′)、方位角(A′)、极化角(P′)满足以下公式:
x′、y′、z′满足
其中
x=cosA·cosE
y=sinA·cosE
z=sinE
上述公式中,ω为自动追踪天线的俯仰角、δ为自动追踪天线的横滚角、为自动追踪天线的航向角,均可通过角运动检测装置进行测量。
主控单元计算出天线当前状态的E′、A′、P′,与理论的E、A、P比较,得到误差信号ΔE=E-E′、ΔA=A-A′、ΔP=P-P′。利用该误差信号,采用自动控制算法,通过比例放大、微分、积分、比例放大来改变天线的方向使其对准卫星。
(4)直射分量的获取
将地面测试设备移动至测试场景周围的开阔地,对卫星信号进行滤波处理,得到信号的样本利用公式 得到信号在该场景下的LOS分量。
(5)目标测量场景卫星信号的获取
将地面测试设备移动至目标测量场景,接收卫星信号,并进行降噪预处理。采用二阶巴特沃斯滤波器,对频带外的噪声进行滤波。其幅频响应为:其中n为滤波器阶数,ωc为截止频率。
(6)卫星信号特征提取
经过滤波后,I、Q两路数据通过的运算得到信号的包络,利用进行归一化处理,并计算出r′概率分布p(r′)。由 计算r′残余概率分布(CCDF),得到卫星信号的一阶统计量。同时数据处理单元生成与发射的长码高斯伪随机序列相同的序列,与接收到的下变频信号进行滑动互相关,其时域公式为:
Rpq(τ)=IFFT(FFT(p)×conj[FFT(q)])
其中Rpq(τ)为复数互相关值,p和q分别为接收到的下变频信号和数据处理单元生成的长码高斯伪随机序列,得到卫星信号的二阶统计量。
(7)信道模型拟合
利用基于非线性最小二乘法的初值遍历自适应算法对卫星信号的一阶参量进行拟合。需要拟合的模型为f(b,d,μ,r)。其中b,d,μ为拟合参数,r为自变量。假设公式yi=f(b,d,μ,ri)-εi成立,yi是卫星信号包络概率分布实测值,f(b,d,μ,ri)是预测值。则εi为拟合后的误差函数。要使拟合误差函数最小,则可得
具体步骤如下:
步骤1:主控单元将环境监视单元得到的环境信息反馈给数据处理单元。数据处理单元根据场景确定算法参数迭代初值范围,以及三个参数的搜索窗。b0∈[bmin,bmax],d0∈[dmin,dmax],μ0∈[μmin,μmax]。若测试场景为空旷地,b0∈[0.1,0.3],d0∈[0.1,0.3],μ0∈[0.1,0.3];若测试场景为阴影遮蔽环境,b0∈[0.1,0.3],d0∈[0.6,0.8],μ0∈[-0.6,-0.4];若测试场景为城市密集区b0∈[0.5,0.7],d0∈[0.1,0.3],μ0∈[-0.6,-0.4]。确定迭代初值b0=bmin,d0=dmin,μ0=μmin。给定拟合误差门限e,最大迭代次数it,梯度最小门限gmin,步长最小门限hmin。初始化阻尼系数λ0
步骤2:根据当前参数,得到函数值
其中yi是由卫星信号得到的实测数据,f(b0,d0,μ0,ri)是由初值得到的模型函数值。若f2(b,d,μ)≤e,则退出迭代;若f2(b,d,μ)>e,则进行步骤3。
步骤3:计算拟合函数的雅可比矩阵
步骤4:将bk,dk,μk带入得到Jk,fk。并计算梯度gk=Jk Tfk。若||gk||≤gmin,则退出迭代;若||gk||>gmin,则进行步骤5。
步骤5:根据公式(Jk TJkkI)h=-gk,其中I为3×3单位矩阵,求解迭代步长h。若||h||≤hmin,则退出迭代;若||h||>hmin,则执行步骤6。
步骤6:通过公式得到新的参数若k<it,则返回步骤2,再次执行;否则执行步骤7。
步骤7:若b0∈[bmin,bmax],d0∈[dmin,dmax],μ0∈[μmin,μmax]。设本初值的迭代最优解为fn 2(b,d,μ)min。令b0=b0+0.01,d0=d0+0.01,μ0=μ0+0.01。返回步骤2,进行下一组初值迭代。若b0>bmax,d0>dmax,μ0>μmax,则执行步骤8.
步骤8.fn 2(b,d,μ)min中的最小值即为迭代最优解。得到所对应的的b,d,μ,带入信道模型公式f(b,d,μ,r),得到符合当前测试场景的信道模型。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。

Claims (10)

1.一种场景自适应的卫星信号接收及处理方法,其特征在于包括如下步骤:
(1)产生测试用例;
(2)测量环境判定,然后自适应调整增益;
(3)测量基础参数,然后自适应调整天线;
(4)获取直射分量;
(5)获取目标测量场景卫星信号;
(6)提取卫星信号特征;
(7)拟合信道模型。
2.根据权利要求1所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的产生测试用例的方法为:
生成长码高斯伪随机序列,通过BPSK方式进行调制,经上变频向卫星发射,接收采集卫星信号的I、Q两路数据,通过IQ数据,测量得到卫星信号的包络、相位、多径效应、相干带宽、时延功率谱密度。
3.根据权利要求2所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的BPSK调制后,经上变频向卫星发射的信号需要使得信号到达地表时功率不小于120dBm。
4.根据权利要求1或2所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的测量环境判定,然后自适应调整增益的方法为:
接收采集卫星信号的I、Q两路数据时,对测试场景进行场景图片采集,然后对场景图片进行分析,若测试场景为空旷地,将信号接收增益调为50dB,若测试场景为阴影遮蔽地区,将信号接收增益调为60dB,若测试场景为城市密集区,将信号接收增益调为55dB。
5.根据权利要求4所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的测量基础参数,然后自适应调整天线的方法为:
计算测量位置的理论卫星仰角E、方位角A、极化角P为
然后得到天线当前状态的实测卫星仰角E′、方位角A′、极化角P′为
其中
x=cosA·cosE
y=sinA·cosE
z=sinE
ω为自动追踪天线的俯仰角、δ为自动追踪天线的横滚角、为自动追踪天线的航向角、为卫星轨位的经度、为测量场景所在的经度、β为测量场景所在的纬度;
计算出天线误差信号ΔE=E-E′、ΔA=A-A′、ΔP=P-P′,采用自动控制算法,通过比例放大、微分、积分、比例放大来改变天线的方向使其对准卫星。
6.根据权利要求5所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的获取直射分量的方法为:
将地面测试设备移动至测试场景周围的开阔地,对卫星信号进行滤波处理,得到信号的样本根据n=max(i)得到信号在当前场景下的LOS分量,i为正整数。
7.根据权利要求6所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的获取目标测量场景卫星信号的方法为:
将地面测试设备移动至目标测量场景,接收卫星信号,并进行降噪预处理,然后采用二阶巴特沃斯滤波器,对频带外的噪声进行滤波,使得幅频响应为其中,m为滤波器阶数,ωc为截止频率,ω为通带频率。
8.根据权利要求7所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的提取卫星信号特征的方法为:
将滤波后的I、Q两路数据通过的运算得到信号的包络,再利用进行归一化处理,计算出r′概率分布p(r′),然后由 计算r′残余概率分布CCDF,得到卫星信号的一阶统计量,然后生成与发射的长码高斯伪随机序列相同的序列,与接收到的下变频信号进行滑动互相关,得到Rpq(τ)=IFFT(FFT(p)×conj[FFT(q)]),进而得到卫星信号的二阶统计量;其中,Rpq(τ)为复数互相关值,p和q分别为接收到的下变频信号和数据处理单元生成的长码高斯伪随机序列。
9.根据权利要求8所述的一种场景自适应的卫星信号接收及处理方法,其特征在于:所述的拟合信道模型的方法为:
步骤1:基于非线性最小二乘法的初值遍历自适应算法对卫星信号的一阶参量进行拟合,根据场景确定算法参数迭代初值范围以及搜索窗,b0∈[bmin,bmax],d0∈[dmin,dmax],μ0∈[μmin,μmax],若测试场景为空旷地,b0∈[0.1,0.3],d0∈[0.1,0.3],μ0∈[0.1,0.3],若测试场景为阴影遮蔽环境,b0∈[0.1,0.3],d0∈[0.6,0.8],μ0∈[-0.6,-0.4],若测试场景为城市密集区b0∈[0.5,0.7],d0∈[0.1,0.3],μ0∈[-0.6,-0.4],然后确定迭代初值b0=bmin,d0=dmin,μ0=μmin,给定拟合误差门限e,最大迭代次数it,梯度最小门限gmin,步长最小门限hmin,初始化阻尼系数λ0
步骤2:计算函数值为
其中,yi为由卫星信号得到的第i个实测数据,ri为根据地i个I、Q两路数据计算得到信号包络,f(b0,d0,μ0,ri)是由初值得到的模型函数值,若f2(b,d,μ)≤e,则退出迭代;若f2(b,d,μ)>e,则转入步骤3;
步骤3:计算拟合函数的雅可比矩阵
步骤4:根据Jk、fk计算得到梯度gk=Jk Tfk,若||gk||≤gmin,则退出迭代;若||gk||>gmin,则转入步骤5,其中,k为正整数,且初值为1;
步骤5:根据公式(Jk TJkkI)h=-gk,I为3×3单位矩阵,求解迭代步长h,若||h||≤hmin,则退出迭代;若||h||>hmin,则执行步骤6;其中,λk为第k次迭代时的阻尼系数,黄色部分为Jk
步骤6:根据得到新的参数若k<it,则返回步骤2,再次执行;否则执行步骤7;
步骤7:若b0∈[bmin,bmax],d0∈[dmin,dmax],μ0∈[μmin,μmax],则本初值的迭代最优解为fn 2(b,d,μ)min,令b0=b0+0.01,d0=d0+0.01,μ0=μ0+0.01,返回步骤2,进行下一组初值迭代;若b0>bmax,d0>dmax,μ0>μmax,则执行步骤8;fn 2是不同迭代初值下第n次拟合的最优解;
步骤8:将fn 2(b,d,μ)min中的最小值作为迭代最优解,并得到对应的b,d,μ,带入信道模型公式f(b,d,μ,r),得到符合当前测试场景的信道模型。
10.一种场景自适应的卫星信号接收及处理系统,其特征在于包括信号发射模块、主控模块、环境监视模块、信号接收模块,信号接收模块包括自动追踪天线、伺服控制单元;其中:
信号发射模块生成长码高斯伪随机序列,调制后经上变频向卫星发射,信号接收模块中的自动追踪天线采集卫星信号的I、Q两路数据,测量得到接收此信号卫星信号的包络、相位、多径效应、相干带宽、时延功率谱密度,主控模块控制环境监视模块对测试场景进行场景图片采集,然后对采集的场景图片进行分析,确定测试场景,进而调整信号接收增益,信号接收模块中的自动追踪天线对测试所处地理位置进行测量与记录,得到测量位置的理论卫星仰角、方位角、极化角,信号接收模块中的伺服控制单元控制自动追踪天线在静止或移动过程中均能实现自动调整对星,主控模块获取直射分量、目标测量场景卫星信号,提取卫星信号特征,进而完成信道模型拟合。
CN201811543803.6A 2018-12-17 2018-12-17 一种场景自适应的卫星信号接收及处理方法 Active CN109725333B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811543803.6A CN109725333B (zh) 2018-12-17 2018-12-17 一种场景自适应的卫星信号接收及处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811543803.6A CN109725333B (zh) 2018-12-17 2018-12-17 一种场景自适应的卫星信号接收及处理方法

Publications (2)

Publication Number Publication Date
CN109725333A true CN109725333A (zh) 2019-05-07
CN109725333B CN109725333B (zh) 2020-09-18

Family

ID=66296119

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811543803.6A Active CN109725333B (zh) 2018-12-17 2018-12-17 一种场景自适应的卫星信号接收及处理方法

Country Status (1)

Country Link
CN (1) CN109725333B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110505004A (zh) * 2019-08-13 2019-11-26 北京邮电大学 一种s波段卫星移动通信的卫星信号获取方法及装置
CN110928705A (zh) * 2019-10-15 2020-03-27 湖南大学 面向高性能计算应用的通信特征模型方法及系统
CN111211862A (zh) * 2019-12-17 2020-05-29 西安空间无线电技术研究所 一种遥感卫星信道码编码增益测试方法
CN112977887A (zh) * 2020-12-14 2021-06-18 深圳市星楷通讯设备有限公司 一种卫星的锁定方法及设备、计算机可读存储介质
CN113037354A (zh) * 2019-12-09 2021-06-25 大唐移动通信设备有限公司 一种系统设计方法及电子设备
CN113298113A (zh) * 2021-04-06 2021-08-24 北京交通大学 基于列车车载卫星定位观测数据的轨道沿线环境分类方法
CN113552565A (zh) * 2021-07-21 2021-10-26 中国矿业大学 针对sar数据高噪声及大梯度变化区域的相位解缠方法
WO2022022336A1 (zh) * 2020-07-28 2022-02-03 华为技术有限公司 业务卫星运行状态诊断方法及相关装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6784831B1 (en) * 2003-05-05 2004-08-31 Tia Mobile, Inc. Method and apparatus for GPS signal receiving that employs a frequency-division-multiplexed phased array communication mechanism
CN1980099A (zh) * 2005-12-01 2007-06-13 中国科学院上海微系统与信息技术研究所 高分辨率的实时多径信道测试方法、数据处理方法及装置
CN103457678A (zh) * 2013-02-22 2013-12-18 南京信息工程大学 多云天气条件下的气象卫星复合信道仿真模型的建立方法
CN104539341A (zh) * 2014-12-26 2015-04-22 西安电子科技大学 无线宽带多天线mimo全双工系统主动式回波自干扰抑制方案
CN105158775A (zh) * 2015-07-18 2015-12-16 西安电子科技大学 卫星定位的方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6784831B1 (en) * 2003-05-05 2004-08-31 Tia Mobile, Inc. Method and apparatus for GPS signal receiving that employs a frequency-division-multiplexed phased array communication mechanism
CN1980099A (zh) * 2005-12-01 2007-06-13 中国科学院上海微系统与信息技术研究所 高分辨率的实时多径信道测试方法、数据处理方法及装置
CN103457678A (zh) * 2013-02-22 2013-12-18 南京信息工程大学 多云天气条件下的气象卫星复合信道仿真模型的建立方法
CN104539341A (zh) * 2014-12-26 2015-04-22 西安电子科技大学 无线宽带多天线mimo全双工系统主动式回波自干扰抑制方案
CN105158775A (zh) * 2015-07-18 2015-12-16 西安电子科技大学 卫星定位的方法和装置

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110505004A (zh) * 2019-08-13 2019-11-26 北京邮电大学 一种s波段卫星移动通信的卫星信号获取方法及装置
CN110928705A (zh) * 2019-10-15 2020-03-27 湖南大学 面向高性能计算应用的通信特征模型方法及系统
CN110928705B (zh) * 2019-10-15 2023-09-22 湖南大学 面向高性能计算应用的通信特征分析方法及系统
CN113037354B (zh) * 2019-12-09 2022-06-21 大唐移动通信设备有限公司 一种系统设计方法及电子设备
CN113037354A (zh) * 2019-12-09 2021-06-25 大唐移动通信设备有限公司 一种系统设计方法及电子设备
CN111211862A (zh) * 2019-12-17 2020-05-29 西安空间无线电技术研究所 一种遥感卫星信道码编码增益测试方法
CN111211862B (zh) * 2019-12-17 2022-07-05 西安空间无线电技术研究所 一种遥感卫星信道码编码增益测试方法
WO2022022336A1 (zh) * 2020-07-28 2022-02-03 华为技术有限公司 业务卫星运行状态诊断方法及相关装置
CN112977887A (zh) * 2020-12-14 2021-06-18 深圳市星楷通讯设备有限公司 一种卫星的锁定方法及设备、计算机可读存储介质
CN113298113A (zh) * 2021-04-06 2021-08-24 北京交通大学 基于列车车载卫星定位观测数据的轨道沿线环境分类方法
CN113298113B (zh) * 2021-04-06 2023-09-26 北京交通大学 基于列车车载卫星定位观测数据的轨道沿线环境分类方法
CN113552565A (zh) * 2021-07-21 2021-10-26 中国矿业大学 针对sar数据高噪声及大梯度变化区域的相位解缠方法
CN113552565B (zh) * 2021-07-21 2023-07-18 中国矿业大学 针对sar数据高噪声及大梯度变化区域的相位解缠方法

Also Published As

Publication number Publication date
CN109725333B (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN109725333A (zh) 一种场景自适应的卫星信号接收及处理方法
CN109217956B (zh) 一种无人机电磁干扰通信环境半物理仿真方法及装置
Rice et al. Wideband channel model for aeronautical telemetry
CN110990505B (zh) 一种基于神经网络的Loran-C ASF修正方法
CN109916372B (zh) 信道状态信息不准确下的无人机基站最佳高度计算方法
Vogel et al. Measurement and modeling of land mobile satellite propagation at UHF and L-band
CN112383370B (zh) 一种卫星时变信道的建模仿真方法及系统
Bhardwaj et al. Geometrical-empirical channel propagation model for human presence at 60 GHz
CN106443721A (zh) 一种快速收敛的卫星导航天线阵抗干扰方法
CN105763274A (zh) 一种卫星移动通信系统中无线传输环境的仿真系统及方法
CN105704071A (zh) 基于信息序列的自适应渐消扩展卡尔曼粒子滤波的多普勒频移估计方法
Pastore et al. Refractivity inversions from point-to-point X-band radar propagation measurements
Ding et al. A time-varying transition channel model for air-ground communication
Tobie et al. Hybrid navigation filters performances between GPS, Galileo and 5G TOA measurements in multipath environment
Reyes-Guerrero Experimental broadband channel characterization in a sea port environment at 5.8 GHz
Das et al. Effects of transionospheric signal decorrelation on Global Navigation Satellite Systems (GNSS) performance studied from irregularity dynamics around the northern crest of the EIA
Biscarini et al. Statistical Modeling of Atmospheric Propagation Channel at W-Band through Sun-Tracking Microwave Radiometric Measurements for Non-Geostationary Satellite Links
Shi et al. Neural network based localization using outdoor lte measurements
Wessinger Improving near Surface Refractivity Estimates in Marine Environments
Rice et al. A wideband channel model for aeronautical telemetry—Part 1: Geometric considerations and experimental configuration
CN117082536B (zh) 一种基于强化学习的空地网络协同覆盖方法
Abdallah et al. Multipath mitigation of 5g signals via reinforcement learning for navigation in urban environments
Al-Saegh et al. Tracking-and scintillation-aware channel model for GEO satellite to land mobile terminals at Ku-band
Zheng et al. A modified s-band satellite channel simulation model
Zhou et al. 5G-based measurements and characterizations of low-altitude tethered balloon multipath channel

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