fortran用高斯消去法，矩阵求逆法、三角分解法、追赶法中的两种解线性方程组.doc

fortran:用高斯消去法，矩阵求逆法、三角分解法、追赶法中的两种解线性方程组 悬赏分：125 | 解决时间：2010-7-15 20:58 | 提问者：flygg888 一物理系统可用下列线性方程组来表示： 从文件中读入m1、m2和θ的值，求a1、a2、N1 和N2的值。其中g取9.8，输入θ时以角度为单位。 要求： （1）分别用两种方法（例如高斯消去法、矩阵求逆法、三角分解法、追赶法等），定义求解线性方程组Ax=b的子程序，要求该子程序能求解任意线性方程组。 （2）在主程序中分别调用上面定义的两个子程序，并对求解结果进行对比分析。 （3）绘制以上两个方法所求得的方程解的数据分布图。 请将答案发送至flygg888@126.com，答案合格者，将继续追加高分！ 问题补充： 请提供全面一点的答案（不要只求逆矩阵），谢谢了 最佳答案 !高斯消去法 subroutine agaus(a,b,n,x,l,js) dimension a(n,n),x(n),b(n),js(n) double precision a,b,x,t l=1 !逻辑变量 do k=1,n-1 d=0.0 do i=k,n do j=k,n if (abs(a(i,j))d) then d=abs(a(i,j)) js(k)=j is=i end if end do end do !把行绝对值最大的元素换到主元位置 if (d+1.0==1.0) then l=0 else !最大元素为0无解 if(js(k)/=k) then do i=1,n t=a(i,k) a(i,k)=a(i,js(k)) a(i,js(k))=t end do !最大元素不在K行，K行 end if if(is/=k) then do j=k,n t=a(k,j) a(k,j)=a(is,j) a(is,j)=t !交换到K列 end do t=b(k) b(k)=b(is) b(is)=t end if !最大元素在主对角线上 end if !消去 if (l==0) then write(*,100) return end if do j=k+1,n a(k,j)=a(k,j)/a(k,k) end do b(k)=b(k)/a(k,k) !求三角矩阵 do i=k+1,n do j=k+1,n a(i,j)=a(i,j)-a(i,k)*a(k,j) end do b(i)=b(i)-a(i,k)*b(k) end do end do if (abs(a(n,n))+1.0==1.0) then l=0 write(*,100) return end if x(n)=b(n)/a(n,n) do i=n-1,1,-1 t=0.0 do j=i+1,n t=t+a(i,j)*x(j) end do x(i)=b(i)-t end do 100 format(1x,fail) js(n)=n do k=n,1,-1 if (js(k)/=k) then t=x(k) x(k)=x(js(k)) x(js(k))=t end if end do return end program main dimension a(4,4),b(4),x(4),js(4) double precision a,b,x real m1,m2,j open(1,file=laiyi.txt) read(1,*)m1,m2,j close(1) n=4 print*,m1,m2,j a(1,1)=m1*cos(3.14159*j/180) a(1,2)=-m1 a(1,3)=-sin(3.14159*j/180) a(1,4)=0 a(2,1)=m1*sin(3.14159*j/180) a(2,2)=0 a(2,3)=cos(3.14159*j