原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一维搜索方法与终止准则实验
实验报告
课程名称：_现代机械设计设计方法__ 项目名称：__一维搜索方法与终止准则____
实验预习部分：
实验目的：
2．实验原理：
3．实验仪器：
4．实验步骤及方法：
实验过程记录部分：
1.
%typbound.m;进退法1
function [lowbound,upbound]=typbound11(x0,step0,startopint,searchdirection)
step=step0;
f0=tryobjfun11(x0,startopint,searchdirection);
x1=x0+step0;
f1=tryobjfun11(x1,startopint,searchdirection);
if f1=f0
while true
step=2*step;
x2=x1+step;
f2=tryobjfun11(x2,startopint,searchdirection);
if f1=f2
lowbound=x0;
upbound=x2;
break;
else
x0=x1;
x1=x2;
f0=f1;
f1=f2;
end
end
else
while true
step=2*step;
x2=x0-step;
f2=tryobjfun11(x2,startopint,searchdirection);
if f0=f2
lowbound=x2;
upbound=x1;
break;
else
x1=x0;
x0=x2;
f1=f0;
f0=f2;
end
end
end
%(2)tryobjfun.进退法2
function f=tryobjfun11(a,startopint,searchdirection)
f=objfun11(startopint+a.*searchdirection);
%(3)确定函数，也就是你要确定搜索区间的目标函数，你也可以改，这仅是一个例子。进退法3
function f=objfun(x)
f=x.^2-6.*x+9
end
2.
function[x,y] = goldmin12(f,xa,xb,s)
% 黄金分割法求解函数最小值
% 输入
% f 待优化函数
% a,b 区间
% s 精度
% 输出
% x 最优解
% y 最优解对应的最小值
g = (sqrt(5)-1)/2; % 黄金分割比，0.618
a = xa;
b = xb;
x2 = a + g*(b-a);
x1 = a + b - x2;
y1 = f(x1);
y2 = f(x2);
while abs(b-a) s
if y1y2
b = x2;
x2 = x1;
x1 = a + b - x2;
end
if y1=y2
a = x1;
x1 = x2;
x2 = a + b - x1;
end
y1=f(x1);
y2 = f(x2);
end
x = x1;
y = f(x);
实验报告成绩（百分制）__________ 实验指导教师签字：__________ 验结果与讨论：
数据处理及实验结果分析：
实验报告成绩（百分制）__________ 实验指导教师签字：__________
文档评论（0）