做尽调需要用到的网站,网站设计合同注意事项,做全屏轮播的网站有哪些,如何做好互联网营销2021年全国研究生数学建模竞赛华为杯
E题 信号干扰下的超宽带(UWB)精确定位问题
原题再现#xff1a; 一、背景 UWB#xff08;Ultra-Wideband#xff09;技术也被称之为“超宽带”#xff0c;又称之为脉冲无线电技术。这是一种无需任何载波#xff0c;通过发送纳秒…2021年全国研究生数学建模竞赛华为杯
E题 信号干扰下的超宽带(UWB)精确定位问题
原题再现 一、背景 UWBUltra-Wideband技术也被称之为“超宽带”又称之为脉冲无线电技术。这是一种无需任何载波通过发送纳秒级脉冲而完成数据传输的短距离范围内无线通信技术并且信号传输过程中的功耗仅仅有几十µW。UWB因其独有的特点使其在军事、物联网等各个领域都有着广阔的应用。其中基于UWB的定位技术具备实时的室内外精确跟踪能力定位精度高可达到厘米级甚至毫米级定位。UWB在室内精确的定位将会对卫星导航起到一个极好的补充作用可在军事及民用领域有广泛应用比如电力、医疗、化工行业、隧道施工、危险区域管控等。UWB更多应用场景请参见[4—6]。 UWB的定位技术有多种方法本文仅考虑基于飞行时间Time of Flight, TOF的测距原理它是UWB定位法中最常见的定位方法之一。TOF测距技术属于双向测距技术其通过计算信号在两个模块的飞行时间再乘以光速求出两个模块之间的距离这个距离肯定有不同程度 在室内定位的应用中UWB技术可以实现厘米级的定位精度一般指2维平面定位,并具有良好的抗多径干扰和衰弱的性能以及具有较强的穿透能力。但由于室内环境复杂多变UWB 通信信号极易受到遮挡虽然UWB技术具有穿透能力但仍然会产生误差在较强干扰时数据会发生异常波动通常是时间延时基本无法完成室内定位甚至会造成严重事故。因此信号干扰下的超宽带UWB精确定位问题成为亟待解决的问题。 二、问题描述 为解决信号干扰下的超宽带UWB精确定位问题我们通过实际场景实测采集到一定数量的数据即利用UWB的定位技术TOF采集到锚点anchor与靶点Tag之间的距离希望通过数学建模或算法方法 无论信号是否干扰都可以给出目标物靶点的精确定位3维坐标。 三、实验场景和数据采集 如图所示在5000mm5000mm3000mm的测试环境中分别在4个角落A0A1A2A3放置UWB锚点anchor锚点向所有方向发送信号。Tag是UWB标签靶点即需要定位的目标只在测试环境范围内。Tag接收到4个UWB锚点anchor的信号无论信号是否干扰Tag一般都可以接收到信号利用TOF技术分别解算出对应的4个距离数据。 实验在实验场景1中采集了Tag在324个不同位置在信号无干扰和信号干扰下的UWB数据即每个位置各测试采集2次一次信号无干扰另一次信号有干扰锚点与靶点间有遮挡注意每次采集数据时由于Tag在同一位置会停留一会儿时间而锚点与Tag之间每0.2—0.3秒之间就会发送、接收信号一次所以在同一位置点UWB会采集到多组数据多组数据都代表同一位置的信息组数的多少视Tag在同一位置的时间而定停留的时间越长组数就越多。数据见文件夹“附件1UWB数据集”。 实验场景1 靶点Tag范围5000mm5000mm3000mm 锚点anchor位置单位mm A0001300、 A1500001700、 A2050001700、A3500050001300 四、数据文件说明 1UWB数据集 “附件1UWB数据集”有2个文件夹和1个文件1个文件Tag坐标信息.txt存放324个不同位置的编号及3维坐标信息2个文件夹中1个存放信号无干扰下正常采集的数据各文件名为x.正常.txtx表示对应的位置编号另1个存放信号有干扰下异常采集的数据各文件名为x.异常.txtx表示对应的位置编号。 2数据文件 Tag在每个位置都采集了2个数据文件1个正常另1个异常共有648个数据文件无论正常、异常数据数据格式都一样每个数据文件开头第1行为采集开始行无实际意义接下来每4行为一组表示UWB采集的一组完整数据一组数据表示一个样品如 T:144235622:RR:0:0:950:950:118:1910 T:144235622:RR:0:1:2630:2630:118:1910 T:144235622:RR:0:2:5120:5120:118:1910 T:144235622:RR:0:3:5770:5770:118:1910 这4行数据的含义分别是: Tag标识时间戳Range Report的缩写Tag ID锚点ID该锚点的测距值(mm)测距值的校验值数据序列号数据编号每个数据之间用“”分隔。实际上就是提供了4个锚点到靶点Tag的距离即 A0到靶点距离为950mm A1到靶点距离为2630mm A2到靶点距离为5120mm A3到靶点距离为5770mm 每个数据文件都有多组数据表示在同一位置连续时间内UWB自动采集到的多组数据。 五、完成任务 试根据上述数据完成如下任务 任务1数据预处理清洗 无论是信号无干扰下采集数据或信号干扰下采集数据Tag在同一坐标点上都采集多组数据见附件1中648个数据文件请用某种方法把每个数据文件相应数值抓取出来并转换成二维表矩阵形式txt、Excel或其他数据格式每一行代表一组数据即一个样品然后对这些数据文件进行预处理清洗删除掉一些“无用”异常、缺失、相同或相似的数据样品。经处理后“正常数据”所有数据文件和“异常数据”所有数据文件最后各保留多少组多少个样品数据并重点列出以下4个数据文件经处理后保留的数据矩阵形式 “正常数据”文件夹中 24.正常.txt、 109.正常.txt “异常数据”文件夹中 1.异常.txt、 100.异常.txt 任务2: 定位模型 利用任务1处理后的数据分别对“正常数据”和“异常数据”设计合适的数学模型或算法估计或预测出Tag的精确位置并说明你所建立的定位模型或算法的有效性同时请利用你的定位模型或算法分别对附件2中提供的前5组信号无干扰数据和后5组信号有干扰数据进行精确定位3维坐标 注意1定位模型必须体现实验场景信息 2请同时给出定位模型的3维xyz精度、2维xy精度以及1维的各自精度。 任务3不同场景应用 我们的训练数据仅采集于同一实验场景实验场景1但定位模型应该能够在不同实际场景上使用我们希望你所建立的定位模型能够应用于不同场景。附件3中10组数据采集于下面实验场景2前5组数据信号无干扰后5组数据信号有干扰请分别用上述建立的定位模型对这10 组数据进行精确定位3维坐标 实验场景2 靶点Tag范围5000mm3000mm3000mm 锚点anchor位置单位mm A0001200、 A1500001600、 A2030001600A3500030001200 任务4: 分类模型 上述定位模型是在已知信号有、无干扰的条件下建立的但UWB在采集数据时并不知道信号有无干扰所以判断信号有无干扰是UWB精确定位问题的重点和难点。利用任务1处理后的数据建立数学模型(或算法)以便区分哪些数据是在信号无干扰下采集的数据哪些数据是在信号干扰下采集的数据并说明你所建立的分类模型(或算法)的有效性同时请用你所建立的分类模型或算法判断附件4中提供的10组数据这10组数据同样采集于实验场景1是来自信号无干扰或信号干扰下采集的 任务5: 运动轨迹定位 运动轨迹定位是UWB重要应用之一利用静态点的定位模型加上靶点自身运动规律希望给出动态靶点的运动轨迹。附件5是对动态靶点采集的数据一段时间内连续采集的多组数据请注意在采集这些数据时会随机出现信号干扰请对这个运动轨迹进行精确定位最终画出这条运动轨迹图数据采集来自实验场景1。
整体求解过程概述(摘要) UWB 定位技术是一类室内外精确跟踪能力强定位精度高的导航技术。近年来在各领域得到了快速发展和广泛应用。然而基于UWB的三维高精度定位算法仍然是一项亟待解决的技术问题。本文将基于这样的背景采用提供的试验数据研究一类高精度的分类与定位算法并应用提供的数据完成了方法检验和技术验证。 针对问题1中的数据预处理问题本文统计了样本数据的分布特性得出了大部分样本数据符合正态分布的结论对于样本中测量值异常的数据采用 3σ 原则对于偏离分布中心3倍标准差的数据进行删除考虑到相同相似样本会加大定位算法设计与计算负担采用欧氏距离判定同一样本中相同相似数据并进行删除进一步给出了面向定位算法设计的简化数据。同时应用统计分析方法检验了预处理前后数据分布特性确定了测量数据存在有效偏差和白噪声干扰统计了3附件一样本对应的324项标准差数据得到了标准差符合正态分布且均值为10标准差为3的结论确定了各组传感器性能一致且明确了测量噪声的分布规律。 针对问题二与问题三首先分别应用传统的最小二乘法、基于PSO的定位算法建立了接近传统意义上的四点算法精度分析表明两种算法的求解精度基本一致基本符合xy 轴求解精度较高而z轴求解精度较低的特性。为了进一步建立高精度定位算法采用机器学习方法建立了一类高精度定位算法以测量长度为输入、xyz 坐标数据为输出设计了三层BP 神经网络通过大量训练完成了定位算法设计精度分析表明基于机器学习的定位算法仍然表现出xy轴求解精度较高而z轴求解精度较低的特性但是x、y轴求解精度达到厘米级而z轴求解精度达到了分米级极大地提高了三点定位算法的精度特性。 针对问题四为了进一步解决正常数据、异常数据的分类问题首先基于传统方法出发拟寻找一类分类指标完成对正常、异常数据的二分类。然而基于统计分析建立的测距偏差分类方法和基于平面几何出发的三圆相交面积判别皆未能充分区分正常、异常数据因此本文认为基于有限的先验信息和较大的噪声水平干扰下无法通过有效的传统方法完成对正常数据、异常数据的识别因此仍采用机器学习方法建立了一类有效的分类方法基于附件1的测试表明识别精度达到97.69%并完了对附件四数据的计算。 针对问题五首先从附件5中提取得到了运动过程中靶点到4个锚点的长度信息再使用问题四中的分类算法验证每组数据是否存在异常干扰对于修正后的数据采用前面问题所建立的机器学习的方法求解坐标最后根据坐标即可得到在空间中的运动轨迹。 本文综合了针对UWB技术三维坐标精确定位技术出发展开了相关研究分析了传统方法的不足并应用机器学习方法形成了一套准确性高、适用性强的分类和定位方法有效的解决了题目中的相关问题。
模型假设 假设测量误差仅仅因TOF技术不足和传感器性能有限而引起误差引起的形式为一类有色偏差和一类白噪声误差的组合。 基于数据统计分析可以说明UWO 给出的测距偏差与实际测量距离存在线性相关关系随着实际测量距离的增加UWO 测量结果会被引入一类偏置误差同时由于传感器性能不同和实验引入的误差测量结果还会引入一类白噪声误差假定白噪声幅值与传感器性能有关且假设传感器的性能并不随实验次数和实验环境变化。
问题分析(部分内容由于篇幅限制) 通过数据处理得到了附件1中的全数据考虑题中规定附件1中数据来源于324个不同的位置测量获得首先通过标准差分析检验数据的分布特点。下文为了叙述和绘图方便约定ZC表示正常YC表示异常。 原始数据标准差分析表明正常数据 A0ZCA3ZC 的标准差存在尖端A1ZC、A2ZC的标准差控制在20以内异常数据中每个锚点都存在标准差处于高水平(超过200)的状态选取两幅典型数据解释该问题这是因为异常数据中同一时间内仅有一个锚点存在异常状态但是同一组数据内先后有两个锚点被干扰干扰对测量数据造成的影响体现为两种情况在原有数据基础上引入了超过400的测量偏差或者引入一段均值为400左右的大幅值(波动水平在400上下有幅值为200的波动)波动误差。这里要指出对于异常数据的定义为在正常数据水平上出现了正偏置干扰测量数据原有的测量噪声与正常数据中的测量噪声一致看待。 题中对正常数据的定义为不存在测量干扰的数据但是通过标准差分析发现A0ZCA3ZC 的标准差存在尖端。提取出部分数据并做图分析可知这是因为原始数据中存在异常点具体表现为图像出现尖端。 首先考虑通过几何关系求解靶点的坐标。在实验场景 1 中靶点Tag范围5000mm5000mm3000mm如图5.1 建立空间直角坐标系其中锚点位置如图所示分别为A0001300、A1500001700、A2050001700、A3500050001300而靶点可以出现在5000mm5000mm3000mm立方体空间中任意一处根据靶点和锚点的空间关系符合公式 式中xyz为靶点坐标xi、yi、zi为锚点的坐标i为锚点的编号i1,2,3,4。空间问题需要至少4个锚点来确定坐标而上述方程为三元二次方程组。在理想情况下三个锚点坐标信息以及这三个锚点到靶点的距离即可求解方程得到两组解通过第4个锚点的信息确定靶点的最终坐标在实际情况下由于测距数据存在误差满足4个方程的解可能无实数解故不能用此方法直接来求解靶点的坐标。 由于通过锚点位置信息和锚点测距值无法求解直接求解靶点坐标故利用锚点坐标和锚点的测距值设计一个期望函数通过函数值来反映靶点的坐标状况。
论文缩略图 全部论文及程序请见下方“ 只会建模 QQ名片” 点击QQ名片即可
程序代码
import osimport reimport numpy as npimport pandas as pdfrom os import listdirfrom collections import Counterdef missing_data_processing(content,j):a []time re.compile(rT:(\w):RR).findall(content)b list(Counter(time))for i in b:if (Counter(time)[i] 4):a.append(i)if (int(j)109):print(len(a))return adef three_sigmal(data):n 3 #n*sigmaprint(data)ymean np.mean(data, axis0) # 均值
ystd np.std(data, axis0)threshold1 ymean- n * ystdthreshold2 ymean n * ystdoutlier [] # 将异常值保存
final []# 标准差
# 删除掉异常值后的数据
for i in range(0, len(data)):#print(np.abs(data[i]-ymean))if (np.abs(data[i]-ymean) n * ystd).any() :outlier.append(data[i])else:final.append(data[i])return outlier#normal_data_path D:/竞赛/数学建模/2021 年 E 题/附件 1UWB数据集/正常数据abnormal_data_path D:/竞赛/数学建模/2021 年 E 题/附件 1UWB 数据集/异常数
据label_path D:/竞赛/数学建模/2021 年 E 题/附件1UWB数据集/Tag坐标信息.txtfilename1 listdir(normal_data_path)filename2 listdir(abnormal_data_path)1. import matplotlib.pyplot as plt2. import numpy as np3. import pandas as pd4. import math5. from sklearn import datasets, linear_model, discriminant_analysis6. from sklearn.model_selection import train_test_split7. from sklearn.preprocessing import StandardScaler #数据预处理8. import csv9. import joblib10.11. def test_ridge(*data):12. X_train, X_test, y_train, y_test data13. ridgeRegression linear_model.Ridge()14. # ridgeRegression linear_model.LinearRegression()15. rdridgeRegression.fit(X_train, y_train)16. joblib.dump(rd, yc_model.pkl)17.18. #print(预测性能得分: %.2f % ridgeRegression.score(X_test, y_test))19. return ridgeRegression.predict(X_test)20.21. #测试不同的α值对预测性能的影响
22. def test_ridge_alpha(*data):23. X_train, X_test, y_train, y_test data24. alphas [0.01, 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000]25. scores []26. for i, alpha in enumerate(alphas):27. ridgeRegression linear_model.Ridge(alphaalpha)28. ridgeRegression.fit(X_train, y_train)29. #scores.append(ridgeRegression.score(X_test, y_test))30. return alphas, scores31.32. def show_plot(alphas, scores):33. figure plt.figure()34. ax figure.add_subplot(1, 1, 1)35. ax.plot(alphas, scores)36. ax.set_xlabel(r$\alpha$)37. ax.set_ylabel(rscore)38. ax.set_xscale(log)39. ax.set_title(Ridge)40. plt.show()41.42.43.
44. if __name__ __main__:45.46. A0[0,0,130]47. A1[500,0,170]48. A2[0,500,170]49. A3[500,500,130]50.51. url ./预测/16维训练异常.csv52. data pd. read_csv(url, sep,,headerNone)53. datanp.array(data)54.55. X_datadata[:,4:16]56.57. Ydata[:,16:19]58.59. # url_test ./预测/场景1/场景1_异常.csv60.61. # url_test ./预测/16维训练异常.csv62. url_test ./预测/16维测试异常.csv63. # url_test ./预测/场景2/场景2_异常.csv64. data_test pd. read_csv(url_test, sep,,headerNone)65. data_testnp.array(data_test)66.67. X_data_testdata_test[:,4:16]68. Y_testdata_test[:,16:19]69.70. X_trainX_data71. X_testX_data_test72. Y_trainY73. Y_testY_test74. print(Y_test)75.76. predicttest_ridge(X_train, X_test, Y_train, Y_test)77. print(predict)78.79. dis((predict[:,0]-Y_test[:,0])**2(predict[:,1]-Y_test[:,1])**2(predict[:,2]-Y_test[:,2])**2)80. dis_meannp.sqrt(dis)81. print(np.mean(dis_mean))82.83.84. dis((predict[:,0]-Y_test[:,0])**2(predict[:,1]-Y_test[:,1])**2)85. dis_meannp.sqrt(dis)86. print(np.mean(dis_mean))全部论文及程序请见下方“ 只会建模 QQ名片” 点击QQ名片即可 文章转载自： http://www.morning.rycd.cn.gov.cn.rycd.cn http://www.morning.jpbpc.cn.gov.cn.jpbpc.cn http://www.morning.hrzhg.cn.gov.cn.hrzhg.cn http://www.morning.kgsws.cn.gov.cn.kgsws.cn http://www.morning.pltbd.cn.gov.cn.pltbd.cn http://www.morning.nxdqz.cn.gov.cn.nxdqz.cn http://www.morning.slkqd.cn.gov.cn.slkqd.cn http://www.morning.zfgh.cn.gov.cn.zfgh.cn http://www.morning.nqlx.cn.gov.cn.nqlx.cn http://www.morning.yrddl.cn.gov.cn.yrddl.cn http://www.morning.rnfwx.cn.gov.cn.rnfwx.cn http://www.morning.xprq.cn.gov.cn.xprq.cn http://www.morning.qxwgx.cn.gov.cn.qxwgx.cn http://www.morning.mypxm.com.gov.cn.mypxm.com http://www.morning.jcxyq.cn.gov.cn.jcxyq.cn http://www.morning.trtdg.cn.gov.cn.trtdg.cn http://www.morning.wdhlc.cn.gov.cn.wdhlc.cn http://www.morning.gpcy.cn.gov.cn.gpcy.cn http://www.morning.wrtw.cn.gov.cn.wrtw.cn http://www.morning.bbmx.cn.gov.cn.bbmx.cn http://www.morning.wpwyx.cn.gov.cn.wpwyx.cn http://www.morning.hhfwj.cn.gov.cn.hhfwj.cn http://www.morning.mhnb.cn.gov.cn.mhnb.cn http://www.morning.egmux.cn.gov.cn.egmux.cn http://www.morning.zqzzn.cn.gov.cn.zqzzn.cn http://www.morning.jrsgs.cn.gov.cn.jrsgs.cn http://www.morning.yzxlkj.com.gov.cn.yzxlkj.com http://www.morning.bmfqg.cn.gov.cn.bmfqg.cn http://www.morning.wschl.cn.gov.cn.wschl.cn http://www.morning.qtsks.cn.gov.cn.qtsks.cn http://www.morning.thntp.cn.gov.cn.thntp.cn http://www.morning.ycwym.cn.gov.cn.ycwym.cn http://www.morning.hjssh.cn.gov.cn.hjssh.cn http://www.morning.thxfn.cn.gov.cn.thxfn.cn http://www.morning.rzmzm.cn.gov.cn.rzmzm.cn http://www.morning.qnklx.cn.gov.cn.qnklx.cn http://www.morning.rkqqf.cn.gov.cn.rkqqf.cn http://www.morning.kwnbd.cn.gov.cn.kwnbd.cn http://www.morning.xqgfy.cn.gov.cn.xqgfy.cn http://www.morning.qmsbr.cn.gov.cn.qmsbr.cn http://www.morning.nlkhr.cn.gov.cn.nlkhr.cn http://www.morning.mhrzd.cn.gov.cn.mhrzd.cn http://www.morning.nkllb.cn.gov.cn.nkllb.cn http://www.morning.rchsr.cn.gov.cn.rchsr.cn http://www.morning.ruyuaixuexi.com.gov.cn.ruyuaixuexi.com http://www.morning.rgmls.cn.gov.cn.rgmls.cn http://www.morning.qbfwb.cn.gov.cn.qbfwb.cn http://www.morning.rjmb.cn.gov.cn.rjmb.cn http://www.morning.mlffg.cn.gov.cn.mlffg.cn http://www.morning.wspyb.cn.gov.cn.wspyb.cn http://www.morning.mbmh.cn.gov.cn.mbmh.cn http://www.morning.nqmdc.cn.gov.cn.nqmdc.cn http://www.morning.xnltz.cn.gov.cn.xnltz.cn http://www.morning.myxps.cn.gov.cn.myxps.cn http://www.morning.hrtfz.cn.gov.cn.hrtfz.cn http://www.morning.brlcj.cn.gov.cn.brlcj.cn http://www.morning.bxsgl.cn.gov.cn.bxsgl.cn http://www.morning.cbchz.cn.gov.cn.cbchz.cn http://www.morning.kaweilu.com.gov.cn.kaweilu.com http://www.morning.kphsp.cn.gov.cn.kphsp.cn http://www.morning.gpryk.cn.gov.cn.gpryk.cn http://www.morning.tmcmj.cn.gov.cn.tmcmj.cn http://www.morning.ygkk.cn.gov.cn.ygkk.cn http://www.morning.dqxnd.cn.gov.cn.dqxnd.cn http://www.morning.xfcjs.cn.gov.cn.xfcjs.cn http://www.morning.dpbgw.cn.gov.cn.dpbgw.cn http://www.morning.jcyyh.cn.gov.cn.jcyyh.cn http://www.morning.ywzqk.cn.gov.cn.ywzqk.cn http://www.morning.tfcwj.cn.gov.cn.tfcwj.cn http://www.morning.bkgfp.cn.gov.cn.bkgfp.cn http://www.morning.wrtxk.cn.gov.cn.wrtxk.cn http://www.morning.rggky.cn.gov.cn.rggky.cn http://www.morning.gzzxlp.com.gov.cn.gzzxlp.com http://www.morning.fbqr.cn.gov.cn.fbqr.cn http://www.morning.trfh.cn.gov.cn.trfh.cn http://www.morning.kcfnp.cn.gov.cn.kcfnp.cn http://www.morning.dxzcr.cn.gov.cn.dxzcr.cn http://www.morning.tktcr.cn.gov.cn.tktcr.cn http://www.morning.mngh.cn.gov.cn.mngh.cn http://www.morning.ntqqm.cn.gov.cn.ntqqm.cn
