本发明属于无线通信技术领域,尤其涉及一种变加速无人机移动中继的通信方法。
背景技术:
近年来,无人机协作通信技术成为了无线通信领域研究的热点。相比于传统的地面通信,无人机协作通信易于实现按需分配,从而具有更高的成本效益;无人机具有高度的移动性,因而部署起来更加灵活快捷;无人机与地面终端的信道链路多为视距链路,能提供更好的信道环境。因此,无人机将在未来的无线通信领域中发挥极为重要的作用,其应用主要包括:(1)作为临时基站;(2)作为移动中继;(3)用于物联网。
目前,大量文献研究了无人机作为移动中继时的无人机中继协作通信系统的信息容量和频谱效率的优化问题。与此同时,无人机具有有限的机载能量,因此节能问题被认为是无人机通信的一个关键问题。当前大部分研究仅考虑匀速或匀加速无人机中继通信,未考虑变加速情况,没有充分考虑无人机移动中继初、末速度及加速度大小的约束,也未能考虑总吞吐量要求(任务量)对无人机中继系统的影响。
技术实现要素:
为克服现有技术的缺点和不足,本发明的目的在于提供一种变加速无人机移动中继的通信方法。
本发明是这样实现的,一种变加速无人机移动中继的通信方法,该方法包括以下步骤:
s1、根据无人机移动中继r在源节点s与目的节点d之间的加速飞行模式与信息收发形式,建立无人机移动中继信息的最优收发时间分配函数;
根据最优收发时间分配函数将无人机移动中继r的最优加速度的求解拆分为两个过程:无人机中继从源节点接收至少q比特数据时的最低能耗以及无人机中继向目的节点发送q比特数据时的最低能耗;
s2、设定算法分别对上述两个过程求解并得到各过程的最优飞行速度,根据两个最优飞行速度得出无人机中继的最优加速度;
s3、无人机中继根据源节点处的初速度和所述最优加速度进行飞行,并且按照所述最优收发时间分配函数进行数据信息的接收和转发。
优选地,在步骤s1中,所述加速飞行模式为:固定高度位置下,无人机移动中继从源节点正上方以初速度从源节点出发,以先加速、后减速的方式直线飞行,并到达目的节点正上方停止。
优选地,在步骤s1中,所述信息收发形式为:所述无人机移动中继采用时分双工模式以及解码转发中继协议,且所述无人机移动中继从源节点飞至目的节点的过程中,源节点通过无人机移动中继将至少q比特数据发送给目的节点。
优选地,在步骤(1)中,所述最优收发时间分配函数为:
{λ*(n),n=1,2,…,n}
式中,将地面节点s到d分为n段(n为偶数),每一段的距离为
本发明克服现有技术的不足,提供一种变加速无人机移动中继通信方法及系统,本发明主要提供了一个变加速的无人机中继系统,也就是无人机通过变加速的形式进行飞行,飞行过程中,无人机通过调整接收和发送消息的时间和飞行加速度,以此在完成给源节点转发需要发送的信息比特数的条件下,最小化无人机中继的能量消耗。相比于现有技术的缺点和不足,本发明具有以下有益效果:本发明根据源节点待传输的信息比特数和源节点与目的地节点之间的距离,对时分双工无人机中继的收发时间及飞行加速度进行联合调整,在满足源节点待发送信息能够完全传输至目的地节点的条件下,实现无人机总飞行能耗的最小化,仿真实验验证了所提算法的有效性。
附图说明
图1是本发明方法涉及的系统示意图;
图2是无人机中继从源节点接收至少q比特数据时的最低能耗问题求解算法收敛性及系统能效下界有效性示意图;
图3是不同收发时间分配函数下无人机中继系统的能效性能比较;
图4是不同任务量对无人机移动中继速度的影响示意图;
图5是无人机中继不同飞行方案的能效性能对比;
图6是无人机匀加速与变加速飞行能效性能对比;
图7是不同距离下任务量与能效性能关系示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提供一种变加速无人机移动中继的通信方法,根据源节点待传输信息比特数和源节点与目的节点之间的距离,对时分双工无人机中继的收发时间及飞行加速度进行联合调整,在满足源节点待发送信息能够完全传输至目的节点的条件下,实现无人机总飞行能耗的最小化,具体的,该方法包括以下步骤:
s1、根据无人机移动中继r在源节点s与目的节点d之间的加速飞行模式与信息收发形式,建立无人机移动中继信息的最优收发时间分配函数,根据最优收发时间分配函数将无人机移动中继r的最优加速度的求解拆分为两个过程:无人机中继从源节点接收至少q比特数据时的最低能耗以及无人机中继向目的节点发送q比特数据时的最低能耗
在步骤s1中,如图1所示,一个由单天线地面源节点s、地面目的节点d和空中固定翼无人机中继节点r组成的三节点无人机移动中继通信系统,固定翼无人机中继以初速度v0从源节点s正上方出发,沿直线飞向目的节点d,以末速度vf到达目的节点d正上方停止,假设节点s与d之间的距离为l,源节点s需要至少将q比特数据发送给目的节点d,由于阻塞,节点s与d之间无法直接进行通信,所有信息都需通过无人机中继r进行转发,无人机中继r工作在固定高度h处,实际情况下,h可以是无人机避免碰撞的最低高度,这样可以保证无人机在工作过程中不需要频繁的进行升降;
这里只考查无人机的飞行阶段,忽略其起飞、降落过程,假设节点s和d分别位于坐标(0,0,0)和(l,0,0)处,无人机从s的正上方出发,这里无人机r工作在时分双工模式,采用的是解码转发(df)中继协议。
假设无人机r与地面节点s、d间的通信链路由视距链路组成,此时无人机通信系统的时变信道遵循自由空间路径衰落模型,同时假设无人机的数据缓冲区足够大。
为了便于数学说明,将地面节点s到d分为n段(n为偶数),每一段的距离为
假设无人机在每段路径中的位置选取在每段的中点,则在第n(n∈{1,2,…,n})段中,无人机与s和d的距离可分别写为
假设无人机中继r运动产生的多普勒效应能被完全消除,则在第n个段中,节点s和d到无人机中继r的信道增益hsr(n)和hrd(n)可分别表示为
其中,β表示单位距离的信道功率损耗,其值取决于天线增益、载波频率等,α表示大尺度衰落因子,由于无人机与地面节点之间的信道由视距链路构成,因此,α取2,此时,在第n段,无人机中继r或目的节点d可能接收到的信号yr(n)或yd(n)可以表示为:
yr(n)=hsr(n)xs(n) zr(n),n=1,…,n(5)
yd(n)=hrd(n)xr(n) zd(n),n=1,…,n(6)
其中,xs(n)表示节点s发送的信号,xr(n)表示节点r发送的信号,zr(n)和zd(n)表示第n段r和d处均值为0,方差为σ2的高斯白噪声;
假设地面节点s和无人机中继r的发射功率分别为ps和pr,此时源节点s到中继节点r的链路以及中继节点r到目的节点d的链路的瞬时信道容量rsr(n)和rrd(n)可分别写为:
其中,b表示信道带宽,
由于无人机中继工作在时分双工模式,为此定义收发时间分配函数λ(n)
当λ(n)=0时,表示第n段仅源节点s向无人机中继r发送数据,当λ(n)=1时,表示第n段仅无人机中继r向目的节点d发送数据;
无人机中继r从s正上方沿直线飞向d的过程中,从s接收到的总信息比特数和向d发送的总信息比特数可分别表示为
其中,
表示无人机r通过第n段的时间,v(n)和a(n)表示第n段无人机r的初速度与加速度;
需要说明的是:在第n段,无人机的初速度和加速度分别是v(n)和a(n),则在第n段末尾,即,在第n 1段初,无人机r的速度可由式(13)给出;此外,无人机r在整个飞行过程中,从源节点s正上方出发的初速度为v0,到达目的节点d正上方的末速度为vf,因此,对于v0和vf有下面的式(14)成立;
v(n 1)=v(n) a(n)t(n),n=1,2,…,n-1(13)
v(1)=v0,v(n 1)=vf(14)
式(13)中的t(n)表示无人机r通过第n段的时间,式(14)中的v(n 1)表示第n段的末速度;
由于在每个段中,无人机r不能自己产生数据,只能转发已经从源节点s接收到的信息,因此存在下面的信息因果约束:
需要说明的是:在第1段中,无人机中继r中没有数据发送给d,因此有λ(1)=0,同样,在第n段中,r若仍从s接收数据,则此段接收的数据将无法发送给d,因此有λ(n)=1,因此系统吞吐量
对于系统吞吐量需要满足式(17)给出的约束条件:
即无人机中继系统至少将q比特数据转发给目的节点d;
一般情况下,无人机中继r的能耗由2个部分组成,1)无人机维持飞行和向前推进所需要的能耗;2)无人机所携带的通信设备运行和无人机本身电路的静态能耗;
需要说明的是:实际情况下,通信和静态电路的能耗远小于无人机飞行推进的能耗,因此,这里忽略通信和静态电路的能耗,式(18)给出了固定翼无人机沿直线飞行的能耗:
其中,v(n)表示第n段无人机的飞行速度,t(n)表示通过第n段的时间,m表示无人机的质量,v0和vf分别表示无人机的初、末速度,g表示重力加速度,c1和c2是无人机的空气动力参数,有c1=ρcd0s/2、c2=2w2/[(πe0ar)ρs],这里ρ表示空气密度,cd0表示零升力阻力系数,s表示参考面积(机翼面积),e0是oswald效率(翼展效率)其典型值在0.7到0.85之间,w表示无人机整体的重量,ar表示无人机翼的纵横比;
由于固定翼无人机在飞行过程中存在一定的起飞速度和降落速度,这里假设无人机的初速度等于末速度,即,此时沿直线飞行的无人机中继系统的能耗可表示为:
这里,t(n)表示无人机r通过第n段的时间,v(n)是无人机r在第n段的初速度,并且v(n)满足式(13)和式(14);
由式(16)和式(19),无人机移动中继系统的能效可表示为:
在确保中继系统至少转发q比特数据的条件下,研究通过优化无人机中继加速度、收发时间分配,实现无人机移动中继通信系统的能效最大化,该优化问题可描述为:
s.t.|a(n)|≤amax,n=1,…,n(21b)
vmin≤v(n)≤vmax,n=1,…,n(21c)
v(n 1)=v(n) a(n)t(n),n=1,2,…,n-1(21f)
v(1)=v0,v(n 1)=vf,v0=vf(21g)
其中,t(i)表示无人机r通过第i段的时间,具体可由式(12)给出,
下面考虑ps=pr情况,然后求解问题(21),即,式(21a)、式(21b)、式(21c)、式(21d)、式(21e)、式(21f)和式(21g)组成的优化问题,给出最优收发时间函数及飞行加速度的求解算法;
由于问题(21)是一个非凸问题,无法直接求解,因此,将问题(21)分解为两个子问题:1)在固定加速度下,求解最优收发时间分配函数,2)将求得的最优收发时间分配函数代入原问题,求解出最优加速度。
上述的第1个子问题可转化为标准lp问题通过凸优化工具直接求解,对于第2个子问题,取其目标函数的下界,通过泰勒公式将非线性约束转化为线性约束,解得近似问题的闭式解,最后利用迭代算法对近似问题的闭式解进行优化;
因为,在给定加速度{a(n),n=1,…,n}下,问题(21)的目标函数的分母为常数,上述第1个子问题可写为:
这里,对于收发时间分配函数λ(i)的求解属于整数规划问题,直接求解复杂度高,因此,考虑将收发时间分配函数λ(i)松弛为[0,1]区间上的连续变量ξ(i),将整数规划变量连续化,此外,由式(13)可知,当k=n时,有
0≤ξ(n)≤1,n=1,2,…,n(23d)
其中,
问题(23),即,式(23a)、式(23b)、式(23c)和式(23d)组成的优化问题,可以用矩阵的形式重新表述为:
s.t.ω2x≤b(24b)
0≤ξ(n)≤1(24d)
其中,ω1=[rrd(2)rrd(3)…rrd(n)],x=[ξ(2)ξ(3)…ξ(n)]t,
c=[rrd(2)rrd(3)…rrd(n)],
b=[rsr(1)rsr(1) rsr(2)…rsr(1) rsr(2) … rsr(n-1)]t;
问题(24),即,式(24a)、式(24b)、式(24c)和式(24d)组成的优化问题,属于标准lp问题,可以通过凸优化工具直接求得解,然后将求得的解采用“舍入化整”的方法,可得最优收发时间分配函数{λ*(n),n=1,…,n},此时,对于最优收发时间分配函数,存在如下定理:
定理1:当pr=ps,v0=vf时,最优收发时间分配函数{λ*(n),n=1,…,n}为
定理1说明:在无人机数据缓冲区足够大时,无人机移动中继的收发方式一定是在更靠近源节点s时全部接收数据,靠近目的节点d时全部发送数据,无人机移动中继工作过程中,不存在先接收后发送,然后再接收的情况,一定是接收数据全部完成后再开始发送;
下面求解第2个子问题,即,将最优收发时间分配函数代入问题(21),再求解出最优加速度;
下面将对问题(21)进一步进行处理,根据式(17)可给出问题(21)目标函数的下界:
在给定收发时间分配函数情况下,问题(21)可转化为:
s.t.|a(n)|≤amax,n=1,…,n(27b)
vmin≤v(n)≤vmax,n=1,…,n(27c)
需要说明的是:当n→∞时,2la(n)/n→0,此时,式(12)给出的t(n)可以通过等价无穷小公式转化为:
将式(28)代入式(19),则式(27a)可转化为
此时,加速度a(n)可表示为:
由式(13)和式(30),约束条件式(27b)和式(27c)可以合并表示为式(31)和式(32),
式(31)和式(32)中的vmin(n)和vmax(n)是中间变量,用于方便后面的数学表示。
由定理1可知,无人机移动中继r在数据缓冲区足够大情况下,在全部接收完成后才会向目的节点d发送消息,此时,约束条件式(27d)可等效为式(33)、式(34)和式(35):
根据最优收发时间分配函数,问题(27),即,式(27a)、式(27b)、式(27c)和式(27d)组成的联合优化问题,可拆分为两个问题,即,分别求解:1)无人机中继r从节点s接收至少q比特数据时的最低能耗,2)无人机中继r向节点d发送q比特数据时的最低能耗。
s2、设定算法分别对上述两个过程求解并得到各过程的最优飞行速度,根据两个最优飞行速度得出无人机中继的最优加速度
在步骤s2中,无人机中继r从节点s接收至少q比特数据时的最低能耗问题可描述为:
vmin(n)≤v(n)≤vmax(n)(36c)
由于式(36b)是非线性约束,直接求解问题(36),即,式(36a)、式(36b)和式(36c)组成的联合优化问题,计算复杂度高,因此需要对问题(36)进行化简;
取满足式(36b)和式(36c)约束的vk(n),k∈{1,2,3,…},这里的vk(n)表示第k次迭代后的v(n)值(具体迭代算法见下面的算法1),对式(36b)在vk(n)处进行泰勒展开,有:
此时,问题(36)可重新描述为:
s.t.vmin(n)≤v(n)≤vmax(n)(38b)
对于问题(38),即,式(38a)、式(38b)和式(38c)组成的联合优化问题,首先忽略约束条件式(38b),然后采用拉格朗日乘子法,其拉格朗日函数为:
这里,θ1为拉格朗日乘子,此时,kkt条件可由式(40)、式(41)、式(42)和式(43)给出;
θ1≥0(43)
需要说明的是:满足式(40)、式(41)、式(42)、式(43)和式(38b)的解即为问题(38)的最优解,下面将对其进行分情况讨论。
case1:当θ1=0时,
由式(40)可得
定理2:若θ1=0,当
case2:当θ1>0时,
由式(40)可得
定理3:若θ1>0,当
同样,对于无人机中继r向节点d发送q比特数据时的最低能耗问题,可描述为
vmin(n)≤v(n)≤vmax(n)(44c)
对于问题(44),即,式(44a)、式(44b)和式(44c)组成的联合优化问题,采用解问题(36)相同的方法,则kkt条件可由式(45)、式(46)、式(47)和式(48)给出;
θ2≥0(48)
这里,θ2为拉格朗日乘子,下面将对其进行分情况讨论;
case1:当θ2=0时,
定理4:若θ2=0,当
case2:当θ2>0时,
定理5:若θ2>0,当
对于问题(38),根据定理2和定理3,给出一种迭代算法获得问题(38)的解,步骤如下:
算法1:求解问题(38)的迭代算法(这里vk(n)表示第k次迭代后的v(n)值)
步骤11、令k=0、index=0,定义一个大于零但接近于零的容差变量ε,初始化
步骤12、令n=2,vk 1(1)=v0;
步骤13、根据式(31)、式(32)和vk 1(n-1)计算得出vmin(n)和vmax(n);
步骤14、假设θ1=0,将
步骤15、假设θ1>0,将
步骤16、令n=n 1,若
步骤17、将求得的vk 1(n)代入式(29)计算系统能耗,将计算出的系统能耗记为ek 1,若|ek 1-ek|≤ε则结束计算;否则令k=k 1,返回步骤12;
步骤18、迭代结束,若index=0,
对于问题(44),采用上述同样的方法,给出一种迭代算法获得问题(44)的解,步骤如下:
算法2:求解问题(44)的迭代算法
步骤21、令k=0、index=0,定义一个大于零但接近于零的容差变量ε,初始化
步骤22、令
步骤23、根据式(31)、式(32)和vk 1(n-1)计算得出vmin(n)和vmax(n);
步骤24、θ2是引入的拉格朗日乘子,假设θ2=0,将
步骤25、假设θ2>0,将
步骤26、令n=n 1,若n>n则执行步骤27,否则返回步骤23;
步骤27、将求得的vk 1(n)代入式(29)计算系统能耗,并将计算出的系统能耗记为ek 1,若|ek 1-ek|≤ε则结束计算;否则令k=k 1,返回步骤22;
步骤28、迭代结束,若index=0,
对于问题(21),即无人机中继的收发时间分配及飞行加速度的联合优化问题的求解,可通过算法3解得:
算法3:无人机中继的收发时间分配及飞行加速度的联合优化问题求解算法
步骤31、通过定理1得到问题(21)的最优收发时间分配函数{λ*(n),n=1,2,…,n};
步骤32、将所得最优收发时间分配函数{λ*(n),n=1,2,…,n}带入问题(21),取其目标函数的下界并通过等价无穷小变化,将问题(21)转化为问题(38)和问题(44);
步骤33、通过算法1解得
步骤34、通过算法2解得
步骤35、令v*(n 1)=vf,将所求得{v*(n),n=1,2,…,n}和初、末速度带入式(30)即可求得最优加速度{a*(n),n=1,2,…,n};
步骤36、输出{λ*(n),n=1,2,…,n}、{v*(n),n=1,2,…,n}和{a*(n),n=1,2,…,n},算法结束。
上述算法3给出的{λ*(n),n=1,2,…,n}和{a*(n),n=1,2,…,n}即为问题(21)的最终解。
s3、无人机中继根据源节点处的初速度和所述最优加速度进行飞行,并且按照所述最优收发时间分配函数进行数据信息的接收和转发。
在步骤例s3中,针对本发明提出的联合实现方法,进行了仿真实验,并且与“吞吐量最大化飞行”、“匀速飞行”、“以最大加速度加速后匀速飞行”、“匀加速飞行”四种方案的系统能效性能进行了比较,实验环境为matlab环境。假设无人机移动中继与地面节点通信链路由los链路构成,且源节点s的发射功率ps与中继r的发射功率pr相等,均为10dbm;无人机中继r及目的节点d处的噪声功率谱密度相同,为-170dbm/hz;单位距离信道增益β=-50db;系统其他参数如表1所示:
表1
图2证明了算法1的收敛性。假定源节点s与目的节点d之间的距离l=1000m,段数n=2000,任务量q=150mbit。由图2可以看出,随着迭代次数的上升,迭代能效值快速趋于稳定,说明算法1单调收敛。将式(26)取得的能效下界与式(20)取得的精确能效值进行对比,图中虚线是由式(26)给出的能效下界,而实线是由式(20)给出的能效准确值,图2说明了所取的能效下界是有效的。
对于不同收发时间分配函数,无人机移动中继的能效也不同。同样假定源节点s与目的节点d之间的距离l=1000m,段数n=2000。假设n=1,2,…,μ段用于无人机中继接收源节点s发送的数据,而n=μ 1,μ 2,…,n段用于向目的节点d转发数据。当μ值比较小时,意味着无人机中继用于接收的段少,系统吞吐量变小;当μ值比较大时,意味着无人机中继用于转发的段比较少,同样吞吐量变小。定理1证明了对于任意任务量的无人机中继系统而言,
图4给出不同任务量q对无人机飞行速度的影响。同样,假设l=1000m,n=2000,对比了q=50mbit、q=300mbit及q=600mbit三种不同任务量对无人机飞行速度的影响。图4中,当任务量比较小时,即q=50mbit时,无人机移动中继在接收段先以较大的加速度加速到最低能耗速度后再匀速飞行,由于存在末速度vf的约束,无人机移动中继在后半段会减速直至目的节点d的上空。当任务量中等,即q=300mbit时,无人机先以最大加速度加速到一定速度后缓慢加速,发送段做接收段的镜像运动。对于任务量较大情况,即q=600mbit时,无人机在开始会以初速度v0(最低速度)匀速飞行一段时间以获得更多的数据,然后做缓慢的加速,为了满足任务量的约束,无人机移动中继在这种情况下全程速度都比较低。图4说明,变加速无人机可以通过改变其飞行速度,满足系统任务量的要求。
图5、图6对比不同飞行模式下的无人机移动中继通信系统的能效性能。同样,假设l=1000m,n=2000。这里对比了“吞吐量最大化飞行”、“匀速飞行”、“以最大加速度加速后匀速飞行”、“匀加速飞行”四种方案下无人机中继系统的能效,具体如图5、图6所示。结果表明:“吞吐量最大化飞行”方案能效最低,原因是无人机低速飞行时能耗较大,“匀速飞行”方案在任务量较小时可以在出发点以较高的速度飞行,忽略了加速(减速)的过程,因此能效要高于本文所提变加速飞行方案,当任务量较大时,匀速飞行方案则需要降低速度以获得吞吐量,能效要低于本文所提变加速飞行方案;“以最大加速度加速飞行”方案和“匀加速飞行”方案无法根据任务量灵活调整加速度,在任务量较低情况下,能效较低。综上,本文所提变加速飞行模式可根据任务量灵活调整飞行速度,获得较高的能效,通过仿真可以看出,其能效要优于其他飞行模式的无人机移动中继。
从图7中可以看出,当l=1000m、任务量q=150mbit时,无人机移动中继能效最大。当l=1500m、任务量q=200mbit时,无人机中继能效最大。随着距离l的增加,无人机移动中继系统取得最优能效的任务量q也增加。因此,从能效的角度,所提算法更加适用于短距离、任务量小的情景或是长距离、任务量大的情景。
以上所述仅为本发明的较佳实施而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
1.一种变加速无人机移动中继的通信方法,其特征在于,该方法包括以下步骤:
s1、根据无人机移动中继r在源节点s与目的节点d之间的加速飞行模式与信息收发形式,建立无人机移动中继信息的最优收发时间分配函数;
根据最优收发时间分配函数将无人机移动中继r的最优加速度的求解拆分为两个过程:无人机中继从源节点接收至少q比特数据时的最低能耗以及无人机中继向目的节点发送q比特数据时的最低能耗;
s2、设定算法分别对上述两个过程求解并得到各过程的最优飞行速度,根据两个最优飞行速度得出无人机中继的最优加速度;
s3、无人机中继根据源节点处的初速度和所述最优加速度进行飞行,并且按照所述最优收发时间分配函数进行数据信息的接收和转发。
2.如权利要求1所述的变加速无人机移动中继的通信方法,其特征在于,在步骤s1中,所述加速飞行模式为:固定高度位置下,无人机移动中继从源节点正上方以初速度从源节点出发,以先加速、后减速的方式直线飞行,并到达目的节点正上方停止。
3.如权利要求1所述的变加速无人机移动中继的通信方法,其特征在于,在步骤s1中,所述信息收发形式为:所述无人机移动中继采用时分双工模式以及解码转发中继协议,且所述无人机移动中继从源节点飞至目的节点的过程中,源节点通过无人机移动中继将至少q比特数据发送给目的节点。
4.如权利要求1所述的变加速无人机移动中继的通信方法,其特征在于,在步骤(1)中,所述最优收发时间分配函数为:
{λ*(n),n=1,2,…,n}
式中,将地面节点s到d分为n段(n为偶数),每一段的距离为