黄金分割法-机械优化设计-C语言程序

黄金分割法的优化设计实验报告学院：机电工程机制自动化11-03班学号：54110326姓名：刘点点1，黄金分割法的程序流程图2，相应流程图的C语言程序下面应用C语言程序运用黄金分割法求一元函数F=x2+2*x的最优解，已知初始区间为-3,5 ,取收敛精度e=10-4。C语言程序如下：#include #include#define f(x) pow(x,2)+2*x #define M 0.618 void main() double y1,y2,x1,x2,x,a,b,e; int n; n=1; printf(请输入收敛精度e=); scanf(%lf,&e); printf(请输入区间左值a=); scanf(%lf,&a); printf(请输入区间右值b=); scanf(%lf,&b); printf(n a b x1 x2 y1 y2n); x1=b-M*(b-a); x2=a+M*(b-a); y1=f(x1); y2=f(x2); printf(%d %.4lf %.4lf %.4lf %.4lf %.4lf %.4lfn,n,a,b,x1,x2,y1,y2); n=n+; do if(y1=e&fabs(y2-y1)/y2)=e); x=(a+b)*0.5; printf(x=%.5lfn,x); getchar(); 3.运营成果:假定经十二次迭代后已满足收敛精度规定，则得x*=1/2(a+b)=1/2(-1.0214-0.9812)=-1.0013，相应的函数极值f(x*)=-0.9999；近似精确值x*=-1,f(x*)=-1,与解析法求得的精确值相似