高等流体力学
班 级：学 号：姓 名：授课老师：日 期：
机设15学硕班 201 张湘楠 毕新胜 2016年7月 1日
1 / 26
一、研究报告内容：
1、λ＝0.618的证明、一维搜索程序作业； 2、单位矩阵程序作业；
3、连杆机构问题 ＋ 自行选择小型机械设计问题或其他工程优化问题； （1）分析优化对象，根据设计问题的要求，选择设计变量，确立约束
条件，建立目标函数，建立优化设计的数学模型并编制问题程序； （2）选择适当的优化方法，简述方法原理，进行优化计算； （3）进行结果分析，并加以说明。 4、写出课程实践心得体会，附列程序文本。
5、为响应学校2014年度教学工作会议的改革要求，探索新的课程考核评价方法，特探索性设立一开放式考核项目，占总成绩的5%。 试用您自己认为合适的方式（书面）表达您在本门课程学习方面的努力、进步与收获。（考评将重点关注您的独创性、简洁性与可验证性）。
二、研究报告要求
1、报告命名规则：学号－姓名－《机械优化设计》课程实践报告.doc 2、报告提交邮址：（收到回复，可视为提交成功）。
追 求：问题的工程性，格式的完美性，报告的完整性。 不 追 求：问题的复杂性，方法的惟一性。 评判准则：独一是好，先交为好；切勿拷贝。
2 / 26
目录：
λ＝0.618的证明、一维搜索程序作业
① 关于??0.618的证明……………………………………………………4 ② 一维搜索的作业
采用matlab进行编程…………………………………………… 5 采用C语言进行编程……………………………………………… 7 单位矩阵程序作业
① 采用matlab的编程………………………………………………… 9 ② 采用c语言进行编程………………………………………………… 9 机械优化工程实例
① 连杆机构…………………………………………………………………11 ② 自选机构…………………………………………………………………16 课程实践心得…………………………………………………………………… 20 附列程序文本…………………………………………………………………… 21 进步，努力，建议………………………………………………………………25
3 / 26
一、λ＝0.618的证明、一维搜索程序作业
①关于??0.618的证明
黄金分割法要求插入点?1，?2的位置相对于区间[a,b]两端具有对称性，即
?1?b??(b?a)
?2?a??(b?a)
其中?为待定常数。
此外，黄金分割法还要求在保留下来的区间内再插入一点所形成的区间新三段，与原来的区间三段具有相同的比例分布。
黄金分割法还要求在保留下来的区间内再插一点所形成的区间新三段，与原来的区间三段有相同的比例分布。
设原区间[a,b]的长度为1，如图一所示，保留下来的区间[?1,?2]长度为?，区间缩短率为?。为了保持相同的分别比例。插入新点?3应在??1-??位置上，?1在原区间的1-?1位置应相当于在保留区间的?位置。 故有：
2a
α1
1-λ λ α2 b
1-???2
?2??-1?0
5-1?0.618 解得??2
1 a
λ(1-λ) λ2 λ 图二 α3
图一 α1
α2 2 4 / 26
编写0.618的程序，并计算下列问题 NO.1: a?0 b?2? f?x??cosx
NO.2: a?0 b?10 f?x???x?2??3
2（1）采用MATLAB进行编程
%%
%fun.m:黄金分割法求极值点 %输入数据
% a –搜索区间下限 % b –搜索区间上限 % e – 精度 %输出数据： % x – 极小值点 %其他常量： % c1,c2,c3 – 区间
% r –黄金分割比例0.618 %%
%定义函数
function x=fun(a,b,e) r=0.618;
c1=b-r*(b-a);c2=a+r*(b-a); y1=f(c1);y2=f(c2);
while (abs((b-a)/b)>e)&&(abs((y2-y1)/y2)>e) if y1>=y2
a=c1;c1=c2;y1=y2; c2=a+r*(b-a);y2=f(c2); else
b=c2;c2=c1;y2=y1; c1=b-r*(b-a);y1=f(c1); end end
x=0.5*(a+b); end
对f函数的确立
10 function y=f(x)
y=(x-2)^2+3; end
11 function y=f(x)
y=cosx end
5 / 26
程序框图
机械优化设计大作业
