Fortran语言基础资料讲解.docx

1、Fortran语言基础资料讲解For t r an 语 言 基 础数值分析程序设计Part I Fortran语言基础COMPAQ VISUAL FORTRAN 6.50编译器的使用0.1 编译器简介高级语言以及汇编语言的程序代码在没有转换成机器代码前，计算机是 无法执行的。编译器的功能是将高级语言的程序代码翻译成计算机可执行的机 器码，也就是生成扩展名为EXE, COM的文件。0.2 Visual Fortran 的使用Visual Fortran 起源于 Microsoft 的 Fortran PowerStation 4.Q 这套工具后来 卖给Digital公司继续开发，第二个版本称为

2、Digital Visual Fortran 5.0， Digital 被Compaq并购之后，接下来的版本 6.0和6.5称为Compaq Visual Fortran。下 面的介绍以Compaq Visual Fortran 6.5作范例。Visual Fortran被组合在一个叫做 Microsoft Visual Studio的图形接口开发 环境中。Visual Studio提供一个统一的使用接口，这个接口包括文字编辑功 能、Project管理功能、调试工具等。而编译器则被组合到 Visual Studio中，VF和VC+使用相同的使用接口。Visual Fortran 6.5除了完全

3、支持Fortran 95的语法外，扩展功能方面提供 了完整的Windows程序开发工具，专业版还含有IMSL数值计算连接库。另 外还可以和VC+直接互相连接使用，也就是把 Fortran和C语言的程序代码 混合编译成一个执行文件。安装好Compaq Visual Fortran后，运行Developer Studio就可以开始编译 Fortra n 程序了。运行 Developer Studio启动 Visual Fortran,默认程序名称为 CompaqVisual Fortran 6.5选择File菜单中的New选项在弹出的对话框中，选择 Project标签。在Project中选择For

4、tran ConsoleApplication。在Project Name的文本框中给定 Project的名字，Location会显示出整个Project的工作目录。点击”0K”选择 An empty project选项，单击 Finish 按钮。这个画面只有在VF6.5中才会出现，他显示Project打开后自动生成的文 件，直接单击OK按钮。这是刚设置完成的Project后的界面，目前还没有任何程序代码。以上是建立Project的方法。Visual Studio的环境是以Project作为编译单 位，*.dsp或*dsw是记录Project文件。打开Project后，还要把程序代码加入 Pr

5、oject中才能编译。下面是添加程序代码文件的方法。再次选择File菜单中的New选项在对话框中，选择 Files标签，选择Fortran Free Format Source File并为文件命名。编写一个程序或打开一个已有的程序选择Build菜单中的Execute选项，VF会编译并执行编译好的程序程序执行的结果最后再一次简单地说明编译程序的过程：（1） 建立一个新的Project （File/New，选择Project选项卡，选择 Fortran ConsoleProgram格式，给定 Project名称）。Project会保存成 * dsw 文件。（2） 生成一个新的程序文件（File/

6、New，选择Files选项卡，选择Fortran FreeFormat SourceFile，给定文件名），或是插入一个已有的程序 文件（选项Project/Add to Project/Files）。程序代码会保存成*f90 或*.for文件。（3） 用Build菜单中的Execute选项来编译并运行程序，或只是单击 Build选项来只作编译，不运行程序。（4） 要写新的程序可以另外建立一个新的 Project,或是直接更换Project 中的文件。千万不要把两个独立的程序文件放在同一个 Project中， 否则导致编译过程出现错误。（5） 下次修改程序时，可以直接使用 File/Open

7、Workspace来打开*.dsw 的Project工程文件。作业：1、熟悉Fortran编译器的使用方法。1 Fortran语言基本概念1.0程序书写格式Fortran程序代码的编写格式有两种：Free Format （自由格式）和FixedFormat （固定格式）。固定格式属于旧式写法，他在编写版面上有很多限制。自由格式是Fortran90之后的新写法，取消了许多旧的限制。 Fortran程序代码扩展名为*.F或*FOR的文件是以Fixed Format编写的程序；以*.F90为扩展名的文件，是 以Free Format编写的程序。今后我们编写的程序建议采用自由格式编写。1.0.1 Fi

8、xed Format （固定格式）在固定格式中，规定了程序代码每一行中每个字段的意义。第 772个字符是可以用来编写程序的字段。每一行的前 5个字符只能是空格或者是数字，数字用来作为“行代码”。每一行的第6个字符只能是空格或者“0”以外的字符。第1个字符：如果是字母C,c或者星号*，这一行文本会被当作说明批 注，不会被编译。第15个字符：如果是数字，就是用来给这一行程序代码取个代号。不 然只能是空格。第6个字符：如果是“0”以外的任何字符，表示这一行程序会接续上一行;第772个字符：Fortran程序代码的编写区域。第73个字符之后：不使用，超过的部分会被忽略，有的编译器会发出错 误信息。C

9、FIXED FORMAT DEMOprogram mai nwrite（*,*） 15000.789 write (*,*)1Hello100 write（*,*） Hello10 stopend注意：程序中的空格，没有任何意义。固定格式是配合早期使用穿孔卡 片输入程序所发明的格式。熟悉固定格式，有助于阅读早期的 Fortran程序代码。1.0.2 Free Format （自由格式）自由格式基本上允许非常自由地编写格式，他没有规定每一行的第几个 字符有什么作用。需要注意的事项只有以下几点：（1） 叹号“！”后面的文本都是注释；（2） 每一行可以编写132个字符；（3） 行号放在每行程序的最前

10、面；（4） 一行程序代码的最后如果是符号“ &”,代表下一行程序会和这一行连 接。如果一行程序代码的开头是符号 & ,代表它会和上一行程序连 接。看一个自由格式编写的程序：! Free Formatprogram mai nwrite（*,*） Hello !这也是注解write （*,*） &Hellowr&te （*,*） Helloend1.1字符集字符集是指编写Fortran程序时，所能使用的所有字符及符号。Fortran所 能使用的字符集有：英文字母：AZ及az（英文字母不区分大小写）；数字：09；22个特殊符号：=+ - * / . ！ “ % & ; ? $ _ （还有一个显示不

11、出来的空格符）注意：Fortran是不区分大小写的语言。所以，INTEGERnteger,inteGER是相同的命令。1.2数据类型整数（INTEGER）:整数又可以分为长整型与短整型。在个人计算机中 长整型占用32 bits（4 bytes）的空间，长整型可以保存的数值范围在- 2147483648+2147483647（也就是-231+1231）之间；短整型占用 16 bits （2 bytes）的空间，保存数值的范围在-32768+32767之间。浮点数（REAL）:浮点数有两种类型：单精度和双精度。单精度浮点数占用32 bits的空间，有效位数为67位。可记录的最大数双精度浮点数占用6

12、4 bits，有效位数为1516位。可记录的最大数值为1.79 10308，最小数值为 2.23 10 308。复数（COMPLEX ）:单精度和双精度。以a bi表示的复数，其中a,b是由两个浮点数来记录的字符型（CHARACTER ）:计算机除了记录数之外，可以记录一段文本 或称作字符串。记录一个字符需要一个字节。逻辑型（LOGICAL ）:只有两种结果：真（TRUE）和假（FALSE） 在二进制中，通常以1代表真，0代表假。使用Fortran编写程序，必须事先声明数据类型。不同数据类型必须经过 转换才能互通。1.3数学表达式Fortran使用的数学运算符号，根据运算优先级顺序排列如下:+

13、 加法-减法* 乘法 /除法* 乘幕（乘方）（）括号越是下面的运算符号，运算的优先级越高。Fortran的数学表达式和手写的差别主要有三点：(1) 乘幕要连用两个星号*,例如43必须写成4*3 ；(2) 乘号不能省略，2a应当写作2*a；(3)2(e f)除法算式(a b)(c d)应当写作(a+b)*(c+d)/(2*(e+f);1.4变量命名在Fortran中，变量命名要注意一些原则：(1)变量的名称以使用英文字母为主，可以含有下划线或数字，但前 缀必须是英文字母；(2)变量名的长度最多31个字符；(3)变量名不要与Fortran的命令名相同，也不要与程序名相同；(4)程序中变量名不区分大

14、小写。(5)1.5 Fortran程序结构例如:program main write (*,*) Hello stop end作业:1、变量命名有哪些要求?2、变量的类型有哪些?2 输入输出及声明WRITE命令语法：WRITE(UNIT, FMT)注释：屏幕输出。UNIT输出位置FMT输出格式。默认格式WRITE(*,*)。默认输出格式unit=6，也就是屏幕输出。program mai nwrite(*,*) Hellostopend program mai n其他的输出格式参考Fortran的帮助。PRINT命令program mai nprint *,Hellostopend progr

15、am mai nPrint的用法和write大致相同，只是print后面不适用括号，而且只有一个星 号。这个星号的意义是不限定输出格式。声明所谓声明是指：在程序代码中，程序员向编译器要求预留一些存放数据的 内存空间。看一个例子：program ex0404in teger aa=3write (*,*) a=,astop endin teger a声明变量为整型。编写程序时，应当根据变量的类型选择适当的 数据类型，数据类型选用不当，可能造成错误。这种错误有时是编译器无法察 觉的。2.3.1整数类型（INTEGER）语法：integer（kind=n） a注释：kind的值决定整型数据类型。ki

16、nd=2短整型，kind=4长整型。默 认为长整型。In teger aA=1.5最终的结果为a=1，这是因为a为整型，会把小数部分无条件舍去。2.3.2浮点数（REAL）语法：real（kind=n） a注释：kind的值决定浮点数的精度类型。kind=4单精度浮点数（有效数 字位数67位），kind=8双精度浮点数（有效数字位数15位）。默认为单精 度浮点数。超大数值的设置：单精度浮点数：1.23E20双精度浮点数：1.232365D10计算机在进行浮点数计算时，需要对位，由于有效数位的限制，可能造 成大数吃小数的情况。2.3.3复数（COMPLEX）语法：complex（kind=n）

17、aa=（x,y）注释：kind的值决定浮点数的精度类型。kind=4单精度浮点数（有效数 字位数67位），kind=8双精度浮点数（有效数字位数15位）。默认为单精 度浮点数。2.3.4字符及字符串（CHARACTER）语法：characte（len=n） a注释：len决定字符串的长度。字符串赋值由双引号定义a=”am a stude nt”program ex0413characte(le n=20) stri ngstri ng = Good mornin g.write(*,*) stringstring(6:) = evening. !重设设定从第6个字符之后的字符串 write(*

18、,*) stri ngend执行结果：Good morning.Good eve ning.注意：Fortran的命令不区分大小写，但是在字符串中是区分大小写的。 关于字符串的有关函数，请查阅相关的帮助。2.2.5 逻辑变量(LOGICAL)语法：logical a注释：逻辑变量的值为“true”和fafse”。逻辑变量输出值为”T”和”F”a=.true.a=.false.逻辑变量用于程序的逻辑判断语句，一般并不输出。READ命令语法：reac(*,*) aprogram ex0417in teger areac(*,*) a !由键盘读入一个整数write(*,*) a !写出读进变量a的

19、内容endRead命令用于实时接受用户从键盘输入数据。Read命令配合输入格式可以 方便读取文件中的数据。有关格式命令参考帮助文件。格式化输入(FORMAT)命令2.5.1格式化输出概论Format命令用来设置输出格式，看一个例子：program ex0420in teger aa=100write(*,100) a !使用行代码100地方设定的格式来输出变数 a100 format(l4)!最前面的100是行代码,把这一行程序代码给一个编号endH Visual StudiMyProjecfelianxicbugI 199iPreise any to coet inus*.程序的执行结果使用

20、4个字符宽输出整数，前面没有多余的空格。Format命令中可以使用很多的格式控制描述。下面列出所有格式命令的 功能(中的选项可以省略)：Aw以w个字符宽来输出字符串；BN 定义文本框中的空位为没有东西，在输入时才需要使用；BZ 定义文本框中的空位为0,在输入时才需要使用；Dw.d以w个字符宽输出指数类型的浮点数，小数部分占用 d个字符宽；EW.dEe以w个字符宽输出指数类型的浮点数，小数部分占用 d个字符宽，指数部分占用e个字符；En w.dEe以指数类型输出浮点数；Esw.dEe以指数类型输出浮点数；Fw.d以w个字符宽输出浮点数，小数部分占用 d个字符宽；Gw.dEe以w个字符宽输出任何种

21、类的数据；Iw.m以w个字符宽输出整数，最少输出 m个数字；Lw以w个字符宽输出T或F的真假值；nX 把输出的位置向右跳过n个位置；/ 代表换行；: 在没有更多数据时结束输出；下面看一个例子：program ex0421in teger areal bcomplex clogical dcharacte(le n=20) ea=10b=12.34c=(1,2)d=.true.e=FORTRANwrite(*,(1X,l5) a !用I来格式化整数 write(*,(1X,F5.2) ) b !用F来格式化浮点数 write(*,(1X,F4.1,F4.1) ) c ! complex 也是浮点

22、数 write(*,(1X,L3) d !用 L 来输出 logical write(*,(1X,A10) e !用 A 来输出字符串endCMfcruSiurt Vkujl Sludic MyPnijAc !12 341 d 2 eI FOftTftANang key to continue上例将输出格式写在Write命令内，下面的例子采用Format命令规定输出格式：PROGRAM ex0422INTEGER AREAL BCOMPLEX CLOGICAL DCHARACTER *(20) EA=10B=12.34C=(1,2)D=.true.E=FORTRAN WRITE(*,100)

23、A !用I来格式化整数 WRITE (*,200) B !用F来格式化浮点数 WRITE (*,300) C ! complex 也是浮点数 WRITE (*,400) D !用 L 来输出 logical WRITE (*,500) E !用A来输出字符串100 FORMAT (1X,I5)200 FORMAT (1X,F5.2)300 FORMAT (1X,F4.1,F4.1)400 FORMAT (1X,L3)500 FORMAT (1X,A10)END详细的输出格式使用方法，参考相关帮助。2.5.2详论格式化输出格式化输出的控制字符非常丰富，但是常用的并不多，所以不需要记住 每一个控制

24、字符。一般来说，“ I、F、E、A、X”是最常用的几个格式，最好能 记住他们的用法。【lw.m】一一以w个字符的宽度输出整数至少输出m个数字write(*, ”(15) ”)100WW100 ! 采用5个字符宽输出，100前补两个空白write(*, ”(13) ”)10000* !设置的输出文本框不足，输出 3个*以示警告write(*, ”(15.3) ”)10WW010 ! 强迫输出声明的其他命令261 IMPLICT 命令Fortran编译器会根据变量名的首字母自动决定变量类型，以 I、J、K、L、M、N开头的变量被视为整型变量。implicit none关闭默认类型功能，所有变量都要

25、事先声明；implicit real(M-P)implicit integer(A,B,C)2.6.2常数的声明方法(PARAMETER)程序中的一些不变量，可以声明为常数，其在程序中是不变的。看一个例子：program ex0429implicit nonereal piparamete(pi=3.14159)write(*,(F4.2) si n(pi/6) end2.6.3设置变量初值程序中的变量可以在声明变量类型时同时赋值program ex0430in teger : a = 1real : b = 2 complex : c = (1,2)characte(len=20) : st

26、r = FORTRAN 90write(*,*) a,b,c,strend在Fortran77中，需要使用Data命令设值初值。264 等价声明(EQUIVALENCE )把两个以上的变量，声明他们使用同一个内存地址，就是等价声明。使用 等价声明可以节省内存、精简程序代码。in teger a, bequivalence(a,b) !声明a,b这两个变量使用同一块内存空间注意：所有的声明都应当放在程序代码的可执行描述之前 。混合运算混合运算是指：在算式中所进行计算的数字类型不同。进行混合运算时， 要先经过数据类型的转换，以免出现错误。program ex0432implicit nonein

27、teger: a=1in teger: b=2real : cc=real(a)/real(b) !经由库函数real把整数转换成浮点数 write(*,(F5.2) cend程序的执行结果：0.5如果不做类型转换，c=a/b=0,出现错误的结果。自定义数据类型利用自定义数据类型可以定义数据结构体，管理不同类型的数据，这对于 数据库建立是非常有用的。program ex0434implicit none!开始建立person这个类型type : pers oncharacte(le n=30) : n ame !人名in teger : age !年龄in teger : height !身高

28、in teger : weight !体重characte(le n=80) : address !地址end type persontype(person) : a!声明一个 person类型的变量write(*,*) NAME:reac(*,*) a%n amewrite(*,*) AGE:reac(*,*) a%agewrite(*,*) HEIGHT:reac(*,*) a%heightwrite(*,*) WEIGHT:reac(*,*) a%weightwrite(*,*) ADDRESS:reac(*,(A80) a%addresswrite (*,100) a%n ame,a%

29、age,a%height,a%weight100 format(Name:,A10/,Age:,l3/,Height:,l3/,Weight:,l3，& Addres:,A50)stopend每个类型为 person的变量中，都有 name, age, height, weight, addres这几个元素可以使用，使用时要加上一个百分号 来取用他们。Perso n%age,pers on %weight。Fortran90中的type功能与C语言中的struct类似。3流程控制与逻辑运算3.1IF语句IF语句基本用法IF语句基本流程：IF (逻辑表达式)THEN程序代码END IFprogr

30、am ex0501implicit nonereal(k in d=4) : speedwrite(*,*) speed:reac(*,*) speedif ( speed 100.0 ) then! speed 100时才会执行下面这一行程序write (*,*) Slow dow n.en d ifstopendIF命令可以与ELSE搭配，逻辑流程为：IF (逻辑表达式)THEN程序代码ELSE程序代码END IFprogram ex0502implicit nonereal(k in d=4) : height !记录身高real(kind=4) : weight !记录体重write(*,*) height:reac(*,*