收藏 分享(赏)

进退法黄金分割法c语言编程.doc

上传人:eukav 文档编号:5998163 上传时间:2019-03-23 格式:DOC 页数:3 大小:20.50KB
下载 相关 举报
进退法黄金分割法c语言编程.doc_第1页
第1页 / 共3页
进退法黄金分割法c语言编程.doc_第2页
第2页 / 共3页
进退法黄金分割法c语言编程.doc_第3页
第3页 / 共3页
亲,该文档总共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;elseh=-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+k2*x+k3;y3=

2、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-a);x=x1;y=k1*x*x+

3、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。

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 网络科技 > 开发文档

本站链接:文库   一言   我酷   合作


客服QQ:2549714901微博号:道客多多官方知乎号:道客多多

经营许可证编号: 粤ICP备2021046453号世界地图

道客多多©版权所有2020-2025营业执照举报