工程结构优化设计—编程作业.docVIP

工程结构优化设计—编程作业 1.用黄金分割法求方程 x3-3x+1=0 在MATLAB中没有专门的函数实现黄金分割法求解线性方程，可通过编写glodf.m函数实现黄金分割法求解，其代码如下: function x=glodf(f,a,b,eps) if(nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if(f1==0) x=a; end if(f2==0) x=b; end if(f1*f20) disp(’两端点函数值乘积大于0!’); return; else t1=a+(b-a)*0.382; t2=a+(b-a)*0.618; f_1=subs(sym(f),findsym(sym(f)),t1); f_2=subs(sym(f),findsym(sym(f)),t2); tol=abs(t1-t2); while(toleps) if(f_1*f_20) a=t1; b=t2; else fa=subs(sym(f),findsym(sym