




已阅读5页,还剩1页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
#include #define N 2float gs(float zN)float f;/f=10*(z0+z1-5)*(z0+z1-5)+(z0-z1)*(z0-z1);f=4+4.5*z0-4*z1+z0*z0+2*z1*z1-2*z0*z1+z0*z0*z0*z0-2*z0*z0*z1;/f=1.5*z0*z0+0.5*z1*z1-z0*z1-2*z0;return(f);float ywyh(float xN,float t,float sN,float zN)float q=0.618,e,AN,f2,a,b,c=0.1,d=0.1;int i,j=0;a=0;b=t;f0=gs(x);for(i=0;iN;i+)zi=xi+b*si;f1=gs(z);if(f1f0)dot=2*t;b=b+t;f0=f1; for(i=0;if1) a=b-t; else break;while(f0f1);elset=0-t;doa=a+t;f1=f0; for(i=0;if0)b=a-t;t=2*t;else break;while(f1f0);doj=j+1;A0=b-q*(b-a);for(i=0;iN;i+)zi=xi+A0*si;f0=gs(z);A1=a+q*(b-a);for(i=0;if1)a=A0;A0=A1;f0=f1; A1=a+q*(b-a); for(i=0;i50) break;elseb=A1;A1=A0;f1=f0;A0=b-q*(b-a); for(i=0;i50) break;if(f1-f0)/f10)e=(f1-f0)/f1;elsee=0-(f1-f0)/f1;while(ec);if(A1-A0)/A00)e=(A1-A0)/A0;elsee=0-(A1-A0)/A0;if(e0)e=(A1-A0)/A0;elsee=0-(A1-A0)/A0;while(ed);if(f0f1)for(i=0;iN;i+)zi=xi+A0*si;return(f0);elsefor(i=0;iN;i+)zi=xi+A1*si;return(f1);void main()float xoN,ssNN,zN,c,d,t,fmin,xoooN, xooN,fiN,sN,xN,xeN,fN,fio,dm,dml,df,dx,dxl;int i,j,k,ks,m;for(i=0;iN;i+)printf(请输入初始点x%d=,i);scanf(%f,&xoi);printf(请输入迭代精度Ef,Ex=);scanf(%f,%f,&c,&d);printf(请输入步长T=);scanf(%f,&t);k=0;ks=0;for(i=0;iN;i+)xoooi=xoi;for(i=0;iN;i+)for(j=0;jN;j+)if(i=j) ssij=1; else ssij=0;dofor(i=0;iN;i+)xooi=xoooi;k=k+1;for(i=0;iN;i+) for(j=0;jN;j+) sj=ssij;fmin=ywyh(xooo,t,s,z);for(j=0;jN;j+)xj=zj;fii=fmin;for(j=0;jN;j+)xoooj=xj;for(i=0;iN;i+)si=xi-xooi;xei=2*xi-xooi;f0=gs(xoo);f1=fmin;f2=gs(xe);fio=f0;dm=fio-fi0;m=0;for(i=0;iN-1;i+)dml=fii-fii+1;if(dmdml)dm=dml;m=i+1;else ;if(f0+f2-2*f1)*(f0-f1-dm)*(f0-f1-dm)0.5*dm*(f0-f2)*(f0-f2)&f2f0)fmin=ywyh(xooo,t,s,z);f1=fmin;for(j=0;jN;j+)xj=zj;for(i=0;iN-1;i+) if(im)for(j=0;jN;j+) ssij=ssij;elsefor(j=0;jN;j+) ssij=ssi+1j;for(i=0;if2)for(i=0;iN;i+)xi=xei;f1=f2;fmin=f2;else ;for(i=0;i0) df=(f0-f1)/f1;elsedf=0-(f0-f1)/f1;if(xoo0-x00)dx=xoo0-x0;elsedx=x0-xoo0;for(i=1;i0) dxl=xooi-xi; else dxl=xi-xooi;if(dxc|dfd);for(i=0;iN;i+)zi=xi; for(i=0;iN;i+) print
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 明确责任的工作目标设定计划
- 如何提升财务团队的协作效率计划
- 2025年鞋用乳液胶粘剂项目合作计划书
- 2025年医用冷疗项目发展计划
- 2025年其它核材料及相关特殊材料合作协议书
- 远程在线教育平台学习免责协议
- 电动汽车充电桩安装施工合同
- Rac-Ganoderic-acid-C2-生命科学试剂-MCE
- 财务顾问聘用协议
- 工作总结写作培训
- 食用油营销整合规划(含文字方案)
- 2025采煤工作面巷道平面布置示意图
- 庆阳煤炭资源开发调研报告
- 苏教版科学五年级下15《升旗的方法》教案
- 桥博常见问题
- 现代工业发酵调控绪论
- 超高性能混凝土项目立项申请(参考模板)
- 电缆桥架招标文件范本(含技术规范书)
- 试车场各种道路施工方案设计
- 贵州省电梯日常维护保养合同范本
- 泰国主要城市中英泰文对照
评论
0/150
提交评论