FORTRAN实验报告
《FORTRAN语言》实验（上机）报告
FORTRAN 90 程序编程规范
FORTRAN 90 程序编程规范 Fortran 90 编程规范，使程序代码高度组织化，更加易读、易懂、易于维护，程序更加高效。使编出的程序更易懂、易于维护。 1 语言选择 数值预报创新系统软件开发应避免使用Fortran77 的某些过时特征以Fortran 90不一致的特征。选择Fortran 90 作为开发语言，并采用Fortran 90 的新功能，如动态内存的分配(dynamic memory allocation)、递归（recursion ）, 模块(modules)、POINTER 、长变量名、自由格式等。 Fortran 77其中某些只是一些冗余的功能，这些功能已经过时，另外，还有一些在Fortran90 中被证明是不好的用法，建议不要使用。 2 Fortran 90 的新特性 2.1.1 建议使用的Fortran 90 新特性 建议使用Fortran 90 提供的模块（module ），并用Use ONLY 指定module 中哪些变量或派生类型定义可用于调用程序。 尽量使用数组下标三元组，这样可优化并减少所需的代码行数。为提高可读性，要在括号内表明数组的维数，例如： 1dArrayA(:) = 1dArrayB(:) + 1dArrayC(:) 2dArray(: , :) = scalar * Another2dArray(: , :) 当访问数组的子集时，例如在有限差分等式中，可以通过使用下标三元组实现。例如：2dArray(: , 2:len2) = scalar *（ & Another2dArray(:, 1:len2 -1) & - Another2dArray(:, 2:len2) & ) 对程序单元（program units ）命名，并使用End program ，End subroutine ，End interface ，End module 等结构再次指定“program unit ”的名称。 在逻辑表达式中使用>、 >=、 ==、 <、 <=、 /=，它们分别代 替.gt.、.ge.、.eq.、.lt.、.le.、.ne. 。新的表示方法更接近标准的数学符号 在变量定义中始终使用“：：”；始终用“DIMENSION ”定义数组形状；始终用（len=）的语法格式声明字符变量的长度。
fortran语法手册
1 FORTRAN77四则运算符 + - * / ** (其中**表示乘方) 在表达式中按优先级次序由低到高为：+或-→*或/→**→函数→() 2 FORTRAN77变量类型 2.1 隐含约定：I-N规则 凡是以字母I，J，K，L，M，N六个字母开头的，即认为是整型变量，其它为实型变量。 如IMPLICIT REAL (I,J) 三种定义的优先级别由低到高顺序为：I-N规则→IMPLICIT语句→类型说明语句，因此，在程序中IMPLICIT语句应放在类型说明语句之前。 2.4 数组的说明与使用 使用I-N规则时用DIMENSION说明数组，也可在定义变量类型同时说明数组，说明格式为：数组名(下标下界,下标上界），也可省略下标下界，此时默认为1，例： DIMENSION IA(0:9),ND(80:99),W(3,2),NUM(-1:0),A(0:2,0:1,0:3) REAL IA(10),ND(80:99)使用隐含DO循环进行数组输入输出操作：例如WRITE(*,10) ('I=',I,'A=',A(I),I=1,10,2) 10FORMAT(1X,5(A2,I2,1X,A2,I4)) 2.5 使用DATA语句给数组赋初值 变量表中可出现变量名,数组名,数组元素名,隐含DO循环，但不许出现任何形式的表达式：例如 DATA A,B,C/-1.0,-1.0,-1.0/ DATA A/-1.0/,B/-1.0/,C/-1.0/ DATA A,B,C/3*-1.0/CHARACTER*6 CHN(10)
DATA CHN/10*' '/INTEGER NUM(1000) DATA (NUM(I),I=1,500)/500*0/,(NUM(I),I=501,1000)/500*1/ 3 FORTRAN77程序书写规则 程序中的变量名，不分大小写； 变量名称是以字母开头再加上1到5位字母或数字构成，即变更名字串中只有前6位有效； 一行只能写一个语句； 程序的第一个语句固定为PROGRAM 程序名称字符串 某行的第１个字符至第５个字符位为标号区，只能书写语句标号或空着或注释内容； 某行的第１个字符为Ｃ或＊号时，则表示该行为注释行，其后面的内容为注释内容； 某行的第６个字符位为非空格和非０字符时，则该行为上一行的续行，一个语句最多可有19个续行； 某行的第７至７２字符位为语句区，语句区内可以任加空格以求美观； 某行的第７３至８０字符位为注释区，８０字符位以后不能有内容。 4 FORTRAN77关系运算符 .GT. 大于 .GE. 天于或等于 .LT. 小于 .LE. 小于或等于 .EQ. 等于 .NE. 不等于 .AND. 逻辑与 .OR. 逻辑或 .NOT. 逻辑非 .EQV. 逻辑等 .NEQV. 逻辑不等 运算符优先级由高到低顺序为:()→**→*或/→+或-→.GT.或.GE.或.LT. 或.LE.或.EQ.或.NE.→.NOT.→.AND.→.OR.→.EQV.或.NEQV 5 FORTRAN77语句
工程分析程序设计上机作业(三)
工程分析程序设计 上机作业（三） 模块化编程（2） 上机目的：进一步掌握内部例程、外部例程、接口块、模块等功能的使用方法。联系例程重载、例程递归的使用方法。 1、 编写一个子例程子程序，SUM(S ，T ，N1，N2)。把整型数N1到N2进行求和，并把求和的结果放置到S ，把整型数N1到N2进行求积，并把求积的结果放置到T 。并用这个子程序来计算 (1234)(345678)(3*4*5*6)(1*2*3)y =++++++++++- 2、 编写一个模块程序，提供以下服务：定义出常量π、e 。定义出子程序，实现求和21n i n =∑、 求阶乘n !。并在主程序中计算如下结果：从键盘上输入整数n 、实型数A 、R 、R 0，求21! n i n n =?? ???∑（实型）和0202n n R R An R e R R π??- ????? ???（实型）。 3、 编写函数子程序GDC 求两个数的最大公约数。求最大公约数的算法如下：把两个数中大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。 调用此函数，求1260，198，72三个数的最大公约数。 选作： 编写一个函数子程序求sinh(x)，以下分三个步骤完成：（1）用递归算法，求出1!1!n n x x x n n n -=?-；（2）自行编程求出234012!3!4!!n x n x x x x e x n ∞=≈+++++???=∑，要求计算精度是610!n x n -<[不能用Fortran 的内部函数EXP(x)]；（3）计算sinh()2x x e e x --=
fortran程序案例题汇编(14道)
1．Fibonacci数列定义如下： F1=1 F2=1 F n=F n-1+F n-2 (n>2) 求Fibonacci数列的前30项。 integer f(30),i f(1)=1 f(2)=2 do i=3,30 f(i)=f(i-1)+f(i-2) enddo print*,f end 2．输入10个学生的总分，求每个学生的名次integer s(10),a(10),i,j do i=1,10 read*,s(i) enddo do i=1,10 a(i)=1 do j=1,10 if(s(i) 3．给定一组数，按照从小到大的顺序输出。integer a(10) integer p do i=1,10 read *,a(i) enddo do j=1,9 p=j do i=j+1,10 if (a(i) 4．输入若干名学生的学号和三门课程（语数英）的成绩，要求从键盘输入一个学生的学号，能打印出该学生的三门功课成绩和总分。 character*6,dimension(:),allocatable::xue integer,dimension(:,:),allocatable::g integer,dimension(:),allocatable::zong integer i,j,n character*6,xh print *,"请输入学生的个数" read *,n allocate(xue(n)) allocate(g(n,3)) allocate(zong(n)) do i=1,n read *,xue(i),(g(i,j),j=1,3) enddo do i=1,n zong(i)=0 do j=1,3 zong(i)=zong(i)+g(i,j) enddo enddo print *,"请输入你要打印的学生的学号" read *,xh do i=1,n if(xue(i)==xh)then print *,(g(i,j),j=1,3),zong(i) exit endif enddo end Fortran 复习题汇总program exer101_1 implicit none realpi,r,t,s C 说明符号常量 parameter(pi=3.1415926) C 打开用于保存结果的文件 open(9,file='mydata.dat') write(*,*)"请输入r和t的值：" read*,r,t s=t/360.0*pi*r**2 C 在屏幕上显示结果 print*,'r=',r,'t=',t,'s=',s C 采用格式说明语句显示结果 write(*,100)'r=',r,'t=',t,'s=',s 100 format(1X,3(A2,F10.5,1X)) C 将结果保存至文件中 write(9,*)'r=',r,'t=',t,'s=',s end program exer101_2 implicit none realpi.r,t,s parameter(pi=3.1415926) !说明符号常量 open(9,file='mydata.dat') !打开用于保存结果的文件write(*,*)"请输入r和t的值：" read*,r,t s=t/360.0*pi*r**2 print*,'r=',r,'t=',t,'s=',s !在屏幕上显示结果 write(*,100)'r=',r,'t=',t,'s=',s !采用格式说明语句显示结果100 format(1X,3(A2,F10.5,1X)) write(9,*)'r=',r,'t=',t,'s=',s !将结果保存至文件中 end 1 2 program ex202 implicit none real F,C print*,"请输入一个华氏温度：" read*,F 实验一 上机日期:2012年2月23日 1.实验目的: (1)了解FORTRAN90与软件开发环境的关系. (2)掌握FORTRAN90上机实操作过程. (3)掌握FORTRAN90软件开发环境启动方法. (4)理解有关工作区,项目,文件创建方法. (6)熟悉软件开发环境图形用户界面. (7)掌握程序输入,编译,运行方法. (8)了解软件开发环境常用菜单,工具按钮,环境窗口基本功能. 2上机过程记录 (1)主要记录输入程序的原代码: Program main Real length,width,radius Real ::unitprice=0.5,totalprice,area Print*,’请输入图形长宽值:’ Read*,length,width Print*”请输入绿地半径:” Read*,radius Area=length*width-3.1415926*radius**2 Totalprice=unitprice*area Print*,’地块总价为:’totalprice’万元’ End (2)调试过程: D:\shiyan01\xm1\chengxul.f90 D:\shiyan01\xm1\chengxul.f90(6) : Error: Syntax error, found CHARACTER_CONSTANT '请输入绿地半径:' when expecting one of: , 1)实例3—求多个半径下的圆周长 ! z3.f90 --Fortran95 ! FUNCTIONS: ! z3 - Entry point of console application. !************************************************************************* ! PROGRAM: z3 ! PURPOSE: Entry point for the console application. !************************************************************************ program z3 ! 求多个半径下的圆周长 ! 主程序 ! PROGRAM Z3 PRINT *, 'R=',1.2,'C=',C(1.2) PRINT *, 'R=',3.4,'C=',C(3.4) PRINT *, 'R=',15.6,'C=',C(15.6) PRINT *, 'R=',567.3,'C=',C(567.3) END program z3 !子程序 FUNCTION C(R) PI=3.1415926 C=2*PI*R RETURN ! Body of z3 end 2)实例4—键盘与显示器输入/输出 a)Fortran 基本操作 b)程序指令 ! ZXZ_I_O.f90 ! FUNCTIONS: ! ZXZ_I_O - Entry point of console application. ! PROGRAM: ZXZ_I_O ! PURPOSE: Entry point for the console application. !***************输入、输出样式种种************************** program ZXZ_I_O implicit none !变量声明的位置 INTEGER(2) i; INTEGER(4) j; INTEGER(4) m; REAL n INTEGER A,B ! Variables PRINT*,'输入整数A'; READ*, A PRINT*,'输入整数B'; READ*, B B=A+B PRINT*,'B=A+B=',B WRITE(*,*) 'A*B=',A*B PRINT* ,'以上为计算机的计算结果，注意B的值' 1、RANDOM_NUMBER Syntax ['sint?ks] n. 语法 CALL RANDOM_NUMBER (harvest结果) Intrinsic Subroutine（固有子程序）：Returns a pseudorandom number greater than or equal to zero and less than one from the uniform distribution. 返回大于或等于0且小于1，服从均匀分布的随机数 2、RNNOA/ DRNNOA （Single/Double precision） Generate pseudorandom numbers from a standard normal distribution using an acceptance/rejection method. 产生服从标准正态分布的随机数 Usage（用法） CALL RNNOA (NR, R) Arguments（参数） NR— Number of random numbers to generate. (Input) 要产生随机数的个数 R— Vector of length NR containing the random standard normal deviates. (Output) 输出长度为NR，随机正态分布的向量 Comments（注解） The routine RNSET can be used to initialize the seed of the random number generator. The routine RNOPT can be used to select the form of the generator. 程序RNSET可以用来初始化随机数发生器的种子 Example In this example, RNNOA is used to generate five pseudorandom deviates from a standard normal distribution. INTEGER ISEED, NOUT, NR REAL R(5) EXTERNAL RNNOA, RNSET, UMACH C CALL UMACH (2, NOUT) NR = 5 ISEED = 123457 CALL RNSET (ISEED) CALL RNNOA (NR, R) WRITE (NOUT,99999) R 99999 FORMAT (' Standard normal random deviates: ', 5F8.4) END Output Standard normal random deviates: 2.0516 1.0833 0.0826 1.2777 -1.2260 工程分析程序设计上机作业（五） 模块化编程 上机目的：掌握内部例程、主程序、外部例程、模块等功能的使用方法。 1、有一个六边形，求其面积。为求面积，作了3条辅助线。如图所示：（提示，三角形面 积 2 a b c s ++ =，a、b、c为三个边长）。要求用内部函 数来计算每个三角形的面积。 =21 program main implicit none real :: A(9)=(/10,30,16,13,21,14,20,36,28/) integer i real s s=area(A(3),A(4),A(7))+area(A(2),A(7),A(8))+ & area(A(1),A(8),A(9))+area(A(5),A(6),A(9)) print*,'The area is',s contains function area(x,y,z) implicit none real x,y,z real p,area p=(x+y+z)/2 area=sqrt(p*(p-x)*(p-y)*(p-z)) end function end 2、编写一个子例程子程序，SUM(S，T，N1，N2)。把整型数N1到N2进行求和，并把求 和的结果放置到S，把整型数N1到N2进行求积，并把求积的结果放置到T。并用这个子程序来计算： y=++++++++++- (1234)(345678)(3*4*5*6)(1*2*3) program main implicit none integer A(8) integer i,b,c,d,e,f,g,m,n,y b=0;d=0;g=1;n=1 do i=1,8 A(i)=i end do call sum(b,c,1,4) call sum(d,e,3,8) call sum(f,g,3,6) call sum(m,n,1,3) y=b+d+g-n print*,'y=(1+2+3+4)+(3+4+5+6+7+8)+(3*4*5*6)-(1*2*3)=',y contains subroutine sum(s,t,N1,N2) integer s,t,N1,N2,i do i=N1,N2 s=s+A(i) t=t*A(i) end do end subroutine end 3、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中 大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。 不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。 “专家系统”实验报告 专业：智能科学与技术 班级：1102 学号：0909113310 姓名：宋佳欣 2014年6月5日 实验1 专家系统与CLIPS 实验内容提要： 专家系统不仅是具有特定的结构，同时，它还有多种知识表示形式，及其相应的知识推理方法。因此，直接采用程序设计语言实现专家系统是比较困难的，而采用适当的专家系统开发工具则可以大大降低系统开发的难度。因此，本章结合麦粒肿系统专断专家系统开发的需要，首先安排了对专家系统及其开发工具CLIPS基本理论知识，以及CLIPS基本理论知识，以及CLIPS基本使用方法的概要介绍，目的是为第6章“麦粒肿诊断专家系统开发”做好准备。 实验重点： 专家系统的一般结构 产生式规则的表示及其正向推理 专家系统开发工具的结构和特点 专家系统开发工具CLIPS的基本语法 专家系统开发工具CLIPS的基本使用 实验目的 掌握专家系统的一般结构 掌握产生式规则的表示及其正向推理 掌握专家系统开发工具的结构和特点 了解专家系统开发工具CLIPS的设计目标与特点 掌握专家系统开发工具CLIPS的基本语法 掌握专家系统开发工具CLIPS的基本使用 实验内容及要求 学习专家系统的一般结构 学习产生式规则的表示及其正向推理 学习专家系统开发工具的结构和特点 学习专家系统开发工具CLIPS的设计目标与特点 学习专家系统开发工具CLIPS的基本语法 学习专家系统开发工具CLIPS的基本使用 实验步骤 本实验包括以下六个具体步骤，分别是：学习专家系统的一般结构、学习产 生式规则的表示及其正向推理、学习专家系统开发工具的结构和特点、学习专家系统开发工具CLIPS的设计目标与特点、学习专家系统开发工具CLIPS的基本语法，以及学习专家系统开发工具CLIPS的基本使用。 1 学习专家系统的一般结构： 专家系统（ES，Expert System）是利用大量专家知识，运用只是推理方法解决特定领域中实际问题的计算机程序系统。专家系统由产生式规则专家系统发展而来，主要包括知识库、推理机、综合数据库、知识获取程序、人机接口，以及解释程序六大部分。其中，知识库、推理机、综合数据库，以及人际接口是目前大多数专家系统的主要内容，而知识获取程序和解释则是所有专家系统都期望具有的两个模块。在专家系统的一般结构中，知识库用于存放各领域专家提供的经验性知识：综合数据库用于存放关于问题求解的初始数据、中间结果、假设、目标，以及最终求解结果（事实）；推理机在一定的控制策略下针对综合数据库中的当前信息，识别和选取知识库中对当前问题求解有用的知识进行推理；知识获取程序完成把专家只是按一定表示形式输入到知识库中，部分代替知识工程师进行知识的自动获取，实现专家系统的字学习，不断完善知识库；解释程序可以根据用户提问，对系统给出的结论、求解过程，以及系统当前求解状态提供说明，便于用户理解系统的问题求解，增加用户对系统求解结果的信任程度，在知识库的完善过程中，便于专家或知识工程师发现和定位知识库中的错误，便于各领域的专业人员或初学者能够从问题的求解过程中得到直观学习；人机接口将专家爱或用户的输入信息翻译为系统可接受的内部形式，同时，也把系统向专家或用户输出的信息转化成人类易于理解的外部形式。 2 学习产生式规则的表示及其正向推理 目前，用于专家系统的知识表示形式主要有谓词逻辑、产生式规则、语义网络、框架，以及剧本。其中，采用产生式规则知识表示形式建立的专家系统最为流行和广泛，这主要因为：首先，产生式规则知识表示容易被人理解；其次，它基于演译推理，保证了推理结果，的准确性；最后，大量产生式规则可以连接成多棵推理树，树的宽度反映了实际问题的范围，树的深度反映了问题的难度。 产生式规则一般表示为：如果A成立，则B成立，简记为A->B,其中A是前提，B是结论。一般来说，不同的产生式规则之间都是有联系的，即谋条产生式规则的前提实力外一条产生式规则的结论。如果把规则库（知识库）中所包含的总目标（某些规则的结论）作为根节点，按规则的前提和结论展开形成树的形式，则形成的树就被称为推理树或知识树，它把规则库中的所有规则都连接了起来。 产生式规则的推理分为正向推理和逆向推理两种，这里主要介绍正向推理。 F O R T R A N 第2章FORTRAN90基础知识： 程序单元的概念: fortran90程序是一种分块形式的程序，整个程序由若干程序单元块组成。每个程序只有一个主程序单元。各单元体的程序体形式上相同。 程序单元可以是主程序、子程序（外部过程或内部过程）、模块MODULE （供其他程序单元引用即把该程序单元内的全部语句复制到程序单元中）或块数据程序单元BLOCK 。 语言元素：指在程序设计中用到的基本成分，如字符集、常量、变量、记号（标号、关键字、名字、常数、运算符和定界符）以及其他的基本结构等。 字符集：英文字母、阿拉伯数字、下划线、21个特殊字符（$和？没有规定用法）。 数据结构： 整型INTEGER （34-2下划线后面是种别参数），n 位字长的计算机的数据表示范围一般为12~211-+---n n ，种别参数由内在询问函数KIND 送回，种别值提供的最小范围由内在函数SELECTED-INT-KIND 送回，表示十进制幂的范围由内在函数RANGE 送回； 实型REAL ，小数形式和指数形式；复型COMPLEX （种别类 型参数值取实部与虚部中较大者）； 字符型CHARACTER，由一对单撇号或双撇号之间的字符序列组成； 逻辑型LOGICAL。 派生数据类型TYPE； 数组INTEGER,DIMENSION(1,50)::A，可直接对数组元素进行运算如3*A,SQRT(A)； 字符子串，在字符串CHARACTER(LEN=80)::ROW中，ROW(1:3)就表示字符串ROW中第1到第3个元素组成的子串。 变量名的命名规则：不能超过31个字符；组成成分是字母、数字、下划线；第一个字符必须是字母。 第3章基本语句： 类型说明语句：类型说明（种别说明），属性说明：：变量名表 尽量避免把两个实数作相等或不相等的比较。淘汰隐式说明IMPLICIT NONE 种别说明：种别参数即对可移植数据精度和范围进行选择的机制 KIND(X) 返回变元X的种别参数值 SELECTED-REAL-KIND(n,m) 产生一个种别值，它表示某一精度和范围。N指十进制有效位数，m指明值范围内以10为底的幂次。 计算圆周率 REAL R,R1,R2,PI ISEED=RTC() N0=0 N=300000 DO I=1,N R1=RAN(ISEED) R2=RAN(ISEED) R=SQRT(R1*R1+R2*R2) IF(R<1.0)N0=N0+1 END DO PI=4.0*N0/N WRITE(*,*)PI END 一）蒙特卡洛计算生日问题 假设有N个人在一起，各自的生日为365天之一，根据概率理论，与很多人的直觉相反，只需23个人便有大于50％的几率人群中至少有2个人生日相同。 INTEGER M(1:10000), NUMBER1(0:364), NUMBER2 REAL X,Y ISEED=RTC() DO J=1, 10000 NUMBER1=0 X=RAN(ISEED) NUMBER1(0)=INT(365*X+1) JJJ=1 DO I=1,364 Y=RAN(ISEED) NUMBER2=INT(365*Y+1) ETR=COUNT(NUMBER1.EQ.NUMBER2) IF (ETR= =1) THEN EXIT ELSE JJJ=JJJ+1 M(J)=JJJ NUMBER1(I)=NUMBER2 END IF END DO END DO DO I=1,10000 IF(M(I).LE.23) SUM=SUM+1 END DO PRINT *,SUM/10000 END 二）MONTE CARLO SIMULATION OF ONE DIMENSIONAL DIFFUSION 蒙特卡罗计算一维扩散问题 INTEGER X,XX(1:1000,1:1000) REAL XXM(1:1000) ! X:INSTANTANEOUS POSITION OF ATOM ! XX(J,I)：X*X ,J:第几天实验,I:第几步跳跃 ! XXM(I): THE MEAN OF XX WRITE(*,*) "实验天数JMAX，实验次数IMAX" READ(*,*) JMAX,IMAX ISEED=RTC() DO J=1,JMAX !第几天实验 X=0 !!! DO I=1,IMAX !第几步跳跃 RN=RAN(ISEED) IF(RN<0.5)THEN X=X+1 ELSE X=X-1 END IF XX(J,I)=X*X END DO END DO OPEN(1,FILE="C:\DIF1.DAT") DO I=1,IMAX XXM=0.0 XXM(I)=1.0*SUM(XX(1:JMAX,I))/JMAX !! WRITE(1,*) I, XXM(I) END DO CLOSE(1) END 三维的！ 三）通过该程序了解FORTRAN语言如何画图（通过像素画图）USE MSFLIB INTEGER XR，YR !在的区域中画一个圆 PARAMETER XR=400,YR=400 INTEGER R, S(1:XR,1:YR) X0=XR/2 ! 圆心位置X0，YO Y0=YR/2 R=MIN(X0-10,Y0-10) !圆半径 S=0 !像素的初始状态（颜色） DO I=1,XR DO J=1,YR IF((I-X0)**2+(J-Y0)**2<=R**2)S(I,J)=10 IER=SETCOLOR(S(I,J)) IER=SETPIXEL(I,J) END DO END DO END Read 的规则： 按行读取，每次读一行，读完后光标自动跳到下一行的开头，空格和逗号代表结束（所以空格和逗号也是读取的一部分的话，需要使用“输入输出格式”） 如果想要将一行数据读入数组，代码为： Read(10,*) s(:,:,:) 不用规定输入输出格式，因为会根据s(:,:,:)确定需要读入的数字的个数，然后fortran会按部就班的读取，甚至文件中当前这一行的数字个数不足以填满s(:,:,:)的时候，read会自动跳到下一行继续搜索数字，直到填满s(:,:,:)为止。 但是需要注意给数组赋值的顺序：read会把它搜索到的第一个数字给s(1,1,1),第二个给s(2,1,1),第三个给s(3,1,1)… 程序9 1: 将read(unit=field,fmt="(A79)",iostat=status)buffer 中的A79改为A2，结果只输出每行的前两个字符，说明read是按整行整行读取的。 中间空了几行之后，空行之后的内容还是能被读取和输出，这说明，空行和空白是不一样的：空行也算是一种文本内容，因此不会终止读取，而空白意味着结束。 !读取文件 program main implicit none character(len=79)::filename="number.txt",buffer integer,parameter::field=10 integer::status=0 logical alive inquire(file=filename,exist=alive) if(alive)then open(unit=field,file=filename) do while(.true.) read(unit=field,fmt="(A79)",iostat=status)buffer if(status/=0)exit write(*,"(A79)")buffer end do else write(*,*)filename,"does't exist." end if pause stop end program main ============================================= 附number.txt =============================== 1234555666 《复合材料力学》课程上机指导书（力学101，力学C101-2） 河北工业大学机械学院力学系 2013年9月 目录 作业1单向板刚度及柔度的计算(2学时) (1) 作业2单向板的应力、应变计算(2学时) (2) 作业3绘制表观工程常数随 的变化规律(3学时) (3) 作业4绘制强度准则的理论曲线(包络线)(3学时) (4) 作业5层合板的刚度计算(3学时) (5) *作业6层合板的强度计算(4学时) (6) 附录作业提交说明……………………………………………. . 7 注：带“*”的题目可作为自愿选做题。 作业1 单向板刚度及柔度的计算 一、要 求 （1）选用FORTRAN 、VB 、MAPLE 或MATLAB 编程计算下列各题； （2）上机报告内容：源程序、题目内容及计算结果； （3）材料工程常数的数值参考教材自己选择； （4）上机学时：2学时。 二、题 目 1、已知单层板材料工程常数1E ，2E ，12G ，计算柔度矩阵[S ]和刚度矩阵[Q ]。（玻璃/环氧树脂单层板材料的MPa 1090.341?=E ，MPa 1030.142?=E ， MPa 1042.0412?=G ，25.021=μ， MPa 1001=σ，MPa 302-=σ，MPa 1012=τ） 2、已知单层板材料工程常数1E ，2E ，12G ，21μ及θ，计算柔度矩阵][S 和刚度矩阵][Q 。（M P a 1090.341?=E ，MPa 1030.142?=E ，MPa 1042.0412?=G ，25.021=μ，?=30θ） 作业2 单向板的应力、应变计算 一、要 求 1、选用FORTRAN 、VB 、MAPLE 或MATLAB 编程计算下列各题； 2、上机报告内容：源程序、题目内容及计算结果； 3、材料工程常数的数值参考教材自己选择； 4、上机学时：2学时。 二、题 目 1、已知单向板的应力x σ、y σ、xy τ，工程常数1E ，2E ，12G ，21μ及θ，求x ε、 y ε、xy γ；1σ、2σ、12τ；1ε、2ε、12γ。 （知?=30θ，应力MPa 160=x σ，MPa 60=y σ，MPa 20=xy τ，工程常数MPa 1090.341?=E ，MPa 1030.142?=E ，MPa 1042.0412?=G ，25.021=μ，?=30θ） 2、已知1σ、2σ、12τ，工程常数1E ，2E ，12G ，21μ及θ，求1ε、2ε、12γ；x ε、y ε、 xy γ；x σ、y σ、xy τ。 （知MPa 1001=σ，MPa 302-=σ，MPa 1012=τ，MPa 1090.341?=E ，MPa 1030.142?=E ，MPa 1042.0412?=G ，25.021=μ，?=30θ） 蒙特卡罗方法 实验一 实验报告 蒙特卡罗方法实验一实验报告 一、实验目的 1、了解蒙特卡罗方法方法的基本思想； 2、掌握蒙特卡罗方法计算面积、体积的方法； 3、掌握由已知分布的随机抽样方法。 二、实验原理 Monte Carlo 方法，又称统计模拟方法或计算机随机模拟方法，是一种基于“随机数”进行数值模拟的方法，一种采用统计抽样理论近似求解物理或数学问题的方法。 倘若待求量可以表述成某些特征量的期望值、某些事件出现的概率或两者的函数形式，那么可采用蒙特卡罗方法求解。在求解某些特征量的期望值或某些事件出现的概率时，必须构建合符实际的数学模型。例如采用蒙特卡罗方法计算某函数所围面积时，构建的数学模型是构造一已知面积的可均匀抽样区域，在该区域投点，由伯努利定理大数定理可知，进入待求区域投点的频率依概率1收敛于该事件出现的概率（面积之比）。 由已知分布的随机抽样方法指的是由已知分布的总体中抽取简单子样。抽样方法有： 直接抽样方法：离散型分布随机抽样方法、连续型分布直接抽样方法；挑选抽样方法;复合抽样方法；随机抽样一般方法：加抽样方法、减抽样方法、乘抽样方法、乘加抽样方法、乘减抽样方法、对称抽样方法、替换抽样方法、多为分布抽样方法、积分抽样方法；随机抽样其他方法：偏倚抽样方法、近似分布抽样方法、近似-修正抽样方法。 三、实验内容 1、安装所需计算工具（MA TLAB 、fortran 、C++等）； 2、编写一伪随机数发生器；（如乘加同余a=1366,c=150889,M=714025、a=9301,c=49297,M=233280；乘同余a=16807,M=232 -1；或采用其它方法） 以下内容选取一个采用自编伪随机数发生器进行计算，其余采用工具软件中自带伪随机数发生器进行计算。 3、求解以下区域的面积、体积： 3.1、给定曲线y =2 – x 2 和曲线y 3 = x 2，曲线的交点为：P 1( – 1，1 )、P 2( 1，1 )。曲线围成平面有限区域，用蒙特卡罗方法计算区域面积； 3.2、计算22 22 11z x y z x y ?≥+? ?≤+--??所围体积 其中{(,,)|11,11,02}x y z x y z Ω=-≤≤-≤≤≤≤。 4、对以下已知分布进行随机抽样： 4.1、()() []2 3 321,0,12 f x x x x =+ -∈； 4.2、()() ()[]11,1,21E f x f x x E k E = ?∈+ 程序1 PROGRAM chaper1272 INTEGER :: A(10)=(/5,7,4,8,12,2,10,3,9,11/) INTEGER :: sum=0,ave !打开一数据文件,设置一个有格式直接存取文件,将10个数分2个记录写入文件。 !数据文件生成2个记录,每个记录行长度相同,记录长度为25 OPEN(1,FILE='input21.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=25) WRITE(1,"(5I5)",REC=1)(A(I)+10,I=1,5) !按格式说明将头5个数写入第1个记录 WRITE(1,"(5I5)",REC=2)(A(I)+10,I=6,10)!按格式说明将后5个数写入第2个记录 READ(1,"(5I5)",REC=2)(A(I),I=6,10) !按格式说明从第2个记录中读取后5个数 READ(1,"(5I5)",REC=1)(A(I),I=1,5) !按格式说明从第1个记录中读取头5个数 DO I=1,10 sum=sum+A(I) ENDDO ave=sum/10 !打开一个最大记录长度为22的有格式顺序存取文件 OPEN(2,FILE='input22.dat',FORM='FORMATTED',ACCESS='D IRECT',RECL=22) WRITE(2,"('10个数之和为：',I5)",REC=1) sum !输出1记录行,记录长度为22 WRITE(2,"('10个数平均值为：',I5)",REC=2)ave !输出1记录行,记录长度为22 WRITE(2,"(A)",REC=3) '程序运行正常结束。' !输出一个记录行,记录长度为22 END 程序2 PROGRAM average REAL sum,ave INTEGER n OPEN(1,file='score.dat') PRINT*,'正在统计平均成绩,请等待。' sum=0.0;n=0 DO READ(1,*,END=100) S sum=sum+s n=n+1 ENDDO 100 ave=sum/n FORTRAN第二次实习作业 2014301610307 -D-马莉1 、根据最新个人所得税计算方法，如下面的个税税率表。工资、薪金所得，适用七级超额累进税率，税率为百分之三（3% ）至百分之四十五（45% ）。3500 元为起征点，即只需要缴纳高于3500 元部分的个税。请根据下表编程，通过给定个人收入，计算其所需邀纳的个税。 program jiaoshui implicit none integer m real n character i 100 print*,'应缴纳所得税是否含税，请回答Y或N：' read*,i print*,'请输入工资所得：' read *,m if(i=='Y')then select case(m) case(0:1500) n=m*3/100-0 case(1501:4500) n=m*10/100-105 case(4501:9000) n=m*20/100-555 case(9001:35000) n=m*25/100-1005 case(35001:55000) n=m*30/100-2775 case(55001:80000) n=m*35/100-5505 case default n=m*45/100-13505 end select elseif(i==’N’)then select case(m) case(0:1455) n=m*3/100-0 case(1456:4155) n=m*10/100-105 case(4156:7755) n=m*20/100-555 case(7756:27255) n=m*25/100-1005 case(27256:41255) n=m*30/100-2775 case(41256:57505) n=m*35/100-5505 case default n=m*45/100-13505 end select else print*,"输入格式不正确，请重新输入：" goto 100 endif print*,n End 2、编写程序求某天在当年中是第多少天，例如2010 年1 月10 号是2010 年的第10 天。要求分别使用IF 、CASE结构来编写。program tianshu2 implicit none integer M1(12),i,y,m,d,sum,n data M1/31,28,31,30,31,30,31,31,30,31,30,31/ sum=0 100 print *,"请输入日期格式例如：2010 10 23" read 10,y,m,d 10 format(i4,1x,i2,1x,i2) if(m>12.or.d>31)thenFortran 上机题汇总
Fortran 90上级实验报告
fortran程序实例
Fortran常用函数
fortran习题参考答案
专家系统实验报告
FORTRAN知识点总结
几个fortran程序
fortran心得
2013秋复合材料力学上机作业
蒙特卡罗实验报告
fortran程序30个
Fortran作业2