黃金分割法機(jī)械優(yōu)化設(shè)計(jì)C語(yǔ)言程序

1、黃金分割法的優(yōu)化設(shè)計(jì)實(shí)驗(yàn)報(bào)告學(xué)院：機(jī)電工程機(jī)制自動(dòng)化 11-03 班學(xué)號(hào)： 541102010326姓名：劉點(diǎn)點(diǎn)1，黃金分割法的程序流程圖開始給定Q、h、總A 0.618b d)兒一/2i）n + 乂 (h £/)J2 J 先)1 、a11b 口2y2 V1tO + / (/? tf) y2-na2)o力 _ /（丹 _&） 兒一心】）Ib2,對(duì)應(yīng)流程圖的C語(yǔ)言程序F面應(yīng)用C語(yǔ)言程序利用黃金分割法求一元函數(shù)F=xA2+2*x的最優(yōu)解，已知初始區(qū)間為-3,5,取收斂精度e=10-4。C語(yǔ)言程序如下:#in clude<stdio.h>#include<math

2、.h>#define f(x) pow(x,2)+2*x#define M 0.618void main()double y1,y2,x1,x2,x,a,b,e;int n;n=1;printf(" 請(qǐng)輸入收斂精度 e="); scanf("%lf",&e);printf(" 請(qǐng)輸入?yún)^(qū)間左值 a="); scanf("%lf",&a);printf(" 請(qǐng)輸入?yún)^(qū)間右值 b="); scanf("%lf",&b);y2n");%.4lfn

3、",n,a,b,x1,x2,y1,y2);printf("n a b x1 x2 y1 x1=b-M*(b-a);x2=a+M*(b-a);y1=f(x1);y2=f(x2);printf("%d %.4lf %.4lf %.4lf %.4lf %.4lf n=n+;doif(y1<y2)b=x2;x2=x1;y2=y1; x1=b-M*(b-a);y1=f(x1);prin tf("%d%.4lf%.4lf%.4lf%.4lf%.4lf %.4lfn", n,a,b,x1,x2,y1,y2);n=n+;elsea=x1;x1=x2;y1

4、=y2;x2=a+M*(b-a);y2=f(x2);prin tf("%d%.4lf%.4lf%.4lf%.4lf%.4lf%.4lfn", n,a,b,x1,x2,y1,y2);n=n+;while(fabs(b-a)/b)>=e&& fabs(y2-y1)/y2)>=e);x=(a+b)*0.5;prin tf("x=%.5lfn",x); getchar();AAA b 二-Hu.- -UJJ 冃丄戸主fl3.運(yùn)行結(jié)果:鳳昭yr呂 紳左 區(qū)囘右-3.006 -3,6000-1.8326xl1114 8326x2yii94

5、4001151-0.66S4y27.&671-0.98?& -0-30&Ba-1151 -0-98?6-0.?876-0-8881-0.&6E4-1.3868-0.8&04-0_?87&-1.3868-0-6654-1 .1114-0.987&-0-9965-0.9410-0-83S8-0.99G5-0_?730-0.9410-1 .1114-0.9410-1.04t3-1.00-1.8009-0.9410-i-0.9812-1.04C3-0.9812-1.0214-1.0000-1.0214-0.9812-1.0061-0.9?&6-1.anna-1.0000'rcsa Anu <cu to cent in ue假定經(jīng)十二次迭代后已滿足收斂精度要求，則得 x*=1/2(a+b)=1/2(