Fortran90几个重要的实验报告.docx

1、Fortran90几个重要的实验报告Fortran90语言程序设计 上机实验报告河南科技大学规划与建筑工程学院班级： 工力112 姓名： 学号： 111409070210 成绩： 2013年 5月 19 日实验一（1） 熟悉FORTRAN90软件开发环境一 实验目的1 了解FORTRAN90与软件开发环境的关系。2 掌握FORTRAN90上机实验基本操作过程。3 掌握FORTRAN90软件开发环境启动方法。4 理解有关工作区、项目、文件的基本概念。5 掌握工作区、项目、文件创建方法。6 熟悉Compaq Visual Fortran软件开发环境图形用户界面。7 掌握程序的输入、编辑、编译、构建

2、、运行方法。8 了解软件开发环境常用菜单、工具按钮、环境窗口基本功能。9 掌握文件、文件夹的移动、复制、删除等操作方法。二 课堂实验内容1 问题描述现有长1500米，宽1000米的地块要拍卖建造住宅小区，地块中央有一半径为100米的圆形公共绿地不属于拍卖范围，每平方米地价为5000元。编写程序计算该地块拍卖总价，并输出。要求地块长宽和绿地半径数据从键盘输入。2 算法设计通过对该问题的分析研究，设计求解算法，并绘制流程图如图1-1所示结束图1-1 计算地块总价程序流程图3. 程序编写根据算法和流程图，设计和编写程序如下： !班级：工力112 !姓名： !时间：201304 !计算加油站月销售额

3、PROGRAM main PARAMETER (pi=3.1415) INTEGER : number REAL ： cubage,volumes,sale REAL ： radius=5.0,heigh=3.0,unitprice=500.0 PRINT*,请输入月销售油罐数量： READ*，number cubage=pi*radius*radius*heigh volumes=cubage*number sale=volumes*unitprice PRINT*, 加油站月销售额：，sale，元。 END4. 运行结果 请输入月销售油罐数量:4加油站月销售额： 471225.0 元。Pr

4、ess any key to continue 三课外练习1 问题描述. 1有一半径为5，高为20的圆柱型物体，求圆柱体积，并输出。编写程序实现之。 2有一半径为5，高为20的圆锥型物体，求圆锥体积，并输出。编写程序实现之。2 算法设计结束3. 程序编写程序2： !班级：工力112 !姓名： !时间：2101304 !计算圆锥体积 PROGRAM main PARAMETER (pi=3.1415) REAL ： r,h,V PRINT*,请输入圆锥底半径和高： READ*，r,h V=pi*r*r*h/3 PRINT*, 圆锥体积：，V END根据流程图，分别编写求解上述两问题的程序：程序1

5、、程序2。程序1： !班级：工力112 !姓名： !时间：201304 !计算圆柱体积 PROGRAM main PARAMETER (pi=3.1415) REAL ： r,h,V PRINT*,请输入圆柱底半径和高： READ*，r,h V=pi*r*r*h PRINT*, 圆柱体积：，V END4. 运行结果请输入圆柱底半径和高：5，20 圆柱体积： 1570.750Press any key to continue请输入圆锥底半径和高：5，20 圆锥体积： 523.5833Press any key to continue实验一（2） 掌握FORTRAN90基础知识一 实验目的1 了解

6、FORTRAN90字符集。2 掌握FORTRAN90实体名称的语法描述规则和具体使用要求。3 了解FORTRAN90关键字意义，以及与名称的区别。4 掌握FORTRAN90程序组织结构和语句排列次序的基本原则。5 掌握FORTRAN90自由书写格式及三个标志符的主要作用。6 理解数据、数据表示、数据类型概念，熟悉FORTRAN90数据类型。7 进一步熟二 课堂实验内容1 问题描述 有一呈等腰梯形的鱼塘，上底边长为150米，下底边长为210米，高为40米，在鱼塘周围建起了围墙，如图2-1所示。每平方米鱼塘产鱼4公斤，每公斤鱼的价格为10元。编写程序计算该鱼塘的产量和产值，以及围墙的长度。要求梯形

7、的上底边长、下底边长、高从键盘输入。 c 图2-1 鱼塘示意图2 算法设计结束3. 程序编写 !班级：工力112 !姓名： !时间：201304 !计算鱼塘产量和产值 PROGRAM main INTEGER ：p=4,q=10 REAL ： a,b,c,h REAL ： area,perimeter，yield，sales PRINT*，请输入梯形鱼塘的上底、下底和高值： READ*，a,b,h area=(a+b)*h/2 yield=area*p sales=yield*q c=SQRT(b-a)/2)*(b-a)/2)+h*h) perimeter=a+b+2*c PRINT*， 鱼塘

8、产量：，yield，公斤。 PRINT*， 鱼塘产值：，sales，元。 PRINT*， 围墙长度：，perimeter，米。 END 4. 运行结果请输入梯形鱼塘的上底、下底和高值：150,210,40 鱼塘产量： 28800.00 公斤。 鱼塘产值： 288000.0 元。 围墙长度： 460.0000 米。Press any key to continue三. 课外练习1 问题描述：有一长5米、宽3米的长方形铁皮，在四个角对称挖了四个半径为0.2米的圆洞， 根据需要对这块铁皮镀一层铜，镀铜的价格是每平方米50元。计算镀铜费用。 2 算法设计开始输入矩形边长及圆洞半径：length、wid

9、th、radius计算铁皮的面积：area计算镀铜费用：rate输出3. 程序编写程序3-1: !班级：工力112 !姓名： !时间：201304 !计算镀铜费用 PROGRAM main INTEGER ： p=50 REAL ： length,width,radius REAL ： area,rate PRINT*,请输入铁皮长宽和圆洞半径： READ*，length,width,radius area=length*width-4*3.1415*radius*2 rate=p*area PRINT*, 镀铜费用：，rate END4. 运行结果请输入铁皮长宽和圆洞半径：5,3,0.2 镀

10、铜费用: 724.8680Press any key to continue实验二（1） 内部数据类型与简单顺序程序设计1 实验目的1 了解内部数据类型及其KIND值参数概念。2 掌握常量、变量、表达式、赋值语句的语法规则和基本功能。3 掌握常用标准函数的基本功能和使用要求。4 掌握表控输入输出的基本功能和使用要求。5 掌握顺序结构概念和顺序结构程序设计方法。进一步熟悉软件开发环境和掌握程序的输入、编辑、编译、构建和运行方法。二 课堂实验内容1 问题描述已知两复数：3.5+4.8i、7.5-3.2i。计算其和、差、积、商，并输出计算结果。编写程序实现之。要求两复数从键盘输入。2 算法设计开始输

11、入两个复数：c1,c2计算两数之和并赋值于a计算两复数之差并赋值于b计算复数之积赋值给c计算复数之商赋值给d输出a,b,c,d结束3. 程序编写根据算法和程序流程图，设计和编写程序如下。 !班级：工力112 !姓名： !时间：201304 !计算两复数的和、差、积、商 PROGRAM main COMPLEX ： c1,c2 complex： a,b,c,d PRINT*，请输入两个复数： READ*，c1,c2 a=c1+c2；b=c1-c2 c=c1*c2；d=c1/c2 PRINT*， 复数之和：，a PRINT*， 复数之差：，b PRINT*， 复数之积：，c PRINT*， 复数之

12、商：，d END !复数变量占有8字节4. 运行结果 请输入两个复数：3.5004.8007.500-3.2 复数之和： (11.00000,1.600000) 复数之差： (-4.000000,8.000000) 复数之积： (41.61000,24.80000) 复数之商： (0.1637840,0.7098812)Press any key to continue三. 课外练习1 问题描述某炮兵阵地，某次使用大炮向敌人阵地猛烈发射炮弹。炮弹的运行轨迹满足方程：f(x)=-0.0000125x2+0.07x-50,计算炮弹能打击的最大高度和最远距离。编写程序实现之。问题分析：（提示）x1设

13、：y=f(x)=ax2+bx+c，轨迹线为一抛物线，如图3-5所示。ymax图3-5 一元二次方程轨迹曲线一元二次方程求根公式为：假定方程只有两个不同实数根。2 算法设计开始请输入a,b,c计算两根x1,x2计算射程length和最大高度输出射程和最大高度：length、h结束!班级:工程力学112!姓名:!时间:20130426晚上!计算计算大炮的最大打击高度和最远距离，并输出PROGRAM mainREAL:a,b,c ,d,x1,x2,xREAL:length ,hPRINT*,请输入方程参数a,b,cREAD*,a,b,cd=b*b-4*a*cx1=(-b+SQRT(d)/(2*a)x

14、2=(-b-sqrt(d)/(2*a)x=0.5*(x1+x2)h=a*x*x+b*x+clength=x2-x1PRINT*,最大高度是:,h=,hPRINT*,最远射程是:,length=,lengthEND3. 程序编写4. 运行结果 请输入方程参数a,b,c-0.0000127,0.07,-50 最大高度是:h= 46.45669 最远射程是:length= 3825.183Press any key to continue实验二（2） 数据有格式输入输出一 实验目的1 了解数据有格式输入输出的主要作用、基本方法和使用规则。2 掌握数据有格式输入输出语句(READ、PRINT、WRIT

15、E)的基本功能和语法规则。3 掌握格式说明语句(FORMAT)的基本功能和语法规则。4 掌握I、F、E、A、L、X、/、等若干常用格式编辑符的格式要求和基本功能。5 掌握数据有格式输入输出的数据组织和格式说明方法。6 进一步掌握顺序结构程序设计方法。7 进一步熟悉软件开发环境。二 课堂实验内容1 问题描述通常需要计算(a+b)n的展开式系数，展开式系数满足杨辉三角关系。计算并输出n=1，2，3，4，5时的展开式系数及系数之和，即下面的杨辉三角。编写程序实现之。字母“a”、字母“b”、二项式阶数n从键盘输入。 1 1n=1 1 1 2n=2 1 2 1 4n=3 1 3 3 1 8n=4 1 4

16、 6 4 1 16n=5 1 5 10 10 5 1 32要求输入输出数据以规范、整洁、直观形式给出。程序中应采用有格式输入输出。二项式阶数n为2位整数，项符号为长度等于1的字符，系数为3位整数，系数和4位整数。2 算法设计开始根据提示输入二项式第一、二项字母和阶数：x、y、n计算n=1时二项式系数c11,c12及系数和s1计算n=2时二项式系数c21,c22,c23及系数和s2计算n=3时二项式系数c31,c32,c33,c34及系数和s3计算n=4时二项式系数c41,c42,c43,c44,c45及系数和s4。按输出格式要求输出数据。计算n=5时二项式系数c51,c52,c53,c54,c

17、55,c56及系数和s5。结束3. 程序编写!班级：工力112!姓名：!时间：20105!计算并输出二项式系数（杨辉三角）PROGRAM yanghui_triangleINTEGER(1) : nINTEGER: x,yINTEGER : c11=1,c12=1,s1INTEGER : c21=1,c22,c23=1,s2INTEGER : c31=1,c32,c33,c34=1,s3INTEGER: c41=1,c42,c43,c44,c45=1,s4INTEGER : c51=1,c52,c53,c54,c55,c56=1,s5WRITE(*,(1X,请输入二项式第一项字母(1个字母):

18、,)READ(*,(A) xWRITE(*,(1X,请输入二项式第二项字母(2个字母):,)READ(*,(A) yWRITE(*,(1X,请输入二项式阶数(2位整数):,)READ(*,(I2) ns1=c11+c12c22=c11+c12; s2=c21+c22+c23c32=c21+c22; c33=c22+c23; s3=c31+c32+c33+c34c42=c31+c32;c43=c32+c33;c44=c33+c34; s4=c41+c42+c43+c44+c45c52=c41+c42; c53=c42+c43; c54=c43+c44; c55=c44+c45; s5=c51+c

19、52+c53+c54+c55+c56WRITE(*,(1X,下面给出(,A1,+,A1,)1到,I2, 次方的展开式系数及系数之和。) x,y,nWRITE(*,(1X,阶数,9X,杨辉三角,7X,系数之和) WRITE(*,(1X,15X,I3) 1WRITE(*,(1X,n=,I2,8X,2(1X,I3),9X,I4) 1,c11,c12,s1WRITE(*,(1X,n=,I2,6X,3(1X,I3),7X,I4) 2,c21,c22,c23,s2WRITE(*,(1X,n=,I2,4X,4(1X,I3),5X,I4) 3,c31,c32,c33,c34,s3WRITE(*,(1X,n=,

20、I2,2X,5(1X,I3),3X,I4) 4,c41,c42,c43,c44,c45,s4WRITE(*,(1X,n=,I2, 6(1X,I3),1X,I4) 5,c51,c52,c53,c54,c55,c56,s5END 请输入二项式第一项字母(1个字母):1 请输入二项式第二项字母(2个字母):1,1 请输入二项式阶数(2位整数):5 下面给出(1+1)1到 5次方的展开式系数及系数之和。 阶数 杨辉三角 系数之和 1 n= 1 1 1 2 n= 2 1 2 1 4 n= 3 1 3 3 1 8 n= 4 1 4 6 4 1 16 n= 5 1 5 10 10 5 1 32Press a

21、ny key to continue4. 运行结果三. 课外练习1 问题描述使用数据有格式输入输出功能，输出由4个字符abcd组成的图形。编写程序实现之。 a b b c c ddddddddddddd c c b b a问题分析：（提示）a、b、c、d四个字符可以是任意可显示字符，这四个字符从键盘输入得到。若输入的四个字符为*、#、&、$，则输出结果为： * # # & & $ & & # # *由相同字符组成的字符串,可用标准函数REPEAT(str,I)得到。如REPEAT($,3)= $。2 算法设计开始结束3. 程序编写4. 运行结果实验三（1） 选择结构序程序设计一 实验目的1 熟

22、悉关系表达式和逻辑表达式。2 掌握逻辑IF、块IF、多支块IF、SELECT语句的语法规则和使用要求。3 掌握嵌套选择结构和嵌套选择语句的使用要求和实现方法。4 掌握选择结构概念和选择结构程序设计方法。二 课堂实验内容1 问题描述计算下面分段函数的值，编写程序实现之。x02 算法设计输入x的值5 程序编写 !班级：工力112 !姓名： !时间：201305 !计算分段函数的值 PROGRAM example REAL(4) ： x,y PRINT*,请输入x的值： READ*，x IF (x1000假真x500y=0.8*xx200真y=0.97*x假x0真y=xy=0.95*x输出y!班级:

23、工程力学112!姓名：!时间：201305!程序用途：计算超市购物优惠金额PROGRAM mainreal(4): x,yPRINT*,请输入购货金额read*,xIF (x=1000) THEN y=0.8*xELSE IF(X=500) THEN y=0.9*x ELSE IF(X=200) then y=0.95*x else IF(X=100) THEN y=0.97*x else y=xEND IFPRINT*,优惠金额为：,yEND结束3. 程序编写请输入购货金额1200 优惠金额为： 960.0000Press any key to continue 请输入购货金额600 优惠金

24、额为： 540.0000Press any key to continue 请输入购货金额250 优惠金额为： 237.5000Press any key to continue 请输入购货金额150 优惠金额为： 145.5000Press any key to continue 请输入购货金额50 优惠金额为： 50.00000Press any key to continue4. 运行结果请输入购货金额1200 优惠金额为： 960.0000Press any key to continue 请输入购货金额600 优惠金额为： 540.0000Press any key to conti

25、nue 请输入购货金额250 优惠金额为： 237.5000Press any key to continue 请输入购货金额150 优惠金额为： 145.5000Press any key to continue 请输入购货金额50 优惠金额为： 50.00000Press any key to continue实验三(2) 循环结构序程序设计一 实验目的1 掌握关系表达式和逻辑表达式。2 掌握循环控制条件的描述。3 掌握三种循环结构的特征及相互之间的等价转换。4 掌握DO、DO WHILE语句的语法规则和使用要求。5 掌握循环结构概念和循环结构程序设计方法。二 课堂实验内容1 问题描述验证哥德巴赫猜想。哥德巴赫提出，一个不小于6的偶数必定能表示为两个素数之和。如：8=3+5。求6至20之间的所有偶数的两个素数之和。编写程序实现之。输出格式为：6=3+38=3+510=3+710=5+5 开始2 算法设计请输入偶数n从3开始验