五点平滑 fortran

下面是一个简单的五点平滑的Fortran代码示例：

program five_point_smooth
   implicit none
   
   integer :: i, n
   real :: x(100), y(100), y_smooth(100)
   
   ! 读入数据
   open(unit = 1, file = "data.txt")
   n = 0
   do
      read(1, *, iostat = i)
      if (i /= 0) exit
      n = n + 1
      read(1, *) x(n), y(n)
   end do
   close(1)
   
   ! 应用五点平滑
   do i = 3, n-2
      y_smooth(i) = (y(i-2) + y(i-1) + y(i) + y(i+1) + y(i+2)) / 5.0
   end do
   
   ! 输出结果
   open(unit = 2, file = "smoothed_data.txt")
   do i = 1, n
      write(2, *) x(i), y_smooth(i)
   end do
   close(2)
   
end program five_point_smooth

在这个例子中，我们首先读入数据文件 "data.txt"，然后对数据应用五点平滑算法，并将结果写入文件 "smoothed_data.txt"。请注意，我们从第三个数据点开始应用平滑算法，并在倒数第二个数据点结束，以确保平滑后的数据点数量与原始数据点数量相同。