《进退法黄金分割法c语言编程》由会员分享,可在线阅读,更多相关《进退法黄金分割法c语言编程(3页珍藏版)》请在技术文库上搜索。
1、进退法#include “stdio.h“ void main() float y,y1,y2,y3,x,x1,x2,x3,h,h0,a,b,k1,k2,k3; scanf(“%f%f%f%f%f“, x=x1;y=k1*x*x+k2*x+k3;y1=y;h=h0;x2=x1+h;x=x2;y=k1*x*x+k2* x+k3;y2=y; if(y1=y2) h=2*h;x3=x2+h;x=x3;y=k1*x*x+k2*x+k3;y3=y; else h=- 1*h;x3=x1;y3=y1;x1=x2;y1=y2;x2=x3;y2=y3;h=2*h;x3=x2+h;x=x3;y=k1 *x*x+。
2、k2*x+k3;y3=y; while(y2=y3)x1=x2;y1=y2;x2=x3;y2=y3;h=2*h;x3=x2+h;x=x3;y=k1*x*x+k2*x+k3;y3=y;if(h0)a=x1;b=x3;printf(“a=%f,b=%f“,a,b);elsea=x3;b=x1;printf(“a=%f,b=%f“,a,b);getchar();getchar(); 黄金分割法#include “stdio.h“ void main() float y,y1,y2,x,x1,x2,a,b,c,k1,k2,k3; scanf(“%f%f%f%f%f%f“, x1=a+0.318*(b-。
3、a);x=x1;y=k1*x*x+k2*x+k3;y1=y;x2=a+0.618*(b-a);x=x2;y=k1*x*x+k2*x+k3;y2=y;for(;)if(y1y2)b=x2;x2=x1;y2=y1;x1=a+0.382*(b- a);x=x1;y=k1*x*x+k2*x+k3;y1=y;elsea=x1;x1=x2;y1=y2;x2=a+0.618*(b- a);x=x2;y=k1*x*x+k2*x+k3;y2=y;if(b-a=c)break;x=0.5*(a+b);y=k1*x*x+k2*x+k3;printf(“x=%f,y=%f“,x,y);getchar();getchar(); 我在 C 语言编的程序定义的 Y=k1X2-k2X+k3,因此可以输入任意的函 数,以 Y=X2-10X+36 为例。 (1)进退法中依次输入 k1 k2 k3 X1 h 的值,取 X1为 2,h 为 1,最后输出两端的值为 3 和 9。 (2)黄金分割法中,依次输入 k1 k2 k3 a b c(c 为最后的精度), 取 a 为 1,b 取为 25,c 取为 0.006,最后输出 X 值为 5.4374,Y 为 11.191319。。