Fortran学习18：数值计算方法1

最新推荐文章于 2026-04-06 09:19:07 发布

原创

最新推荐文章于 2026-04-06 09:19:07 发布 · 2.2k 阅读

标签

#学习 #开发语言

本文介绍了使用Fortran实现的三种非线性方程求解方法：二分法、割线法及牛顿法，并提供了具体代码示例。通过对比不同方法的优缺点，帮助读者理解数值计算中的求根技巧。

指针之后的问题，包括module、调用其他语言以及GUI，我后面有时间再做吧，感觉对于我来说没啥大的感觉，因为module他其实和c++的class是异曲同工的，没有多少需要用大力气学习的地方，至于调用其他语言，那就更算了，没见过谁Fortran调用c和VB的，感觉没啥用，要调用，也是python调用c++，或者后面看看Python能不能调用Fortran，至于GUI，算了算了，书上的gui画面难看死了，还不如抽时间去学一学Qt，或者是python的canvas、pyqt5和mayavi，再不济学学matlab作图Fortran的作图，就了解一下，知道他还能做gui就行了。

现在主要还是数值计算方法，这个非常重要，虽然Fortran他是又imsl库来解决的，可是你要是不知道数值计算如何进行的，那后面正反演方程你就不会写，更做不出来东西。

一、求解非线性函数：

主要是三种方法：1、二分法(Bisection)2、割线法(Secant)3、牛顿法(Newton,也是最重要的)

1、Bisection：

基本思想：

1、先猜两个值a,b，要让f(a)*f(b)<0，即根据零点定理，再[a,b]之间必有一点c使得f(c)=0

2、取 $c=\frac{a+b}{2}$ ,如果f(c)=0，就找到了正解，如果不等于零，计算f(c)，判断f(a)*f(c)以及f(c)*f(b)，看哪个小于零，在根据零点定理，重复1、2，直到f(c)=0或者是c的取值达到精度。

来一个 $f(x)=sin(x)-x$

real function func(input)
    implicit none
    real ::input
    func=input-sin(input)
    return
end function
program main
    implicit none
    REAL,PARAMETER::lim=0.0001
    REAL i,j,mid
    REAL,external::func
    ! [-1,1]
    i=1.0
    j=-1.0
    do while((func(i)==0).or.(i-j<lim))
        if(func(i)==0)then
        write(*,*)i
        end if
        exit
        if(func(j)==0)then
        write(*,*)j
        end if
        exit
        if(func(i)*func(j)<0)then
            mid=(i+j)/2
            if(func(mid)*func(i)<0)then
                i=mid
            else
                j=mid
            end if
        end if
    end do 
    write(*,*)mid
stop "over"
end program

这个程序的输出很简单，因为f(0)=0了，所以输出是这样的：