首页 机械优化设计阻尼牛顿法

机械优化设计阻尼牛顿法

举报
开通vip

机械优化设计阻尼牛顿法目录第1章选择方法及思路.............................11.1概述...................................................................................................................................1..1.1.11.优化设计.......................................................................

机械优化设计阻尼牛顿法
目录第1章选择方法及思路.............................11.1概述...................................................................................................................................1..1.1.11.优化设计......................................................................................................1.1.21优化设计的 思想 教师资格思想品德鉴定表下载浅论红楼梦的主题思想员工思想动态调查问卷论语教育思想学生思想教育讲话稿 ....................................................................................................1.1.31优化设计的步骤....................................................................................................1.2优化设计的方法...................................................................................................................................1.1.2.1分类.....................................................................................................................1.1.2.22常用的优化方法....................................................................................................第2章阻尼牛顿法计算应.用........................42.1阻尼牛顿法的计算步骤..................................................42.25阻尼牛顿法的程序框图..............................................................................................................2.3实例解析...................................................................................................................................5.2.46阻尼牛顿法的程序编程..............................................................................................................39第章 总结 初级经济法重点总结下载党员个人总结TXt高中句型全总结.doc高中句型全总结.doc理论力学知识点总结pdf ................................................................................................................11.11.1.11.1.2,,,.,,.,,,1.1.2123451.21.2.11122123124121.2.211210.6181:??=?(?1:??)?2f(x)f(x)0.61821nX=[x1,x2,⋯xn]Tminf(X)XX*,minf(X)=f(X*)X*f(X*)2(Powell)1)2)An×nnS1S2S1TAS2=0S1S2A3)(Powell)(Powell)4)n5)6)f(X)X(k)ψ(X)ψ(X)Xψ*X(k+1)22.11)x0k02)f(xk)2f(xk)(2f(xk))1dk(2f(xk))1f(xk)3)xk1xkdkkkkd4)xk1xxxkk1n22.22-12.3()4x2x2fx1,x2x1212x03,0T4(x2)+2(x32x)10f(x0)1128x4x2112x0212x248x50)41442f(x0)114848x02f(x0)1448124112124721247484848511210x112412811f(x)036x131x031816234(x2)3+2(x2x)f(x1)11229791x81413128112(x422)97921x)111219216231f(x)1f(x1)1.940.97274362979147552.4x1882.37712389//.cpp:Definestheentrypointfortheconsoleapplication.#incxl2ude2010.4082x#include#include#includedoublefun1(doubleq1,doubleq2){return(pow((q1-2),4)+pow((q1-2*q2),2));//f(x1,x2)=(x1-2)*(x1-2)*(x1-2)*(x1-2)+(x1-2*x2)*(x1-2*x2)}doublefun2(doubleg,doublex,doubley,doubler1,doubler2){return(pow((x+g*y-2),4)+pow((x+g*y-2*(r1+g*r2)),2));//}voidmain(){doubleA[2][1],B[2][2],C[2][1],D[2][1],X[2][1];doubleE[2][1]={4,3};//x0intt=0,i=0,j=0;doubleE0,x1,x2,x3,h(0.1);doubley1,y2,y3,m;doublea,b,k=0.618,a1,a2,f1,f2;printf("");//std::cin>>E0;do{D[0][0]=E[0][0];D[1][0]=E[1][0];A[0][0]=4*(D[0][0]-2)*(D[0][0]-2)*(D[0][0]-2)+2*D[0][0]-4*D[1][0];A[1][0]=-4*(D[0][0]-2*D[1][0]);//A[0][0],A[1][0]B[0][0]=1.0/(12.0*(D[0][0]-2)*(D[0][0]-2));B[0][1]=1.0/(24.0*(D[0][0]-2)*(D[0][0]-2));B[1][0]=1.0/(24.0*(D[0][0]-2)*(D[0][0]-2));B[1][1]=(6.0*(D[0][0]-2)*(D[0][0]-2)+1)/(48.0*(D[0][0]-2)*(D[0][0]-2));//B[0][0],B[0][1],B[1][0],B[1][1]C[0][0]=-(B[0][0]*A[0][0]+B[0][1]*A[1][0]);C[1][0]=-(B[1][0]*A[0][0]+B[1][1]*A[1][0]);//C[0][0],C[1][0]dk//2x1=0;x2=x1+h;y1=fun2(x1,D[0][0],C[0][0],D[1][0],C[1][0]);y2=fun2(x2,D[0][0],C[0][0],D[1][0],C[1][0]);if(y2>y1){h=-h;x3=x1,y3=y1;x1=x2,y1=y2;x2=x3,y2=y3;}x3=x2+h;y3=fun2(x3,D[0][0],C[0][0],D[1][0],C[1][0]);while(y3=f2){a=a1;a1=a2;f1=f2;a2=a+k*(b-a);f2=fun2(a2,D[0][0],C[0][0],D[1][0],C[1][0]);}else{b=a2;a2=a1;f2=f1;a1=b-k*(b-a);f1=fun2(a1,D[0][0],C[0][0],D[1][0],C[1][0]);}j++;}while(fabs((b-a)/b)>=E0&&fabs((f2-f1)/f2)>=E0);m=0.5*(a+b);//mE[0][0]=D[0][0]+m*C[0][0];E[1][0]=D[1][0]+m*C[1][0];printf("%d%15f10%15f10\n",t,E[0][0],E[1][0],fun1(E[0][0],E[1][0]));t++;}while(fabs(E[0][0]-D[0][0])>=E0&&fabs(E[1][0]-D[1][0])>=E0);X[0][0]=E[0][0];X[1][0]=E[1][0];printf("%d\n",t);printf("(x1,x2)=(%f10,%f10)\n",X[0][0],X[1][0]);printf("f(x1,x2)=%f10\n",fun1(X[0][0],X[1][0]));1.940.973
本文档为【机械优化设计阻尼牛顿法】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
个人认证用户
王淇
热爱文库,热爱新浪。
格式:pdf
大小:232KB
软件:PDF阅读器
页数:10
分类:
上传时间:2023-03-11
浏览量:0