电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

机械优化设计黄金分割法 外推法

机械优化设计黄金分割法 外推法_第1页
1/7
机械优化设计黄金分割法 外推法_第2页
2/7
机械优化设计黄金分割法 外推法_第3页
3/7
郑州大学机械优化设计部分程序1、外推法2、黄金分割法3、二次插值法4、坐标轮换法5、随机方向法6、四杆机构优化设计1、外推法源程序:#include#include#define R 0、01double fun(double x){ double m;m=x*x-10*x+36;return m;}void main(){double h0=R,y1,y2,y3,x1,x2,x3,h;x1=0;h=h0;x2=h;y1=fun(x1);y2=fun(x2);if(y2>y1){h=-h;x3=x1;y3=y1;x1=x2;y1=y2;x2=x3;y2=y3;}x3=x2+h;y3=fun(x3);while(y3#include#define f(x) x*x*x*x-5*x*x*x+4*x*x-6*x+60 double hj(double *a,double *b,double e,int *n) { double x1,x2,s; if(fabs((*b-*a)/(*b))<=e) s=f((*b+*a)/2); else { x1=*b-0、618*(*b-*a); x2=*a+0、618*(*b-*a); if(f(x1)>f(x2)) *a=x1; else *b=x2; *n=*n+1; s=hj(a,b,e,n); } return s; } void main() {double s,a,b,e,m; int n=0; printf("输入 a,b 值与精度 e 值\n");scanf("%lf %lf %lf",&a,&b,&e); s=hj(&a,&b,e,&n); m=(a+b)/2; printf("a=%lf,b=%lf,s=%lf,m=%lf,n=%d\n",a,b,s,m,n); }运行过程及结果:输入 a,b 值与精度 e 值-350、0001a=3、279466,b=3、279793,s=22、659008,m=3、279629,n=213、二次插值法源程序:#include#includeint main(void){ double a1,a2,a3,ap,y1,y2,y3,yp,c1,c2,m; double j[3]; int i,h=1; void finding(double a[3]); finding(j); a1=j[0]; a2=j[1]; a3=j[2]; m=0、001; double f(double x); y1=f(a1); y2=f(a2); y3=f(a3); for(i=1;1>=1;i++) { c1=(y3-y1)/(a3-a1); c2=((y2-y1)/(a2-a1)-c1)/(a2-a3); ap=0、5*(a1+a3-c1/c2); yp=f(ap); if(fabs((y2-yp)/y2)0) { if(y2>=yp){ a1=a2;y1=y2; a2=ap;y2=yp;} else{ a3=ap;y3=yp;} } else if(y2>=yp){ a3=a2;y3=y2; a2=ap;y2=yp;} else{a1=ap;y1=yp;} } double x,y; if(y2<=yp){ x=a2;y=y2;} else{ x=ap;y=yp;} printf("a*=%f\n",x); printf("y*=%f\n",y); return 0;}double ...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

机械优化设计黄金分割法 外推法

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部