爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8024|回复: 7

[程序设计] matlab 五点三次平滑算法

[复制链接]

新浪微博达人勋

发表于 2013-8-27 13:49:43 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
原文地址:matlab 五点三次平滑算法作者:deman
http://www.ilovematlab.cn/thread-71818-1-1.html
这里提供一个函数mean5_3(五点三次平滑算法)对数据进行平滑处理:

load V1.mat
subplot 211; plot(V1); ylim([2000 7000]); grid;
y=mean5_3(V1,5);
subplot 212; plot(y); ylim([2000 7000]); grid;
function y=mean5_3(x,m)
% x为被处理的数据
% m 为循环次数
n=length(x);
a=x;
for k=1: m
b(1) = (69*a(1) +4*(a(2) +a(4)) -6*a(3) -a(5)) /70;
b(2) = (2* (a(1) +a(5)) +27*a(2) +12*a(3) -8*a(4)) /35;
for j=3:n-2
b (j) = (-3*(a(j-2) +a(j+2)) +12*(a(j-1) +a(j+1)) +17*a(j)) /35;
end
b (n-1) = (2*(a(n) +a(n-4)) +27*a(n-1) +12*a(n-2) -8*a(n-3)) /35;
b (n) = (69*a(n) +4* (a(n-1) +a(n-3)) -6*a(n-2) -a(n-4)) /70;
a=b;
end
y =a;
处理后得图如下,如果还想更平滑,可以修改调用函数时的参数m。

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-8-27 14:11:06 | 显示全部楼层
没有图
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-8-27 14:11:35 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-11-25 14:22:26 | 显示全部楼层
为什么我出不来第二个图?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-2 10:31:29 | 显示全部楼层
想问楼主,只进行一次平滑,效果咋样?我用fortran编了一个,没啥变化
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-12-4 11:53:52 | 显示全部楼层
谢谢楼主分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-12-7 22:43:04 | 显示全部楼层
多谢楼主!正好我现在需要!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-4-19 16:13:57 | 显示全部楼层
多谢楼主分享
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表