代码中的参数设置（如噪声协方差矩阵）可以根据实际应用场景进行调整，以适应不同的环境条件。
文章目录
代码介绍
本课题围绕三维空间中的目标定位与轨迹跟踪展开研究，采用基于到达时间（TOA, Time of Arrival）的测距定位方法，并结合扩展卡尔曼滤波（EKF）对目标运动轨迹进行估计与优化。针对实际应用中定位结果易受测量噪声、时钟误差以及环境扰动影响的问题，本文构建了一个包含多个三维锚节点的定位模型，通过对目标与各锚节点之间的传播时间进行测量，完成目标位置的初始估计，并进一步利用 EKF 对连续时刻的位置信息进行滤波处理，从而提升轨迹解算的平滑性与定位精度。
在实现过程中，课题首先建立目标在三维空间中的运动轨迹，并根据锚节点坐标计算其与目标之间的真实距离，随后引入测量噪声生成含误差的 TOA 数据，利用伪逆法求解目标位置。为了进一步抑制观测误差带来的累积影响，程序在定位结果基础上引入扩展卡尔曼滤波，对目标状态进行递推估计与修正。最终通过三维轨迹图、轴向误差曲线以及 RMSE 统计结果，对滤波前后的定位性能进行对比分析，验证了所提方法在动态目标定位场景中的有效性。
本课题具有较强的工程应用价值，可为无线传感网络、室内定位、无人机导航、移动机器人跟踪等场景提供一种可实现的三维定位与轨迹估计方案，同时也为后续研究更复杂环境下的高精度定位算法提供了基础。
运行结果
定位示意图与轨迹输出：
程序结构：
三维误差曲线：
命令行误差特性输出：
源代码
程序结构：
部分源代码如下：
matlab
% TOA测距定位，三维任意（>3)个锚节点，对一个未知点定位、带EKF的轨迹解算
% 2025-03-13/Ver1
clear;clc;close all;
rng(0);
%% 主程序
c = 3e8; %信号传输速度，即光速
range_err = 1e-9; %时钟与时间计算误差
point1 = [1,1,1]; %待求点坐标真值
% 生成目标的运动
positions = repmat(point1,21,1)+[0:0.2:4;0:-0.2:-4;zeros(1,21)]';
for i1 = 1:size(positions,1)
point1 = positions(i1,:);
n = 9; %定义锚节点数量
baseP = 2*[sin(1:n)+0.01*[1:n]+1;cos(4*(1:n))+0.01*[1:n]+1;cos(2*(1:n))+0.01*[1:n]+1]';
R_real = sqrt(diag((point1-baseP)*(point1'-baseP')));
TOA = R_real/c+range_err*randn; %含噪声的传播时间
R_calcu = TOA*c;
p_out(i1,:) = position_3dim(R_calcu,baseP);
end
%% EKF部分
%% 滤波模型初始化课题扩展方向
滤波算法优化 ：
将 EKF 扩展为 UKF、CKF 或粒子滤波（PF），提升强非线性条件下的定位精度与稳定性。
自适应滤波研究 ：
根据环境变化动态调整噪声协方差矩阵，实现复杂环境下的鲁棒定位。
复杂环境误差抑制 ：
针对 NLOS、多径效应等问题，研究鲁棒估计与误差补偿方法。
智能化定位方法 ：
引入深度学习方法，对 TOA 测距误差进行预测与补偿，提升定位精度。
定位性能理论分析 ：
增加 CRLB、GDOP 等理论指标分析，研究系统定位精度上限与锚节点布局影响。
如需帮助，或有导航、定位滤波相关的代码定制需求，请点击下方卡片联系作者