一种基于轨道角动量波束模式扩展的宽度测量方法

1.本发明涉及无线电测量技术领域，尤其涉及一种基于轨道角动量波束模式扩展的宽度测量方法。


背景技术：

2.基于无线电方法的测量技术主要应用在测距方面，按照其工作原理可以分为三种：第一种是脉冲测距，也称为时间测距，记录脉冲信号发射点到接收点的传播时延，再根据电波的传播速度计算两点之间的距离；第二种相位测距，根据电波从发射点到接收点的相位差来测量两点之间的距离；第三种是频率测距，利用时变频率源随时间变化的特性，记录发射点到接收点之间的频率变化量来测量两点之间的距离。
3.这三种方法分别利用了电波的相位、频率以及固有传播速度的特性进行测距，但这些电波特性难以应用到物体的尺度测量上；随着对电波特性的深入研究，发现了轨道角动量(orbital angular momentum，oam)这一全新的维度。当携带轨道角动量的电波在传播过程中受到干扰时，其轨道角动量的模式随之发生扰动，产生模式间串扰，能量泄露到其他模式；根据模式之间的串扰与干扰源之间对应关系，即可实现对干扰源的检测或测量，因此，如何通过轨道角动量来进行检测是现阶段需要考虑的问题。


技术实现要素：

4.本发明的目的在于克服现有技术的缺点，提供了一种基于轨道角动量波束模式扩展的宽度测量方法，可以对物体尺度进行直接测算。
5.本发明的目的通过以下技术方案来实现：一种基于轨道角动量波束模式扩展的宽度测量方法，包括以下步骤：
6.步骤s1，确定待检测物的距离s与最大测量宽度w
max
。
7.步骤s2，确定轨道角动量波束的模式l，计算公式如下：
[0008][0009]
其中f-1
表示不同类型发射源的轨道角动量模式l与波束扩散夹角对应关系。
[0010]
步骤s3，采用抛物方程法，计算不同遮挡宽度下的模式扩展场e，具体实现过程为：
[0011]
步骤s31，根据电波频率f、测量区间尺寸s、w
max
划分计算网格，设置x、y、z方向的网格步长分别为δx、δy、δz；
[0012]
步骤s32，根据轨道角动量波束的模式l、计算网格构建初始口径场u0；
[0013]
步骤s33，根据待测量尺寸w对传播环境进行建模；
[0014]
步骤s34，采用分布傅里叶变换法对模式扩展场进行迭代求解，计算公式如下：
[0015]
[0016]
其中分别表示傅里叶变换与逆傅里叶变换，k
x
表示传播常数在x方向上的分量。
[0017]
步骤s4，将笛卡尔坐标系下的场值分布转换成柱坐标下的场分布，计算公式为：
[0018][0019][0020]
步骤s5，根据oam波束不同模式之间的正交性，计算各个模式的能量强度，计算公式如下：
[0021][0022]
其中n表示计算的扩展模式阶数，i表示虚数符号。单个模式的功率可以表示为：
[0023][0024]
最后可以得到模式扩展场的功率谱，计算公式如下：
[0025][0026]
步骤s6：使用最小二乘法对功率谱与遮挡宽度曲线进行多项式拟合，具体实现过程为：
[0027]
步骤s61：给定函数w＝f(p)，设p为powern的主模，在p1,p2,p3,
…
,pm处的函数值为w1,w2,w3,
…
,wm，其中n表示扩展模式阶数，m表示给定拟合点的个数。设多项式f(p)＝a0 a1p a2p2
……
akpk，其中k表示多项式的阶数，使得
[0028][0029]
为了求得符合条件的系数ai，对上式每个系数求偏导数，即
[0030]
[0031]
对以上方程组化简整理后，用矩阵表示为
[0032][0033]
步骤s62：求解上述矩阵中的各个系数ai，即可得到关于p的多项式拟合函数w(p)。
[0034]
步骤s7：根据函数w(p)实现待检测物体尺寸宽度的快速测量。
[0035]
本发明具有以下优点：一种基于轨道角动量波束模式扩展的宽度测量方法，根据电波的轨道角动量这一全新的维度特性，拓展了无线电方法在测量技术方面的应用，可以实现对物体尺度的直接快速测量。
附图说明
[0036]
图1为本发明提出的基于轨道角动量波束模式扩展的宽度测量方法的流程图；
[0037]
图2为本发明实验例中未受到遮挡情况下模式为2的轨道角动量功率谱图；
[0038]
图3为本发明实验例中受到1m宽度的遮挡后的轨道角动量模式扩展场分布；
[0039]
图4为本发明实验例中受到遮挡后模式扩散场从-5到 5阶的轨道角动量功率谱分布；
[0040]
图5为本发明实验例中在各个宽度遮挡情况下，主模式的功率谱与遮挡宽度之间的对应关系散点图，与通过最小二乘法拟合出的多项式曲线结果；
[0041]
图6为本发明实验例中在0～1之间随机选取的20个测试值，表示遮挡物的宽度尺寸；
[0042]
图7为本发明实验例中所选的20个遮挡宽度测试值所测算出的主模功率谱分布，与拟合的多项式曲线值对比结果。
具体实施方式
[0043]
为使本技术实施例的目的、技术方案和优点更加清楚，下面将结合本技术实施例中附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本技术实施例的组件可以以各种不同的配置来布置和设计。因此，以下结合附图中提供的本技术的实施例的详细描述并非旨在限制要求保护的本技术的保护范围，而是仅仅表示本技术的选定实施例。基于本技术的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本技术保护的范围。下面结合附图对本发明做进一步的描述。
[0044]
如图1所示，本发明涉及一种基于轨道角动量波束模式扩展的宽度测量方法，包括以下步骤：
[0045]
步骤s1，确定待检测物的距离s与最大测量宽度w
max
。
[0046]
步骤s2，确定轨道角动量波束的模式l，计算公式如下：
[0047][0048]
其中f-1
表示不同类型发射源的轨道角动量模式l与波束扩散夹角对应关系。
[0049]
步骤s3，采用抛物方程法，计算不同遮挡宽度下的模式扩展场e，具体实现过程为：
[0050]
步骤s31，根据电波频率f、测量区间尺寸s、w
max
划分计算网格，设置x、y、z方向的网格步长分别为δx、δy、δz；
[0051]
步骤s32，根据轨道角动量波束的模式l、计算网格构建初始口径场u0；
[0052]
具体为，根据轨道角动量波束的模式l，通过公式构建初始口径场u0，其中e0表示电磁波幅值，i表示虚数单位，表示方位角。
[0053]
步骤s33，根据待测量尺寸w对传播环境进行建模；
[0054]
步骤s34，采用分布傅里叶变换法对模式扩展场进行迭代求解，计算公式如下：
[0055][0056]
其中分别表示傅里叶变换与逆傅里叶变换，k
x
表示传播常数在x方向上的分量。
[0057]
步骤s4，将笛卡尔坐标系下的场值分布转换成柱坐标下的场分布，计算公式为：
[0058][0059][0060]
步骤s5，根据oam波束不同模式之间的正交性，计算各个模式的能量强度，计算公式如下：
[0061][0062]
其中n表示计算的扩展模式阶数，i表示虚数符号。单个模式的功率可以表示为：
[0063][0064]
最后可以得到模式扩展场的功率谱，计算公式如下：
[0065][0066]
步骤s6：使用最小二乘法对功率谱与遮挡宽度曲线进行多项式拟合，具体实现过程为：
[0067]
步骤s61：给定函数w＝f(p)，设p为powern的主模，在p1,p2,p3,
…
,pm处的函数值为w1,w2,w3,
…
,wm，其中n表示扩展模式阶数，m表示给定拟合点的个数。设多项式f(p)＝a0 a1p
a2p2
……
akpk，其中k表示多项式的阶数使得
[0068][0069]
为了求得符合条件的系数ai，对上式每个系数求偏导数，即：
[0070][0071]
对以上方程组化简整理后，用矩阵表示为：
[0072][0073]
步骤s62：求解上述矩阵中的各个系数ai，即可得到关于p的多项式拟合函数w(p)。
[0074]
步骤s7：根据函数w(p)实现待检测物体尺寸宽度的快速测量。
[0075]
具体为，根据函数w(p)中，遮挡宽度w与功率谱p的对应关系，实现待检测物体尺寸宽度的快速测量。
[0076]
进一步地，具体为，根据函数w(p)＝a0 a1p a2p2
……
akpk，其中a0,a1,
…
,ak各系数在步骤s62中已给出，再带入功率谱p的值，即可求得对应的遮挡宽度w，实现待检测物体尺寸宽度的快速测量。
[0077]
本发明的实施例具体说明如下：
[0078]
采用本发明提出的基于轨道角动量波束模式扩展的宽度测量方法，具体包括以下步骤：
[0079]
步骤s100，选定所要测量的最大宽度为w
max
，测量距离为s；
[0080]
步骤s200，根据最大测量宽度w
max
以及测量距离s，确定所采用的轨道角动量波束模式(如图2所示)；
[0081]
其中图2表示模式为2的轨道角动量波束；
[0082]
步骤s300，根据电波频率f、测量区间尺寸s、w
max
划分计算网格，设置x、y、z方向的网格步长分别为δx、δy、δz，具体参数设置见表1；
[0083]
序号参数名设定值1f3ghz2s10m3w
max
1m4x0.05m5y0.05m6z1m
[0084]
步骤s400，采用抛物方程法计算各个尺度遮挡下的模式扩散场ed(如图3所示)；
[0085]
其中图3表示遮挡物尺寸为1m时的模式扩散场；
[0086]
步骤s500，将笛卡尔坐标系下的场值分布ed转换成柱坐标下的场分布ec；
[0087]
步骤s600，计算模式扩散场ec的功率谱(如图4所示)；
[0088]
其中图4表示遮挡物尺寸为1m时的模式扩散场，从-5阶到 5阶模式的功率谱；
[0089]
步骤s700，采用最小二乘法的多项式拟合法，对离散数据进行曲线拟合(如图5所示)；
[0090]
其中图5表示若干组离散数据以及拟合曲线，根据离散点与曲线的贴合程度可以看出拟合效果很好；
[0091]
步骤s800，在0～1的宽度范围内随机取20个测试数据(如图6所示)；计算该20个测试宽度主模式的对应功率谱，测试结果如图7；
[0092]
如图7所示的测试结果来看，20个测试数据与曲线上的值吻合非常好，说明该测量方案是可行的，遮挡物的宽度测算结果与真实值的均方根误差为0.0012，说明测量该测量方案精度较高。
[0093]
以上所述仅是本发明的优选实施方式，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、修改和环境，并能够在本文所述构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附权利要求的保护范围内。