一维信号光滑去毛刺处理

2021年11月21日 5点热度 0条评论

  以下是威斯康星大学交通系博士申请遇到的一道题,有兴趣的可以交流

Problem:

 Implement theKalman Filter (Use identity matrix for statepropagation)or other smoothing filter of your choice on the enclosed “speed24.csv” .Compare the filtered and the original data in a time-series plot (x: timeinterval, y: speed data). Report ME, MAE, and MAPE of the filtering results. Ifthe filter contains any parameters, e.g.smoothing window,smoothing coefficient, experiment to find the best fittingparameters that produces the best trend of the original data.( 平滑窗口,平滑参数,试验找到最好的拟合参数以得到原始数据的最好趋势。)

以下是matlab实现的算法

[plain] 
view plain
 copy
 

  1. function xdense = linear_interp(x, t, tdense)  
  2. xdense = zeros(size(tdense));  
  3.   
  4. for i = 1:length(tdense)  
  5.     ind_b = find(tdense(i) >= t, 1, 'last');  
  6.     if isempty(ind_b)  
  7.         ind_b = 1;  
  8.         ind_e = 1;  
  9.         lambda = 1;  
  10.     elseif ind_b == length(t)  
  11.         ind_e = ind_b;  
  12.         lambda = 0;  
  13.     else  
  14.         ind_e = ind_b + 1;  
  15.         lambda = (tdense(i) - t(ind_b)) / (t(ind_e) - t(ind_b));  
  16.     end;  
  17.       
  18.     xdense(i) = (1 - lambda) * x(ind_b) + lambda * x(ind_e);  
  19. end;  

[plain] 
view plain
 copy
 

  1. %--------------------------------------------------------------------------  
  2. % apply fourier analysis on each window size, estimate the dominant  
  3. % frequency, and remove high freqs...  
  4.   
  5. % partition the signals and perform fourier analysis so that we get the  
  6. % dominant freq, then smooth the signal using sigma that changes with the  
  7. % dominant freq...  
  8. %--------------------------------------------------------------------------  
  9.   
  10. clear all;  
  11. %clc;  
  12. %close all;  
  13.   
  14. % 可以调整的参数