现代设计黄金分割法复合形法实验报告word文档良心出品
- 格式:docx
- 大小:94.43 KB
- 文档页数:35
篇一: 3dmax 实训报告计算机 3dsmax 实习报告 室外三维建模所谓三维设计就是利用电脑进行设计与创作, 以产生真实的立体场景与动画。
虽然在这个学 期对3dsmax 的接触才知道3dsmax 是一个广泛应用于游戏开发、后期制造、影视特效及专业 视觉设计领域的一款功能强大的三维设计软件,它是集专业建模、动画、渲染一休的三维解 决方案。
特别是3damax 以其强大的功能卓越的表现力被广大电脑设计人所睛睐,成为当今较热门的设计软件。
在学习3dsmax 的这个学期里,3dsmax 创作的每一个细节都在刺激着我的神经,变幻无穷的 3d 建模,它的每一个创造都给我带来无比的震憾的惊喜与灵感 ?? 今天我要用自己所学到的知识来创造一个初级建模(室外建模)虽然这是我第一次曾未有的 偿试;不过相信抱着对 3d 的渴望,相信3dsmax 可以帮助我将难度复杂的室外模型简单地实 现出来。
下面是此次实习室外建模的几大重要路径:1 创建模型2 修改3 装饰4 环境渲染现在我来跟大家分享下这个初级建模的基本制造: (所谓初级建模: 是指利用几何体创造面板, 二维图形创建面扳中的现有模型来进行的建模操作,包括标准几何体的创建、扩展几何体创 建和二维图形的创建) 。
首先:一 创建模型1击创建命令面板中单击“图形”按钮,“创建”选择“线”按钮。
利用创建二维图形在顶视图描绘出建筑所在定的位置以及大小。
2 单击创建面板,单击“几何体”按钮, 进入“标准基本体”创建面板单击“长方体” 按钮,在“顶视图”刚创建好二维图形(平面图),照大小位置拉出几何体,其 它照样拉好之后,如右图:3 选择命令面板,单击“修改”按钮,打开修改面板设置长方体的参数,在这里主要设置高度即可, 为了方便修理在每一个几何体都给予一个名称。
因此先来设置房 1 的高度参数,(自 己认为合适即可) 。
4 单击“选择按钮” ,选择顶视图的的房 2 几何体; 照前步骤一样来设置几何体的高度。
机械优化设计报告姓名:刘洋学号: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、实验结果(总结/方案)黄金分割法求解极值实例。
研究报告黄金分割在生活中的应用东北育才学校马艺宸一.黄金分割的定义之比等于整体与较大部分之比,其比值为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. 网页设计在网页设计中,黄金分割法也是一种比较常用的设计原则。
例如,网页布局、按钮大小、文本位置等都可以按照黄金分割设计,这样可以让网页看起来更加优美和协调。
总结黄金分割法是一种非常美学的比例关系,被广泛应用于各个领域。
黄金分割法比例的应用可以让设计更加美观和协调,符合人们的审美标准,从而产生更好的视觉效果和感官体验。
现代设计方法实验报告篇一:现代设计方法实验报告《现代机械设计方法学》实验报告班级: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.文献资料法:欧几里得的<<几何原本>>.<<算盘书>>等等调查结果:经过调查发现以下几点:1.早在公元前6世纪古希腊的毕达哥拉斯学派就研究过正五边形和正十边形的作图,说明那时他们已经触及甚至掌握了黄金分割.公元前4世纪,古希腊数学家欧多克索斯第一个系统研究了这一问题,并建立起比例理论., 公元前300年前后欧几里得撰写<<帕乔利>>时吸收了欧多克索斯的研究成果,进一步系统论述了黄金分割,成为最早的有关黄金分割的理论的论著.到19世纪黄金分割这一名称才逐渐通行.黄金分割有许多有趣的性质,人类对它的实际应用也很广泛.最著名的例子是优选学中的黄金分割法或0.618法,是由美国数学家基弗于1953年首先提出的,70年代在中国推广.2.这个数学在自然界中和人们生活中到处可见:人们的肚脐是人体总长的黄金分割点,人的膝盖是肚脐到脚的黄金分割点.大多数门窗的宽长之比也是0.618-.建筑师们对数学0.618特别偏爱,无论是古埃及的金字塔,还是巴黎的圣母院,都与0.618有关3.0.618也广泛应用于战争中,在战略战役中一个极为迷人而神秘的数字,而且它还有着一个很动听的名字----黄金分割率,它是古希腊著名哲学家`数学家毕达哥拉斯于2500多年前发现的.古往今来,这数字一直被后人奉为科学与美学的金科玉律/在艺术史上,几乎所有的杰出作品都不谋而合地验证这一著名的黄金分割率. 在武器装备上,我们也能很容易的发现黄金分割率无处不在.在大炮射击中,如果某种间瞄火炮的最大射程为12公里,最小射程为4公里,则其最佳射击距离在9公里左右,为最大射程的2/3,与0.618十分接近.在防御战斗中,第一道防线的兵力通常为总数的2/3,第二道防线的兵力通常为总数的1/3.4.黄金分割与人的关系相当密切.近年来,在研究黄金分割与人体的关系时,发现了人体结构中有14个”黄金点”.12个”黄金矩形”和两个”黄金指数”. 黄金指数(1)反映鼻口关系的鼻唇指数(2)反映眼口关系的目唇指数0.618.作为一个人体健美的标准尺度之一,是无可非议的,但不能忽略其存在着”模糊特性”,它与其它美学参数一样,都有一个允许变化的幅度,受种族`地区`个体差异的制约.5.医学与0.618有着千丝万缕的联系,它可解释人为什么在环境22至24摄氏度时感觉最舒适.因为人的体温37与0.618的乘积为22.8.而且这一温度中肌体的新陈代谢`生理节奏和生理功能均处于最佳状态.科学家们还发现,当外界环境温度为人体温度的0.618时,人感到最舒服,现代医学研究还表明,0.618与养生之道息息相关.6.在大自然中,植物叶子,千姿百态,生机盎然,给大自然带来了绿色世界.尽管叶子形态随种而不同,但细心观察还是会有发现的.有些植物的花瓣及主干上枝条的生长,也符合这种规律.你从植物茎的顶端向下看,发现上下层中相邻的两片叶子之间约成137.5度.可计算得到360-137.5=222.5 137.5/222.5也约为.618结果: 通过研究探索发现黄金分割不仅应用于生活实际,还应用于数学.艺术和美术中,黄金数还存在于大自然中.从这次研究中,我们学会了很多,也懂了许多课本中没有的知识,明白了数学美是不同于其它的美,这种美是独特的`内在的,它具有严格的比例美`艺术美`和谐美.通过这次研究,最大的收获不仅是了解黄金分割点,重要的是学会了一种审美的角度,一种审美的观点,这种审美观源于大千世界中,源于事物中的存在的黄金分割比.”美是到处都有的,不是缺少美,而是缺少发现”.如果我们能积极地去寻找,不论是什么难题都可以克服,只要有恒心,就能完成.数学,如果正确看待它,不但拥有真理,而且也具有至高无上的美,是一种冷而严肃的美.这种美不是投合我们天性微弱的方面,这种美没有绘画或音乐那样华丽的服饰,它可以纯净到崇高的地步,能够达到严格的只有伟大的艺术能显示的那种完美的境界.数学,对我来说,是那样富有魅力,在生活中只要我们善于观察,善于思考,将所学的知识与生活结合起来将回感到生活的乐趣.生活中处处都应用着数学的知识.就像黄金数一样.。
研究性学习设计方案模板研究课题名称:黄金分割在生活中广泛应用设计者姓名所在学校所教年级七年级研究学科数学联系电话电子邮件一、课题背景、意义及介绍1、背景说明(怎么会想到本课题的):生活中并不缺少美,只是缺少发现。
黄金分割正是人们从生活中发现的美。
黄金分割是一种数学比例关系。
由公元前六世纪古希腊数学家毕达哥拉斯发现,有严格的比例性、艺术性、和谐性,蕴含着丰富的美学价值。
应用时一般取1.618,就像圆周率在应用时取3.14一般。
这个神奇的比例关系被证实于很多学科领域和日常生活的各个方面。
2、课题的意义(为什么要进行本课题的研究):21世纪的数学教学的理念是“人人学有价值的数学,人人都获得必需的数学,不同的人在数学上得到不同的发展”。
而课程标准中也指出:数学学习应该从学生的生活经验和已有知识背景出发,让他们在自主探索和合作交流中真正理解和掌握基本的数学知识。
黄金分割是数学的经典之一,为了让学生明白黄金分割的真面目,了解它在生活实际中的应用,也为了让学生能更深入的了解美丽和谐的概念,让学生激起对数学的兴趣,故设计了这个课题。
3、课题介绍本课题重点解决以下问题:1、黄金分割率由来。
2、黄金分割率的特点。
3、黄金分割率与美感。
4、黄金分割在生活中的应用。
5、正确认识黄金分割率。
采用查阅资料文献、网络搜索相关资料、实际测量法,根据所收集资料和调查的结果进行分析等方式,让学生体会到了调查研究的重要性。
在研究的过程中,学生亲身体现收集资料的成功与失败,获得了在课堂上从没有过的情感体现和社会经历,学会了组员间的相处和互助,培养了团队精神。
同时激发了学生学习数学的热情,又开拓了视野,增长了才智,这些都将成为学生成长过程中的宝贵财富,必将终身受益。
二、研究性学习的教学目的和方法(可按新课程标准的三维目标(或布鲁姆目标分类法)进行研究性学习的教学目的和方法的阐述)1、知识与技能(1)了解黄金分割、黄金分割点、黄金分割数的概念;(2)体验黄金分割在生活中的广泛应用。
《现代设计理论与方法》实验报告、实验目的机械优化设计是一门实践性较强的课程,学生通过实际上机计算可以达到以下目的: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-①宀兰。