FORTRAN实验报告2016LI - 360文档中心

F90实验报告

实验1 FORTRAN 90的编程环境与简单程序设计

一、实验目的

1.熟悉FORTRAN90的编程环境——Fortran PowerStation 4.0，并能够在此环境下编写FORTRAN90的程序。

2.掌握FORTRAN90程序中所用到的类型说明语句、赋值语句、读、写语句等语句的语句格式，理解并掌握这些语句的功能。

3.掌握简单程序设计方法，能够用简单的语句编写一些简单的应用程序。

4.掌握FORTRAN90源程序的结构特点与书写规范。

二、实验要求

1.根据题目，设计各题算法并画出流程图，按照流程图编写程序。

2.对每个程序认真进行语法检查与逻辑检查。

3.上机调试程序。按各个程序输入计算机，并且编译和连接，如果在编译、连接中有错误信息，则要根据这些信息找到程序中的错误所在，然后修改程序，直至正确。

4.运行程序，记录数据输入和数据输出的情况。

5.写出实验报告。

三、实验内容

1.开机；进入Fortran PowerStation4.0集成工作环境，了解工作环境及其菜单的使用，输入如下FORTRAN90源程序，并以A1.F90存盘。

PROGRAM EX1

IMPLICIT NONE

PRINT *，'2+3=',2+3

PRINT *，'This is a program！'

END PROGRAM EX1

编译、链接、运行后，结果是：

2.输入并调试以下程序，查看运行结果。

PROGRAM EX2

IMPLICIT NONE

INTEGER::A,B,C

READ *,A,B

C=A+B

PRINT *,'A+B=',C

FORTRAN实验报告——牛顿插值

III / 4
NBDF=DIAG(N).R(N) DO I=N,2,-1 NBDF=NBDF*(S+I-1)/I+DIAG(N).R(I-1) ENDDO NBDF=Y(N)+S*NBDF END FUNCTION NBDF END MODULE INTERPOLATION 测试程序: PROGRAM TEST_1 !TEST For Newton Interpolatoin Formula with Divided Difference USE INTERPOLATION REAL*8::XSTAR=0.23D0,YSTAR CHARACTER(20)::FNAME='DATA01.TXT' CALL INIT1(FNAME) !初始化结点 PRINT '(2A13)','X','Y' DO I=0,N PRINT '(2F13.7)',X(I),Y(I) END DO PRINT* YSTAR=NIFDD(XSTAR) !插值函数 PRINT 100,'XSTAR','YSTAR','ERROR',XSTAR,YSTAR,YSTAR-(EXP(XSTAR)-EXP(-XSTAR))/2 CALL DESTROY() 100 FORMAT(3A13/,3F13.7)

工程分析程序设计实验报告

工程分析程序设计实验报告（FORTRAN 95）

顺序结构程序设计

实验内容 1

1.问题描述

有一个边长为8.5m的八边形草地。计算八边形草地的面积，并输出。要求边长和边数从键盘输入。边长为单精度实数，面积为双精度实数，边数为100以内的整数。

2.流程图

3.程序代码

Program main

Parameter (pi=3.1415926)

Integer (1)::n

real::a

Real (8):: b, p, area1,area

print*,'请输入多边形边长和边数：'

Read*, a，n

b=a/2/sin (pi/n)

p=a/2+b

area1=sqrt(p*(p-a)*(p-b)**2)

Area=area1*n

print*,n,'边形草地面积：', area

End

4.运行结果

实验内容 2

1．问题描述

贷款月利率为R，从银行贷款D元，每月偿还P元，计算还清贷款月数，并输出。编写程序实现。要求R、D和P从键盘输入。计算还清贷款月数通过以下公式计算:

⁄

M=(lgP−lg(P−D∗R))lg(1+R)

其中：M为还款月数结果为实数，最终结果要求为整数。

2.流程图

3.程序代码

Program main

Integer M1

Real R,D,P,M

print*,'请输入贷款月利率、从银行贷款钱数与每月偿还钱数：'

Read*,R,D,P

M= (log10 (P)-log10 (P-D*R))/log10 (1+R)

M1=int(M)+1

print*,’还清贷款月数为:：’,M1

End

4.运行结果

选择结构程序设计

Fortran实验报告样例(1)

n n 1)1(312111+-+⋅⋅⋅-+-实验报告样例

南京信息工程大学Fortran 语言程序设计实验(实习)报告实验(实习)名称循环程序设计实验(实习)日期得分指导老师

系专业班级姓名学号

一.实验目的:

1. 掌握DO 、DO WHILE 语句的语法规则和使用要求。

2. 掌握循环结构概念和循环结构程序设计方法。

3. 掌握强制性终止循环语句EXIT 和CYCLE 的基本功能和使用规则，以及相互区别。

4. 掌握嵌套循环结构概念和应用。

二.实验内容1:

1．问题描述

给定5名学生的考试成绩S,评判每个学生的成绩等级,并输出。编写程序实现之。成绩按以下标准评定等级：

优:90≤S ≤100;良:80≤S<90;中:70≤S<80; 及格:60≤S<70;不及格:S<60。

2．问题分析

已知：5个学生成绩：s 1、s 2、s 3、s 4、s 5，从键盘输入，数据类型为整型。

求：5个学生成绩等级，g1、g2、g3、g4、g5，数据类型为字符型。

求解该问题的基本思想是：从键盘依次读取学生成绩数据，每读一个学生成绩，就按照成绩等级评定标准评定等级，并输出该学生等级。

由于有确定的5个学生，所以读取学生成绩和评定成绩等级共5次。该问题显然是一个重复次数确定的重复处理问题，可通过“计数型”循环结构和DO 循环语句实现。

用整型变量cj 表示学生成绩，用字符型变量grade 表示成绩等级。用整型变量I 作为循环控制变量。

通过分析，设计算法，绘制框图，如图6-2所示。

3. 实验步骤:

1）算法流程图：

fortran实验报告

FORTRAN实验报告

实验⽇期2012年3⽉6⽇

⼀、实验⽬的

1、了解FORTRAN90软件开发环境

2、学会启动Fortran，掌握程序的编辑、编译、构建、运⾏⽅法。

3、了解软件开发环境常⽤菜单、⼯具按钮、环境窗⼝基本功能。掌握⽂件、⽂件夹的移动、复制、删除等操作⽅法。⼆、实验内容

1、运⾏Microsoft Developer Studio软件开发环境。

2、理解有关⼯作区（Workspaces）、项⽬（Projects）、⽂件（Files）的基本概念。

3、掌握⼯作区（Workspaces）、项⽬（Projects）、⽂件（Files）创建⽅法。

4、掌握程序的编辑（Edit）、编译（Compile）、构建（Build）、运⾏（Run）⽅法

三、实验步骤

1、建⽴⽂件夹

2、启动Microsoft Developer Studio软件开发环境，并创建⼯作区、项⽬

3、创建⽂件

4、编译项⽬内源程序⽂件

5、构建可执⾏程序⽂件。

6、运⾏可执⾏程序⽂件

四、实验结果与分析

实验结果：

FORTRAN实验报告

实验⽇期2013年3⽉13⽇

⼀、实验⽬的

1、掌握源程序的书写格式；

2、掌握FORTRAN 90中基本数据类型及其描述、运算；

3、掌握基本语句的使⽤

⼆、实验内容

编写简单程序代码

1、整数类型

program ex0404 integer a

a=3

write(*,*) "a=",a

stop

End

program ex0406 integer(kind=4) write write=2+2*4-3

【报告】彩虹实验报告

【关键字】报告

彩虹实验报告

篇一：XX1344075王彩虹实验四

南京信息工程大学实验（实习）报告实验（实习）名称蒙古高压与我国气温异常的关系分析实验日期指导教师马利系计软院专业软件工程年级13班次3 姓名王彩虹学号XX1344075

一、实验目的

（1）掌握气象要素相互影响研究的基本方法。

（2）进一步巩固Fortran顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。

（3）了解站点资料的数据结构。

二、实验内容

1 问题描述

利用1951-XX年1月蒙古高压标准化强度指数（data/h-p.dat）和1951-XX年1月中国160站气温资料(data/t1601.dat)，根据3.5节计算1月蒙古高压强度与中国160站气温的相关关系。要求以“*.dat”和“*.grd”两种格式保存1951-XX年1月蒙古高压强度与我国气温的同期相关数据。

2 问题分析

已知：1951-XX年1月蒙古高压标准化强度指数序列和1951-XX年中国160站逐月气温站点资料。计算：1951-XX年1月蒙古高压强度与中国160站气温同期相关系数。

三、实验要求

（1）分析问题，理解所给出的方法与技术，理清问题思路，分解任务，设计算法。

（2）输出数据采用有格式输入输出，使输出数据规范、醒目、简洁。

（3）用注释的方法指出程序中函数调用的起始和结束位置，并为函数进行注释。

（4）巩固顺序结构和循环结构的程序设计方法，重点掌握文件、变量、函数、数组的使用方法、数据的有格式输入输出方法。

FORTRAN模拟计算器程序

FORTRAN实验报告

给水排水1001班

武俊槟

1004030134

模拟计算器程序

一、实验目的：

设计一个程序来模拟一个简单的手持计算器，程序支持算数运算+，-，*，/，=以及C（清除），A（全清除）等操作。

二、实验要求：

在运行程序时，显示一个窗口，等待用户输入，用户可以从键盘输入要计算的表达式，输入的表达式显示在串口中，用户键入‘=’符号后，窗口显示结果。

三、实验步骤：

1

2、程序代码见附录。

3、测试数据：经检验程序无错

四、实验心得：

经过这次编程，我熟练掌握了VF课程所教授的知识，理论联系实际，这次作业提高了我的编程水平，发现了自己平时所忽略的知识点，同时也增强了团队合作意识。

在基本确定编程思路画出大致流程草图后，我们对编程做出了安排，大家分工合作。在编程过程中，我们都或多或少的遇到一些困难，有时可以通过自己的查阅相关资料得到解决，而有时却不得甚解，此时，团队的力量便展现了出来，大家共同探讨遇到的问题。通过大家的共同努力，一起解决在编程中所遇到的问题如：如何实现清屏，输出格式的美观化等。尤其是实现清除命令的过程中，我们最后通过文件操作解决了这个问题。

最后再将大家做的程序进行整合，对程序进行debug，查找程序中的错误，大家共同探讨出现的错误，不断地将程序进行改进，最后实现编程要求并美观化。

编程不仅仅但是编程，它最大乐趣是在编程的过程中不断对自己VF知识的总结和完善，同时也考验一个团队的合作意识。通过这次编程，我觉得我的收获是巨大的，编程锻炼了我在遇到问题时的解决能力及思维方式，加强了我的团队合作意识，使我感觉到了团队的力量，并提高了个人分析能力综合素质。

FORTRAN实验报告

《FORTRAN 语言》实验（上机）报告班级·学号姓名实验日期任课教师实验名称实验一顺序结构程序设计验证型一、实验目的（1）熟悉Fortran PowerStation4.0的集成开发环境和上机步骤；（2）掌握源程序的书写格式； (3) 掌握FORTRAN 90中基本数据类型及其描述、运算；(4) 掌握基本语句的使用。二、实验内容：（1）编写程序，将华氏温度转化为摄氏温度，公式为：。要求：输入一个华氏温度)32(9

5-=F C 值F ，按照公式将其转化为摄氏温度值C 。（2）编写程序，其功能是：输入a,b,x 的值，计算并输出y 的值。 b x e y ax 5ln 2+=当x=1.234,a=-5,b=0.9876时，y= 。三、实验平台 Fortran PowerStation 4.0四、程序清单

五、调试和测试结果

六、教师批语与成绩评定：

年月日

《FORTRAN语言》实验（上机）报告

班级·学号姓名实验日期任课教师

实验名称实验二选择结构程序设计验证型

一、实验目的

（1）掌握逻辑表达式的书写方法和逻辑型数据的使用；

（2）学会正确使用逻辑运算符和逻辑表达式；

(3)掌握块IF结构、块CASE结构、逻辑IF语句和算术IF语句。

二、实验内容：

（1）编写程序，输入年、月，求该月的天数。

用year、month分别表示年、月，day表示每月的天数。注意：

①每年的1，3，5，7，8，10，12月，每月31天；4，6，9，11月每月30天；2月闰年为29天，平年为28天。

②年份能被4整除，但不能被100整除，或者能被400整除的年均为闰年。

土建类fortran实践.doc

FORTRAN 语言课程设计任务书

一、实践目的

通过本课程设计，培养程序设计能力以及综合解决实际问题的能力。通过自己分析问题、寻求算法、编写、调试程序的过程，掌握FORTRAN 程序设计与调试方法，提高灵活运用所学知识解决问题的能力。

二、设计任务（选择一题）

学号为3的倍数选第一题，3的倍数多1选第二题，3的倍数多2选第三题，

1．线性方程组求解问题。

一物理系统可用下列线性方程组来表示：

⎥⎥⎥⎥⎦

⎤⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡θ-θ-θθθ--θg m g m N N a a m m m m 2121212

1

11001cos 000sin 00cos 0sin 0sin cos 从文件中读入m1、m2和θ的值，求a1、a2、N1 和N2的值。其中g 取9.8，输入θ时以角度为单位。

要求：

（1）分别用两种方法（例如高斯消去法、矩阵求逆法、三角分解法、追赶法等），定义求解线性方程组Ax=b 的子程序，要求该子程序能求解任意线性方程组。

（2）在主程序中分别调用上面定义的两个子程序，并对求解结果进行对比分析。

（3）绘制以上两个方法所求得的方程解的数据分布图。

2．线性病态方程组问题。

下面是一个线性病态方程组：

⎥⎥⎥⎦

⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1）求方程的解。

（2）将方程右边向量元素b 3改为0.53，再求解，并比较b 3的变化和解的相对变化。

河北大学工商学院Fortran程序设计实验报告

年级2010级学号2010483036姓名王丽丽成绩

专业土木工程实验地点B3-216 指导教师袁满

实验项目输入输出和声明语句实验日期2012-3-2

一、实验目的：练习fortran语言中的输入输出和声明语句的使用

二、实验要求：阅读书上例题并调试运行

三、实验结果及分析

program e41

implicit none

integer a

real b

complex c

logical d

character*(20) e

a=10

b=12.34

c=(1,2)

d=.true.

e="fortran"

write(*,100) a 运行结果：

100 format(1x,i5)

write(*,200) b

200 format(1x,f5.2)

write(*,300) c

300 format(1x,f4.1,f4.1)

write(*,400) d

400 format(1x,l3)

write(*,500) e

500 format(1x,a10)

End

年级2010级学号2010483036 姓名王丽丽成绩

专业土木工程实验地点B3-216 指导教师袁满

实验项目输入输出和声明语句实验日期2012-3-9

一、实验目的：练习fortran语言中的输入输出和声明语句的使用

二、实验要求：练习使用输入输出和声明语句

三、实验结果及分析

（1）program t0402 运行结果：

real r,mj,zc

write(*,*)"请输入圆的半径"

read(*,*) r

pai=3.14

mj=r*r*pai

zc=2*pai*r

write(*,*) "圆的面积是：",mj

Fortran实验报告

Fortran实验报告ex05

4代码

program ex0502

integer(kind=16) :: n,m

real(kind=4) :: x,t,a,b,c,d,e,i

write(*,*)'input n'

read(*,*) n

if(n<=0) stop

i=4.0

t=i/real(n)

x=-2.0

a=0

b=x**2+(-x)**2+sin(x)+sin(-x)

do m=1,n-1

a=a+(x+m*t)**2+sin(x+m*t)

end do

b=(b+2*a)*t/2.0

e=16.0

d=3.0

c=e/d

write(*,*)'real answer=',c

write(*,*)'s=',b

write(*,*)(b-c)/c

end

5 结果

6 分析

当N足够大时，由于分成小梯形的份数过多，每一分的面积过小，在前面一部分梯形面积求和后，下一个及之后所有梯形的面积都被当做小到不对结果产生影响而被舍去。这是误差的来原及精度不随N增大而增大的原因。

解决方案有二，一是更改变量精度，使可累加的梯形面积的下限降低。二是先将一部分梯形面积相加，再将其求和。代码：

program ex0501

integer(kind=16) :: n,m,l

real(kind=4) :: x,t,i,e,b,c,d,f

real(kind=4),allocatable :: a(:),s(:)

write(*,*)'input n'

read(*,*) n

write(*,*)'input l' ！l为每l个梯形先求面积和

read(*,*)l

数值计算方法实验报告(例)

实验报告

一、实验目的

二、实验内容

三、实验环境

四.实验方法

五、实验过程

1实验步骤

2 关键代码及其解释

3 调试过程

六、实验总结

1．遇到的问题及解决过程

2．产生的错误及原因分析

3．体会和收获。

七、程序源代码：

八、教师评语

实验报告

一.试验目的：练习用数值方法求解给定的非线性方程。二.实验内容：求解人口方程: )1(5

.43e 1004.156-+

=λλλ

e

要求误差小于410-。

三.实验环境：PC 计算机，FORTRAN 、C 、C ++、VB 任选一种。四.实验方法：牛顿法

牛顿法简述：牛顿法是一种特殊的迭代法，其迭代公式为：

,2,1,0,)

()

(1='-

=+k x f x f x x k k k k ，

当数列{}k x 收敛时，其极限值x 即为方程的解。

定理：给定方程],[,0)(b a x x f ∈=

1）设0)()(<b f a f ；

2）)(x f ''在],[b a 上不变号，且],[,0)(b a x x f ∈≠'； 3）选取],[0b a x ∈，满足0)()(00>''x f x f ；

则牛顿法产生的序列{}k x 收敛于0)(=x f 在],[b a 内的唯一解x 。五.实验过程:

1．编程: 用C 语言编出牛顿法的源程序。

2. 开机, 打开C 语言编译程序，键入所编程序源代码.

3. 调试程序, 修改错误至能正确运行.

六．实验总结:

（1）牛顿法收敛速度快，但初值不容易确定，往往由于初值取得

不当而使迭代不收敛或收敛慢，但若能保证)()(1+>K K x f x f （称为下山条件），则有可能收敛。把新的近似值看作初值的话会比原来的取得好，有可能落入局部收敛的邻域。

FORTRAN课程设计报告

中南大学

本科生课程设计（实践）设计报告（程序设计基础）

题目常微分方程数值解学生伟

指导教师红雨

学院土木建筑学院

专业班级土建类0908

学生学号1208090814

计算机基础教学实验中心

2010年07月07日

FORTRAN 课程设计报告

摘要

本题是通过长期的观察统计，得到一个食饵一捕食者系统的数学模型如下:

捕食者（鲨鱼）：a=0.5，入 2=0.02, X 2（0）= 2; 绘制食饵和鲨鱼数量随时间的变化情况图。

再由已知的初始条件，通过一阶线性微分方程组得到食饵和鲨鱼数量随

时间变化关系曲线。这一过程首先要通过四阶龙格库塔法得到多组鲨鱼和食饵随时间变化的数据，再通过

FORTRAN 编译程序,

QuickWin 的图形操作得到变化曲线。

关键词：四阶龙格一库塔法、FORTRAN 、QuickWin 图形操作

设计思想（一）总体设想

1、生成数据

用四阶龙格库塔法，取不同的时刻对应鲨鱼和食饵数量，得到多组数据。

2、 Q uickWin 图形操作画出基本图形界面

模仿课本及例题做出一个基本的图形界面，包括界面的大小、颜色、名称等，以为后面绘图作准备。

dx 1 dt

x i (r i

1

X 2)

dX

2

X 2( dt

2

X 1）

，食饵: r i = 1,入 i = 0.1, X i (0)= 25;

3、设计主程序连接数据与QuickWin程序，画出相应曲线

将得到的多组数据与绘图程序连接起来，并最终绘制鲨鱼、食饵随时间变化曲线。

（二）系统模块结构简介

常微分方程的数

值解

设计主程序连接数

据文件和图形界面，

绘制曲线

用Fortran实现《湿度查算表》的电子化

作者：徐元魁

来源：《科学与财富》2016年第34期

（民航华东空管局气象中心上海 200335）

摘要：空气湿度湿表征大气物理状态的一个要素，是民航气象观测基本观测项目之一。在应急人工观测时，通过人工查阅《湿度查算表》的方式得到相对湿度和露点温度值。人工查表的过程比较繁琐，如遇复杂天气，值班人员还要忙于通报天气，很容易造成迟发报。本文将以中国气象局编写的《湿度查算表》为依据，介绍用Fortran语言作为工具进行程序设计，实现电子化《湿度查算表》的思路与方法。经检验，程序计算湿度结果与人工查算结果相吻合。

关键词：湿度查算表，Fortran，二分法

一、引言

空气湿度湿表征大气物理状态的一个要素，是民航气象观测基本观测项目之一。当自动观测设备发生故障时，气象观测员需要在室外进行应急人工观测。其中湿度项目的观测采用传统的干湿表法，将百叶箱内读取到的干湿球温度，通过人工查阅《湿度查算表》的方式得到相对湿度和露点温度值。在查表过程中，由于需要对湿球温度进行订正，因此需要来回翻阅多次，不仅费时，而且很容易出错，如遇复杂天气，值班人员还要忙于通报天气，很容易造成迟发报。如能由计算机计算代替人工查算，将会大大提高观测员的工作效率，帮助观测员节约时间，从而有更多精力进行报文的编制与检查等工作。

本文将以中国气象局编写的《湿度查算表》为依据，介绍用Fortran语言作为工具进行程序设计，实现电子化《湿度查算表》的思路与方法。

二、《湿度查算表》的查算方法简介

《湿度查算表》主要由表1湿球结冰部分，表2湿球未结冰部分以及附表的湿球温度订正值组成。表1和表2中，每栏居中的数值为干球温度，其余每行分别给出了订正参数（n），湿球温度（tw），水汽压（e），相对湿度（U）和露点温度（td）的一组数值。附表为不同型号干湿表的湿球温度订正值，不同干湿表经过各自的湿球温度订正值订正后，就可以从表1或表2查取湿度。比如，浦东机场气象台使用通风速度为0.4m/s的球状干湿表，应使用附表3进行湿球温度值的订正。

Fortran课程设计报告_范文

中南大学

课程设计报告书

课程名称：Fortran语言课程设计

班级：土木工程2005级06班

学号：1201050624

姓名：成丕富

指导教师：蔡旭晖

二○○六年六月

一、课程设计任务书

(3)

二、需求分析

(5)

三、概要设计

(5)

四、系统实现

(6)

五、详细设计

(10)

六、源程序代码

(12)

七、总结

(28)

八、参考文献

(28)

《FORTRAN 程序设计课程设计》任务书

一、实践目的

通过本课程设计，培养利用程序设计技术进行计算机与专业知识的交叉应用，综合解决实际问题的能力。通过自己分析问题、寻求算法、编写、调试程序、解决问题的过程，掌握多模块结构的程序设计与调试方法，以及程序模块间的数据传递方法，建立提出问题、分析问题、解决问题的思维方式，促进灵活运用所学知识解决问题能力的提高。

二、设计任务

制作一个信息管理系统，从以下题目中任选其一。

1．上网查询各省会城市2005年生产总值(GDP)、农林牧渔业总产值、工业总产值，将数据。以统一格式保存在一个纯文本文件tongji.dat 中。编程实现：

（1）从文件tongji.dat 中读取数据，计算全国省会2005年GDP 的平均值，统计GDP 高于平均值的城市数；

（2）按第一产业（即农林牧渔业总产值与工业总产值之和）总产值的升序排序，将排序后的数据存入文件fenxi.dat 中。

（3）统计工业总产值在1000亿元以下、1000~2000亿元、2000~3000亿元、3000亿元以上的城市数，并绘制分布饼图和折线图。

（4）计算各城市2005年较2004年相比的GDP 增长幅度，计算并输出自己所在省的省会城市GDP 增长幅度在各省会城市中的排名。

fortran数值计算基础

数值计算基础

实验一直接法解线性方程组的 (2)

实验二插值方法 (11)

实验三数值积分 (5)

实验四常微分方程的数值解 (7)

实验五迭代法解线性方程组与非线性方程 (9)

实验一直接法解线性方程组

一、实验目的

掌握全选主元消去法与高斯-塞德尔法解线性方程组。

二、实验内容

分别写出Guass 列选主元消去法与追赶法的算法，编写程序上机调试出结果，要求所编程序适用于任何一解线性方程组问题，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。

1、用Guass 列选主元消去法求解方程组

⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--5.58.37.33.47.11.85.16.93.51.53.25.2321x x x

2、用追赶法求解方程组

⎥⎥

⎥

⎥⎥⎥⎦⎤⎢

⎢⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-----000010210000210000210000210000

254321x x x x x 三、实验仪器设备与材料

主流微型计算机

四、实验原理

1、Guass 列选主元消去法对于AX =B

1）、消元过程：将（A|B ）进行变换为)~|~(B A ，其中A ~

是上三角矩阵。即：

⎪⎪⎪⎪

⎪⎭

⎫

⎝⎛→

⎪⎪⎪⎪⎪⎭⎫ ⎝⎛n nn n

n

n nn

n n n

n b a b a b a a b a a a b a a a b a a a

0010122111221

222221111211 k 从1到n-1

a 、列选主元

选取第k 列中绝对值最大元素ik n