选主元素法实验报告.docVIP

选主元素高斯消去法求解非线性方程组的实验报告 一、实验题目 利用选主元素高斯消去法求解以下线性方程组： 二、实验引言 1.实验目的 (掌握用直接法，即选主元素高斯消去法求解线性方程组的基本思想和步骤，熟。悉fortran语言精确求解该方程组； (并能够利用选主元的高斯消去法解任意阶数的线性方程组。 2.实验意义 消去法虽程序复杂，对高阶矩阵易受到计算机容量的限制，适用于中小型方程组，但该方法精确度高且工作量小 三、算法设计 1.选主元消去法基本思路 （1）消元过程。对于k=1,2,...,n-1，进行下列运算： ①选主元，确定r，t使得若art=0，则系数矩阵为奇异的，停止计算；否则进行下一步。 ②交换A中的r、t两行及t、k两列，并记下交换的号码t、k。 ③对于i=k+1，k+2，…n；j=k+1，k+2，…n+1，计算 （2）回代过程。对于k=n，n-1，…，1，计算 （3）恢复解的顺序 算法流程图 主程序 选主元素 四、程序设计 program main integer,parameter::n=3 real::a(3,3)=(/10.0,-1.0,-1.0,-1.0,10.0,-1.0,-2.0,-2.0,5.0/) real::b(3)=(/7.2,8.3,4.2/) integer::k,i,j,l real::m,t do k=1,n m=abs(a(k,k)) l=k do i=k+1,n if(abs(a(i,k))m) then d=a(i,k) l=i end if if (d/=0) then write(*,*) ! end if if(l/=k) then t=a(l,j) a(l,j)=a(k,j) a(k,j)=t t=b(l) b(l)=b(k) b(k)=t end if end do 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(k,j)*a(i,k) end do b(i)=b(i)-b(k)*a(i,k) end do end do do i=n,1,-1 do j=i+1,n b(i)=b(i)-a(i,j)*b(j) end do end do write(*,*) b(1),b(2),b(3) end 五、结果及讨论 1.实验结果 1.100000 1.200000 1.300000 2.讨论分析 （1）误差 运用选主元素法，减小了舍入误差对解的影响。 （2）收敛性 具有较高的收敛性，收敛速度较快。 六、算法评价 优点：方法精确，较约当消去法进一步节省了计算量。 缺点：程序复杂，对高阶矩阵易受到计算机容量的限制，适用于中小型方程组。 附约当消去法程序： program jordan implicit none integer::k,i,j integer,parameter::m=3,n=3 real::a(m,n)=(/10,-1,-1,-1,10,-1,-2,-2,5/) real::b(m)=(/7.2,8.3,4.2/) open (unit=10,file=1.txt) do k=1,n do j=k+1,n a(k,j)=a(k,j)/a(k,k) end dob(k)=b(k)/a(k,k) do i=1,n if(i/=k) then 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 if end do end do write(10,*) b(1),b(2),b(3) end 程序运行结果： 1.100000 1.200000 1.300000