基于子集带符号残差检验的非视距鉴定方法、装置和介质制造方法及图纸

本发明专利技术公开基于子集带符号残差检验的非视距鉴定方法、装置和介质。本发明专利技术的包括对定位测量数据进行分组、对分组数据进行中间估计、利用中间估计值计算完整的初始子集带符号残差矩阵、基于初始子集带符号残差矩阵检测非视距候选假设、从候选假设中确立最终的非视距基站组合。本发明专利技术利用子集带符号残差(剩余)矩阵，能够有效降低前端估计对残差值的影响,提升了现有算法的鉴定准确率；本发明专利技术提出的ML检验，既统一了多个候选假设情形下的选取流程和准则，又能依据严格计算得出的最大似然度量进行抉择，提高了候选假设的选取精准度。提高了候选假设的选取精准度。提高了候选假设的选取精准度。

【技术实现步骤摘要】
基于子集带符号残差检验的非视距鉴定方法、装置和介质


[0001]本专利技术属于无线电定位
，具体涉及基于子集带符号残差检验的非视距鉴定方法、装置和介质。

技术介绍

[0002]非视距(Non
‑
Line
‑
Of
‑
Sight，称NLOS)现象常见于无线电定位环境中，信号传播由于障碍物的遮挡，使得实际的传播时延大于真实距离所需时延，从而引发较大的定位误差，对蜂窝系统而言，这种现象尤其高发于城市小区中，NLOS问题已成为限制定位精度进一步提升的瓶颈。
[0003]Chan和Tsui等基于所设计的视距(Line
‑
Of
‑
Sight，称LOS)环境近似最大似然估计法(即AML法)，提出了利用AML中间估计位置残差概率分布特性进行非视距(Non
‑
Line
‑
Of
‑
Sight，称NLOS)基站数据鉴定的方法(即Chan法)。当所有基站测量值都为LOS时，文中所定义的归一化残差为中心的χ2分布，而如果含有NLOS测量值，则该归一化残差服从非中心的χ2分布。通过设定适当的χ2分布检验阈值，挑选符合要求的候选假设(基站下标集合的某子集，对应为某个假定的LOS基站组合)。当候选假设不止一个时，以检验量超出阈值范围最少数量对应的假设确立最佳LOS基站组合。当基站数目处于系统最小要求值时，文中还提出了δ检验法，用以提高算法的适应性。
[0004]与Chan法处理过程类似，刘霞提出了一种改进的NLOS鉴定方法(即Liu法)，仍然基于AML中间估计下归一化残差的χ2分布检验，利用重新定义的距离残差，Liu法不需要对每种假设的子假设作检验，且不用反复计算近似克拉美劳下界(CRLB)等参数，因而大幅降低Chan法的复杂度。并且，距离残差仅依赖中间估计，其精度相较位置残差更优。Liu法在处理最小数目基站情形时，仍沿用了Chan法中的δ检验法。
[0005]Chan法复杂度较高，应用范围都非常有限，Liu法的复杂度有所降低，但因仍基于χ2分布，受限于归一化独立高斯变量的形式，计算上有较多的冗余；另外，不论Chan法和Liu法，其最小数目基站下的检验均为δ检验，而δ检验使得算法性能在最小数目基站情形下降幅很大，因而为改善此时的性能，更加简单高效地统一鉴定方法也需要进一步探索。

技术实现思路

[0006]本专利技术的目的在于提升了非视距鉴定方法的鉴定性能和适用范围，统一候选假设择优和最小数目基站下的鉴定方法，改善最小基站数目下检验法的鉴定效果，提高非视距鉴定的精度。
[0007]为实现上述目的，本专利技术提供了基于子集带符号残差检验的非视距鉴定方法，包括：
[0008]步骤S100,对定位测量数据进行分组的步骤；
[0009]该步骤中，从M个定位测量数据中任取其中大于等于D
min
个的定位测量数据作为一个定位测量数据组；所述M个定位测量数据是参与移动站定位且位置已知的基站接收来自
移动站的同一信号而获得的数据，D
min
表示计算移动站位置所需定位测量数据的最少个数，
[0010]其中，M≥D
min
，第m个基站的位置坐标记为x
m
＝(x
m
,y
m
)
T
，1≤m≤M，第m个基站获得的定位测量数据记为定位测量数据中含有均值为零的高斯随机噪声，记随机噪声向量n＝[n
1 n2…
n
M
]T
，n
m
为第m个定位测量数据中的随机噪声，随机噪声向量的协方差矩阵为Q；定位测量数据组中定位测量数据下标的集合记为集合S
k
，共有N个满足该条件的下标集合S
k
，即1≤k≤N，按高维到低维的顺序排列，依次为包含所有定位测量数据下标的集合S1＝{1,2,...,M}，包含M
‑
1个所有定位测量数据下标的集合S2＝{1,2,...,M
‑
1}，S3＝{1,2,...,M
‑
2,M}，
…
，依次类推，最后是包含D
min
个定位测量数据下标的集合为S
N
＝{1,2,...,D
min
}；
[0011]步骤S200，按照分组用每个分组对应的定位测量数据进行中间估计的步骤；
[0012]该步骤中，对每个集合S
k
所对应的定位测量数据组，分别进行单次视距估计，得到N个中间估计值N个中间估计值其中，1≤k≤N，LOS(*)表示依据定位测量数据组*进行单次视距估计；
[0013]步骤S300，利用得到的中间估计值计算完整的初始子集带符号残差矩阵E0的步骤；
[0014]该步骤中，根据每个集合S
k
的中间估计值求其M维子集带符号残差向量该向量由M个元素构成，即：
[0015][0016]其中，第m个元素定义为：
[0017][0018]其中，1≤m≤M，m∈S
k
表示集合S
k
中元素m对应的定位测量数据的带符号残差；
[0019]用上述N个按行组合构成一个N
×
M维的矩阵，即初始子集带符号残差矩阵E0为：
[0020][0021]步骤S400，基于初始子集带符号残差矩阵E0检测非视距候选假设的步骤；
[0022]该步骤中，按照从低维到高维的顺序检测当前维度D上是否存在满足非视距检验规则的集合H
k
；其中，集合H
k
为集合S
k
在全集{1,2,...,M}中的补集，即H2＝{M}，H3＝{M
‑
1}，
…
，H
N
＝{M,M
‑
1,...,D
min
+1}；包括步骤：
[0023]步骤H
D
，对于当前维度D,0≤D＜M
‑
D
min
，基于初始子集带符号残差矩阵E0检测个
D维集合H
k
是否为满足非视距检验规则的集合，若D维集合H
k
是满足非视距检验规则的集合，则将D维集合H
k
标记为非视距候选假设，并且在检测完当前维度所有集合之后转入步骤S500，若不存在满足非视距检验规则的集合，则升高一个维度重复上述检测，若直到D＝M
‑
D
min
‑
1时仍不存在满足非视距检验规则的集合，则进入步骤
[0024]所述非视距检验规则为：根据给定的置信水平α以及由置信水平α确定的阈值λ
m
，或者根据给定的阈值λ
m
以及由阈值λ
m
确定的置信水平α，0＜α＜1，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于子集带符号残差检验的非视距鉴定方法，其特征在于包括：步骤S100,对定位测量数据进行分组的步骤；该步骤中，从M个定位测量数据中任取其中大于等于D
min
个的定位测量数据作为一个定位测量数据组；所述M个定位测量数据是参与移动站定位且位置已知的基站接收来自移动站的同一信号而获得的数据，D
min
表示计算移动站位置所需定位测量数据的最少个数，其中，M≥D
min
，第m个基站的位置坐标记为x
m
＝(x
m
,y
m
)
T
，1≤m≤M，第m个基站获得的定位测量数据记为定位测量数据中含有均值为零的高斯随机噪声，记随机噪声向量n＝[n
1 n2ꢀ…ꢀ
n
M
]
T
，n
m
为第m个定位测量数据中的随机噪声，随机噪声向量的协方差矩阵为Q；定位测量数据组中定位测量数据下标的集合记为集合S
k
，共有N个满足该条件的下标集合S
k
，即1≤k≤N，按高维到低维的顺序排列，依次为包含所有定位测量数据下标的集合S1＝{1,2,...,M}，包含M
‑
1个定位测量数据下标的集合S2＝{1,2,...,M
‑
1}，S3＝{1,2,...,M
‑
2,M}，
…
，依此类推，最后是包含D
min
个定位测量数据下标的集合为S
N
＝{1,2,...,D
min
}；步骤S200，按照分组用每个分组对应的定位测量数据进行中间估计的步骤；该步骤中，对每个集合S
k
所对应的定位测量数据组，分别进行单次视距估计，得到N个中间估计值间估计值其中，1≤k≤N，LOS(*)表示依据定位测量数据组*进行单次视距估计；步骤S300，利用得到的中间估计值计算完整的初始子集带符号残差矩阵E0的步骤；该步骤中，根据每个集合S
k
的中间估计值求其M维子集带符号残差向量该向量由M个元素构成，即：其中，第m个元素定义为：其中，1≤m≤M，m∈S
k
表示集合S
k
中元素m对应的定位测量数据的带符号残差；用上述N个按行组合构成一个N
×
M维的矩阵，即初始子集带符号残差矩阵E0为：步骤S400，基于初始子集带符号残差矩阵E0检测非视距候选假设的步骤；该步骤中，按照从低维到高维的顺序检测当前维度D上是否存在满足非视距检验规则的集合H
k
；其中，集合H
k
为集合S
k
在全集{1,2,...,M}中的补集，即H2＝{M}，H3＝{M
‑
1}，
…
，H
N
＝{M,M
‑
1,...,D
min
+1}；包括步骤：
步骤H
D
，对于当前维度D,0≤D＜M
‑
D
min
，基于初始子集带符号残差矩阵E0检测个D维集合H
k
是否为满足非视距检验规则的集合，若D维集合H
k
是满足非视距检验规则的集合，则将D维集合H
k
标记为非视距候选假设，并且在检测完当前维度所有集合之后转入步骤S500，若不存在满足非视距检验规则的集合，则升高一个维度重复上述检测，若检测完D＝M
‑
D
min
+1维度仍不存在满足非视距检验规则的集合，则进入步骤所述非视距检验规则为：根据给定的置信水平α以及由置信水平α确定的阈值λ
m
，或者根据给定的阈值λ
m
以及由阈值λ
m
确定的置信水平α，0＜α＜1，λ
m
＞0，记落在阈值区间ρ
m
＝[
‑
λ
m
,λ
m
]之外的带符号残差剩余矩阵E
k
元素个数为l
k
，若满足l
k
≤floor((1
‑
α)
·
(M
‑
D)
·
N
k
)，则判定集合H
k
满足非视距假设，即在置信水平α下集合H
k
对应的所有定位测量数据为非视距；所述floor(*)表示对*的向下取整，所述N
k
为带符号残差剩余矩阵E
k
的行数，所述带符号残差剩余矩阵E
k
按如下步骤得到：步骤E1，从初始子集带符号残差矩阵E0中标识出第m个定位测量数据参与中间估计的所有行号，构成集合W
m
,1≤m≤M；步骤E2，根据集合H
k
中的D个元素h1，h2，
…
，h
D
，逐一查找对应的集合而后求其并集步骤E3，从初始子集带符号残差矩阵E0中删除所有并集V
k
中的行，即删除行号a∈V
k
的行，剩余行列构成一个子集带符号残差剩余矩阵E
k
；所述阈值λ
m
近似满足下式，近似满足下式，其中，f
m
(n
m
)为第m个定位测量数据中单一维度高斯噪声的概率密度函数，q
mm
为随机噪声向量n＝[n
1 n2ꢀ…ꢀ
n
M
]
T
的协方差矩阵Q中的第m个主对角线元素；步骤当前维度D＝M
‑
D
min
，将个M
‑
D
min
维集合H
k
全部标记为非视距候选假设，转入步骤S500；步骤S500，从非视距候选假设中确定最终非视距基站组合的步骤；该步骤中，若步骤S400检测出仅有一个非视距候选假设，则输出该非视距候选假设下标对应的基站为最终的非视距基站组合；若存在多个非视距候选假设，则针对所有非视距候选假设的集合H
k
，根据其对应的子集带符号残差剩余矩阵E
k
，计算其对应的似然度量L
k
，其中，n
u
为子集带符号残差剩余矩阵E
k
第u行的非零元素所构成的列向量，代表子集带符号残差剩余矩阵E
k
第u行所对应的基站组合在去除非视距基站影响后的带符号残差值向量，也即近似的剩余噪声向量；Q
u
为去除非视距数据所对应的行与列后，依n
u
中数据顺序整理而成的新的协方差矩阵；A
u
为n
u
的向量维度，也即的Q
u
行数或列数，选取似然度量L
k
最大
的非视距候选假设，将对应的集合H
k
中对应的基站作为最终的非视距基站组合，检验结束。2.如权利要求1所述的基于子集带符号残差检验的非视距鉴定方法，其特征在于，所述定位测量数据是TOA数据或者TDOA数据，并且D
min
＝3。3.如权利要求1所述的基于子集带符号残差检验的非视距鉴定方法，其特征在于，随机噪声向量为零均值独立的高斯白噪声时，随机噪声向量的协方差矩阵Q退化为对角矩阵，所有非对角元素均取值为零，即步骤S300中还对子集带符号残差向量进行归一化，归一化后的子集带符号残差向量为步骤S400采用归一化后的子集带符号残差向量进行计算，通过查表得到不同置信水平α下的阈值λ
m
取值，或者根据不同的阈值区间ρ
m
＝[
‑
λ
m
,λ
m
]的阈值λ
m
直接查出对应的置信水平α。4.如权利要求3所述的基于子集带符号残差检验的非视距鉴定方法，其特征在于，步骤S500中，求似然度量L
k
最大值采用求子集带符号残差剩余矩阵E
k
的各元素平方残差和J
k
最小值替代，即：选取平方残差和J
k
最小的非视距候选假设，将对应的集合H
k
中对应的基站作为最终的非视距基站组合。5.一种基于子集带符号残差检验的非视距鉴定装置，其特征在于包括定位测量数据分组模块、中间估计模块、初始子集带符号残差矩阵计算模块、非视距候选假设检测模块和最终非视距基站组合确定模块；其中，定位测量数据分组模块，用于从M个定位测量数据中任取其中大于等于D
min
个的定位测量...

【专利技术属性】
技术研发人员：江成能，王红，江新超，房利达，王松，
申请(专利权)人：广州红达投资有限公司，
类型：发明
国别省市：