fr共轭梯度法计算例题

最优化方法第四次作业

f(x)?4x1?4x2?4x1x2?12x2。题目：利用FR-共轭梯度法求解无约束优化问题min2x?R22初始点x(0)?(?0.5,1)T.

??g0,k?0;d???k?1?,k?0.??gk??k?1d

Tgg?k?1?kkTgk?1gk?1?k?

一、程序

function [x,val,k]=frcg(fun,gfun,x0) %功能：用FR共轭梯度法求解无约束问题min f（x） %输入：x0是初始点，fun,gfun分别是求目标函数和梯度 %输出：x,val分别是近似最优点和最优值，k是迭代次数 maxk=5000; rho=0.6; sigma=0.4; k=0;

epsilon=1e-4; n=length(x0); while(k

g=feval(gfun,x0);%计算梯度

itern=k-(n+1)*floor(k/(n+1)); itern=itern+1; %计算搜索方向 if(itern==1) d=-g; else

beta=(g'*g)/(g0'*g0); d=-g+beta*d0; g

1 FLUENT变量梯度宏

C_R_G C_P_G C_U_G C_V_G C_W_G C_T_G C_H_G C_YI_G C_R_RG C_P_RG C_U_RG C_V_RG C_W_RG C_T_RG C_H_RG C_YI_RG FLUENT计算中单元变量有两种梯度：非限制梯度（后缀_G）和重建梯度（后缀_RG），帮助文档建议如果需要通过单元变量重建计算单元面上的变量则使用重建梯度计算更精确。 说明：

（1） 只有当求解器正在求解某一变量时才能读取该变量的梯度。例如，定义一个能量源

项时，可以读取温度梯度（C_T_G），但是不能读取速度梯度（如C_U_G）。这是因为求解器不断释放它不需要的变量存储。如果你要保留这些变量梯度，可以使用TUI命令: solve/set/expert 设置“Keep temporary solver memory from being freed?” [Yes]。当然这样做，计算需要更多内存。

（2） C_R_G 只有density-based solver中才能使用； C_P_G只有pressure-based solver中

才能使用。

用挣值法控制成本

（一）三个费用值挣值法是通过分析项目成本目标实施与项目成本目标期望之间的差异，从而判断项目实施的费用、进度绩效的一种方法。挣值法主要运用三个成本值进行分析，它们分别是已完成工作预算成本、计划完成工作预算费用和已完成工作实际成本。

1 ．已完成工作预算成本

已完成工作预算成本为 BCWP ，是指在某一时间已经完成的工作（或部分工作），以批准认可的预算为标准所需要的成本总额。由于业主正是根据这个值为承包商完成的工作量支付相应的成本，也就是承包商获得（挣得）的金额，故称挣得值或挣值。 BCWP =已完成工程量 ×预算成本单价

2 ．计划完成工作预算成本

计划完成工作预算成本，简称 BCWS ，即根据进度计划，在某一时刻应当完成的工作（或部分工作），以预算为标准计算所需要的成本总额。一般来说，除非合同有变更， BCWS 在工作实施过程中应保持不变。 BCWS =计划工程量×预算成本单价

3 ．已完成工作实际成本

已完成工作实际成本，简称 ACWP ，即到某一时刻为止，已完成的工作（或部分工作）所实际花费的成本金额。

（二）挣值法的计算公式在三个成本值的基础上，可以确定挣值法的四个评价指标，它们也都是时间的函数。

1 ．成本偏

一、读识方格网图

方格网图由设计单位（一般在1: 500的地形图上）将场地划分为边长

a=10?40m 的若干

方格,与测量的纵横坐标相对应 ，在各方格角点规定的位置上标注角点的自然地面标高 （H ）和设

计

标高（Hn ）,如图1-3所示.

图1-3

方格网法计算土方工程量图

二、场地平整土方计算

考虑的因素：

① 满足生产工艺和运输的要求；

② 尽量利用地形，减少挖填方数量；

③ 争取在场区内挖填平衡，降低运输费;

+0 39

2 *0 03

3 0 19

4

-Q 朴

5

羽一阴

4? 67

K

44対 ?.81

44

4|昶 00

\

-MM

-270 04

+2SW

\

诃砧

7 +0S0

、

L E -0 05

9

w

弗 59

4? 35

4J &5

43

4』|7 4J 77

44 67

-35 29

+ 142 71

*

-I56G7

T.97

n +Q.7I

13 + 0.44

+o.oe

42.94 地血杯矗/

、復廿札简

42 *0 43(51

4S 1? 43 47 J3 “

43 7? 44 1? 43 7?

42.58

④有一定泄水坡度，满足排水要求 .

⑤场地设计标高一般在设计文件上规定，如无规定:

A. 小型场地一一挖填平衡法；

B. 大型场地一一最佳平面设计法（用最小二乘法，使挖填平衡且总土方量最小）。

1、初步标高（按

xxxx大学实验报告

学生姓名 ＿xxx＿ 学 号＿xxxxxxx＿ 年级班级＿xxxxxxx＿ 实验项目＿xxxxxxxx＿ 实验时间＿xxxxxxxxx＿

实验二

一、实验目的:

1. 充分熟悉复指数函数find、sigshift、sigfold函数的使用； 2. 熟悉序列的加、减、乘、除、移位、折叠的计算； 3. 能够画出结果的图形。 二、实验步骤:

1. 用help查找find、sigshift、sigfold函数的使用情况； 2. 编辑并生成函数sigadd.m(序列相加)

function [y,n] = sigadd(x1,n1,x2,n2) % 实现 y(n) = x1(n)+x2(n) % [y,n] = sigadd(x1,n1,x2,n2)

% y = 在包含n1 和 n2 的n点上求序列和 % x1 = 在 n1上的第一序列

% x2 = 在 n2上的第二序列(n2可与 n1不等)

n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度

仅供个人参考

西京学院数学软件实验任务书

课程名称 For personal use only in study and research; 07 not for commercial use 学号 姓名 commercial use 李亚强 For personal use only in study and research; not for 数学软件实验 班级 数0901 实验课题 线性方程组的最速下降法与共轭梯度法 不得用于商业用途

仅供个人参考

For personal use only in study and research; 熟悉线性方程组的最速下降法与共轭梯度法 not for commercial use 实验目的 运用Matlab/C/C++/Java/Maple/Mathematica等其中实验要求 一种语言完成 不得用于商业用途

仅供个人参考

线性方程组的最速下降法 线性方程组的共轭梯度法 实验内容 成绩 教师 不得用于商业用途

仅供个人参考

实验五实验报告

一、实验名称：最速下降法与共轭梯度法解线性方程组。 二、实验目的：进一步熟悉理解掌握最速下降法与共轭梯度法解法思路，提高matlab编程能力。

地基沉降量计算

地基变形在其表面形成的垂直变形量称为建筑物的沉降量。

在外荷载作用下地基土层被压缩达到稳定时基础底面的沉降量称为地基最终沉降量。

一、分层总和法计算地基最终沉降量

计算地基的最终沉降量，目前最常用的就是分层总和法。

（一）基本原理

该方法只考虑地基的垂向变形，没有考虑侧向变形，地基的变形同室内侧限压缩试验中的情况基本一致，属一维压缩问题。地基的最终沉降量可用室内压缩试验确定的参数（ei、Es、a）进行计算，有：

变换后得：

或

式中：S--地基最终沉降量（mm）；

e1--地基受荷前（自重应力作用下）的孔隙比；

e2--地基受荷（自重与附加应力作用下）沉降稳定后的孔隙比； H--土层的厚度。 计算沉降量时，在地基可能受荷变形的压缩层范围内，根据土的特性、应力状态以及地下水位进行分层。然后按式（4-9）或（4-10）计算各分层的沉降量Si。最后将各分层的沉降量总和起来即为地基的最终沉降量：

（二）计算步骤

1）划分土层

如图4-7所示，各天然土层界面和地下水位必须作为分层界面；各分层厚度必须满足Hi≤0.4B（B为基底宽度）。 2）计算基底附加压力p0 3）计算各分层界面的自重应力σ4）确定压缩层厚度 满足σz=0.2σ

szsz

xxxx大学实验报告

学生姓名 ＿xxx＿ 学 号＿xxxxxxx＿ 年级班级＿xxxxxxx＿ 实验项目＿xxxxxxxx＿ 实验时间＿xxxxxxxxx＿

实验二

一、实验目的:

1. 充分熟悉复指数函数find、sigshift、sigfold函数的使用； 2. 熟悉序列的加、减、乘、除、移位、折叠的计算； 3. 能够画出结果的图形。 二、实验步骤:

1. 用help查找find、sigshift、sigfold函数的使用情况； 2. 编辑并生成函数sigadd.m(序列相加)

function [y,n] = sigadd(x1,n1,x2,n2) % 实现 y(n) = x1(n)+x2(n) % [y,n] = sigadd(x1,n1,x2,n2)

% y = 在包含n1 和 n2 的n点上求序列和 % x1 = 在 n1上的第一序列

% x2 = 在 n2上的第二序列(n2可与 n1不等)

n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度

第一章法调整经济关系的一般理论

一、养马专业户邓某发现路边有一匹无人看管的马，把马带回家。回家后到处打听都没有发现附近有谁丢了马，于是就将该马与自己的马一起饲养，没过几天，马配了种，怀上了小马驹。不久，邻村的何某找到邓某，说明是自己丢了马，经确认丢失的母马确实为何某所有。何某要求把马领走，但邓某认为母马可以领走，但母马产下的马驹应归自己所有，因为是由于自己的精心饲养和照料才使母马怀上马驹，但何某认为自己的母马生下的马驹理所当然属于自己所有，双方争执不下。

1．邓某照料马的行为法律上属于什么行为？请说明理由。

2．邓某要求日后产的马驹归自己所有的主张能否成立？请说明理由。

3．假如在邓某照料马时，马突然发狂，踢伤邓某，邓某某为此支付医疗费200元，应如何处理？

4．假如在邓某照料马期间，马突然遭雷击死亡，几天后何某某找到邓某要求赔偿，如何处理？

答案：答：(1)邓某某照料马的行为构成无因管理。根据法律规定，没有法定或者约定的义务，为避免他人利益受损失进行管理或者服务的，构成无因管理。因此，邓某某为他人照料丢失的马的行为构成无因管理，作为无因管理人，其有权要求受益人偿付因此而支付的必要费用。

(2)邓某某要求日后产下的马驹

功能梯度材料

摘要: 功能梯度材料(FGM)是20世纪80年代出现的一种新型功能材料，它具有耐高温、耐磨损等许多优良性能，由于其结构和性能的优异特性,已成为材料领域研究的热点。本文主要介绍梯度功能材料的概念和特性，着重介绍梯度功能材料的制备方法，对国内外功能梯度材料的研究进展进行了综述,并重点介绍梯度功能材料现阶段的应用及发展趋势，对其前景的展望。 关键词：功能梯度材料、制备方法、应用前景 1．

引言

随着国防高科技的发展,超耐热材料在航天和航空工业中扮演着极其重要的角色。由于航空器在大气圈内高速飞行,机头尖端和机器发动机燃烧室中的内壁温度高达2000Κ以上,使材料处于极其恶劣的服役条件下,一方面,机体材料的外表面和燃烧室的内壁要具有超耐热性能，另一方面，机体材料的内表面和燃烧室的外壁要分别具有低温和超低温性能。而传统的金属材料和超耐热材料都难以满足这种服役条件。为此,人们试图将金属材料与超耐热材料复合,或在金属材料表面涂覆超耐热涂层以达到上述要求,但是,由于两种性质不同的材料的热膨胀系数相差太大,导致这两种材料接合界面的热应力太大, 因此使复合材料或涂层脱落、龟裂而失效。[1]

1987年, 日本科学技术厅航空宇宙技术研究所和东北大学