把爬行轨迹做分形处理是个很扎实的切入点,不过盒计数法在离散轨迹采样里很容易受网格尺度和边界效应干扰。从实验数据的角度看,你提到D≈1.27,这个数值落在常见室内受限游走区间,但仅凭旧手机IMU数据,采样频率和漂移误差必须先过一遍标定。en fait,把原始加速度二次积分得到的位移序列直接喂给分形算法,累积误差会呈t²放大,维数估算很容易被系统性高估0.1到0.15。
定量实验的第一步永远是控制噪声源。手机MEMS加速度计的零偏不稳定性通常在0.5–2 mg/√Hz量级,如果不做互补滤波融合陀螺仪数据,低频段的1/f^β谱线大概率是传感器本底噪声叠加运动伪影。早年我在实验室做反应物扩散轨迹追踪时,也踩过类似的坑。后来改用激光位移传感器配合高精度编码器,采样提到200Hz以上,再用小波阈值去噪重构相空间,得到的幂律指数才稳定下来。你手头的数据如果原始采样低于50Hz,奈奎斯特频率会直接截断0.3Hz以上的有效信息,这时候做对数坐标下的线性拟合,截断频率的选择就值得商榷了。
你提到混沌吸引子投影,这个建模方向很有意思。拉瓦锡当年建立定量化学体系时,面对看似无序的气体逸散,核心思路也是引入精确的边界条件与质量守恒约束。把轨迹映射到状态空间,如果确实存在低维吸引子,应该计算Lyapunov指数谱或关联维数。单纯看D=1.27还不够,需要观察轨迹在重构空间里的折叠拓扑。插座和猫砂盆作为“势阱”,或许可以用朗之万方程里的多稳态势场来参数化,噪声强度与势垒高度的比值直接决定了状态跃迁频率。这比单纯归因于神经退耦更便于设计可重复的对照实验。
如果有条件,建议用开源的惯性导航算法重处理一遍原始CSV,先把坐标系对齐和重力补偿做扎实。顺便问一句,你用的盒计数是固定网格还是动态自适应?β拟合时的加权最小二乘是怎么设置的?周末打算炖锅红酒炖牛肉,文火收汁的时候正好跑一遍蒙特卡洛模拟,看看能不能复现那个谱线特征。等你更新干净数据集。