现代设计黄金分割法复合形法实验报告word文档良心出品
3DMAX实习报告总结word文档良心出品
篇一： 3dmax 实训报告计算机 3dsmax 实习报告 室外三维建模所谓三维设计就是利用电脑进行设计与创作， 以产生真实的立体场景与动画。
虽然在这个学 期对3dsmax 的接触才知道3dsmax 是一个广泛应用于游戏开发、后期制造、影视特效及专业 视觉设计领域的一款功能强大的三维设计软件，它是集专业建模、动画、渲染一休的三维解 决方案。
特别是3damax 以其强大的功能卓越的表现力被广大电脑设计人所睛睐，成为当今较热门的设计软件。
在学习3dsmax 的这个学期里，3dsmax 创作的每一个细节都在刺激着我的神经，变幻无穷的 3d 建模，它的每一个创造都给我带来无比的震憾的惊喜与灵感 ?? 今天我要用自己所学到的知识来创造一个初级建模（室外建模）虽然这是我第一次曾未有的 偿试；不过相信抱着对 3d 的渴望，相信3dsmax 可以帮助我将难度复杂的室外模型简单地实 现出来。
下面是此次实习室外建模的几大重要路径：1 创建模型2 修改3 装饰4 环境渲染现在我来跟大家分享下这个初级建模的基本制造： （所谓初级建模： 是指利用几何体创造面板， 二维图形创建面扳中的现有模型来进行的建模操作，包括标准几何体的创建、扩展几何体创 建和二维图形的创建） 。
首先：一 创建模型1击创建命令面板中单击“图形”按钮，“创建”选择“线”按钮。
利用创建二维图形在顶视图描绘出建筑所在定的位置以及大小。
2 单击创建面板，单击“几何体”按钮， 进入“标准基本体”创建面板单击“长方体” 按钮，在“顶视图”刚创建好二维图形（平面图），照大小位置拉出几何体，其 它照样拉好之后，如右图：3 选择命令面板，单击“修改”按钮，打开修改面板设置长方体的参数，在这里主要设置高度即可， 为了方便修理在每一个几何体都给予一个名称。
因此先来设置房 1 的高度参数，（自 己认为合适即可） 。
4 单击“选择按钮” ，选择顶视图的的房 2 几何体； 照前步骤一样来设置几何体的高度。
最优化课程设计--黄金分割法及其算法实现(3
机械优化设计报告姓名：刘洋学号：S12080203054院系：机械工程学院专业：机械设计及理论2012年 12月 4日摘要最优化理论和方法日益受到重视，已经渗透到生产、管理、商业、军事、决策等各个领域，而最优化模型与方法广泛应用于工业、农业、交通运输、商业、国防、建筑、同学、政府机关等各个部门及各个领域。
伴随着计算机技术的高速发展，最优化理论与方法的迅速进步为解决实际最优化问题的软件也在飞速发展。
其中，MATLAB软件已经成为最优化领域应用最广的软件之一。
有了MATLAB 这个强大的计算平台，既可以利用MATLAB优化工具箱（OptimizationToolbox）中的函数，又可以通过算法变成实现相应的最优化计算。
关键词：优化、黄金分割法、最速下降法、MATLAB、算法AbstractOptimization theory and methods and more attention, have penetrated into the production, management, business, military, decision-making and other fields, and optimization models and methods widely used in industry, agriculture, transportation, commerce, defense, construction, students, government various departments and agencies and other fields. With the rapid development of computer technology,optimization theory and methods for the rapid progress of the optimization problem to solve practical software is also developing rapidly. Which, MATLAB software has become the most optimization software is one of the most widely used. With this powerful computing platform MATLAB, either using MATLAB optimization toolbox (OptimizationToolbox) in the function, but also can achieve the appropriate algorithm to optimize into the calculation.Key words: Optimization、Golden section method、steepest descent method、MATLAB、algorithm目录摘要 (2)第一章绪论 (5)第二章黄金分割法的基本思想与原理 (6)2.1 黄金分割法的基本思路 (6)2.2 算法流程图 (7)2.3 用matlab编写源程序 (7)2.4 黄金分割法应用举例 (8)第三章最速下降法的基本思想与原理 (9)3.1 最速下降法的基本思路 (9)3.2 算法流程图 (11)3.3 用matlab编写源程序 (11)3.4 最速下降法应用举例 (13)第四章惩罚函数法的基本思想与原理 (13)4.1 惩罚函数法的基本思路 (13)4.2 算法流程图 (14)4.3 用matlab编写源程序 (14)4.4 最速下降法应用举例 (16)第五章总结 (17)参考文献 (18)第1章绪论在人类活动中，要办好一件事（指规划、设计等），都期望得到最满意、最好的结果或效果。
最优化方法(黄金分割和进退法)实验报告
一维搜索方法的MATLAB 实现姓名： 班级：信息与计算科学 学号： 实验时间: 2014/6/21 一、实验目的：通过上机利用Matlab 数学软件进行一维搜索，并学会对具体问题进行分析。
并且熟悉Matlab 软件的实用方法，并且做到学习与使用并存，增加学习的实际动手性，不再让学习局限于书本和纸上，而是利用计算机学习来增加我们的学习兴趣。
二、实验背景： 黄金分割法它是一种基于区间收缩的极小点搜索算法，当用进退法确定搜索区间后，我们只知道极小点包含于搜索区间内，但是具体哪个点，无法得知。
1、算法原理黄金分割法的思想很直接，既然极小点包含于搜索区间内，那么可以不断的缩小搜索区间，就可以使搜索区间的端点逼近到极小点。
2、算法步骤用黄金分割法求无约束问题min (),f x x R ∈的基本步骤如下：（1）选定初始区间11[,]a b 及精度0ε>，计算试探点： 11110.382*()a b a λ=+- 11110.618*()a b a μ=+-。
（2）若k k b a ε-<，则停止计算。
否则当()()k k f f λμ>时转步骤（3）。
当()()k k f f λμ≤转步骤（4）。
（3） 11111110.382*()k k k k k k k k k k a b b a b a λλμμ+++++++=⎧⎪=⎪⎨=⎪⎪=+-⎩转步骤（5）（4）转步骤（5）（5）令1k k =+，转步骤（2）。
算法的MATLAB 实现function xmin=golden(f,a,b,e) k=0;x1=a+0.382*(b-a); x2=a+0.618*(b-a); while b-a>e f1=subs(f,x1); f2=subs(f,x2); if f1>f2 a=x1; x1=x2; f1=f2;x2=a+0.618*(b-a); else b=x2; x2=x1; f2=f1;x1=a+0.382*(b-a); end k=k+1; endxmin=(a+b)/2; fmin=subs(f,xmin)fprintf('k=\n'); disp(k);3、实验结果(总结/方案)黄金分割法求解极值实例。
(完整word版)黄金分割在生活中的应用
研究报告黄金分割在生活中的应用东北育才学校马艺宸一．黄金分割的定义之比等于整体与较大部分之比，其比值为1∶0.618或1。
618∶1，即长段为全段的0.618.0.618被公认为最具有审美意义的比例数字.上述比例是最能引起人的美感的比例，因此被称为黄金分割。
二．黄金分割在生活中的应用（一）艺术中的黄金分割1。
人体上的黄金分割。
最完美的人体：肚脐到脚底的距离/头顶到脚的距离=0.618。
最漂亮的脸庞：眉毛到脖子的距离/头顶到脖子的距离=0。
618。
达·芬奇的《蒙娜丽莎》、拉斐尔笔下温和俊秀的圣母像，都有意无意地用上了这个比值。
人们公认的最完美的脸型——“鹅蛋"形，脸宽与脸长的比值约为0.618，如果计算一下翩翩欲仙的芭蕾演员的优美身段，可以得知,他们的腿长与身长的比值也大约是0.618,组成了人体的美.2.中国最古老的古琴，处处透着黄金分割的神奇，琴背两池，左龙右凤。
控制琴弦发音的枢纽有三：轸，凫掌，凤嗉.琴有五弦，音有八度,琴节为徽。
“以琴长全体三分损一，又三分益一，而转相增减”，全弦共有十三徽。
把这些排列到一起,二池，三纽，五弦,八音，十三徽。
多么奇妙的排列，恰是费波那奇数，而两个相邻费波那奇数比率则越来越接近黄金分割率，是有意还是巧合？看来,中国古人对黄金分割的领悟与运用，与西方确有异曲同工之妙.3.1483年左右，达芬奇画的一副未完成的油画，包围着圣杰罗姆躯体的黑线，就是一个黄金分割的矩形，当时达芬奇似乎有意利用这一黄金分割的比值.“检阅”是法国印象派画家舍勒特的一副油画，它的画杠结构比例也正是0.618的比值。
英国在画家斐拉克曼的名著《希腊的神话和传说》一书中，工绘有96幅美人图。
每一幅画上的美人都妩媚无比婀娜多姿.如果仔细量一下她们的比例也都也雅典娜相似。
4。
音乐家发现，二胡演奏中,“千金"分弦的比符合0。
618∶1时，奏出来的音调最和谐、最悦耳。
5。
希腊古城雅典有一座用大理石砌成的神妙，神庙大殿中央的女神像是用象牙和黄金雕成的。
黄金分割法及其应用
黄金分割法及其应用黄金分割法及其应用黄金分割法，又称为黄金比例、黄金分割比等，是一种比例关系，源自于古希腊文化。
它指的是，将一条线段分割为两部分，使其中一部分与另一部分之和的比等于整条线段与其中一部分的比。
这个比例值被称为“黄金分割比”，通常表示为1:φ（phi），φ是一个无理数，约等于1.6180339887。
应用黄金分割法在设计、艺术、建筑等领域广泛应用，被认为是一种非常美学的比例关系。
以下是一些常见的应用方法：1. 黄金矩形黄金矩形是一种矩形，其长和宽按照黄金分割比例进行分割。
这种矩形具有一种非常美学的形态，被广泛应用于设计和艺术领域。
例如，著名的维特鲁威斯男爵的画作中，经常使用黄金矩形比例来构图。
2. 身体比例黄金分割法在人体比例上也有应用。
例如，人体的身高和臂展、腿长等比例，都可以按照黄金分割比例进行分割。
这种比例关系在雕塑和肖像绘画中经常被使用，可以使得作品更加真实生动，具有感染力。
3. 建筑设计建筑中的黄金分割法也常常应用。
例如，建筑的外观比例、窗户的位置和尺寸等都可以按照黄金分割比例进行分配。
这种比例关系能够创造一种和谐而宁静的感觉，符合人们的审美标准。
4. 广告设计广告设计中常常也会使用黄金分割法。
例如，在广告中，图片、文字和背景的比例、位置、大小等都可以进行合理的黄金分割设计，从而产生更好的视觉效果。
5. 网页设计在网页设计中，黄金分割法也是一种比较常用的设计原则。
例如，网页布局、按钮大小、文本位置等都可以按照黄金分割设计，这样可以让网页看起来更加优美和协调。
总结黄金分割法是一种非常美学的比例关系，被广泛应用于各个领域。
黄金分割法比例的应用可以让设计更加美观和协调，符合人们的审美标准，从而产生更好的视觉效果和感官体验。
基于黄金分割法的产品造型设计探究
产 品造型设计是设计师 使用者和产品三者建立联系的纽带 ， 通 过产品外观形态 的设计 ，设 计 师利用产品唤起使 用者 内心 的情感 ， 使得 产品更易于被使用者接受和喜 爱。 来源于 自然界而在 生活 中广泛存在的黄金分割比例 ，已经成 了根深 蒂固的审美法则。而且科学研究表明 ， 黄金分割 比例是最和谐 最协 调的比例 ， 深深影响着人的视觉 系统 ， 给人带来直观的美 的感受 。在产品设 计中， 为塑造具有视觉 吸引力的产品外观 ， 提升产品外观造 型的美感价值 ， 黄金分割 比例以其先 天 的美感法则及数的规律性将成为简 易可行而又卓有成效的设计 法则 。
Ａ
Ｃ
Ａ Ｏ： Ａ Ｂ ＝ Ｂ Ｃ： Ａ Ｃ≈ ０ ． ６１ ８
Ｂ
以及在整个过程 中应 该遵 循怎样的设计原则 ， 并探究 出一种较
完善且易操作的方法 ， 引导设计师利用黄金分割法达到我们对
产 品 外 观 进 行 设 计 比例 的美学价 值 黄金分割 比例的美感来源于人类 的认知偏好 。 人类通过 观
产品造型是产 品给使用者最直观的 感受 ， 因此它的重要性
特殊 美学属性产生 的心理反应 。他将实验对象 限定在人造物一 可见一斑。产 品造 型的设计方法也有很 多，需 要综合考虑 内部
一
产 品中 ，并从 中测量数以千计 的各种矩 形物体 ，结果发现大 构件 、外部环境 、 人 机工程 、 材料 、工艺、成本等众 多因素 ，因
中 ，在鱼 类 的身 形构造 中，总之 ， 在 自然 界里 处处 都 可以 看 以充分的刺激 ，黄金分割比例深深影响 了人的视觉 系统 ，使人 到 黄金 分割 比例 的踪 迹 ；巨石阵 、金 字塔 、泰姬 陵这 些人 造 产生 了与生俱来的视觉美感 。 建筑 ，其 中也同 样蕴 含 了黄金 分割 的数 理 关系 。细 细想来 ， 黄金分割 比例 的美学价值是 经过 实验证明了的。１ ９ 世纪德
黄金分割法在大分子设计实践中的应用
黄金分割法在聚硅烷树脂粘着剂中的应用
黄金分割法在聚硅烷树脂粘着剂中的应用
实例2：黄金分割法在AM/AMPS共聚物合成中的应用
采用过硫酸铵和亚硫酸氢钠氧化还原体系引发 丙烯酰胺( AM) 和2-丙烯酰胺-2-甲基丙磺酸 ( AMPS) 水溶液进行二元共聚，得到 AM/AMPS共聚物。黄金分割法评定确定最佳 合成条件为: 引发剂0．002888%，单体20%， AMPS∶ AM=20∶80( 质量百分比) ，反应温 度45．52℃。 实验结果测定的最佳实验条件与黄金分割法评 定的一致。
/view/1816.htm
实例1：黄金分割法在聚硅烷树脂粘着剂中的应用
制备一种油漆，将聚硅烷溶于二甲苯。为了提高粘着 剂的粘着力，现需要掺杂固体填料，如石墨、MoS2、 滑石。
如何确定最佳添加剂含量？
运用黄金分割法逐步限定含量范围黄金分割法
Billiau-Loreau, Myriam; Delacroix, Alain; Porte, Catherine.Optimization ofan adhesive formulation by experimental designs. Process Control and Quality (1999), 11(4), 313-321.
黄金分割法在摩擦材料配方设计中的应用
原理：
1.黄金分割序列黄 金分割率体现了部 分与部分的比值等 于部分与全体的比 值的最佳比例关系， 继续分割线段，可 得黄金分割序列
黄金分割法在摩擦材料配方设计中的应用
2、摩擦材料配方设计
黄金分割法在摩擦材料配方设计中的应用
3、摩擦材料配方优化
借助灰色相关度分析，可分析每个组分对摩擦性能的正或负贡献
现代设计方法实验报告
现代设计方法实验报告篇一：现代设计方法实验报告《现代机械设计方法学》实验报告班级：08机设（4）班学号：XX 姓名：李成成绩：景德镇陶瓷学院机电学院实验一、有限元分析（一）目的：1、初步掌握有限元软件分析力学问题的过程，包括几何建模、网格划分等前处理功能，掌握各种计算结果的阅读。
2、掌握材料数据、载荷、约束的添加方法。
（二）要求：学生独立完成一个算例的有限元分析，并阅读其计算结果，提交一个算例的分析报告。
（三）计算实例 1、问题的描述为了考察铆钉在冲压时，发生多大的变形，对铆钉进行分析。
铆钉圆柱高：10mm 铆钉圆柱外径：6mm 铆钉下端球径：15mm 弹性模量：2.06E11 泊松比：0.3铆钉材料的应力应变关系如下：1、有限元模型。
3、应力云图，可选主应力或σx、σy、τxy、Von Mises应力、Tresca应力之一输出结果图片，指明你所选的应力的最大值及其位置。
(本文来自：小草范文网:现代设计方法实验报告) （三）思考题：1、如果要提高边界处计算精度，一般应如何处理？在边界处划分网格2、有限元网格划分时应注意哪些问题？选取的时候要将编号显示出来，这样就可以更好的选择，网格尽可能的小，这样结果就越准确。
实验二、优化实验（一）目的：初步掌握利用ANSYS软件或MATLAB软件对问题进行分析。
（二）要求：学生独立完成一个算例的分析，并给出算例的计算结果。
（三）算例1. 实际问题梁的形状优化，优化目的是使梁的体积最小，同时要求梁上的最大应力不超过30000psi，梁的最大挠度不大于0.5in,沿长度方向梁的厚度可以变化，但梁端头的厚度为定值t，采用对称建模。
使用两种方法进行优化，两种方法优化结果。
篇二：现代设计方法实验报告准考证号：1370姓名：倪帅彪主考院校：河南科技大学专业名称：080302机械制造及自动化（独立本科段）现代设计方法实验报告实验一AutoCAD使用的基本知识一、实验目的与要求：（1）掌握AutoCAD的安装和起动（2）了解AutoCAD操作界面组成二、实验设备：AutoCAD安装软件、多媒体电脑等。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《现代设计理论与方法》实验报告、实验目的机械优化设计是一门实践性较强的课程，学生通过实际上机计算可以达到以下目的：1. 加深对机械优化设计方法的基本理论和算法步骤的理解;2. 培养学生独立编制或调试计算机程序的能力；3. 掌握常用优化方法程序的使用方法；4 .培养学生灵活运用优化设计方法解决工程实际问题的能力。
、实验项目、学时分配及对每个实验项目的要求1.明确黄金分割法基本原理、计算步骤及程序框图; 吐入「土 2•编制或调试黄金分割法应用程序； 1 黄金分割法 2八' " 3 •用测试题对所编程序进行测试；4•撰写实验报告。
1.明确复合形法基本原理、计算步骤及程序框图 等；2 复合形法 4 2•编制或调试复合形法应用程序；3 •用测试题对所编程序进行测试；4•撰写实验报告。
二、测试题1. 黄金分割法程序测试题1)rn"何二?-10r+36,取坷=0 ,卜皿1, 沪程序如下:#in clude<stdio.h> #in clude<c oni o.h> #in clude<math.h> #defi ne e 0.00001 序实验项目 学时号实验要求#define tt 0.01float function(float x)float y=pow(x,2)-10*x+36;//return(y); void finding(float a[3],float f[3])float t=tt,a1,f1,ia;int i;f[0]=function(a[0]);for(i=0;;i++)a[1]=a[0]+t;f[1]=function(a[1]);if(f[1]<f[0]) break;if(fabs(f[1]-f[0])>=e)t=-t;a[0]=a[1];f[0]=f[1];else{if(ia==1) return;t=t/2;ia=1;for(i=0;;i++)a[2]=a[1]+t;f[2]=function(a[2]);if(f[2]>f[1]) break;t=2*t;a[0]=0;// 初始区间的下界值 求解的一维函数a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];if(a[0]>a[2])a1=a[0];f1=f[0];a[0]=a[2];f[0]=f[2];a[2]=a1;f[2]=f1;return;}float gold(float *ff)float a1[3],f1[3],a[4],f[4];float aa;int i;finding(a1,f1);a[0]=a1[0];f[0]=f1[0];a[3]=a1[2];f[3]=f1[2];a[1]=a[0]+0.382*(a[3]-a[0]);a[2]=a[0]+0.618*(a[3]-a[0]);f[1]=function(a[1]);f[2]=function(a[2]);for(i=0;;i++)if(f[1]>=f[2])a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];a[2]=a[0]+0.618*(a[3]-a[0]);f[2]=function(a[2]);else{a[3]=a[2];f[3]=f[2];a[2]=a[1];f[2]=f[1];a[1]=a[0]+0.382*(a[3]-a[0]);f[1]=function(a[1]);if((a[3]-a[0])<e)aa=(a[1]+a[2])/2;*ff=function(aa);break;return(aa);void main()float xx, ff;xx=gold(&ff);printf("\n The Optimal Design Result Is:\n"); printf("\n\tx*=%f\n\tf*=%f", xx, ff);getch();运行结果：2) mil SI*-5J?+4J?-fix+fiO 取舟=0 折二Ml A HT*程序如下：#in cludevstdio.h> #in clude<c oni o.h> #in clude<math.h> #defi ne e0.00001 #defi ne tt 0.01float fun cti on( float x)求解的一维函数float y=po w(x,4)-5* pow(x,3)+4* po w(x,2)-6*x+60;//return(y);void fin di ng(float a[3],float f[3])float t=tt,a1,f1,ia;int i;a[0]=0;// 初始区间的下界值f[0]=fu nctio n(a[0]);for(i=0;;i++)a[1]=a[0]+t;f[1]=fu nctio n(a[1]);if(f[1]<f[0]) break;{if(fabs(f[1]-f[0])>=e)t=-t;a[0]=a[1];f[0]=f[1];else{if(ia==1) return;t=t/2;ia=1;for(i=0;;i++)a[2]=a[1]+t;f[2]=function(a[2]);if(f[2]>f[1]) break;t=2*t;a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];if(a[0]>a[2])a1=a[0];f1=f[0];a[0]=a[2];f[0]=f[2];a[2]=a1;f[2]=f1;return;}float gold(float *ff)float a1[3],f1[3],a[4],f[4];float aa;int i;finding(a1,f1);a[0]=a1[0];f[0]=f1[0];a[3]=a1[2];f[3]=f1[2];a[1]=a[0]+0.382*(a[3]-a[0]);a[2]=a[0]+0.618*(a[3]-a[0]);f[1]=function(a[1]);f[2]=function(a[2]);for(i=0;;i++)if(f[1]>=f[2])a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];a[2]=a[0]+0.618*(a[3]-a[0]);f[2]=function(a[2]);else{a[3]=a[2];f[3]=f[2];a[2]=a[1];f[2]=f[1];a[1]=a[0]+0.382*(a[3]-a[0]);f[1]=function(a[1]);if((a[3]-a[0])<e)aa=(a[1]+a[2])/2;*ff=function(aa);break;return(aa);void main(){float xx, ff;{float t=tt,a1,f1,ia; xx=gold(&ff);prin tf("\n The Op timal Desig n Result ls:\n"); prin tf("\n\tx*=%f\n\tf*=%f", xx, ff); getch();}运行结果如下:3)(x+I)(x-2)^,其中讪，取坷 7, A(UH ,程序如下: #in clude<stdio.h>#in clude<c oni o.h>#in clude<math.h>#defi ne e 0.00001#defi ne tt 0.01float fun cti on( float x)float y=(x+1)* po w((x-2),2);//return(y); void fin di ng(float a[3],float f[3]) 求解的一维函数int i;{float t=tt,a1,f1,ia;a[0]=0;// 初始区间的下界值f[0]=function(a[0]);for(i=0;;i++)a[1]=a[0]+t;f[1]=function(a[1]); if(f[1]<f[0]) break;if(fabs(f[1]-f[0])>=e)t=-t;a[0]=a[1];f[0]=f[1];else{if(ia==1) return;t=t/2;ia=1;for(i=0;;i++)a[2]=a[1]+t;f[2]=function(a[2]);if(f[2]>f[1]) break;t=2*t;a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];if(a[0]>a[2])a1=a[0];f1=f[0];a[0]=a[2];f[0]=f[2];a[2]=a1;f[2]=f1;}return;}float gold(float *ff)float a1[3],f1[3],a[4],f[4];float aa;int i;finding(a1,f1);a[0]=a1[0];f[0]=f1[0];a[3]=a1[2];f[3]=f1[2];a[1]=a[0]+0.382*(a[3]-a[0]);a[2]=a[0]+0.618*(a[3]-a[0]);f[1]=function(a[1]);f[2]=function(a[2]);for(i=0;;i++)if(f[1]>=f[2])a[0]=a[1];f[0]=f[1];a[1]=a[2];f[1]=f[2];a[2]=a[0]+0.618*(a[3]-a[0]);f[2]=function(a[2]);else{a[3]=a[2];f[3]=f[2];a[2]=a[1];f[2]=f[1];a[1]=a[0]+0.382*(a[3]-a[0]);f[1]=function(a[1]);if((a[3]-a[0])<e)aa=(a[1]+a[2])/2;*ff=function(aa);break;return(aa);void mai n()float XX, ff;xx=gold(&ff);prin tf("\n The Op timal Desig n Result ls:\n");prin tf("\n\tx*=%f\n\tf*=%f", xx, ff);getch();}运行结果如下:2.复合形法程序测试题1) =齣W = 2p-环D取：1|6[-5 6) ^曰一5 8] * = 4 f = l『程序如下:{}#in clude "math.h"#i nclude "stdio.h"#in clude "stdlib.h"#define E1 0.001#define ep 0.00001#define n 2#define k 4double af;int i,j;double X0[n],XX[n],X[k][n],FF[k];double a[n],b[n];double rm=2657863.0;double F(double C[n])double F;F=pow(C[0]-2,2)+pow(C[1]-1,2);return F;int cons(double D[n])if((D[1]-pow(D[0],2)>=0)&&((2-D[0]-D[1])>=0)) return 1;elsereturn 0;void bou()a[0]=-5,b[0]=6; a[1]=-5,b[1]=8;{}double r()double r1,r2,r3,rr;r1=pow(2,35);r2=pow(2,36);r3=pow(2,37);rm=5*rm; if(rm>=r3){rm=rm-r3;}if(rm>=r2){rm=rm-r2;}if(rm>=r1){rm=rm-r1;}rr=rm/r1;return rr;void produce(double A[n],double B[n])int jj;double S;s1: for(i=0;i<n;i++)S=r();XX[i]=A[i]+S*(B[i]-A[i]);if(cons(XX)==0){goto s1;}for(i=0;i<n;i++)X[0][i]=XX[i];for(j=1;j<k;j++)for(i=0;i<n;i++)S=r();for(j=1;j<k;j++){X[j][i]=A[i]+S*(B[i]-A[i]);}for(i=0;i<n;i++)X0[i]=0;for(jj=1;jj<j+1;jj++)X0[i]+=X[jj][i];X0[i]=(1/j)*(X0[i]);if(cons(X0)==0)goto s1;for(i=0;i<n;i++){XX[i]=X[j][i];}while(cons(XX)==0)for(i=0;i<n;i++)X[j][i]=X0[i]+0.5*(X[j][i]-X0[i]); XX[i]=X[j][i];main()}}double EE,Xc[n],Xh[n],Xg[n],Xl[n],Xr[n],Xs[n],w; int l,lp,lp1;bou();s111:produce(a,b);s222:for(j=0;j<k;j++)for(i=0;i<n;i++)XX[i]=X[j][i];FF[j]=F(XX);for(l=0;l<k-1;l++)for(lp=0;lp<k-1;lp++)lp1=lp+1;if(FF[lp]<FF[lp1])w=FF[lp];FF[lp]=FF[lp1];FF[lp1]=w;for(i=0;i<n;i++)XX[i]=X[lp][i];X[lp][i]=X[lp1][i];X[lp1][i]=XX[i];for(i=0;i<n;i++)Xh[i]=X[0][i];Xg[i]=X[l][i];Xl[i]=X[k-1][i];for(i=0;i<n;i++)}{Xs[i]=0;for(j=0;j<k;j++)Xs[i]+=X[j][i];Xs[i]=1/(k+0.0)*Xs[i];EE=0;for(j=0;j<k;j++)EE+=pow((FF[j]-F(Xs)),2);EE=pow((1/(k+0.0)*EE),0.5); if(EE<=E1)goto s333;for(i=0;i<n;i++)Xc[i]=0;for(j=1;j<k;j++)Xc[i]+=X[j][i];Xc[i]=1/(k-1.0)*Xc[i];if(cons(Xc)==1)}}af=1.3;ss:for(i=0;i<n;i++)Xr[i]=Xc[i]+af*(Xc[i]-Xh[i]);if(cons(Xr)==1)if(F(Xr)>=F(Xh)) if(af<=ep)for(i=0;i<n;i++)Xh[i]=Xg[i];af=1.3;goto ss;else {af=1/2.0*af;goto ss;}elsefor(i=0;i<n;i++)X[0][i]=Xr[i];goto s222;else {af=1/2.0*af;goto ss;}}}Q —善— h 却王N —才二岀—£1H+J f D +&—*}%+叱唱心亠一----F一= D -瓦d-j各nqwo盂舄90010国0畧dopISQVEP匕旦WEP戈芒s m p H二二_X T =.」&一匸七(++g v =?D 」o 」兰X)ZI-=UV&H(U壬一」tlls0106宀兰_XH S q 兰o x 丄亘①S-①宀兰。