最小二乘法fortran应用。ppt课件.pptVIP

最小二乘法;主要内容;AX = b;超定方程组：方程的个数超过未知量的个数，一般来说 是没有解的。即对于任意一组数(x1,x2,…,xn);使得：;即：;可以证明如果如果A是列满秩的，则该方程存在唯一解，且使得;例 1：用最小二乘法求下列超定方程组的近似解;解：;故得方程组;2. 程序编写; if( m n ) then C = transpose(a) D = matmul( c, a ) E = matmul( c, b ) X = D .ix. E ! invert(D) * E else if( m == n ) then X = A .ix. B else write(*,*) No roots! end if ; do i = 1, n write(*,*) X(i) end do stop end program least_square;如何拟合线性y=ax+b;program least_square use IMSL integer, parameter :: m = 6 integer, parameter :: n = 2 integer :: i double precision :: A(m,n), B(m), X(n)