合工大机械优化设计课程实践报告

合工大机械优化设计课程实践报告

ID:47668405

大小:729.00 KB

页数:18页

时间:2020-01-27

合工大机械优化设计课程实践报告_第1页
合工大机械优化设计课程实践报告_第2页
合工大机械优化设计课程实践报告_第3页
合工大机械优化设计课程实践报告_第4页
合工大机械优化设计课程实践报告_第5页
资源描述:

《合工大机械优化设计课程实践报告》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、.word格式.合肥工业大学《机械优化设计》课程实践研究报告班级:机械设计制造及其自动化12-3班学号:姓名:授课教师:王卫荣日期:2015年11月14日.专业.专注..word格式.目录一、一维搜索程序作业.......................................31.λ=0.618的证明............................................32.编写0.618法程序并计算.....................................4二、单位矩阵程序作业..........

2、............................6三、连杆机构问题和自选工程优化问题..........................71.连杆机构问题.............................................72.自选工程优化问题.........................................14四、课程实践心得体会........................................18.专业.专注..word格式.一、一维搜索程序作业1.λ=0.618的证明黄金分割法

3、,又称作0.618法,适用于[a,b]区间上的任何单谷函数求极小值问题。黄金分割法是建立在区间消去法原理基础上的试探方法,即在搜索区间[a,b]内适当插入两点α1、α2,并计算其函数值。α1、α2将区间分成三段。应用函数的单谷性质,通过函数值大小的比较,删去其中一段,使搜索区间得以缩短。然后再在保留下来的区间上做同样的位置,如此迭代下去,使搜索区间无限缩小,从而得到极小点的数值近似值。黄金分割法要求插入点α1、α2的位置相对于区间[a,b]两端点具有对称性,即图1-1黄金分割法α1=b–λ(b–a)α2=a+λ(b–a)(3-1)其中,λ为待

4、定常数。下面证明λ=0.618。除对称性要求外,黄金分割法还要求保留下来的区间内再插入一点所形成的区间新三段,与原来区间的三段具有相同的比例分布。设原有区间[a,b]长度为1如图1-1所示,保留下来的区间[a,b]长度为λ,区间缩短率为λ。为了保持相同的比例分布,新插入点α3应在λ(1–λ)位置上,α1在元区间的1–λ位置应相当于在保留区间的λ²位置。故有1–λ=λ²即λ²+λ–1=0取方程正数解得.专业.专注..word格式.若保留下来的区间为[α1,b],根据插入点的对称性,也能推得同样的λ的值。2.编写0.618法程序并计算(1)0.6

5、18法程序:#include#includefloatm=0.618;floatfun(floatt){floaty;y=cos(t);returny;}main(){floata,b,eps;printf("min=");scanf("%f",&a);%输入函数下限%printf("max=");scanf("%f",&b);%输入函数上限%floatt1,t2,t,f1,f2,min;printf("eps=");scanf("%f",&eps);%输入精度%while((b-a)/b>=eps){t

6、1=a+(1-m)*(b-a);t2=a+m*(b-a);f1=fun(t1);f2=fun(t2);if(f1>=f2){a=t1;t1=t2;f1=f2;t2=a+m*(b-a);f2=fun(t2);}else{b=t2;t2=t1;f2=f1;t1=a+(1-m)*(b-a);f1=fun(t1);}}t=(a+b)/2;.专业.专注..word格式.min=fun(t);printf("最优点t=%f",t);%输出最优点t%printf("最优值f=%f",min);}%输出最优值f%(2)计算:1.a=0,b=2π,Y=

7、cos(t)2.a=0,b=10,y=(t-2)*(t-2)+3.专业.专注..word格式.二、单位矩阵程序作业编写生成单位矩阵的程序:程序文本#includevoidmain(void){inta[100][100];intN,i,j;printf("请输入所要输出矩阵的阶数(最多100阶):");scanf("%d",&N);printf("输出的矩阵阶数为%d",N);printf("N");/*****制作表头*****/for(i=0;i

8、i+1);printf("");for(i=0;i

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。