一种基于卡尔曼滤波的直升机旋翼遮挡预测方法技术

本发明专利技术公开了一种基于卡尔曼滤波的直升机旋翼遮挡预测的方法，其包括：当直升机接收到被遮挡的信号时，对信号功率进行滑动平均运算，基于预设的功率阈值，统计遮挡周期内遮挡周期长度以及遮挡时间长度，根据旋翼的初始遮挡状态，计算旋翼的遮挡周期以及遮挡率，通过卡尔曼滤波得到下一个周期的遮挡周期以及遮挡率的预测值，然后计算下一周期的缝隙时间长度的预测值，即完成了旋翼遮挡预测。利用该方法可以用于直升机返向链路中，实现突发传输，并且能够有效地减少旋翼遮挡对信号传输的影响，提高信道利用率。

A prediction method for helicopter rotor shelter based on Calman filtering

The invention discloses a method of helicopter rotor occlusion prediction based on Calman filter. It includes: when the helicopter receives the shielded signal, the signal power is slid averaging, based on the predetermined power threshold, the length of the occlusion cycle and the length of the occlusion time are calculated in the period of the occlusion period, according to the rotor. The occlusion period and the occlusion rate of the rotor are calculated in the initial occlusion state. The occlusion cycle of the next cycle and the prediction value of the occlusion rate are obtained by Calman filter. Then the prediction value of the length of the slot time in the next period is calculated, that is, the rotor occlusion prediction is completed. This method can be used in the back link of the helicopter to realize the burst transmission, and can effectively reduce the influence of the rotor occlusion on the signal transmission and improve the channel utilization.

【技术实现步骤摘要】
一种基于卡尔曼滤波的直升机旋翼遮挡预测方法
本专利技术涉及直升机抗旋翼遮挡的卫星通信
，具体涉及一种直升机旋翼遮挡预测方法。
技术介绍
目前，直升机在反恐怖、抢险救灾、处理边防突发事件等非战争军事行动中起到越来越重要的作用。直升机卫星通信系统是由直升机和地面固定站通过同步卫星构成点对点通信系统。由于卫星通信电波传播方式是直射波，要求在无遮挡的条件下通信，而直升机飞行过程中，旋翼桨叶会周期性的遮挡天线，造成通信信号的周期性衰落，影响正常通信。直升机卫星通信的返向链路是指直升机经卫星到地面基站之间的通信链路。在返向链路中，常采用旋翼同步突发，直升机在发送数据时，首先检测旋翼缝隙，预测突发时间窗口，然后周期性的在旋翼同步缝隙中将通信信息突发出去。旋翼遮挡预测一般采用能量和信噪比估计的方式，但该方式预测精度往往较差，如何精确的预测旋翼遮挡，保证信息无损失的传输，提高信道利用率显得十分必要。
技术实现思路
本专利技术的专利技术目的在于：针对上述存在的问题，提供一种基于卡尔曼滤波的直升机旋翼遮挡预测方法，本专利技术通过预测下一周期的遮挡周期长度以及遮挡时间长度，来减少旋翼遮挡对信号传输的影响，从而提高信道的利用率以及信号传输的可靠性。为了实现上述专利技术目的，本专利技术采用的技术方案是：计算接收信号的信号功率，并对功率进行滑动平均运算处理；基于预置的功率阈值，判断滑动平均运算处理后的信号功率是否大于(或满足大于等于也可)功率阈值，若是，则判定信号处于缝隙时间内；否则信号处于遮挡时间内；并统计旋翼的遮挡周期长度len_p以及旋翼的遮挡时间长度len_s；根据旋翼的初始遮挡周期To、旋翼遮挡周期长度len_o，计算旋翼的当前遮挡周期T以及遮挡率τ，所述遮挡周期所述遮挡率τ＝len_s/len_p；将遮挡周期T以及遮挡率τ作为状态变量，通过卡尔曼滤波预测旋翼下一周期的遮挡周期以及遮挡率根据预测的旋翼的遮挡周期和遮挡率计算旋翼的缝隙时间长度的预测值。进一步的，计算接收信号的信号功率，对功率进行滑动平均运算处理具体为：计算接收信号的每个采样点的信号功率x(n)，其中n为采样点标识符；基于预设的滑动平均长度len，对各采样点的信号功率x(n)进行滑动平均运算，得到滑动平均后的信号功率s(n)，其中s(n)＝sum(x(n:n+len-1))/len，函数sum(x(n:n+len-1))表示第n个采样的信号功率x(n)到第n+len-1个采样点的信号功率x(n+len-1)的累加和。进一步的，若第n个采样点的滑动平均后的信号功率s(n)大于功率阈值，则当前采样点处于缝隙时间内；否则处于遮挡时间内；并统计第一次处于遮挡状态的第一个采样点到下一次处于遮挡状态的第一个采样点的采样点数，得到旋翼的遮挡周期长度len_p；统计一个遮挡周期内处于遮挡状态的第一个采样点到最后一个采样点的采样点数，得到旋翼的遮挡时间长度len_s。其中，根据预测的遮挡周期以及遮挡率计算缝隙时间长度的预测值具体为：根据公式得到旋翼的遮挡周期长度的预测值根据公式得到旋翼的遮挡时间长度的预测值根据公式得到旋翼的缝隙时间长度的预测值。综上所述，由于采用了上述技术方案，本专利技术的有益效果是：本专利技术的方法针对直升机返向链路通信，提出了一种基于卡尔曼滤波的旋翼遮挡预测方法，该方法能够有效的预测出缝隙时间，从而保证信息无损失发送，并且该方法实现简单，已在硬件中得到实现应用。附图说明图1为本专利技术的处理流程示意图；图2为直升机旋翼遮挡模型；图3为实施例1的测试处理过程示意图；图4为实施例1的旋翼遮挡模型；图5为实施例1的遮挡周期内数据点数对比图图6为实施例1的遮挡率对比图。具体实施方式为使本专利技术的目的、技术方案和优点更加清楚，下面结合实施方式和附图，对本专利技术作进一步地详细描述。本专利技术通过对接收到的信号功率进行滑动平均运算，基于预置阈值统计遮挡周期内遮挡周期长度以及遮挡时间长度，计算旋翼的遮挡周期以及遮挡率，通过卡尔曼滤波得到下一个周期的遮挡周期以及遮挡率的预测值，然后计算下一周期的缝隙时间长度的预测值，即可实现遮挡预测，参见图1，其具体步骤为：步骤A：计算接收的信号的功率，并对功率进行滑动平均运算，得到信号平均功率；步骤B：基于预置的功率阈值，当信号平均功率小于该功率阈值时，则判断信号处于遮挡时间内，否则处于缝隙时间内，并统计遮挡周期长度以及遮挡时间长度；步骤C：根据旋翼的初始遮挡状态，计算旋翼的遮挡周期以及遮挡率；步骤D：有别于现有的预测方式(将计算得到的遮挡周期以及遮挡率直接作为下一时刻的预测值)，本专利技术通过卡尔曼滤波预测下一周期的遮挡周期以及遮挡率。步骤E：根据预测的遮挡周期以及遮挡率计算缝隙时间长度的预测值。其中，各步骤具体实现过程分别为：所述步骤A中，直升机旋翼遮挡会对信号的功率造成衰减，因此可以根据信号功率判断信号是否处于遮挡状态。当直升机接收到信号后，计算信号功率，对该信号功率进行滑动平均运算，将滑动平均后信号功率记为s(n)，则s(n)＝sum(x(n:n+len-1))/len,n＝1,2,…，其中，x(n)为第n个采样点的信号功率，len为滑动平均的长度，sum(x(n:n+len-1))表示第n个采样的信号功率x(n)到第n+len-1个采样点的信号功率x(n+len-1)的累加和。所述步骤B中，根据步骤A中得到的滑动平均后的信号功率，设置功率阈值TAP，当功率小于该功率阈值时，判断此时处于遮挡状态，否则当功率大于该功率阈值时，判断此时处于非遮挡状态，即：统计第一次处于遮挡状态的第一个点到下一次处于遮挡状态的第一个点的长度(即采样点个数)，则可得到遮挡周期长度len_p，统计一个遮挡周期内处于遮挡状态的第一个点到最后一个点的长度(采样点个数)，则可得到遮挡时间长度len_s，如图2所示。所述步骤C中，将遮挡时间长度除以遮挡周期长度得到遮挡率τ，即τ＝len_s/len_p，并利用遮挡周期长度计算遮挡周期T：其中To表示初始的遮挡周期，len_o表示初始的遮挡周期长度。所述步骤D中，在原有的旋翼遮挡预测的基础上，即将计算得到的遮挡周期以及遮挡率作为下一时刻的预测值，加入了卡尔曼滤波理论，根据卡尔曼滤波原理，将遮挡周期T以及遮挡率τ作为状态变量X，即则此时卡尔曼滤波的状态方程以及量测方程为：其中卡尔曼滤波预测递推公式为：状态预测：状态预测误差的协方差矩阵：Pk/k-1＝Pk-1+Qk-1；卡尔曼增益：Kk＝Pk/k-1(Pk/k-1+Rk)-1；状态估计：Pk＝(I-Kk)Pk/k-1；状态估计误差的协方差矩阵：状态预测的输出即为下一个周期内的遮挡周期以及遮挡率的预测值。所述步骤E中，根据预测得到的遮挡周期以及遮挡率，计算遮挡周期长度的预测值即以及遮挡时间长度的预测值即因此可以得到缝隙时间长度的预测值即实施例1：本实施例所使用仿真平台为MATLABR2011b，其测试过程如图3所示，其中遮挡模型如图4所示，遮挡周期To为50ms，一个遮挡周期内数据长度为40000，遮挡率τ为0.25，一个遮挡周期内遮挡时间长度为10000，遮挡衰减最大为20db，信噪比为3db。首先，随机产生1，-1数据作为发送数据，经遮挡模型处理并加入高斯噪声，然后计算加噪后数据的功率，并进行滑动本文档来自技高网...

【技术保护点】
1.一种基于卡尔曼滤波的直升机旋翼遮挡预测方法，其特征在于，包括下列步骤：计算接收信号的信号功率，并对功率进行滑动平均运算处理；基于预置的功率阈值，判断滑动平均运算处理后的信号功率是否大于功率阈值，若是，则判定信号处于缝隙时间内；否则信号处于遮挡时间内；并统计旋翼的遮挡周期长度len_p以及旋翼的遮挡时间长度len_s；根据旋翼的初始遮挡周期To、旋翼遮挡周期长度len_o，计算旋翼的当前遮挡周期T以及遮挡率τ，所述遮挡周期

【技术特征摘要】
1.一种基于卡尔曼滤波的直升机旋翼遮挡预测方法，其特征在于，包括下列步骤：计算接收信号的信号功率，并对功率进行滑动平均运算处理；基于预置的功率阈值，判断滑动平均运算处理后的信号功率是否大于功率阈值，若是，则判定信号处于缝隙时间内；否则信号处于遮挡时间内；并统计旋翼的遮挡周期长度len_p以及旋翼的遮挡时间长度len_s；根据旋翼的初始遮挡周期To、旋翼遮挡周期长度len_o，计算旋翼的当前遮挡周期T以及遮挡率τ，所述遮挡周期所述遮挡率τ＝len_s/len_p；将遮挡周期T以及遮挡率τ作为状态变量，通过卡尔曼滤波预测旋翼下一周期的遮挡周期以及遮挡率根据预测的旋翼的遮挡周期和遮挡率计算旋翼的缝隙时间长度的预测值。2.如权利要求1所述的方法，其特征在于，计算接收信号的信号功率，对功率进行滑动平均运算处理具体为：计算接收信号的每个采样点的信号功率x(n)，其中n为采样点标识符；基于预设的滑动平均长度len，对各采样点的信号功率x...

【专利技术属性】
技术研发人员：戴彬彬，何春，李浩，朱立东，王剑，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51