计算机程序设计基础—FORTRAN

中南大学

本科生课程设计(实践)任务书、设计报告（计算机程序设计基础—FORTRAN）

题目线性病态方程组问题

学生姓名刘鑫

指导教师陈语林

学院土木工程学院

专业班级土建类1007班

学生学号1208100724

计算机基础教学实验中心

2011年7 月1日

一、任务内容

利用文件读写、数组、子程序、循环等内容解决线性病态方程组求解、求条件数问题，并绘制方程解的数据分布图。

二、系统开发设计思路

1、系统结构图

三、系统功能及系统详细设计

1、系统功能

通过调用子程序，能求出任意线性病态方程组的解以及系数矩阵的范数并绘出解的数据分布图

求解：全选主元高斯消去法，利用全选主元、归一化、消去、回代过程求出任意方程组的解。再结合矩阵求积子程序，利用循环提高解的精确度从而求出线性病态方程组的解。

求逆矩阵：利用高斯-约当法求逆矩阵。

绘图：利用QuickWin绘图，并设置背景颜色、绘图颜色、标注相应数据等。2、模块

求解求逆矩阵高斯消去法绘图：利用QinckWin 绘图

四、遇到的问题与解决的方法

1、读写过程中无法将文件插入到主程序当中。

解决方法：在主程序当中进行文件操作，先产生文件后调用。 2、用高斯消去法求解得到的解误差很大

解决方法：设定一个EPS 作为精确度，再用三中的求解子程序求解。 3、在求解范数中输出的为7.48E+28 解决方法：请教老师后发现子程序中缺少定义导致相除是分母为0，后添加定义。

五、总结

通过两周的计算机实践，我的计算机基础以及程序编写能力均有很大的提高，特别是对Fortran 这门语言有了更深的了解，使我在今后的学习生活中能更加熟练地使用这门实用的技术。

六、参考文献

[1] 刘卫国，蔡旭辉.FORTRAN90程序设计教程.北京：京邮电大学出版社，2007

[2] 徐士良.FORTRAN 常用算法程序集.北京：清华大学出版社，1995

FORTRAN程序设计复习题及答案

FORTRAN程序设计复习题一、选择题 B (1)下列各FORTRAN表达式中合法的是 A) S+T*2P >= B) .NOT. (A*B+C) C) A2+B2/(C+D) <= D) (A+B).NOT．A*B．GT．(.NOT.只跟一个表达式) C (2)数学式(3/5)ex+y的FORTRAN表达式是 A) 3*EXP(X+Y)/5 B) 3*E* *(X+Y)/ C) (3/5)*EXP(X+Y)D) EXP(X+Y) D (3)下列FORTRAN77表达式中不合法的是 A) A．GT．B．EQV．C．GT．D B) A．AND．B．AND．C．AND．D C) .NOT．(X．LE．D) A．LT．B．LT．C．LT．D D(4)下列叙述中不正确的是 A) FORTRAN子程序可以单独编译 B) 对一个FORTRAN源程序进行编译和连接无误后可生成可执行文件 C) 即使编译和连接都正确无误，FORTRAN程序运行时仍可能出错 D) FORTRAN连接的主要任务是把函数库中的函数翻译成机器指令(正确描述：主要任务为连接目标文件) B (5)在下列FORTRAN77运算符中，优先级最高的是 A) .AND. B) .NOT. C) .OR. D) .EQ. B (6)FORTRAN表达式"6/5+9/2**3/2"的值为 A) 33 B) 1 C) 5 D) 3 A (7)下列FORTRAN77表达式中，合法的是： A) .AND.. B) 10.0 C) D) 提示：A）相当于 .AND.（.NOT.（）） D (8)关于编译一个FORTRAN源程序文件，下列说法中错误的是 A) 允许编译只有一个主程序而没有子程序的源文件 B) 允许编译有多个子程序的源文件 C) 允许编译只有一个子程序而没有主程序的源文件 D) 允许编译有多个主程序的源文件 C (9)在FORTRAN77源程序中，续行标志符必须放在 A) 第1列 B) 第1-6列C) 第6列D) 第5列 D (10)下列关于"SUBROUTIN E MAP(X，Y)"语句行的叙述中，不正确的是 A) 这是子程序的第一个语句 B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据提示：子程序无返回值，自定义函数才有） A (11)FORTRAN表达式"2/4+"的值是 A) B) 1 C) D) 0 提示：2/4默认等于整型，=》 D (12)FORTRAN表达式"MOD，"的值是 A) B)0.0 C) D) A (13下列FORTRAN运算符中，优先级最低的是 A）逻辑运算符.AND. B）算术运算符*

(完整)Fortran经典编程语言笔记(你值得拥有)

FORTRAN笔记 2014.10.07 目录第七讲_FORTRAN的基本知识.ppt (2) FORTRAN语言程序设计初步 (2) FORTRAN源程序的书写格式（以77为例） (2) 变量 (2) 变量类型 (2) 算术运算符和运算优先级 (3) 赋值语句 (3) 参数语句（PARAMETER语句） (3) END语句 (3) PAUSE语句 (3) 逻辑运算和选择结构 (4) 关系表达式 (4) FORTRAN中数组的定义及使用 (4) 其他 (5) 1. fortran语言定义CHARACTER*6 TTL(14,3),CNAM(400)是什么意思? (5) 2. fortran里character*10 是什么意思 (5) 3. Fortran中kind是什么函数？ (5)

第七讲_FORTRAN的基本知识.ppt FORTRAN语言程序设计初步 FORTRAN是Formula Translation的缩写，意为“公式翻译”，它是为科学、工程问题或企事业管理中的那些能够用数学公式表达的问题而设计的，其数值计算的功能较强。常用的是FORTRAN77和FORTRAN90两种标准。 1、一个程序由若干个程序单位组成。主程序和每一个子程序分别是一个独立的程序单位。 2、每一个程序单位都是以“END”结束的。 3、一个程序单位包括若干行。 1）语句行。由一个FORTRAN语句组成。 2）非语句行，即注释行。 4、FORTRAN程序中的语句可以没有标号，也可以有标号，根据需要而定。标号的作用是标志一个语句以便被其他语句引用。 5、一个程序单位中各类语句的位置是有一定规定的。 6、FORTRAN源程序必须按一定的格式书写。 FORTRAN源程序的书写格式（以77为例）每一行有80列，分别如下： 1、第1-5列为标号区。一行中第一列为“C”或“*”，该行即被认为是注释行。 2、第6列为“续行标志区”，如果在一行的第6列上写一个非空格和非零的字符，则该行作为其上一行的续行。 3、第7-72列为语句区。 4、第73-80列，注释区。变量变量名：一个变量需要用一个名字（变量名）来识别。在同一个程序单位中不能用同一个变量名代表两个不同的变量。 FORTRAN的变量名按以下规则选定： 1）第一个字符必须是字母，即变量名必须以字母开头； 2）在一个字母后面可以跟1-5为数字或字母。如果选定的变量名超过6个字符，则只有前面6个字符有效。注：在变量名中大写与小写字母是等价的。变量类型整型变量Integer、实型变量Real、双精度变量Double Precision、复型变量Complex、逻辑型变量Logical和字符型变量Character。 1、隐含约定（I-N规则） FORTRAN规定：在程序中的变量名，凡以字母I，J，K，L，M，N六个字母开头的，即认为该变量为整型变量。在程序中，凡是变量名以字母I,J,K,L,M,N,i,j,k,l,m,n开头的变量被默认为整型变量，以其他字母开头的变量被默认为实型变量。 2、用类型说明语句确定变量类型 1）INTEGER语句（整型说明语句） 2）REAL语句（实型说明语句） 3）DOUBLE PRECISION语句（双精度说明语句） 4）COMPLEX语句（复型说明语句） 5）LOGICAL语句（逻辑型说明语句）

fortran语言语法

FORTRAN是世界上最早出现的高级编程语言，是工程界最常用的编程语言，它在科学计算中（如航空航天、地质勘探、天气预报和建筑工程等领域）发挥着极其重要的作用。经过40多年的发展，伴随着FORTRAN语言多次版本的更新及相应开发系统的出现，其功能不断完善，最新版本的开发系统几乎具备了VC、VB的所有特点，如图形界面编程、数据库等。目前，工科院校开设的计算机编程语言课首选仍然是FORTRAN :< 说实话，从科技发展的趋势来说这不是好事。您可以设想一下，如果需要用鹅毛笔抄写大量的古籍是什么感受！强烈建议阅读《发掘C#特性赋予科学计算项目以威力》 1 FORTRAN77四则运算符 + - * / ** (其中**表示乘方) 在表达式中按优先级次序由低到高为： +或-→*或/→**→函数→() 2 FORTRAN77变量类型隐含约定：I-N规则凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。用类型说明语句确定变量类型：可以改变I-N规则

用IMPLICIT语句将某一字母开头的全部变量指定为所需类型如IMPLICIT REAL (I,J) 三种定义的优先级别由低到高顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。数组的说明与使用使用I-N规则时用DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例：DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3)

fortran基本函数

FORTRAN 90标准函数（一） (2012-07-03 17:14:57) 转载▼ 分类：学习标签： fortran 函数教育符号约定： ●I代表整型;R代表实型;C代表复型;CH代表字符型;S代表字符串;L代表逻辑型;A代表数组;P代表指针;T代表派生类型;AT为任意类型。 ●s:P表示s类型为P类型(任意kind值)。s:P(k)表示s类型为P类型(kind值=k)。 ●[…]表示可选参数。 ●*表示常用函数。

注：三角函数名前有C、D的函数为复数、双精度型函数。注：指数函数名、平方根函数名、对数函数名前有C、D的函数为复数、双精度型函数。表4 参数查询函数

atan2函数的值域是多少？我从网上找到一个fortran函数的日志，说此值域是-π~π，但正常反正切函数的值域应该是-π/2~π/2。对atan2函数不够了解，所以不知道你的答案对不对，我个人认为不对。我是用正常的反正切函数atan(v/u)来算的： FORTRAN： if (u>0..and.v>0.) dir=270-atan(v/u)*180/pi if (u<0..and.v>0.) dir=90-atan(v/u)*180/pi if (u<0..and.v<0.) dir=90-atan(v/u)*180/pi if (u>0..and.v<0.) dir=270-atan(v/u)*180/pi if (u==0..and.v>0.) dir=180 if (u==0..and.v<0.) dir=0 if (u>0..and.v==0.) dir=270 if (u<0..and.v==0.) dir=90 if (u==0..and.v==0.) dir=999 其中uv等于零的五种情况要单独挑出来，不然程序会有瑕疵。atan函数换成atand函数的话直接是度数，不用*180/pi 我四个象限和轴都试了，应该没错。最需要注意的问题，一个是函数值域，另一个是uv矢量方向和风向是反着的，并且风向角度数是从正Y轴开始顺时针算，和三角函数里度数从正X轴开始逆时针算不一样。

计算机程序设计语言(FORTRAN语言)

计算机程序设计语言(FORTRAN语言) (总分：36.00，做题时间：90分钟) 一、 (总题数：36，分数：36.00) 1.编译程序能将高级语言编写的源程序转换成( )。 A．解释程序 B．汇编程序 C．映象程序 D．目标程序（分数：1.00） A. B. C. D. √ 解析： 2.一个完整的FORTRAN源程序( )。 A．至少包括一个主程序 B．至少包括一个主程序和一个子程序 C．由一个主程序与一个以上的子程序组成 D．由一个主程序与一个子程序组成（分数：1.00） A. √ B. C. D. 解析： 3.语句函数定义语句在程序内合法的位置是( )。 A．在程序块开头语句之后，END语句之前 B．在程序块开头语句之后，可执行语句之前 C．在说明语句之后，END语句之前 D．在说明语句之后，可执行语句之前（分数：1.00） A. B. C. D. √ 解析： 4.下列关于函数子程序虚实参数的错误说法是( )。 A．可以没有形参数 B．虚实结合的数组长度可以不同 C．实参表与虚参表类型可以不同 D．函数名可以作为虚参

（分数：1.00） A. B. C. √ D. 解析： 5.下列叙述中正确的是( )。 A．FORTRAN程序块中，无名公用语句只能有一个B．FORTRAN子程序中，至少应有一个RETURN语句C．FORTRAN程序块中，最后一行必须是END语句D．FORTRAN程序块中，必须有变量说明语句（分数：1.00） A. B. C. √ D. 解析： 6.运行下面的程序时得不到所需的结果，其主要原因是( )。INTEGER X(11) DATA X/9，8，7，6，5，4，3，2，1，0，-1/ DO 10 1=1，X(1) ，-1 WRITE(*，*)1.0/SQRT(25.0-REAL(X(1) )* * 2) 10 CONTINUE END A．没有给X(11) 赋初值 B．发生除以零的情况 C．发生负数开平方的情况 D．循环参数设置错误（分数：1.00） A. B. C. D. √ 解析： 7.下列数据中，不符合FORTRAN常量表示法的是( )。 A．-25.6 B．2.758D3 C．'FOOT"=' D．TRUE （分数：1.00） A. B. C. D. √ 解析：

计算机技术基础(Fortran)试题库

一、选择题 1、下列关于"SUBROUTINE MAP(X，Y)"语句行的叙述中，不正确的是( ) A) 这是子程序的第一个语句B) 字符串"MAP"是子程序名 C) 变量X是子程序的形参D) 子程序执行后，MAP将返回整型数据 2、FORTRAN表达式"2/4+"的值是( ) A) B) 1 C) D) 0 3、阅读下列FORTRAN程序： PI=3. WRITE(*,' ')PI END 程序运行后输出结果是( ) A) B) 3.1415 C) D) 4、圆的直径存放在整型变量K之中，下列计算圆面积的表达式中正确的是( ) A) *K*K/4 B) *(K*K/4) C) *(K/2)**2 D) *(K/2)*(K/2) 5、要判断"月收入M在2000元以上(含2000元)且5000元以下)不含5000元)"的职工，应该使用的逻辑表达式是( ) A) M .GE. 2000 .AND. M .LE. 5000 B) M .GE. 2000 .OR. M .LE. 5000 C) M .GE. 2000 .AND. M .LT. 5000 D) M .GE. 2000 .OR. M .LT. 5000 6、设FORTRAN子程序如下： SUBROUTINE SUB(I，X) B=I+1 X=X+B END 若主程序中的变量M，B已按I-N规则赋值，则调用上述子程序的语句中正确的是( ) A) CALL SUB(M，M) B) CALL SUB(M+4，B) C) CALL SUB(B，B) D) CALL SUB(M，B+4) 7、下列不是FORTRAN常量的是( ) A) , B) +00 C) 2/3 D) 'Very good!' 8、阅读下列FORTRAN程序： F= DO 20 K=3, 7, 3 DO 20 L=K-1, K

fortran语言内部函数

附录 FORTRAN 90标准函数符号约定： ●I代表整型;R代表实型;C代表复型;CH代表字符型;S代表字符串;L代表逻辑型;A代表数组;P代表指针;T代表派生类型;AT为任意类型。 ●s:P表示s类型为P类型(任意kind值)。s:P(k)表示s类型为P类型(kind值=k)。 ●[…]表示可选参数。 ●*表示常用函数。表1 数值和类型转换函数函数名说明 ABS(x)*求x的绝对值∣x∣。x:I、R,结果类型同x; x:C,结果:R AIMAG(x)求x的实部。x:C,结果:R AINT(x[,kind])*对x取整,并转换为实数(kind)。x:R, kind:I,结果:R(kind) AMAX0(x1,x2,x3,…)*求x1,x2,x3,…中最大值。x I:I,结果:R AMIN0(x1,x2,x3,…)*求x1,x2,x3,…中最小值。x I:I,结果:R ANINT(x[,kind])*对x四舍五入取整,并转换为实数(kind)。x:R, kind:I,结果:R(kind) CEILING(x)*求大于等于x的最小整数。x:R,结果:I CMPLX(x[,y][,kind]))将参数转换为x、(x,0.0)或(x,y)。x:I、R、C, y:I、R,kind:I,结果:C(kind) CONJG(x)求x的共轭复数。x:C,结果:C DBLE(x)*将x转换为双精度实数。x:I、R、C,结果:R(8) DCMPLX(x[,y])将参数转换为x、(x,0.0)或(x,y)。x:I、R、C, y:I、R,结果:C(8) DFLOAT(x)将x转换为双精度实数。x:I,结果:R(8) DIM(x,y)*求x-y和0中最大值,即MAX(x-y,0)。x:I、R, y的类型同x,结果类型同x DPROD(x,y)求x和y的乘积,并转换为双精度实数。x:R, y:R,结果:R(8)

Fortran基本用法小结档

Fortran基本用法小结目录：一、说明二、概述三、数据类型及基本输入输出四、流程控制五、循环六、数组七、函数八、文件一、说明本文多数内容是我读彭国伦《Fortran 95 程序设计》的笔记。只读到第九章，主要是3~9章，都是些最基本的用法（原书共16章）。这里主要摘录了我看书过程中总结的一些Fortran和C不同的地方，主要是语法方面。希望这份笔记能够给学过C但没有接触过Fortran的同学带去一些帮助。要想弄得更清楚些，推荐看一下原书，觉得作者真的写得很好，很清楚；如果有C语言的基础，看完前九章应该很快的，花一两天就行了。觉得如果耐心看完本文，基本功能应该也可以顺利用起来了。另外，由于我之前没有用过Fortran，这次为了赶文档看书又看得很粗浅，大多数东西看过之后都没来得及仔细想，只是按着作者的意思去理解。所以这份笔记还处于纸上谈兵的层次。如果有不妥的地方，希望大家指正。谢谢！文中蓝色的部分是程序代码，!后面的内容为注释。二、概述 1、名词解释 Fortran=Formula Translator/Translation 一看就知道有什么特色了：可以把接近数学语言的文本翻译成机械语言。的确，从一开始，IBM设计它的时候就是为了方便数值计算和科学数据处理。设计强大的数组操作就是为了实现这一目标。Fortran奠定了高级语言发展的基础。现在Fortran在科研和机械方面应用很广。 2、Fortran的主要版本及差别按其发展历史，Fortran编译器的版本其实很多。现在在广泛使用的是Fortran 77和Fortran90。Fortran 90在Fortran 77基础上添加了不少使用的功能，并且改良了77编程的版面格式，所以编程时推荐使用90。鉴于很多现成的程序只有77版本，有必要知道77的一些基本常识，至少保证能够看懂77程序。以下是77和90的一些格式上的区别。

Simple算法_fortran语言编写

Cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc Cccc This computer program was copied from the graduate student course program Cccc of the University of Minnesota. Part of it was re-formulated to meet the Cccc personal computer environment. Some inappropriate expressions were also Cccc corrected. The program is used only for teaching purpose. No part of it Cccc may be published. You may use it as a frame to re-develop your own code Cccc for research purpose. XJTU Instructor, 1995.11 **************************************************************************** *----------------------------MAIN PROGRAM----------------------------------- **************************************************************************** LOGICAL LSTOP COMMON/CNTL/LSTOP **************************************************************************** OPEN(08,FILE='teresul') CALL SETUP0 CALL GRID CALL SETUP1 CALL START 10 CALL DENSE CALL BOUND CALL OUTPUT IF(.NOT.LSTOP) GO TO 15 CLOSE(08) STOP 15 CALL SETUP2 GO TO 10 END *--------------------------------------------------------------------------- SUBROUTINE DIFLOW **************************************************************************** COMMON/COEF/FLOW,DIFF,ACOF **************************************************************************** ACOF=DIFF IF(FLOW .EQ.0.0)RETURN TEMP=DIFF-ABS(FLOW)*0.1 ACOF=0. IF(TEMP .LE. 0. ) RETURN TEMP=TEMP/DIFF ACOF=DIFF*TEMP**5 RETURN END *-------------------------------------------------------------------------- SUBROUTINE SOLVE ****************************************************************************

Fortran基本用法小结

【以下文字转载自D_Atmosphere 讨论区】【原文由superjyq 所发表】我的Fortran基本用法小结高级语言和算法组heavensky 目录：一、说明二、概述三、数据类型及基本输入输出四、流程控制五、循环六、数组七、函数八、文件一、说明本文多数内容是我读彭国伦《Fortran 95 程序设计》的笔记。只读到第九章，主要是3~9章，都是些最基本的用法（原书共16章）。这里主要摘录了我看书过程中总结的一些Fortran和C不同的地方，主要是语法方面。希望这份笔记能够给学过C但没有接触过Fortran的同学带去一些帮助。要想弄得更清楚些，推荐看一下原书，觉得作者真的写得很好，很清楚；如果有C语言的基础，看完前九章应该很快的，花一两天就行了。觉得如果耐心看完本文，基本功能应该也可以顺利用起来了。另外，由于我之前没有用过Fortran，这次为了赶文档看书又看得很粗浅，大多数东西看过之后都没来得及仔细想，只是按着作者的意思去理解。所以这份笔记还处于纸上谈兵的层次。如果有不妥的地方，希望大家指正。谢谢！文中蓝色的部分是程序代码，!后面的内容为注释。二、概述 1、名词解释 Fortran=Formula Translator/Translation 一看就知道有什么特色了：可以把接近数学语言的文本翻译成机械语言。的确，从一开始，IBM设计它的时候就是为了方便数值计算和科学数据处理。设计强大的数组操作就是为了实现这一目标。Fortran奠定了高级语言发展的基础。现在Fortran在科研和机械方面应用很广。 2、Fortran的主要版本及差别按其发展历史，Fortran编译器的版本其实很多。现在在广泛使用的是Fortran 77和Fortran90。Fortran 90在Fortran 77基础上添加了不少使用的功能，并且改良了77编程的版面格式，所

fortran文件操作

fortran文件操作 From: 《Fortran 95 程序设计》彭国伦 1.文件读取的概念：读取：“顺序读取”和“直接读取”。保存：“文本文件”和“二进制文件”。二进制文件：读取快，因为保存格式和数据在内存中的保存方法一样，同时也节省空间。－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－ 2. 文件的操作：－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－－ open的使用：使用open命令打开文件之后，就可以对文件来做输入输出。 example: program ex0901 impicit none open(unit=10, file='hello.txt') ! 打开hello.txt文件, unit指定文件代码，file 指定文件名称。 write(10, *) "hello" !在代码为10的文件中写入hello stop end program ex0901 open中有很多参数可以使用，详细如下： OPEN(UNIT=number, FILE='filename', FORM='...', STATUS='...', ACCESS='...', RECL=length, ERR=label, IOSTAT=iostat, BLANK='...', POSITION='...', ACTION=action, PAD='...', DELIM='...') UNIT='number'： number必须是一个正整数，它可以使用变量或是常量来赋值。number最好避开1,2,5,6。因为2,6是默认的输出位置，也就是屏幕。1，5则是默认的输入位置，键盘。 FILE='filename'：指定要打开的文件名称，文件名要符合系统规定。windows下不区分大小写，unix下则会区分大小写，最好不要使用中文文件名。 FORM='FORMATTED' OR 'UNFORMATTED' FORM字段只有两个值可以设置： FORM='FORMATTED' “文本文件”格式来保存 FORM='UNFORMATTED' “二进制文件”格式保存这一栏不给定时候的默认值是: FORM='FORMATTED' STATUS='NEW' or 'OLD' or 'SCRATCH' or 'UNKNOWN' 用来说明打开一个新的文件或已经存在的旧文件。 STATUS='NEW' 打开一个原本不存在的新文件 STATUS='OLD' 打开一个原来已经存在的文件 STATUS='REPLACE' 若文件已经存在则重新创建一次，原来的内容消失；若不存在则会创建新文件。 STATUS='SCRATCH' 表示要打开一个暂存文盘，这个时候可以不需要指定文件名称，也就是FILE这个一栏可以忽略。因为程序本身会自动取一个文件名，至于文件名是啥也不重要，因为暂存盘会在程序结束后自动删除。 STATUS='UNKNOWN' 由各编译器自定义。通常会同REPLACE的效果。！这一栏不给定时，默认为STATUS='UNKNOWN'。 ACCESS='SEQUENTIAL' or 'DIRECT' 设置读写文件的方法：

Fortran语言

#LOCAL REAL Uds,Uqs,angles,vs,w,o,ang,ta,tb,t0,Ud1,Uq1,E #LOCAL REAL tc,t,mt,pi #LOCAL INTEGER sign,M,nl,n t=TIME tc=$TC nl=DNINT(tc/delt) pi=3.1415926 E=$Udc w=2*pi*50 if(t.lt.delt) then sign=0 M=0 Uds=0 Uqs=0 angles=0 vs=0 o=0 ang=0 ta=0 tb=0 t0=0 Ud1=0 Uq1=0 n=1 end if if(mod(n,nl).lt.1E-7) then !初始化变量 M=M+1 Uds=$Ud Uqs=$Uq angles=$angle !变换到alphabeta坐标直角坐标

Ud1=cos(angles)*Uds+sin(angles)*Uqs Uq1=sin(angles)*Uds-cos(angles)*Uqs ! Ud1=$Ud ! Uq1=$Uq !生成矢量的模和角度 vs=sqrt(Ud1**2+Uq1**2) if(Ud1.gt.0.and.Uq1.ge.0) then ang=atan(Uq1/Ud1) else ang=atan(Uq1/Ud1)+2*pi end if if(Ud1.lt.0) then ang=atan(Uq1/Ud1)+pi end if if(Ud1.eq.0) then if(Uq1.ge.0) then ang=pi/2 else ang=3*pi/2 end if end if !判断矢量的所属,计算角度 if(ang.ge.0.and.ang.lt.pi/3) then sign=1 o=ang ! o=w*t end if if(ang.ge.pi/3.and.ang.lt.2*pi/3) then sign=2 o=ang-pi/3 ! o=w*t-pi/3 end if if(ang.ge.2*pi/3.and.ang.lt.3*pi/3) then sign=3 o=ang-2*pi/3 ! o=w*t-2*pi/3

Fortran95第9章基本计算(一)赋值语句

第9章基本计算(一)赋值语句从计算机的实际状态来看，给出表达式并不意味着计算的开始，真正能够驱动计算的是可执行语句，而最直接驱动表达式的计算过程的就是赋值语句，完成一个赋值步骤，就意味着机器的状态的局部或全局发生了一个根本的变化。从语法的角度来看，一个表达式还只是一个比较复杂的复合词汇，还不能构成一个完整的语句，而赋值语句则是一个在语法意义上的完整的语句，表达式在赋值语句中扮演关键的语法作用。从计算问题的角度来看，赋值可以说是基本的计算步骤，考虑一下我们人工计算的过程就可以发现，任何一个完整的计算步骤，都可以说就是一个赋值步骤，特别是公式演算和数值计算，完成一个局部运算过程的标志，往往就是求出某个中间变量的数值或表达式。 FORTRAN语言作为一种以公式翻译为初衷的高级语言，它的赋值语句的一般形式就是一个数学等式，当然不是那种公式恒等变换得到的恒等式，而是要把运算的结果赋予一个变量，因此它的一般句法形式为： variable = expression variable => expression 可以看到一个赋值语句作为一个完整的语句，分为三个部分： ● 被赋值的变量 ● 赋值符号 ● 计算(表达式) 其中被赋值的变量，通过赋值，就拥有了明确的取值形式，如果该变量含有下标，片断下标或子串范围，则赋值的前提是它们都已经获得具体的取值。基于表达式结果的不同种类，赋值语句分为： ●固有赋值：固有赋值是把任意类型的值赋予一个非指针变量，或把一个指针变量赋予一个与之相关联的目标。 ●自定义赋值：自定义赋值得以构成的前提是存在一个可访问的子例行程序，它包含一个具有ASSIGNMENT形式的赋值界面，其属性与自定义赋值语句里的变量和表达式的属性保持一致。 ●指针赋值：指针赋值把一个指针变量关联到它的目标对象，或者说把一个目标对象赋予到一个指针变量。在形式上它使用符号=>。 ●过滤数组赋值：这个赋值过程是对满足一定条件的数组元素进行赋值，而不是对数组整体进行赋值。 ●并行指标数组赋值。这种赋值形式给出了一种有效的并行机制，能够大规模地对多重指标变量进行赋值。下面我们将分节讨论这5类赋值形式，首先给出一些例子如下：

Fortran95程序设计习题答案

Fortran95程序设计习题答案第四章 1.program main implicit none write(*,*) "Have a good time." write(*,*) "That's not bad." write(*,*) '"Mary" isn''t my name.' end program 2.program main real, parameter :: PI=3 implicit none.14159 real radius write(*,*) "请输入半径长" read(*,*) radius write(*,"(' 面积='f8. 3)") radius*radius*PI end program 3.program main implicit none real grades write(*,*) "请输入成绩" read(*,*) grades write(*,"(' 调整后成绩为 'f8.3)") SQRT(grades)*10.0 end program 4.integer a,b real ra,rb a=2 b=3 ra=2.0 rb=3.0 write(*,*) b/a ! 输出1, 因为使用整数计算, 小数部分会无条件舍去 write(*,*) rb/ra ! 输出 1.5 5.program main implicit none type distance real meter, inch, cm end type type(distance) :: d write(*,*) "请输入长度:" read(*,*) d%meter d%cm = d%meter*100 d%inch = d%cm/ 2.54 write(*,"(f8.3'米 ='f8.3'厘米 ='f8.3'英寸')") d%meter, d%cm, d%inch end program 第五章 1.program main implicit none integer money real tax write(*,*) "请输入月收入" read(*,*) money if ( money<1000 ) then tax = 0.03 else if ( money<5000) then tax = 0.1 else tax = 0.15 end if write(*,"(' 税金为 'I8)") nint(money*tax) end program 2.program main implicit none integer day character(len=20) :: tv write(*,*) "请输入星期几" read(*,*) day select case(day) case(1,4) tv = "新闻" case(2,5) tv = "电视剧" case(3,6) tv = "卡通" case(7) tv = "电影" case default write(*,*) "错误的输入" stop end select write(*,*) tv end program 3.program main implicit none integer age, money real tax write(*,*) "请输入年龄"

《FORTRAN 95程序设计》学习笔记

《FORTRAN 95程序设计》学习笔记 66RPG gg ★目录★ 《FORTRAN 95程序设计》学习笔记 (1) 基础知识(基础、字符串、FORMAT、隐式、TYPE) (1) 流程与控制(if、select、do) (4) 数组(声明、隐式循环、整体操作、可变数组) (5) 函数与子程序(子程序、函数、全局变量) (6) MODULE与面向对象(重载操作符、虚函数) (9) 文件相关(OPEN、WRITE、READ) (10) 指针（指向变量、数组、函数） (11) Visual Fortran 编译器（DLL，VB调用） (12) 数值算法与IMSL（数值算法插件） (14) 常用库函数（数学、数组、零碎、子程序） (15) 基础知识(基础、字符串、FORMAT、隐式、TYPE) ★【小玩意】二进制观察器：装在M.. Visual Studio\DF98\bin，有一个Bitviewer，可以观察变量储存方式 ★【语法】续行：行结尾或行开头使用& 符号；注释：使用! 符号 ★【语法】数学表达式：+ ；- ；* ；/ ；( ；) ；**乘幂 ★【语法】程序结束：STOP (Ruby的exit) ★【语法】输出：write(*,*)，完整写法：write(unit=*,fmt=*) ?建议：少用print，尽量用write ★【语法】声明 ?整型：integer(kind=4) a ；其中kind是使用的bytes数，4 or 2 ◆其他写法：integer*4 a; integer(4) a ?浮点：real(kind=4) a ；有效数位6位（12345678存为1.234567E7），如果是kind8 则为15位有效数字 ◆此外：1E10：单精10^10，1D10：双精10^10 ?复数：complex :: a=(2,3)

fortran初步学习资料

第一章: Fortran语言程序设计初步 Fortran语言的发展概况本节介绍Fortran的起源与发展历史，讲述Fortran由产生到形成标准FortranIV、Fortran77，并进一步形成新标准Fortran90/95的发展历程。 1.1.1 Fortran的历史 a)a)FortranI FortranIV Fortran是目前国际上广泛流行的一种高级语言，适用于科学计算。Fortran 是英文FORmula TRANslatio n的缩写，意为“公式翻译”。它是为科学、工程问题中的那些能够用数学公式表达的问题而设计的语言，主要用于数值计算。这种语言简单易学，因为可以像抄写数学教科书里的公式一样书写数学公式，它比英文书写的自然语言更接近数学语言。Fortran语言是第一个真正推广的高级语言。至今它已有四十多年历史，但仍历久不衰，始终是数值计算领域所使用的主要语言。Fortran语言问世以来，根据需要几经发展，先后推出形成了很多版本。第一代Fortran语言是在1954年提出来的，称为FortranI。它于1956年在IBM 704计算机上得以实现。在此之前编写计算机程序是极为繁琐的，程序员需要详细了解为之编写代码的计算机的指令、寄存器和中央处理器(CPU)等方面的知识。源程序本身是用数学符号(八进制码)编写的，后来采用了助记符，即所谓机器码或汇编码，这些编码由汇编程序转换为指令字。在50年代书写和调试一个程序要很长时间，因为用这种方式编写程序显然是很不方便的，尽管它能使CPU高效地工作。正是这些原因，促使由John Backus率领的IBM公司的一个小组研究开发最早的高级程序设计语言Fortran。其目的是开发一种容易理解、简单易学又能几乎像汇编一样高效运行的语言，他们取得了极大的成功。Fortran 语言作为第一种高级语言不仅是一次创新，也是一次革命。它使程序员摆脱了使用汇编语言的冗长乏味的负担，而且它使得不再只是计算机专家才能编写计算机程序，任何一名科学家或工程技术人员，只要稍加努力学习和使用Fortran，就能按自己的意图编写出用于科学计算的程序。经过不断发展，FortranI形成了很多不同版本，其中最为流行的是1958年出现的FortranII，它对FortranI进行了很多扩充(如引进了子程序)，FortranII 在很多机器上得以实现。其后出现的FortranIII未在任何计算机上实现。1962年出现的FortranIV对原来的Fortran作了一些改变，使得FortranII源程序在FortranIV编译程序下不能全部直接使用，导致了语言不兼容的问题。这样就形成了当时同时使用FortranII和FortranIV两种程序设计语言的局面。

FORTRAN经典入门程序20例

对于FORTRAN的初学者。这些例子可作为小练习。 1.例题：计算工资问题。每小时工资为RATE,如果工作超过40小时，加班呢部分工资是正常时间工资的1.5倍。 C Payroll with overtime program payroll real rate, hours, pay read (*,*) rate, hours if (hours>40) then pay=40*rate+(hours-40)*1.5*rate else pay=hours*rate END IF print *,"rate=", rate print *, "hours=", hours print *,"pay=",pay end 2.学生成绩问题。大于80为A级。大于60小于80为B级。小于60为C级。IF的嵌套。注意空格可以看清楚else if ,end if,print的内容. PROGRAM GRADESTUDENT READ *,GRADE IF (GRADE .GE. 80) THEN PRINT *,GRADE,"GRADE=>A" ELSE IF (GRADE.LT.60) THEN PRINT*,GRADE,"GRADE=>C" ELSE PRINT*,GRADE,"GRADE=>B" END IF END IF END 3.三个数按从小到大排序。 PROGRAM MAXMIN REAL A,B,C,T READ *,A,B,C IF (A.GT.B) THEN T=A A=B B=T ELSE END IF IF (B.GT.C) THEN T=B B=C

PROGRAM EQUATION READ (*,*) X IF (X .GE. 0.0) Y=3*X+6 IF (X .LT. 0.0) Y=-X**2+2*X-8 PRINT *,"X=",X, "Y=",Y END 6.CONTINUE语句。为了使循环的起止范围清晰，使用终端语句例：求5！ program ex1 FACT=1 do 10, N=1,5,1 ##DO 语句标号，循环不变量=初值，终值，步长FACT=FACT*N 10continue print*,FACT end 求1！+2！+3！+···+N！ READ *,N FACT=1

IntelVisualFortran安装使用详细讲解

I n t e l V i s u a l F o r t r a n安装+使用详细讲解好久都没有用Fortran了，这两周由于帮老师做事，就写了两个插值程序。刚刚去论坛里面搜了一下关于IVF的帖子，发现不是很多，个人感觉某些不怎么详细，从我第一次接触IVF 已经有三年了，下面我就把我对于IVF的认识和理解跟大家分享咯，求各位同仁不喜勿喷... 简介:IntelVisualFortran(IVF)和CompaqVisualFortran(CVF)类似，都是用来编写Fortran 语言的一个编译器，IVF看名字就知道，它是由Intel公司开发，而CVF则是由HP公司开发，之前看到一个帖子说，HP公司的CVF开发团队加盟到了Intel公司，或者可以说Intel 公司收购了HP公司的这个团队，而从2006年开始就不再开发新的CVF系列版本了，这么多年以来CVF也确实没有、......之后的版本出现，转而出现的确实IVF2010，IVF2011,IVF2013的新版本的出现。这就说明IVF是以后的主流趋势，而CVF只能成为我们的回忆，值得一提的是，IVF肯定是在CVF的基础上发展起来的，因此我们用到IVF的时候，会有一种熟悉的感觉。 Q1：为什么要用IntelVisualFortran？我第一次接触到IVF是当时我的电脑(Win7)不能正常使用CVF版本，会出现各种意想不到的问题(反正我是遇到过得，搞得编程的心情都没有了)，论坛里面也有一些什么破解的方法，都不管用，于是就在“小木虫”论坛里面看到了IVF，到目前为止，我用IVF在Win7的32位和64位的电脑上没有出现任何不兼容的问题（CVF是不支持64位机子的）。我自己的理解是新版本的出现，先不从稳定不稳定来说，肯定是有它的优势所在，Fortran标准的版本也在逐步的修改，比如将面向对象这些加入进来，以前版本的编辑器或许就不能支持了。因此，用IVF编写Fortran程序是一种与时俱进的趋势。 Q2：如何获取IntelVisualFortran？在百度或者家园论坛里面搜一下就会有很多相关下载的帖子，我自己也搜集了一些，下面就把我目前用到的版本给大家分享链接:提取密码:5p02 链接:c:\iknow\docshare\data\cur_work\_blank提取密码:ady4 第一个是IVF的安装包第二个就是IVF的licence文件（在安装过程中需要用到） Q3：如何在电脑中安装IntelVisualFortran？（建议电脑是XP系统的坛友，还是用CVF，电脑是Win7/8的坛友就可以安装IVF） IVF的安装比CVF的安装稍微麻烦一点的就是在安装过程这一步，在安装上面下载的IVF安装包之前，需要提前安装一个东西，那就是MicrosoftVisualStudio****(VS)这个集成开发环境(IDE)，可以简单的理解为，IVF就是VS里面的一个用于编写Fortran的插件，下面我也提供一个我的VS2010版本的下载。(市场上最新的是VS2013) VisualStudio2010简体中文专业版（MSDN原版下载内置KEY）.iso 链接:c:\iknow\docshare\data\cur_work\_blank提取密码:w7yv