工程分析程序设计 上机作业(三)

工程分析程序设计上机作业（三）

模块化编程（2）

上机目的：进一步掌握内部例程、外部例程、接口块、模块等功能的使用方法。练习例程重载、例程递归的使用方法。

1、利用例程重载编写一个子程序求余数ACR(A,B)。要求，对两个整型数和两个实型数都有

效（实型数相除的余数：两个实型数相除后的商仍然取一个整数，但剩余值小于除数）[不能用Fortran的内部函数MOD(x,y)]。

MODULE MOD

IMPLICIT NONE

INTERFACE ACR

MODULE PROCEDURE ACRREALS,ACRINTEGERS

END INTERFACE

CONTAINS

FUNCTION ACRREALS(A,B)

REAL,INTENT(IN)::A,B

INTEGER C,ACRREALS

C=B/A

ACRREALS=B-A*C

END FUNCTION

FUNCTION ACRINTEGERS(A,B)

INTEGER,INTENT(IN)::A,B

INTEGER C,ACRINTEGERS

C=B/A

ACRINTEGERS=B-A*C

END FUNCTION

END MODULE

PROGRAM MAIN

USE MOD

IMPLICIT NONE

REAL::A=4.1,B=6.9

INTEGER::C=3,D=10

PRINT*,'ACR(A,B)=',ACR(A,B)

PRINT*,'ACR(C,D)=',ACR(C,D)

END PROGRAM

2、编写函数子程序GDC求两个数的最大公约数。求最大公约数的算法如下：把两个数中

大的那个数作为被除数，两数相除得到一个余数。把余数去除除数得到新一轮的余数。

不断重复这一过程直到余数为0，这时的除数就是两个数的最大公约数。

调用此函数，求1260，198，72三个数的最大公约数。

PROGRAM MAIN

IMPLICIT NONE

INTEGER::C=1260,D=198,E=72

INTEGER::F,G,H

F=GDC(C,D)

G=GDC(F,E)

PRINT*,G

CONTAINS

FUNCTION GDC(A,B)

INTEGER,INTENT(IN)::A,B

INTEGER GDC,Z,T,X,Y

X=A

Y=B

IF(X

Z=X

X=Y

Y=Z

ENDIF

DO WHILE(Y/=0)

T=Y

Y=MOD(X,Y)

X=T

END DO

GDC=T

END FUNCTION

END PROGRAM MAIN

3、 编写一个函数子程序求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()2

x x

e e x --=。 PROGRAM MAIN

IMPLICIT NONE

REAL F,I,X,G,N

PRINT*,'X='

READ*,X

I=1

F=X

G=0

DO WHILE(F>=1.0E-6)

I=I+1

CALL FACT(F,I,X)

END DO

DO N=1,I,1

CALL FACT(F,N,X)

G=G+F

END DO

PRINT*,'(2)≈',G

PRINT*,'(3)≈',(G-1/G)/2

CONTAINS

RECURSIVE SUBROUTINE FACT(F,N,X)

REAL F,N,X

IF(N==1)THEN

F=X

ELSE

CALL FACT(F,N-1,X)

F=F*X/N

END IF

END SUBROUTINE

END PROGRAM

4、 编写一程序用Euler 法求解微分方程22dy y x dx

=-，当x=0时，y=1.0。试求出x=0.1,0.2,0.3,0.4….,1.0时的y 值。算法如下：

Euler 法求解))(,('x y x f y =，定解条件：00,

y y x x ==。取向前差分， 令11'()n n y y y h

+=-， n n x x h -=+1 得)(),(*21h O y x f h y y n n n n ++=+

n y y y x f y y x f y y x x ??????????==2111000

0),(,),( PROGRAM MAIN

IMPLICIT NONE

REAL X,Y

X=0.0

Y=1.0

DO WHILE(X<1.0)

Y=Y+0.1*(Y**2-X**2)

X=X+0.1

PRINT*,'X=',X,',Y=',Y

END DO

END PROGRAM

软件工程案例分析

一、 阅读下列系统需求陈述，回答问题1、问题2、问题3和问题4。 某银行准备开发一个网上信用卡管理系统CCMS，该系统的基本功能为： （1）信用卡申请。非信用卡客户填写信用卡申请表，说明所要申请的信用卡类型及申请者的基本信息，提交CCMS登录。如果信用卡申请被银行接受，客户会收到银行的确认函，并告知用户信用卡的有效期及信贷限额；否则银行会发送一封拒绝函给该客户。客户收到确认函后，需再次登录CCMS ，用信用卡号和密码激活该信用卡。激活操作结束后，CCMS将激活通知发送给客户，告知客户其信用卡是否被成功地激活。 （2）月报表生成。在每个月第一天的零点，CCMS为每个信用卡客户创建一份月报表，对该客户上月的信用卡交易情况及交易额进行统计。信用卡客户可以登录CCMS查看月报表，也可以要求CCMS提供打印出的月报表。 （3）信用卡客户信息管理。信用卡客户的个人信息可以在 CCMS中进行在线的管理。每个信用卡客户可以在线查询其个人信息。 （4）信用卡交易记录。信用卡客户使用信息卡进行的每一笔交易都会记录在CCMS中。 （5）交易信息查询。信用卡客户可以登录CCMS查询并核实其信用卡交易记录及交易额。在系统的需求分析阶段，使用用例对系统需求建模。表1—1和表1—2给出了其中两个用例的概要描述。 ［问题1］） 将表1—1和表1—2中的（1）~（10）填充完整。 ［问题2］ 除了表1—1和表1—2给出的用例外，从上述系统陈述中还可以获取哪些由信用卡客户发起的用例？（给出用例名称即可）

［问题3］ 用400字以内文字，简要说明用例获取的基本步骤。 ［问题4］ 用例除了使用表1—1和表1—2所示的形式描述外，还可以使用UML的用例图来表示。分别用50字以内文字，解释UML用例图中扩展用例和抽象用例的内涵。 二、 阅读以下关于工作流系统性能分析的叙述，回答问题1、问题2和问题3。 某企业正在创建一个工作流管理系统，目前正处于过程定义阶段，即创建工作流模型阶段。对于这些工作流模型，除了要考虑工作流的正确性外，工作流的性能也是十分重要的。工作流性能主要反映工作流定量方面的特性，例如，任务的完成时间、单位时间内处理的任务数量、资源的利用率以及在预定的标准时间内完成任务的百分比等等。 图2—1所示的是一个简单的工作流模型（其中单位时间为1小时），它表示这样一个执行过程：每小时将会有20个任务达到c1，这20个任务首先经过处理taskl，再经过处理task2，最终将结果传递到c3。处理taskl和处理task2相互独立。 图2-1 假设性能评价模型符合M/M/1排队模型，在计算性能指标的过程中可以使用下列公式进行计 算：,其中ρ表示资源利用率，表示单位时间内到达的任务数，表示该资源单位时间内能够完成的任务数。 ［问题1］ 计算图2—1所示的工作流模型的下列性能指标： （1）每个资源的利用率； （2）每个处理中的平均任务数L； （3）平均系统时间S； （4）每个处理的平均等待时间W。 ［问题2］

(完整版)哈工大工程热力学习题答案——杨玉顺版

第二章 热力学第一定律 思 考 题 1. 热量和热力学能有什么区别？有什么联系？ 答：热量和热力学能是有明显区别的两个概念：热量指的是热力系通过界面与外界进行的热能交换量，是与热力过程有关的过程量。热力系经历不同的过程与外界交换的热量是不同的；而热力学能指的是热力系内部大量微观粒子本身所具有的能量的总合，是与热力过程无关而与热力系所处的热力状态有关的状态量。简言之，热量是热能的传输量，热力学能是能量？的储存量。二者的联系可由热力学第一定律表达式 d d q u p v δ=+ 看出；热量的传输除了可能引起做功或者消耗功外还会引起热力学能的变化。 2. 如果将能量方程写为 d d q u p v δ=+ 或 d d q h v p δ=- 那么它们的适用范围如何？ 答：二式均适用于任意工质组成的闭口系所进行的无摩擦的内部平衡过程。因为 u h pv =-，()du d h pv dh pdv vdp =-=-- 对闭口系将 du 代入第一式得 q dh pdv vdp pdv δ=--+ 即 q dh vdp δ=-。 3. 能量方程 δq u p v =+d d （变大） 与焓的微分式 ()d d d h u pv =+（变大） 很相像，为什么热量 q 不是状态参数，而焓 h 是状态参数？ 答：尽管能量方程 q du pdv δ=+ 与焓的微分式 ()d d d h u pv =+（变大）似乎相象，但两者 的数学本质不同，前者不是全微分的形式，而后者是全微分的形式。是否状态参数的数学检验就是，看该参数的循环积分是否为零。对焓的微分式来说，其循环积分：()dh du d pv =+???蜒? 因为 0du =??，()0d pv =?? 所以 0dh =??， 因此焓是状态参数。 而对于能量方程来说，其循环积分： q du pdv δ=+???蜒?

数值分析上机作业

数值分析上机实验报告 选题：曲线拟合的最小二乘法 指导老师： 专业： 学号： 姓名：

课题八曲线拟合的最小二乘法 一、问题提出 从随机的数据中找出其规律性，给出其近似表达式的问题，在生产实践和科学实验中大量存在，通常利用数据的最小二乘法求得拟合曲线。 在某冶炼过程中，根据统计数据的含碳量与时间关系，试求含碳量y 与时间t 的拟合曲线。 二、要求 1、用最小二乘法进行曲线拟合； 2、近似解析表达式为()33221t a t a t a t ++=?； 3、打印出拟合函数()t ?，并打印出()j t ?与()j t y 的误差，12,,2,1 =j ； 4、另外选取一个近似表达式，尝试拟合效果的比较； 5、*绘制出曲线拟合图*。 三、目的和意义 1、掌握曲线拟合的最小二乘法； 2、最小二乘法亦可用于解超定线代数方程组； 3、探索拟合函数的选择与拟合精度间的关系。 四、计算公式 对于给定的测量数据(x i ,f i )(i=1,2,…，n )，设函数分布为 ∑==m j j j x a x y 0)()(? 特别的，取)(x j ?为多项式 j j x x =)(? (j=0, 1,…，m )

则根据最小二乘法原理，可以构造泛函 ∑∑==-=n i m j i j j i m x a f a a a H 1 10))((),,,(? 令 0=??k a H (k=0, 1,…，m ) 则可以得到法方程 ???? ??????? ?=????????????????????????),(),(),(),(),(),(),(),(),(),(),(),(1010101111000100m m m m m m m m f f f a a a ????????????????????? 求该解方程组，则可以得到解m a a a ,,,10 ，因此可得到数据的最小二乘解 ∑=≈m j j j x a x f 0)()(? 曲线拟合：实际工作中，变量间未必都有线性关系，如服药后血药浓度与时间的关系；疾病疗效与疗程长短的关系；毒物剂量与致死率的关系等常呈曲线关系。曲线拟合是指选择适当的曲线类型来拟合观测数据，并用拟合的曲线方程分析两变量间的关系。 五、结构程序设计 在程序结构方面主要是按照顺序结构进行设计，在进行曲线的拟合时，为了进行比较，在程序设计中，直接调用了最小二乘法的拟合函数polyfit ，并且依次调用了plot 、figure 、hold on 函数进行图象的绘制，最后调用了一个绝对值函数abs 用于计算拟合函数与原有数据的误差，进行拟合效果的比较。

系统工程的案例分析

金 马 河 2 1 3 离堆

鱼嘴分水工程：把岷江水一分为二 “鱼嘴”是都江堰分水工程，因其形如鱼嘴而得名，它昂头于岷江江心，把岷江分成内外二江。西边叫外江，俗称“金马河”，是岷江正流，主要用于排洪；东边沿山脚的叫内江，是人工引水渠道，主要用于灌溉； 鱼嘴的设置极为巧妙，它利用地形、地势，巧妙地完成分流引水的任务，而且在洪、枯水季节不同水位条件下，起着自动调节水量的作用。 鱼嘴所分的水量有一定的比例,春天，岷江水流量小；灌区正值春耕，需要灌溉,这时岷江主流直入内江,水量约占六成,外江约占四成，以保证灌溉用水；洪水季节，二者比例又自动颠倒过来，内江四成，外江六成，使灌区不受水潦灾害。 在二王庙壁上刻的治水《三字经》中说的“分四六，平潦l áo 旱（潦l áo 即“涝水”,在今陕西省户县境,北入渭河）”,就是指鱼嘴这一天然调节分流比例的功能。“八水绕长安”胜景：黑、涝、新、沣、浐、灞、临、曹运明渠。 我们的祖先十分聪明，在流量小、用水紧张时，为了不让外江40％的流量白白浪费，采用杩m à搓c h a （同杈，三脚木架。中国四川都江堰的活动拦水坝，就是用杩杈和满装卵石的竹笼做成的）截流的办法，把外江水截入内江，整就使内江灌区春耕用水更加可靠。1974 年，在鱼嘴西岸的外江河口

建成一座钢筋混凝土结构的电动制闸，代替过去临时杩搓工程，截流排洪，更加灵活可靠。 “飞沙堰”溢洪道：把多余的洪水和流沙排入外江 在鱼嘴以下的长堤，即分内、外二江的堤叫金刚堤。堤下段与内江左岸虎头岸相对的地方，有一低平的地段，这里春、秋、冬、三季是人们往返于离堆公园与索桥之间的行道的坦途，洪水季节这里浪花飞溅，是内江的泄洪道。 泄洪道，唐朝名“侍郎堰”、“金堤”，后又名“减水河”，它具有泄洪徘砂的显著功能，故又叫它“飞沙堰”。 飞沙堰是都江堰三大件之一，看上去十分平凡，但功能非常实用，可以说是确保成都平原不受水灾的关键要害。 飞沙堰的作用主要是当内江 的水量超过宝瓶口流量上限时，多 余的水便从飞沙堰自行溢出；如遇 特大洪水的非常情况，它还会自行 溃堤，让大量江水回归岷江正流。 另一作用是“飞沙”，岷江从万山 丛中急驰而来，挟着大量泥沙，石 块，如果让它们顺内江而下，就会 淤塞宝瓶口和灌区。飞沙堰真是善 解人意、排人所难，将上游带来的 泥沙和卵石，甚至重达千斤的巨石，从这里抛入外江（主要是巧妙地利用离心力作用），确保内江通畅，确有鬼斧神功之妙。 “深淘滩，低作堰”是都江堰的治水名言，淘滩是指飞沙堰一段、内江一段河道要深淘，深淘的标准是古人在河底深处预埋的“卧铁”。岁修淘滩要淘到卧铁为止，才算恰到好处，才能保证灌区用水。低作堰就是说飞沙堰有一定高度，高了进水多，低了进水少，都不合适。 古时飞沙堰，是用竹笼卵石堆砌的临时工程；如今已改用混凝土浇铸，以保一劳永逸的功效。 内江 外江 金 刚 堤

工程热力学作业.

1-1 一立方形刚性容器，每边长1m ，将其中气体的压力抽至1000Pa ，问其真空度为多少毫米汞柱？容器每面受力多少牛顿？已知大气压力为0.1MPa 。 解：p = 1 000 Pa = 0.001 MPa 真空度mmHg Pa MPa MPa MPa p p p b V 56.74299000099.0001.01.0===-=-= 容器每面受力F ＝p V A = 9 900 Pa×1m 2 =9.9×104 N 1-2 试确定表压力为0.01 MPa 时U 形管压力计中液柱的高度差。(1)U 形管中装水，其密度为1 000 kg/m 3；(2) U 形管中装酒精，其密度为789 kg/m 3。 解： 因为表压力可以表示为p g =ρgΔz ，所以有 g p z g ρ= ? 既有（1）mm m s m m kg Pa g p z g 72.101901972.1/80665.9/10001001.02 36==??=?＝水ρ （2） mm m s m m kg Pa g p z g 34.129729734.1/80665.9/7891001.02 36==??=?＝酒精 ρ 1-7 从工程单位制热力性质查得，水蒸气在500℃、100at 时的比体积和比焓分别为v =0.03347m 3/kg 、h =806.6kcal/kg 。在国际单位制中，这时水蒸气的压力和比热力学能各为多少？ 解： 水蒸气压力p =100at×9.80665×104Pa/at = 9.80665×106Pa=9.80665MPa 比热力学能u=h-pv=806.6kcal ×4.1868kJ/kcal)/kg-9806.65kPa ×0.03347m 3/kg = 3377.073kJ-328.228kJ =3048.845kJ 2-1 冬季，工厂某车间要使室内维持一适宜温度。在这一温度下，透过墙壁和玻璃等处，室内向室外每一小时传出0.7×106kcal 的热量。车间各工作机器消耗的动力为是500PS(认为机器工作时将全部动力转变为热能)。另外，室内经常点着50 盏100W 的电灯，要使该车间的温度保持不变，问每小时需供给多少kJ 的热量？ 解：要使车间保持温度不变，必须使车间内每小时产生的热量等散失的热量 Q = Q 机＋Q 灯＋Q 散＋Q 补 = 0 Q 机 = 500PSh = 500×2.647796×103 kJ = 1.32×106 kJ Q 灯 = 50×100W×3600s = 1.8×107J = 1.8×104 kJ Q 散 = -0.7×106kcal =- 0.7×106×4.1868kJ = -2.93×106 kJ Q 补 = -Q 机-Q 灯＋Q 散 = -1.32×106 kJ-1.8×104 kJ+2.93×106 kJ = 1.592×106 kJ

数值分析上机作业

昆明理工大学工科研究生《数值分析》上机实验 学院：材料科学与工程学院 专业：材料物理与化学 学号：2011230024 姓名: 郑录 任课教师：胡杰

P277-E1 1.已知矩阵A= 10787 7565 86109 75910 ?? ?? ?? ?? ?? ??,B= 23456 44567 03678 00289 00010 ?? ?? ?? ?? ?? ?? ?? ?? ,错误！未找到引用源。 = 11/21/31/41/51/6 1/21/31/41/51/61/7 1/31/41/51/61/71/8 1/41/51/61/71/81/9 1/51/61/71/81/91/10 1/61/71/81/91/101/11?????????????????? （1）用MA TLAB函数“eig”求矩阵全部特征值。 （2）用基本QR算法求全部特征值（可用MA TLAB函数“qr”实现矩阵的QR分解）。解：MA TLAB程序如下： 求矩阵A的特征值： clear; A=[10 7 8 7;7 5 6 5;8 6 10 9;7 5 9 10]; E=eig(A) 输出结果： 求矩阵B的特征值： clear; B=[2 3 4 5 6;4 4 5 6 7;0 3 6 7 8;0 0 2 8 9;0 0 0 1 0]; E=eig(B) 输出结果：

求矩阵错误！未找到引用源。的特征值： clear; 错误！未找到引用源。=[1 1/2 1/3 1/4 1/5 1/6; 1/2 1/3 1/4 1/5 1/6 1/7; 1/3 1/4 1/5 1/6 1/7 1/8; 1/4 1/5 1/6 1/7 1/8 1/9;1/5 1/6 1/7 1/8 1/9 1/10; 1/6 1/7 1/8 1/9 1/10 1/11]; E=eig(错误！未找到引用源。) 输出结果： （2）A= 10 7877565861097 5 9 10 第一步：A0=hess(A);[Q0,R0]=qr(A0);A1=R0*Q0 返回得到： 第二部：[Q1,R1]=qr(A1);A2=R1*Q1

数值分析上机题课后作业全部-东南大学

2015.1.9 上机作业题报告 USER

1.Chapter1 1.1题目 设S N = 1 j 2?1 N j =2 ，其精确值为 )1 1123(21+--N N 。 （1）编制按从大到小的顺序1 1 1311212 22-+??+-+-=N S N ，计算S N 的通用程序。 （2）编制按从小到大的顺序1 21 1)1(111222-+ ??+--+-= N N S N ，计算S N 的通用程序。 （3）按两种顺序分别计算64210,10,10S S S ,并指出有效位数。（编制程序时用单精度） （4）通过本次上机题，你明白了什么？ 1.2程序 1.3运行结果

1.4结果分析 按从大到小的顺序，有效位数分别为：6，4，3。 按从小到大的顺序，有效位数分别为：5，6，6。 可以看出，不同的算法造成的误差限是不同的，好的算法可以让结果更加精确。当采用从大到小的顺序累加的算法时，误差限随着N 的增大而增大，可见在累加的过程中，误差在放大，造成结果的误差较大。因此，采取从小到大的顺序累加得到的结果更加精确。 2.Chapter2 2.1题目 （1）给定初值0x 及容许误差ε，编制牛顿法解方程f(x)=0的通用程序。 （2）给定方程03 )(3 =-=x x x f ,易知其有三个根3,0,3321= *=*-=*x x x ○1由牛顿方法的局部收敛性可知存在,0>δ当),(0δδ+-∈x 时，Newton 迭代序列收敛于根x2*。试确定尽可能大的δ。 ○2试取若干初始值，观察当),1(),1,(),,(),,1(),1,(0+∞+-----∞∈δδδδx 时Newton 序列的收敛性以及收敛于哪一个根。 （3）通过本上机题，你明白了什么？ 2.2程序

《系统工程案例分析》课程实践作业题

《系统工程案例分析》课程实践作业题 1、大学图书馆使用效率调查研究 图书馆是高校的信息中心，为读者用户提供设施资源、文献资源及信息检索资源等服务，其使用效率的高低直接反映了学员的自我学习能力。为分析图书馆的使用效率，需对图书馆设施、文献及信息检索资源现状等进行调查。 要求利用系统工程的相关理论和方法，设计一种方便可行且准确度较高的方法，对本校图书馆使用效率进行调查研究，并调查结果进行分析。

2、医院眼科病床的合理安排的问题 医院就医排队是大家都非常熟悉的现象，例如，患者到门诊就诊、到收费处划价、到药房取药、到注射室打针、等待住院等，往往都需要排队等待接受某种服务。 考虑某医院眼科病床的合理安排的问题。 该医院眼科门诊每天开放，住院部共有病床79张。该医院眼科手术主要分四大类：白内障、视网膜疾病、青光眼和外伤。白内障手术较简单，而且没有急症。目前该院是每周一、三做白内障手术，此类病人的术前准备时间只需1、2天。做两只眼的病人比做一只眼的要多一些，大约占到60%。如果要做双眼是周一先做一只，周三再做另一只。 外伤疾病通常属于急症，病床有空时立即安排住院，住院后第二天便会安排手术。其他眼科疾病比较复杂，有各种不同情况，但大致住院以后2-3天内就可以接受手术，主要是术后的观察时间较长。这类疾病手术时间可根据需要安排，一般不安排在周一、周三，并不考虑在急症范围内。 该医院考虑病床安排时可不考虑手术条件的限制，但在通常情况下白内障手术与其他眼科手术（急症除外）不安排在同一天做。当前该住院部对全体非急症病人是按照FCFS（First come, First serve）规则安排住院，但等待住院病人队列却越来越长，通过数学建模的方法来合理安排住院部的病床，完成下列问题，以提高对医院资源的有效利用。 问题一：试分析确定合理的评价指标体系，用以评价该问题的病床安排模型的优劣。 问题二：试就该住院部当前的情况，建立合理的病床安排模型，以根据已知的第二天拟出院病人数来确定第二天应该安排哪些病人住院。并对你们的模型利用问题一中的指标体系作出评价。 问题三：作为病人，自然希望尽早知道自己大约何时能住院。能否根据当时住院病人及等待住院病人的统计情况，在病人门诊时即告知其大致入住时间区间。 问题四：若该住院部周六、周日不安排手术，请你们重新回答问题二，医院的手术时间安排是否应作出相应调整? 问题五：有人从便于管理的角度提出建议，在一般情形下，医院病床安排可采取使各类病人占用病床的比例大致固定的方案，试就此方案，建立使得所有病人在系统内的平均逗留时间（含等待入院及住院时间）最短的病床比例分配模型。

工程热力学习题解答

1. 热量和热力学能有什么区别？有什么联系？ 答：热量和热力学能是有明显区别的两个概念：热量指的是热力系通过界面与外界进行的热能交换量，是与热力过程有关的过程量。热力系经历不同的过程与外界交换的热量是不同的；而热力学能指的是热力系内部大量微观粒子本身所具有的能量的总合，是与热力过程无关而与热力系所处的热力状态有关的状态量。简言之，热量是热能的传输量，热力学能是能量？的储存量。二者的联系可由热力学第一定律表达式 d d q u p v δ=+ 看出；热量的传输除了可能引起做功或者消耗功外还会引起热力学能的变化。 2. 如果将能量方程写为 d d q u p v δ=+ 或 d d q h v p δ=- 那么它们的适用范围如何？ 答：二式均适用于任意工质组成的闭口系所进行的无摩擦的内部平衡过程。因为 u h p v =-，()du d h pv dh pdv vdp =-=-- 对闭口系将 du 代入第一式得 q dh pdv vdp pdv δ=--+ 即 q dh vdp δ=-。 3. 能量方程 δq u p v =+d d （变大） 与焓的微分式 ()d d d h u pv =+（变大） 很相像，为什么热量 q 不是状态参数，而焓 h 是状态参数？ 答：尽管能量方程 q du pdv δ=+ 与焓的微分式 ()d d d h u pv =+（变大）似乎相象，但两者的数学本 质不同，前者不是全微分的形式，而后者是全微分的形式。是否状态参数的数学检验就是，看该参数的循环积分是否为零。对焓的微分式来说，其循环积分：()dh du d pv =+??? 因为 0du =?，()0d pv =? 所以 0dh =?， 因此焓是状态参数。 而对于能量方程来说，其循环积分： q du pdv δ=+??? 虽然： 0du =? 但是： 0pdv ≠? 所以： 0q δ≠? 因此热量q 不是状态参数。 4. 用隔板将绝热刚性容器分成A 、B 两部分（图2-13），A 部分装有1 kg 气体，B 部分为高度真空。将隔板抽去后，气体热力学能是否会发生变化？能不能用 d d q u p v δ=+ 来分析这一过程？

东南大学数值分析上机作业汇总

东南大学数值分析上机作业 汇总 -标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

数值分析上机报告 院系： 学号： 姓名：

目录 作业1、舍入误差与有效数 (1) 1、函数文件cxdd.m (1) 2、函数文件cddx.m (1) 3、两种方法有效位数对比 (1) 4、心得 (2) 作业2、Newton迭代法 (2) 1、通用程序函数文件 (3) 2、局部收敛性 (4) （1）最大δ值文件 (4) （2）验证局部收敛性 (4) 3、心得 (6) 作业3、列主元素Gauss消去法 (7) 1、列主元Gauss消去法的通用程序 (7) 2、解题中线性方程组 (7) 3、心得 (9) 作业4、三次样条插值函数 (10) 1、第一型三次样条插值函数通用程序： (10) 2、数据输入及计算结果 (12)

作业1、舍入误差与有效数 设∑ =-=N j N j S 2 2 11 ，其精确值为?? ? ??---1112321N N . （1）编制按从小到大的顺序1 1 131121222-? ??+-+-=N S N ，计算N S 的通用程序； （2）编制按从大到小的顺序()1 21 11111222-???+--+-=N N S N ，计算N S 的通用程序； （3）按两种顺序分别计算642101010,,S S S ,并指出有效位数； （4）通过本上机你明白了什么？ 程序： 1、函数文件cxdd.m function S=cxdd(N) S=0; i=2.0; while (i<=N) S=S+1.0/(i*i-1); i=i+1; end script 运行结果（省略>>）： S=cxdd(80) S= 0.737577 2、函数文件cddx.m function S=cddx (N) S=0; for i=N:-1:2 S=S+1/(i*i-1); end script 运行结果（省略>>）： S=cddx(80) S= 0.737577 3、两种方法有效位数对比

数值分析作业思考题汇总

￥ 数值分析思考题1 1、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。 2、相对误差在什么情况下可以用下式代替 3、查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。 4、取 ，计算 ，下列方法中哪种最好为什么（1）(3 3-，（2）(2 7-，（3） ()3 1 3+ ，（4） ()6 1 1 ，（5）99- , 数值实验 数值实验综述：线性代数方程组的解法是一切科学计算的基础与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，因此也称为精确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时，Gauss消去法是目前最基本和常用的方法。如若系数矩阵具有某种特殊形式，则为了尽可能地减少计算量与存储量，需采用其他专门的方法来求解。 Gauss消去等同于矩阵的三角分解，但它存在潜在的不稳定性，故需要选主元素。对正定对称矩阵，采用平方根方法无需选主元。方程组的性态与方程组的条件数有关，对于病态的方程组必须采用特殊的方法进行求解。 数值计算方法上机题目1 1、实验1. 病态问题 实验目的： 算法有“优”与“劣”之分，问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感，反之属于好问题。希望读者通过本实验对此有一个初步的体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 $ r e x x e x x ** * ** - == 141 . ≈)61

Matlab作业3(数值分析)答案

Matlab作业3（数值分析） 机电工程学院（院、系）专业班组 学号姓名实验日期教师评定 1.计算多项式乘法(x2+2x+2)(x2+5x+4)。 答： 2. （1）将(x-6)(x-3)(x-8)展开为系数多项式的形式。（2）求解在x=8时多项 式(x-1)(x-2) (x-3)(x-4)的值。 答：（1） （2）

3. y=sin(x)，x从0到2π，?x=0.02π，求y的最大值、最小值、均值和标准差。 4.设ｘ＝[0.00.30.8 1.1 1.6 2.3]'，y=[0.500.82 1.14 1.25 1.35 1.40]'，试求二次多项式拟合系数，并据此计算x1=[0.9 1.2]时对应的y1。解：x=[0.0 0.3 0.8 1.1 1.6 2.3]'; %输入变量数据x y=[0.50 0.82 1.14 1.25 1.35 1.40]'; %输入变量数据y p=polyfit(x,y,2) %对x,y用二次多项式拟合,得到系数p x1=[0.9 1.2]; %输入点x1 y1=polyval(p,x1) %估计x1处对应的y1 p = -0.2387 0.9191 0.5318 y1 = a) 1.2909

5.实验数据处理：已知某压力传感器的测试数据如下表 p为压力值，u为电压值，试用多项式 d cp bp ap p u+ + + =2 3 ) ( 来拟 合其特性函数，求出a,b,c,d，并把拟合曲线和各个测试数据点画在同一幅图上。解： >> p=[0.0,1.1,2.1,2.8,4.2,5.0,6.1,6.9,8.1,9.0,9.9]; u=[10,11,13,14,17,18,22,24,29,34,39]; x=polyfit(p,u,3) %得多项式系数 t=linspace(0,10,100); y=polyval(x,t); %求多项式得值 plot(p,u,'*',t,y,'r') %画拟和曲线 x = 0.0195 -0.0412 1.4469 9.8267

(完整版)数值计算方法上机实习题答案

1． 设?+=1 05dx x x I n n ， （1） 由递推公式n I I n n 1 51+-=-，从0I 的几个近似值出发，计算20I ； 解：易得：0I =ln6-ln5=0.1823, 程序为： I=0.182; for n=1:20 I=(-5)*I+1/n; end I 输出结果为：20I = -3.0666e+010 （2） 粗糙估计20I ，用n I I n n 51 5111+- =--，计算0I ； 因为 0095.05 6 0079.01020 201 020 ≈<<≈??dx x I dx x 所以取0087.0)0095.00079.0(2 1 20=+= I 程序为：I=0.0087; for n=1:20 I=(-1/5)*I+1/(5*n); end I 0I = 0.0083 （3） 分析结果的可靠性及产生此现象的原因(重点分析原因)。 首先分析两种递推式的误差；设第一递推式中开始时的误差为000I I E '-=，递推过程的舍入误差不计。并记n n n I I E '-=，则有01)5(5E E E n n n -==-=-Λ。因为=20E 20020)5(I E >>-，所此递推式不可靠。而在第二种递推式中n n E E E )5 1(5110-==-=Λ，误差在缩小， 所以此递推式是可靠的。出现以上运行结果的主要原因是在构造递推式过程中，考虑误差是否得到控制， 即算法是否数值稳定。 2． 求方程0210=-+x e x 的近似根，要求4 1105-+?<-k k x x ，并比较计算量。 （1） 在[0，1]上用二分法； 程序：a=0;b=1.0; while abs(b-a)>5*1e-4 c=(b+a)/2;

工程热力学课后作业答案(第七章)第五版

7-1当水的温度t=80℃，压力分别为、、、及1MPa时，各处于什么状态并求出该状态下的焓值。 解：查表知道t=80℃时饱和压力为。 因此在、、、及1MPa时状态分别为过热、未饱和、未饱和，未饱和、未饱和。焓值分别为kg，kJ/kg，335 kJ/kg，kJ/kg，kJ/kg。 7－2已知湿蒸汽的压力p=1MPa干度x=。试分别用水蒸气表和h-s图求出hx，vx，ux，sx。解：查表得：h``＝2777kJ/kg h`= kJ/kg v``＝kg v`＝m3/kg u``= h``－pv``= kJ/kg u`＝h`－pv`= kJ/kg s``= kJ/ s`＝kJ/ hx＝xh``+(1-x)h`= kJ/kg vx＝xv``+(1-x)v`= m3/kg ux＝xu``+(1-x)u`=2400 kJ/kg sx＝xs``+(1-x)s`= kJ/ 7-3在V＝60L的容器中装有湿饱和蒸汽，经测定其温度t＝210℃，干饱和蒸汽的含量mv=，试求此湿蒸汽的干度、比容及焓值。 解：t＝210℃的饱和汽和饱和水的比容分别为： v``＝kg v`＝m3/kg h``＝kg h`= kJ/kg 湿饱和蒸汽的质量： 解之得： x= 比容：vx＝xv``+(1-x)v`= m3/kg 焓：hx＝xh``+(1-x)h`=1904kJ/kg 7－4将2kg水盛于容积为的抽空了的密闭刚性容器中，然后加热至200℃试求容器中（1）压力；（2）焓；（3）蒸汽的质量和体积。 解：（1）查200℃的饱和参数 h``＝kg h`= kJ/kg v``＝kg v`＝kg 饱和压力。 刚性容器中水的比容： ＝m3/kg

安全系统工程在线作业答案

安全系统工程在线作业答案 -标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

第一阶段在线作业单项选择题 第1题安全是 A、没有危险的状态 B、没有事故的状态 C、达到可接受的伤亡和损失的状态 D、舒适的状态 第2题安全系统的认识论产生于（） A、20世纪初 B、20世纪50年代后 C、20世纪末 D、21世纪初 第3题生产事故的特性不包括那一项 A、因果性 B、随机性 C、潜伏性 D、可预防性 E、严重性 第4题海因里希的因果连锁论所建立的多米诺骨牌模型中最后一块骨牌所代表 的是 A、人的不安全行为 B、物的不安全状态 C、伤害 D、事故 第5题在系统安全分析中，FMECA指的是 A、事故树法 B、故障类型及影响分析法 C、原因—后果分析法 D、故障类型、影响及致命度分析法 第6题海因里希对5000多起伤害事故案例进行了详细调查研究后得出海因里希法则，事故后果为

严重伤害、轻微伤害和无伤害的事故件数之比为 A、1:29:300 B、1:10:300 C、1:10:100 D、1:100:500 第7题安全系统工程的研究对象 A、人子系统 B、机器子系统 C、环境子系统 D、人—机—环 第8题下列属于安全系统工程研究内容的是 （） A、地质灾害 B、社会治安 C、安全评价 D、疾病防治 第9题 ()是系统安全的主要观点 A、安全是绝对的 B、安全只是系统运行阶段的要考虑的工作 C、在系统的各个阶段都要进行危险源辨识、评价和控制 D、事故是系统的危险源 第10题事故和隐患是 A、完全相同的 B、后者是前者的可能性 C、后者是前者的必然条件 D、前者是后者的必然条件 第11题安全检查表方法不可以用于 A、方案设计

数值分析上机作业1-1

数值计算方法上机题目1 1、实验1. 病态问题 实验目的： 算法有“优”与“劣”之分，问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感，反之属于好问题。希望读者通过本实验对此有一个初步的体会。 数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 问题提出： 考虑一个高次的代数多项式 ∏=-= ---=20 1)()20)...(2)(1()(k k x x x x x p （E1-1） 显然该多项式的全部根为l ，2，…，20，共计20个，且每个根都是单重的（也称为简 单的）。现考虑该多项式方程的一个扰动 0)(19 =+x x p ε (E1-2) 其中ε是一个非常小的数。这相当于是对（E1-1）中19 x 的系数作一个小的扰动。我们希望比较（E1-1）和（E1-2）根的差别，从而分析方程（E1-1）的解对扰动的敏感性。 实验内容： 为了实现方便，我们先介绍两个 Matlab 函数：“roots ”和“poly ”，输入函数 u ＝roots （a ） 其中若变量a 存储1+n 维的向量，则该函数的输出u 为一个n 维的向量。设a 的元素依次为121,...,,+n a a a ，则输出u 的各分量是多项式方程 0...1121=++++-n n n n a x a x a x a 的全部根，而函数 b=poly(v) 的输出b 是一个n ＋1维变量，它是以n 维变量v 的各分量为根的多项式的系数。可见“roots ”和“Poly ”是两个互逆的运算函数. ve=zeros(1,21); ve(2)=ess; roots(poly(1:20))+ve) 上述简单的Matlab 程序便得到（E1-2）的全部根，程序中的“ess ”即是（E1-2）中的ε。 实验要求： （1）选择充分小的ess ，反复进行上述实验，记录结果的变化并分析它们。如果扰动项的系数ε很小，我们自然感觉（E1-1）和(E1-2)的解应当相差很小。计算中你有什么出乎意料的发现？表明有些解关于如此的扰动敏感性如何？ （2）将方程（E1-2）中的扰动项改成18 x ε或其他形式，实验中又有怎样的现象出现？

《系统工程案例分析》课程实践作业题

系统工程案例分析》课程实践作业题 1、大学图书馆使用效率调查研究 图书馆是高校的信息中心，为读者用户提供设施资源、文献资源及信息检索资源等服务，其使用效率的高低直接反映了学员的自我学习能力。为分析图书馆的使用效率，需对图书馆设施、文献及信息检索资源现状等进行调查。 要求利用系统工程的相关理论和方法，设计一种方便可行且准确度较高的方法，对本校图书馆使用效率进行调查研究，并调查结果进行分析。

2、医院眼科病床的合理安排的问题 医院就医排队是大家都非常熟悉的现象，例如，患者到门诊就诊、到收费处划价、到药房取药、到注射室打针、等待住院等，往往都需要排队等待接受某种服务。 考虑某医院眼科病床的合理安排的问题。 该医院眼科门诊每天开放，住院部共有病床79张。该医院眼科手术主要分四大类：白内障、视网膜疾病、青光眼和外伤。白内障手术较简单，而且没有急症。目前该院是每周一、三做白内障手术，此类病人的术前准备时间只需1、2天。做两只眼的病人比做一只眼的要多一些，大约占到60%。如果要做双眼是周一先做一只，周三再做另一只。 外伤疾病通常属于急症，病床有空时立即安排住院，住院后第二天便会安排手术。其他眼科疾病比较复杂，有各种不同情况，但大致住院以后2-3 天内就可以接受手术，主要是术后的观察时间较长。这类疾病手术时间可根据需要安排，一般不安排在周一、周三，并不考虑在急症范围内。 该医院考虑病床安排时可不考虑手术条件的限制，但在通常情况下白内障手术与其他眼科手术（急症除外）不安排在同一天做。当前该住院部对全体非急症病人是按照FCFS（First come, First serve ）规则安排住院，但等待住院病人队列却越来越长，通过数学建模的方法来合理安排住院部的病床，完成下列问题，以提高对医院资源的有效利用。 问题一：试分析确定合理的评价指标体系，用以评价该问题的病床安排模型的优劣。问题二：试就该住院部当前的情况，建立合理的病床安排模型，以根据已知的第二天拟出院病人数来确定第二天应该安排哪些病人住院。并对你们的模型利用问题一中的指标体系作出评价。 问题三：作为病人，自然希望尽早知道自己大约何时能住院。能否根据当时住院病人及等待住院病人的统计情况，在病人门诊时即告知其大致入住时间区间。 问题四：若该住院部周六、周日不安排手术，请你们重新回答问题二，医院的手术时间安排是否应作出相应调整? 问题五：有人从便于管理的角度提出建议，在一般情形下，医院病床安排可采取使各类病人占用病床的比例大致固定的方案，试就此方案，建立使得所有病人在系统内的平均逗留时间（含等待入院及住院时间）最短的病床比例分配模型。

完整word版2017工程造价案例分析作业题一及答案

》作业题（一）)《工程造价案例分析(一）（课程代码：03942 一、案例分析题项目为某甲级医院门诊楼，某施工单位拟投标此楼的土建工程。造价师根据该A已知万元，假定管理费按施工企业的定额和招标文件，分析得知此楼需人、材、机合计为2500计，工期4.5%计，规费按6.85%计，综合税率按3.51%人、材、机之和的18%计，利润按为1年，不考虑风险。问题：列式计算A项目的建筑工程费用。（计算过程及计算结果保留两位小数，以万元为单位）万元解：土建工程人、材、机之和：2500 万元管理费：2500*18%=450 ）*4.5%=132.75万元利润：（2500+450*6.85%=211.16万元规费：（2500+450+132.75） *3.51%=115.61万元税金：（2500+450+132.75+211.16）建筑工程费用： 2500+450+132.75+11.16+115.61=3209.52万元 二、案例分析题背景：某办公楼的招标人于2016年10月8日向具备承担该项目能力的A、B、 C、D、E共5家投标单位发出投标邀请书，其中说明，10月12~18日9至16时在该招标人总工程师室领取招标文件，11月8日14时为投标截止时间。该5家投标单位均接受邀请，并按规定时间提交了投标文件。但投标单位A在送出投标文件后发现报价估算有较严重的失误，遂赶在投标截止时间前10分钟递交了一份书面声明，撤回已提交的投标文件。 开标时，由招标人委托的市公证处人员检查投标文件的密封情况，确认无误后，由工作人员当众拆封。由于投标单位A已撤回投标文件，故招标人宣布有B、C、D、E共4家投标单位投标，

工程热力学习题集及答案(1)

工程热力学习题集及答案 一、填空题 1．能源按使用程度和技术可分为 常规 能源和 新 能源。 2．孤立系是与外界无任何 能量 和 物质 交换的热力系。 3．单位质量的广延量参数具有 强度量 参数的性质，称为比参数。 4．测得容器的真空度48V p KPa =，大气压力MPa p b 102.0=，则容器内的绝对压力为 54kpa 。 5．只有 准平衡 过程且过程中无任何 耗散 效应的过程是可逆过程。 6．饱和水线和饱和蒸汽线将压容图和温熵图分成三个区域，位于三区和二线上的水和水蒸气呈现五种状态：未饱和水 饱和水 湿蒸气、 干饱和蒸汽 和 过热蒸汽 。 7．在湿空气温度一定条件下，露点温度越高说明湿空气中水蒸气分压力越 高 、水蒸气含量越 多 ，湿空气越潮湿。（填高、低和多、少） 8．克劳修斯积分/Q T δ? 等于零 为可逆循环。 9．熵流是由 与外界热交换 引起的。 10．多原子理想气体的定值比热容V c = g 72R 。 11．能源按其有无加工、转换可分为 一次 能源和 二次 能源。 12．绝热系是与外界无 热量 交换的热力系。 13．状态公理指出，对于简单可压缩系，只要给定 两 个相互独立的状态参数就可以确定它的平衡状态。 14．测得容器的表压力75g p KPa =，大气压力MPa p b 098.0=，则容器内的绝对压力为 173a KP 。 15．如果系统完成某一热力过程后，再沿原来路径逆向进行时，能使系统和外 界都返回原来状态而不留下任何变化，则这一过程称为可逆过程。 16．卡诺循环是由两个 定温 和两个 绝热可逆 过程所构成。 17．相对湿度越 小 ，湿空气越干燥，吸收水分的能力越 大 。（填大、小） 18．克劳修斯积分/Q T δ? 小于零 为不可逆循环。 19．熵产是由 不可逆因素 引起的。 20．双原子理想气体的定值比热容p c = 72g R 。 21．基本热力学状态参数有：（ 压力）、（温度 ）、（体积）。 22．理想气体的热力学能是温度的（单值 ）函数。 23．热力平衡的充要条件是：（系统内部及系统与外界之间各种不平衡的热力势差为零 ）。 24．不可逆绝热过程中，由于不可逆因素导致的熵增量，叫做（熵产）。 25．卡诺循环由（两个可逆定温和两个可逆绝热 ）热力学过程组成。 26．熵增原理指出了热力过程进行的（方向 ）、（限度）、（条件）。 31．当热力系与外界既没有能量交换也没有物质交换时，该热力系为_孤立系_。 32．在国际单位制中温度的单位是_开尔文_。