[MATLAB]三/五/七/九点二次滑动平均

dwcai

三/五/七/九点二次平滑MATLAB程序，起到低通滤波器的作用，克服了滑动平均削弱过多波幅的缺点。话不多说，直接上代码。如有问题可按程序里的联系方式联系我，希望能帮到你！

[attach]8481[attach]

1. % Name:     runningmean.m
   
2. % Function: 计算滑动平均(气象上的三/五/七/九..点二次平滑)的MATLAB程序(函数)
   
3. % Author:   dwcai
   
4. % Email:    1510336267@qq.com
   
5. % Createtime: 20190724
   
6. % Usage:    [rm] = runningmean(X, n, [endpoint])
   
7. % Version:  V0.1
   
8. % Parameters: rm: 滑动平均序列
   
9. %             X:  数据序列
   
10. %             n:  平滑窗口长度(3, 5, 7, 9)
    
11. %             endpoint: 处理端点['ave','drop'],ave算术平均[默认],drop去掉端点值
    
12. % Reference:
    
13. % https://wenku.baidu.com/view/556d418f3086bceb19e8b8f67c1cfad6185fe953.html
    
14. % Test code:
    
15. %--------------以下为测试代码(复制到此函数相同路径的 m 脚本中)------------------
    
16. % ts = randi(10, 50, 1);
    
17. % r3 = runningmean(ts, 3);   % 默认端点值处理为算数平均法
    
18. % r5 = runningmean(ts, 5, 'drop');   % 设置为抛除端点值
    
19. % r7 = runningmean(ts, 7, 'ave');
    
20. % r9 = runningmean(ts, 9);
    
21. % plot(ts)
    
22. % hold on
    
23. % plot(r3)
    
24. % plot(r5)
    
25. % plot(r7)
    
26. % plot(r9)
    
27. % hold off
    
28. %------------------------------测试结束-----------------------------------
    
29. 
    
30. function [rm] = runningmean(X, n, endpoint)
    
31.    if nargin < 2
    
32.     error('输入参数不足!');
    
33.    elseif nargin == 2
    
34.     endpoint = 'ave';    % 默认首尾采用算术平均法
    
35.    end
    
36.    if min(size(X)) > 1
    
37.     error('不能处理维度>=2的数据序列!');
    
38.    end
    
39.    if max(size(X)) <= n
    
40.     error('输入数据个数必须大于滑动点数!');
    
41.    end
    
42.    if isodd(n) == false
    
43.     error('滑动点数n必须为奇数[3,5,7,9].');
    
44.    else
    
45.     half = (n - 1) / 2;  % 求半区间
    
46.    end
    
47.    rm = ones(size(X))*nan;
    
48.    ln = length(X);
    
49.    % 计算1+half:ln-half区间的 n 点二次滑动平均值
    
50.    startp = 1 + half;
    
51.    endp = ln - half;
    
52.    switch n
    
53.     case 3   % 三点二次平滑
    
54.         for i = startp:endp
    
55.             rm(i) = (X(i-1)+2*X(i)+X(i+1))/4;
    
56.         end
    
57.     case 5   % 五点二次平滑
    
58.         for i = startp:endp
    
59.             rm(i) = (-3*X(i-2)+12*X(i-1)+17*X(i)+12*X(i+1)-3*X(i+2))/35;
    
60.         end
    
61.     case 7   % 七点二次平滑
    
62.         for i = startp:endp
    
63.             rm(i) = (-2*X(i-3)+3*X(i-2)+6*X(i-1)+7*X(i)+6*X(i+1)+3*X(i+2)-2*X(i+3))/21;
    
64.         end
    
65.     case 9   % 九点二次平滑
    
66.         for i = startp:endp
    
67.             rm(i) = (-21*X(i-4)+14*X(i-3)+39*X(i-2)+54*X(i-1)+59*X(i)+54*X(i+1)+39*X(i+2)+14*X(i+3)-21*X(i+4))/231;
    
68.         end
    
69.      %case 11... 待扩展
    
70.    end
    
71.    % 处理边界值
    
72.    switch endpoint
    
73.      case 'ave'      % 算术均值平滑
    
74.          % 前段端点
    
75.          for i = 1:startp
    
76.                rm(i) = mean(X(1:i+1));
    
77.          end
    
78.           % 后端端点
    
79.           for i = endp:ln
    
80.                rm(i) = mean(X(i-1:end));
    
81.           end
    
82.      case 'drop'         % 抛弃端点值, 则返回值的首尾出现 NaN 值
    
83.           warning('端点值被抛弃, 返回值的首尾以NaN填充.');
    
84.           disp('端点值被抛弃, 返回值的首尾以NaN填充.');
    
85.    end
    
86. end
    
87. 
    
88. % 奇数判断函数
    
89. function [odd] = isodd(n)
    
90. if rem(n,2) == 1
    
91.     odd = true;
    
92.     else
    
93.     odd = false;
    
94.     end
    
95. end
    

复制代码

weizhudeng

多谢分享，好东西

fyzhu027

很好，学习下。

zihan6688

很不错，学习一下

小小观测员

{:5_235:}

807727227

感谢分享，最近正好需要用到9点平滑