MATLAB中非均匀采样速率转均匀采样速率的方法求助

针对非均匀采样转均匀采样时丢失关键特征、高采样率成本高，以及低通滤波无法消除下降沿尖峰的问题，以下是具体的MATLAB实现方案：

一、非均匀转均匀采样：保留关键特征，控制计算成本 #

放弃过度平滑的样条拟合，改用线性插值配合关键特征点强制采样，既保证均匀性，又不丢失下降沿尖峰这类关键信息：

% 假设你的非均匀时间、速度数据为 t_nonuniform, v_nonuniform
% 1. 设定目标均匀采样的时间范围和步长
t_start = t_nonuniform(1);
t_end = t_nonuniform(end);
% 基于原数据平均采样间隔设置步长，平衡精度与计算量
avg_interval = mean(diff(t_nonuniform));
t_uniform = t_start:avg_interval:t_end;

% 2. 线性插值（保留突变特征，避免样条的过度平滑）
v_uniform = interp1(t_nonuniform, v_nonuniform, t_uniform, 'linear', 'extrap');

% 可选优化：强制保留下降沿尖峰点
% 检测速度突变点（一阶差分超过阈值，阈值可根据数据调整）
diff_v = abs(diff(v_nonuniform));
threshold = 0.1 * max(diff_v);
peak_indices = find(diff_v > threshold) + 1;
peak_times = t_nonuniform(peak_indices);
% 合并尖峰时间到均匀采样序列，去重排序
t_uniform_with_peaks = sort(unique([t_uniform, peak_times]));
v_uniform_with_peaks = interp1(t_nonuniform, v_nonuniform, t_uniform_with_peaks, 'linear', 'extrap');

二、移动平均滤波：消除高频尖峰 #

使用移动平均时，窗口大小是关键——太小无法消去尖峰，太大则会模糊正常信号。推荐结合中值滤波先抑制脉冲尖峰，再用移动平均平滑：

% 方法1：基础移动平均
window_size = 5; % 根据采样步长调整，比如对应0.005秒的时间窗口
v_filtered = movmean(v_uniform_with_peaks, window_size);

% 方法2：中值滤波+移动平均（更适合脉冲类尖峰）
% 先3点中值滤波消除孤立尖峰
v_med = medfilt1(v_uniform_with_peaks, 3);
% 再用3点移动平均平滑
v_final = movmean(v_med, 3);

关键注意点 #

• 采样步长：无需盲目提升采样率，以原数据最小采样间隔对应的奈奎斯特频率为上限即可，避免不必要的计算成本
• 窗口大小：如果尖峰是窄脉冲，窗口大小需覆盖尖峰的持续时间，同时不影响正常下降沿的斜率，可通过对比滤波前后的信号调整
• 插值方法：线性插值是平衡特征保留与平滑度的最优选择，若需要更贴近原数据的细节，也可尝试'nearest'最近邻插值，但会引入锯齿

内容的提问来源于stack exchange，提问作者Yuki.kuroshita