|||
非常抱歉,各位,原来上传的关于高斯权函数平滑的Fortran程序有个致命的错误,现在重新上传修正bug后的,大家如果有下载的,要重新下载下,里面的时间引数搞错了,谢谢。
高斯权函数的卷积平滑法有时亦简称为高斯平滑法,在天文地球动力学的观测资料处理中应用较多,其主要原理是将测量数据y(ti), i = 1,2,...,N与一个高斯曲线分布的权函数作为卷积计算,求得测量数据的平滑值。由于高斯权函数平滑法简单,计算机易于实现,且既可用于等间隔又可用于不等间隔数据,因此应用范围较广。
代码附件:
gaussian_weighting_function_smoothing_Linux64.zip
gaussian_weighting_function_smoothing_Win64.zip
其中,a为高斯权函数的带宽。Matlab和Fortran用同样的测试数据。
a=0.1时,
a=2时,
带宽a的取值起着控制平滑性质的作用。当a值取值较大时,较多的测量数据在计权平均中起作用,而得到平滑度较强的平滑值;当a取值较小时,意味着较少的测量数据在计权平均中起作用,而得到平滑度较弱的平滑值。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-5-16 06:17
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社