Fortran Coder

标题: fortran读取文件中实数的精度问题 [打印本页]

作者: matliu 时间: 2015-12-30 15:07
标题: fortran读取文件中实数的精度问题
我写了一个利用fortran程序读取文件数据的程序，发现数据的实数小数点后的数据总是不对呢？请大家帮忙看看啊，先谢谢了！

程序如下：

[Fortran] syntaxhighlighter_viewsource syntaxhighlighter_copycode



PROGRAM t1

IMPLICIT NONE  



INTEGER :: number_nodes

REAL, DIMENSION(6) :: info_nodes



INTEGER :: ierror

INTEGER :: i,j



! open the file

OPEN (UNIT=8, FILE = 'nodes1.dat', STATUS='OLD', ACTION = 'READ', IOSTAT=ierror)

! read  integers

READ (8,*) number_nodes



! read nodes

DO i = 1,number_nodes

   READ(8,*) info_nodes(i)

END DO



WRITE(*,*) 'number of nodes is', number_nodes



DO i = 1,number_nodes

   WRITE(*,*) 'node ', i , 'is ', info_nodes(i)

END DO



CLOSE(UNIT=8)

END PROGRAM t1

其中文件 nodes1.dat的格式为：
6
0.0
0.2
0.4
0.6
0.8
1.0
利用ubuntu下的gfortran编译，得到的结果为：
number of nodes is          6
node          1 is 0.00000000
node          2 is 0.200000003
node          3 is 0.400000006
node          4 is 0.600000024
node          5 is 0.800000012
node          6 is 1.00000000
其中node2---nodes5后面为什么总有些莫名奇妙小数呢？
刚学fortran不就，请大家不吝赐教.

作者: fcode 时间: 2015-12-30 16:25
这是正常的。并不是读取问题，而是浮点数问题。
单精度浮点只有6-7位有效数字。

gfortran的表控格式，会用9位有效数字来输出单精度，这样就会出现显示误差。
你可以指定输出格式来避免显示误差。但你无法避免内存中的存储误差。

作者: matliu 时间: 2015-12-30 22:50
我明白了，谢谢！

作者: matliu 时间: 2015-12-30 22:53
比如我读取的是0.2，真实数值也是0.2吧，是显示有问题！可以这样理解吧。

作者: 楚香饭 时间: 2015-12-30 23:02
本帖最后由楚香饭于 2015-12-30 23:03 编辑

不，真实数值也不是 0.2，IEEE浮点数无法精确的表达 0.2
就好像十进制无法精确的表达三分之一一样。它只能用 0.3333333 近似

作者: pasuka 时间: 2015-12-31 09:42

matliu 发表于 2015-12-30 22:53
比如我读取的是0.2，真实数值也是0.2吧，是显示有问题！可以这样理解吧。

一般而言，双精度的浮点数机器误差大约是2.2e-16，所以在计算机里面的双精度0.2实际范围在0.2-2.2e-16到0.2+2.2e-16之间

作者: matliu 时间: 2015-12-31 15:25
明白，谢谢，高手！

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)