zhoufcumt的个人博客分享 http://blog.sciencenet.cn/u/zhoufcumt

博文

修正致命bug:Fortran 90 实现高斯权函数平滑法

已有 5497 次阅读 2013-11-2 10:19 |个人分类:Fortran|系统分类:科研笔记| 高斯权函数平滑法

  非常抱歉,各位,原来上传的关于高斯权函数平滑的Fortran程序有个致命的错误,现在重新上传修正bug后的,大家如果有下载的,要重新下载下,里面的时间引数搞错了,谢谢。

高斯权函数的卷积平滑法有时亦简称为高斯平滑法,在天文地球动力学的观测资料处理中应用较多,其主要原理是将测量数据y(ti), i = 1,2,...,N与一个高斯曲线分布的权函数作为卷积计算,求得测量数据的平滑值。由于高斯权函数平滑法简单,计算机易于实现,且既可用于等间隔又可用于不等间隔数据,因此应用范围较广。

代码附件:

gaussian_smooth_Matlab.zip

gaussian_weighting_function_smoothing_Linux64.zip

gaussian_weighting_function_smoothing_Win64.zip

其中,a为高斯权函数的带宽。Matlab和Fortran用同样的测试数据。

a=0.1时,


a=2时,


带宽a的取值起着控制平滑性质的作用。当a值取值较大时,较多的测量数据在计权平均中起作用,而得到平滑度较强的平滑值;当a取值较小时,意味着较少的测量数据在计权平均中起作用,而得到平滑度较弱的平滑值。




https://blog.sciencenet.cn/blog-350278-738254.html

上一篇:Fortran 90: a simple example of array pointer
下一篇:Fortran 90 实现 Vondrak 滤波
收藏 IP: 222.66.117.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-16 06:17

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部