FORTRAN语言课程研究设计.docVIP

FORTRAN语言课程设计 摘要： 科技的日新月异使得计算机领域不断取得新的研究成果。计算机在代替和延伸脑力劳动方面发挥越来越重要的作用，不仅在工业方面而且在日常生活和科研中也越来越离不开计算机。特别是在天体运动方面需要运用到计算机处理大量的数据。这次我选的实践课题是用Jacobi迭代Gauss-Seidel迭代语言接近数学公式的自然描述，在计算机里具有很高的执行效率的最大特性语言是一种极具发展潜力的语言，在数值计算中，Fortran语言仍然不可替代。Fortran90标准引入了数组计算等非常利于矩阵运算的功能。在数组运算时，Fortran能够自动进行并行运算，这是很多编程语言不具备的。运用ortran语言，你能够运用很多现成的函数软件包，所以非常便利。Jacobi迭代Gauss-Seidel迭代Jacobi迭代Gauss-Seidel迭代b。一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，在五个不同的点对小行星作了五次观察，测得轨道上五个点的坐标数据（单位：万公里）如下表所示： P1 P2 P3 P4 P5 X坐标 53605 58460 62859 66662 68894 Y坐标 6026 11179 16954 23492 30890 由开普勒第一定律知，小行星轨道为一椭圆，椭圆的一般方程可表示为： a1x2+2a2xy+a3y2+2a4x+2a5y+1=0 分别将五个点的数据代入椭圆一般方程中，得到线性方程组，求出待定系数a1,a2,a3,a4,a5。绘制行星绕太阳运行的轨道。 从题目中观察到一共需要五组数据联立成方程组才能将a1,a2,a3,a4,a5求出，自然想到运用Jacobi迭代Gauss-Seidel迭代a1,a2,a3,a4,a5和系数。同时为了满足程序的实用性我将Gauss-Seidel迭代Gauss-Seidel迭代Gauss-Seidel迭代Gauss-Seidel迭代法与高斯消元法相似，只是Gauss-Seidel迭代法每迭代一次只需计算一次矩阵与向量的乘法。其运算公式如下： Gauss-Seidel迭代法： x(0)=(x1(0), x2(0), x3(0),…xn(0))T xi(k+1)=( bi- ijxj(k+1)-ijxj(k)) ii (i=1,2, …,n;k=0,1, …) 运用Gauss-Seidel迭代法的算法源程序编码如下： dimension a(n,n),b(n),x(n) real a,b,x,t,s,p,q do i=1,n if(abs((a(i,i))+1.0)==1.0)then l=0 write(*,100) return endif enddo l=100 do i=1,n x(i)=0.0 enddo do p=0.0 l=l-1 do i=1,n t=x(i) s=0.0 do j=1,n if(j/=i)s=s+a(i,j)*x(j) enddo x(i)=(b(i)-s)/a(i,i) q=abs(x(i)-t)/(1+abs(x(i))) if(qp)p=q enddo if((peps).or.(l==0))exit enddo if(l==0)write(*,100) return 100 format(1x,fail) End 运用Gauss-Seidel迭代法的程序框图如下： 要注意的是在这两种方法中也要记录迭代过程的迭代次数，特别地，在运用Gauss-Seidel迭代法时一定要将k的值慢慢加大，防止出现死循环。 总的算法设计好之后还要注意画图程序的设计，利用QuickWin来画图不仅要确定好坐标轴的位置，同时要注意x、y轴坐标的刻度，在描点是也要控制好精度以便于确定点的疏密程度。具体的画图源程序编码如下： use msflib type(qwinfo) winfo type(wxycoord) wxy integer(2) status character*5 c real(8) x,y winfo.type=qwin$max status=setwsizeqq(qwin$framewindow,winfo) status=setwsizeqq(0,winfo) status=setwindow(.true.,-100,-5000,10000,5000) status=setbkcolorrgb(#000000) call clearscreen($gclearscreen) status=setcolorrgb(#00ffff) call moveto_w(-100d0,0d0,wxy) status=lineto_w