机械优化设计大作业2011 - 副本
《机械优化设计》大作业讲解学习
高等流体力学班级：机设15学硕班学号： 2015200813 姓名：张湘楠授课老师：毕新胜日期： 2016年7月 1日一、研究报告内容：1、λ＝0.618的证明、一维搜索程序作业；2、单位矩阵程序作业；3、连杆机构问题＋自行选择小型机械设计问题或其他工程优化问题；（1）分析优化对象，根据设计问题的要求，选择设计变量，确立约束条件，建立目标函数，建立优化设计的数学模型并编制问题程序；（2）选择适当的优化方法，简述方法原理，进行优化计算；（3）进行结果分析，并加以说明。
4、写出课程实践心得体会，附列程序文本。
5、为响应学校2014年度教学工作会议的改革要求，探索新的课程考核评价方法，特探索性设立一开放式考核项目，占总成绩的5%。
试用您自己认为合适的方式（书面）表达您在本门课程学习方面的努力、进步与收获。
（考评将重点关注您的独创性、简洁性与可验证性）。
二、研究报告要求1、报告命名规则：学号－姓名－《机械优化设计》课程实践报告.doc2、报告提交邮址：weirongw@（收到回复，可视为提交成功）。
追求：问题的工程性，格式的完美性，报告的完整性。
不追求：问题的复杂性，方法的惟一性。
评判准则：独一是好，先交为好；切勿拷贝。
目录：λ＝0.618的证明、一维搜索程序作业① 关于618.0=λ的证明……………………………………………………4 ② 一维搜索的作业采用matlab 进行编程…………………………………………… 5 采用C 语言进行编程……………………………………………… 7 单位矩阵程序作业① 采用matlab 的编程………………………………………………… 9 ② 采用c 语言进行编程………………………………………………… 9 机械优化工程实例① 连杆机构...........................................................................11 ② 自选机构...........................................................................16 课程实践心得.............................................................................. 20 附列程序文本.............................................................................. 21 进步，努力，建议 (25)一、λ＝0.618的证明、一维搜索程序作业①关于618.0=λ的证明黄金分割法要求插入点1α，2α的位置相对于区间],[b a 两端具有对称性，即)(1a b b --=λα)(2a b a -+=λα其中λ为待定常数。
机械优化设计范例(共9张PPT)
设计变量
现设 甲矿运往东站x万吨
乙矿运往东站y万吨
则甲矿运往西站200-x万吨
乙矿运往西站260-y万吨 令x=x1,y=x2
所以：X43;1.5(200-x1)+0.8x2+1.6(260-x2) =716-0.5x1-0.8x2（万元）
所以：Min f(X)= 716-0.5x1-0.8x2
约束条件
- x1 ≤0 X1-200 ≤0 -x2 ≤0 x2 - 260 ≤ 0
x1+x2-280≤ 0 100-x1-x2≤0
求解结果
x2 280 260
100
Z
(20,260)
x1=20 x2=260
Minf（X）= 498万元
100
200 280
x1
所以： 乙矿运往西站260-y万吨
Mx2in-f（26X0）≤ =0 498万元 则令甲x=矿x1运,y=往x2西站200-x万吨
最少的运费为498万元 x令1x+=xx21-2,y8=0x≤20
己 x1知+x甲2-、28乙0≤两0煤矿每年的产量分别为200万吨和260万吨，需经过东、西两个车站运外地。 M甲i煤nf（矿X运）往=东49站8和万西元车站的运费价格分别为1元／吨和1．5元／吨，乙煤矿运往东车站和西车站的运费价格分别为0．8元／吨和1．6元／吨 所。以：Min f(X)= 716-0. 煤乙矿应 运怎往样东编站制y万调吨运方案才能使总运费最少？ 己x1知+x甲2-、28乙0≤两0煤矿每年的产量分别为200万吨和260万吨，需经过东、西两个车站运外地。 xM1i+nfx（2-X2）80=≤ 4098万元 现甲设煤矿甲运矿往运东往站东和站西x万车吨站的运费价格分别为1元／吨和1．5元／吨，乙煤矿运往东车站和西车站的运费价格分别为0．8元／吨和1．6元／吨 所 。以：X = [ x1, x2 ]T
机械优化设计优化作业习题
1—61、⎪⎪⎩⎪⎪⎨⎧≥≤-≤+≤+-⋅--=0,31232424min 2121212121x x x x x x x x t s x x f 2、72220:m in 321321≤++≤⋅-=x x x t s x x x f3、机床主轴是机床中重要零件之一，一般为多支承空心阶梯轴。
为了便于使用材料力学公式进行结构分析，常将阶梯轴简化成以当量直径表示的等截面轴。
图1所示的为一根简化的机床主轴。
要求以主轴的自重为目标，对该主轴进行优化设计。
已知条件：主轴材料为45#，内径d=30mm ，外力F=15000N ，许用挠度y 0=0.05mm ，材料的弹性模量E=210GPa ，许用应力[σ]=180MPa 。
300≤ l ≤650， 60≤ D ≤110， 90≤ a ≤150。
l 、d 、a 的量纲均为毫米。
7-121、⎪⎩⎪⎨⎧≥=++≥++⋅++=0,20521532min 21321321321x x x x x x x x t s x x x f2、⎪⎪⎩⎪⎪⎨⎧≥≤+≤+-≤+⋅---+=0,32222625.0min 2121212121212221x x x x x x x x t s x x x x x x f 3、由两根实心圆杆组成对称的两杆桁架，其顶点承受负载为2p ，两支座之间的水平距离为2L ，杆的比重为ρ，弹性模量为E ，屈服强度为σ。
求在桁架不被破坏的情况下使桁架重量最轻的桁架高度h 及圆杆直径d 。
7-121、⎪⎪⎩⎪⎪⎨⎧≥=+-=++=++⋅++=0,102420521532min 21321321321321x x x x x x x x x x x t s x x x f 2、⎩⎨⎧-≥-≤--⋅++++=105.1)12424(min 21212122122211x x x x x x t s x x x x x e f x 3、已知：制造一体积为100m3，长度不小于5m ，不带上盖的箱盒，试确定箱盒的长x1，宽x2，高x3，使箱盒用料最省。
研究生《工程优化设计》大作业,包含小型机械的优化!绝对的原创!
工程优化设计计算机编程实验报告一 、优化设计方法基本原理优化设计就是：根据给定的设计要求和现有的技术条件，应用专业理论和优化方法，在电子计算机上从满足给定的设计要求的许多可行方案中，按照给定的目标自动地选出最优的设计方案。
优化过程是：寻找约束空间下给定函数取极大值（以max 表示)或极小(以min 表示)的过程，优化方法也称数学规划。
现代优化设计方法就是：在计算机上进行的半自动或自动优化，以选出在现有工程条件下的最佳设计方案的一种现代设计方法。
机械优化设计就是：把机械设计与优化设计理论及方法相结合，借助电子计算机，自动寻找实现预期目标的最优设计方案和最佳设计参数。
1.一维优化方法求一元函数 f(x)的极小点和极小值问题就是一维最优化问题。
求解一维优化问题的方法称为一维优化方法。
实际优化问题中一维问题是很少的，大多数问题都是多维的，但一维问题是优化中最简单、最基本的方法，它是解决多维问题的基础。
一维优化方法分为两类：一类是直接法：按某种规律取若干点计算其函数值，然后通过函数值的直接比较来最后确定最优解。
确定初始区间的进退法、黄金分割法等；一类是间接法：要利用函数的导数，故称解析法，牛顿法和二次插值法。
本次实验用到黄金分割法和牛顿法，故重点介绍这两种方法。
1.1黄金分割法黄金分割法是通过不断缩短单峰区间的长度来搜索极小点的一种有效方法，它是搜索区间按比例λ缩小，通过计算和比较()f x 的函数值，以确定取舍区间，因按黄金分割原理：0.618λ=，故此法又称为0.618法。
在搜索区间[a,b]中取两点c 、d ，然后比较c 、d 两点的函数()f c 和()f d ，有三种情况：（1）如()()f c f d >，根据函数的单峰性，极值点必在[,]c b 区间。
（2）如()()f c f d <，极值点必在[,]a d 区间。
（3）如()()f c f d =，极值点在[,]c d 区间。
《机械优化设计》大作业
1、报告命名规则：学号－姓名－《机械优化设计》课程实践报告.doc
2、报告提交邮址：weirongw@（收到回复，可视为提交成功）。
追 求：问题的工程性，格式的完美性，报告的完整性。
不 追 求：问题的复杂性，方法的惟一性。
评判准则：独一是好，先交为好；切勿拷贝。
printf(“区间为(%lf,%lf),精度为%lf\n”,a,b,size);
printf(“序号\t a1\t\t a2\t\t f(a1)\t\t f(a2)\n”);
for(i=0;i<64;i++)
printf(“-“);
printf(“\n”);
i=0;
while((b-a)>size)/*****用精度控制循环次数*****/
如果要计算y=cosx；黄金分割法，则需要将图所示的f脚本中的函数写成如②所示
fun函数表示对matlab的主程序语言。
函数NO.2运算结果：
函数NO.1运算结果：
（2）采用C语言进行编程
#include<stdio.h>
#defineK 0.618
doublef(double);/*****函数值计算函数声明*****/
c1=b-r*(b-a);c2=a+r*(b-a);
y1=f(c1);y2=f(c2);
while (abs((b-a)/b)>e)&&(abs((y2-y1)/y2)>e)
if y1>=y2
a=c1;c1=c2;y1=y2;
c2=a+r*(b-a);y2=f(c2);
else
b=c2;c2=c1;y2=y1;
《机械优化设计》大作业
for(i=0;i<64;i++)
printf(“-“);
printf(“\n”);
i=0;
while((b-a)>size)/*****用精度控制循环次数*****/
{
i++;
a1=b-K*(b-a);/*****按0.618法插入两点*****/
a2=a+K*(b-a);
printf(“%2d:\t%f\ta%f\t%f\t%f\n”,I,a1,a2,f(a1),f(a2));
/*****输出每次计算后a1,a2,f(a1),f(a2)的值*****/
if(f(a1)>=f(a2))
a=a1;
else
b=a2;
}
printf(“所求极小值点为:x=%lf\t极小值f(x)=%f\n”,a,f(b));
（2）选择适当的优化方法，简述方法原理，进行优化计算；
（3）进行结果分析，并加以说明。
4、写出课程实践心得体会，附列程序文本。
5、为响应学校2014年度教学工作会议的改革要求，探索新的课程考核评价方法，特探索性设立一开放式考核项目，占总成绩的5%。
试用您自己认为合适的方式（书面）表达您在本门课程学习方面的努力、进步与收获。（考评将重点关注您的独创性、简洁性与可验证性）。
end
end
x=0.5*(a+b);
end
对f函数的确立
10function y=f(x)
y=(x-2)^2+3;
end
11function y=f(x)
y=cosx
机械优化设计方案三个案例
机械优化设计案例11. 题目对一对单级圆柱齿轮减速器，以体积最小为目标进行优化设计。
2.已知条件已知数输入功p=58kw ，输入转速n 1=1000r/min ，齿数比u=5，齿轮的许用应力[δ]H =550Mpa ，许用弯曲应力[δ]F =400Mpa 。
3.建立优化模型3.1问题分析及设计变量的确定由已知条件得求在满足零件刚度和强度条件下，使减速器体积最小的各项设计参数。
由于齿轮和轴的尺寸（即壳体内的零件）是决定减速器体积的依据，故可按它们的体积之和最小的原则建立目标函数。
单机圆柱齿轮减速器的齿轮和轴的体积可近似的表示为：]3228)6.110(05.005.2)10(8.0[25.087)(25.0))((25.0)(25.0)(25.0222122212221222212212122221222120222222222121z z z z z z z z z z z g g z z d d l d d m u m z b bd m u m z b b d b u z m b d b z m d d d d l c d d D c b d d b d d b v +++---+---+-=++++-----+-=πππππππ式中符号意义由结构图给出，其计算公式为b c d m u m z d d d mu m z D m z d m z d z z g g 2.0)6.110(25.0,6.110,21022122211=--==-===由上式知，齿数比给定之后，体积取决于b 、z 1 、m 、l 、d z1 和d z2 六个参数，则设计变量可取为T z z T d d l m z b x x x x x x x ][][211654321==3.2目标函数为min)32286.18.092.0858575.4(785398.0)(2625262425246316321251261231232123221→++++-+-+-+=x x x x x x x x x x x x x x x x x x x x x x x x x x f3.3约束条件的建立1）为避免发生根切，应有min z z ≥17=，得017)(21≤-=x x g2 )齿宽应满足max min ϕϕ≤≤d b，min ϕ和max ϕ为齿宽系数d ϕ的最大值和最小值，一般取min ϕ=0.9，max ϕ=1.4，得04.1)()(0)(9.0)(32133212≤-=≤-=x x x x g x x x x g3）动力传递的齿轮模数应大于2mm ，得 02)(34≤-=x x g4）为了限制大齿轮的直径不至过大，小齿轮的直径不能大于max 1d ，得0300)(325≤-=x x x g 5）齿轮轴直径的范围：max min z z z d d d ≤≤得0200)(0130)(0150)(0100)(69685756≤-=≤-=≤-=≤-=x x g x x g x x g x x g 6）轴的支撑距离l 按结构关系，应满足条件：l 2min 5.02z d b +∆+≥（可取min ∆=20），得0405.0)(46110≤--+=x x x x g7）齿轮的接触应力和弯曲应力应不大于许用值，得400)10394.010177.02824.0(7098)(0400)10854.0106666.0169.0(7098)(0550)(1468250)(224222321132242223211213211≤-⨯-⨯+=≤-⨯-⨯+=≤-=---x x x x x x g x x x x x x g x x x x g8）齿轮轴的最大挠度max δ不大于许用值][δ，得0003.0)04.117)(445324414≤-=x x x x x x g 9）齿轮轴的弯曲应力w δ不大于许用值w ][δ，得5.5106)1085.2(1)(05.5104.2)1085.2(1)(1223246361612232463515≤-⨯+⨯=≤-⨯+⨯=x x x x x g x x x x x g4.优化方法的选择由于该问题有6个设计变量，16个约束条件的优化设计问题，采用传统的优化设计方法比较繁琐，比较复杂，所以选用Matlab 优化工具箱中的fmincon 函数来求解此非线性优化问题，避免了较为繁重的计算过程。
机械优化设计课后习题答案(可编辑修改word版)
3 1 2 32第一章习题答案1-1某厂每日（8h 制）产量不低于 1800 件。
计划聘请两种不同的检验员，一级检验员的标准为：速度为 25 件／h，正确率为 98％，计时工资为 4 元／h；二级检验员标准为：速度为 15 件／h，正确率为 95％，计时工资 3 元／h。
检验员每错检一件，工厂损失 2 元。
现有可供聘请检验人数为：一级 8 人和二级 10 人。
为使总检验费用最省，该厂应聘请一级、二级检验员各多少人？解：（1）确定设计变量；根据该优化问题给定的条件与要求，取设计变量为X = ⎡x1⎤=⎡一级检验员⎤；⎢x ⎥⎢⎥（2）建立数学模型的目标函数；取检验费用为目标函数，即：⎣2 ⎦ ⎣二级检验员⎦f(X) = 8*4*x1+ 8*3*x2+ 2（8*25*0.02x1+8*15*0.05x2）=40x1+ 36x2（3）本问题的最优化设计数学模型：min f(X) = 40x1+ 36x2X∈R3·s.t. g1(X) =1800-8*25x1+8*15x2≤0g2(X) =x1-8≤0 g3(X)=x2-10≤0g4(X) = -x1≤0g5(X) = -x2≤01-2已知一拉伸弹簧受拉力F，剪切弹性模量G，材料重度r，许用剪切应力[]，许用最大变形量[] 。
欲选择一组设计变量X = [x1x2x ]T= [d D n]T使弹簧重量最轻，同时满足下列限制条件：弹簧圈数n ≥ 3 ，簧丝直径d ≥ 0.5 ，弹簧中径10 ≤D2≤ 50 。
试建立该优化问题的数学模型。
注：弹簧的应力与变形计算公式如下8FD 1 D 8F D3=k 2,k=1+,c=2(旋绕比），=n2解：（1）确定设计变量；s d 3s2c d Gd 4⎡x1 ⎤⎡d ⎤根据该优化问题给定的条件与要求，取设计变量为X = ⎢x ⎥=⎢D ⎥；⎢ 2 ⎥⎢2 ⎥（2）建立数学模型的目标函数；取弹簧重量为目标函数，即：⎢⎣x3⎥⎦ ⎢⎣n ⎥⎦f(X) = 2rx 2x x1 2 34（3） 本问题的最优化设计数学模型：min f (X ) =2 rx 2 x x4X ∈R 3·1 1 8Fx 3x x 高h s.t. g 1(X ) =0.5-x 1 ≤0 g 2(X ) =10-x 2 ≤0 g 3(X )=x 2-50 ≤0g 4(X ) =3-x 3 ≤0g 5(X ) =(1+x 1 2x 2 ) 8Fx 2 - []≤0 x 3g 6(X ) = 2 3 - []≤0 Gx 41-3 某厂生产一个容积为 8000 cm 3 的平底、无盖的圆柱形容器，要求设计此容器消耗原材料最少，试写出这一优化问题的数学模型。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
宁波工程学院机械工程学院
机械优化设计大作业
班级
姓名
学号
教师
机械优化设计大作业
1.题目
行星减速器结构优化设计
NGW型行星减速器应用非常广泛。
1.1结构特点
(1)体积小、重量轻、结构紧凑、传递功率大、承载能力高；
(2)传动效率高，工作高；
(3)传动比大。
1.2用途和使用条件
某行星齿轮减速器主要用于石油钻采设备的减速，其高速轴转速为1300r/min；工作环境温度为-20℃～60℃，可正、反两向运转。
按该减速器最小体积准则，确定行星减速器的主要参数。
2.已知条件
传动比u=4.64，输入扭矩T=1175.4N.m，齿轮材料均选用38SiMnMo钢，表面淬火硬度HRC 45~55，行星轮个数为3。
要求传动比相对误差02
∆u。
.0
≤
弹性影响系数Z E=189.8MPa1/2；载荷系数k=1.05；
齿轮接触疲劳强度极限[σ]H=1250MPa；
齿轮弯曲疲劳强度极限[σ]F=1000MPa；
=2.97；应力校正系数Y Sa=1.52；
齿轮的齿形系数Y
Fa
小齿轮齿数z取值范围17--25；模数m取值范围2—6。
注：
优化目标为太阳轮齿数、齿宽和模数，初始点[24，52，5]T
3.数学模型的建立
建立数学模型见图1，即用数学语言来描述最优化问题，模型中的数学关系式反映了最优化问题所要达到的目标和各种约
束条件。
3.1设计变量的确定
影响行星齿轮减速器体积的独立参数为中心轮齿数、齿宽、模数及行星齿轮的个数,将他们列为设计变量，即：
x=[x
1 x
2
x
3
x
4
]T=[z
1
b m c]T [1]
式中：z
1¯￣
太阳轮齿数；b―齿宽（mm）；m—模数（mm）;行星轮的个数。
通常情况下，行星轮个数根据机构类型以事先选定，由已知条件c=3。
这样，设计变量为：
x=[x
1 x
2
x
3
]T=[z
1
b m]T [1]
3.2目标函数的确定
为了方便，行星齿轮减速器的重量可取太阳轮和3个行星轮体积之和来代替，即：
V=π/4（d
12+Cd
2
2）b
式中：d
1--太阳轮1的分度圆直径，mm；d
2
--行星轮2的分度圆
直径，mm。
将d
1=mz
1，
d
2
=mz
2
，z
2
=z
1
（u－2）/2代入（3）式整理，目标函数
则为：
F（x）=0.19635m2z
1
2b［4+（u－2）2c］[1]
式中u--减速器传动比;c--行星轮个数
由已知条件c=3,u=4.64,因此目标函数可简化为：
F（x）=4.891x
32x
1
2x
2
3.3约束条件的建立
3.3.1限制齿宽系数b/m的范围5≤b/m≤17，得：
g
（x）=5x3-x2≤0[1]
1
g
（x）=x2-17≤0[1]
2
不发生跟切，得：
3.3.2保证太阳轮z
1
g
（x）=17-x1≤0[1]
3
3.3.3限制齿宽最小值，得：
g
（x）=10-x2≤0[1]
4
3.3.4限制模数最小值，得：
g
（x）=2-x3≤0[1]
5
3.3.5按齿面接触疲劳强度条件，有：
g
（x）=750937.3/﹙x1x2x31/2)－[σ]H≤0[1]
6
式中：[σ]H--齿轮接触疲劳强度极限。
3.3.6按齿根弯曲疲劳强度条件，有：
g
（x）=1482000Y Fa Y Sa/﹙x1x2x32)－[σ]F≤0[1]
7
式中：[σ]F--齿轮弯曲疲劳强度极限；Y Fa--齿轮的齿形系数；Y Sa--应力校正系数。
4.优化方法的选择
解决优化问题的方法有很多中，该优化问题涉及到的数据计算比较庞大、处理比较麻烦，含有较多的非线性约束，所以我选用Matlab优化工具箱中的fmincon函数来求解该类非线性优化问题，避免了复杂且费时费力的计算过程。
5.数学模型的求解
5.1由已知数据可得该优化数学模型如下：
min F（x）=4.891x12x2 x32
Subject to:
g
（x）=5x3-x2≤0
1
g
（x）=x2-17x3≤0
2
g
（x）=17-x1≤0
3
g
（x）=10-x2≤0
4
g
（x）=2-x3≤0
5
g
（x）=750937.3/﹙x1x2x31/2)-1250≤0
6
g
（x）=6690340.8/﹙x1x2x32)-1000≤0
7
5.2基于Matlab对数学模型进行求解
在Matlab中首先编写目标函数的M文件myfun.m,返回x 处的函数值f：
function f = myfun(x)
f=4.891*x(1)^2*x(2)*x(3)^2;
由于约束条件中有非线性约束，所以需要编写一描述非线性约束条件的M文件mycon.m：
function[c,ceq]=confun(x)
c(1)=5*x(3)-x(2);
c(2)=x(2)-17*x(3);
c(3)=750937.3/(x(1)*x(2)*sqrt(x(3)))-1250;
c(4)=6690340.8/(x(1)*x(2)*x(3)^2)-1000;
ceq=[];
最后在command window里输入：
options=optimset('largescale','off'); %采用标准算法
x0 = [24; 52; 5]; %给定初值
lb=[17 ;10 ; 2];ub=[25 ;60; 6];%给定变量的上、下限约束
[x,fval,exitflag,output,lambda]=fmincon(@myfun,x0,[],[],[],[],lb,...
ub,@mycon,options) %调用优化过程
5.3最优解和结果分析
运行结果如下图所示：
由图可知，优化的结果为齿数z1=17.0000，齿宽b=27.8872mm，模数m=3.7566mm，总体积V min=5.5628×105mm3。
不过优化结果中的齿数z1必须为整数，齿宽b应圆整为27或者28；对于模数m，必须标准化为3.5或4。
经过计算比较取：z1=17，b=27mm，m=4mm为最优解。
6.结果对比分析
若按初始值减速器的体积V大约为7.3247×105mm3，而优化后的体积V则为6.1063×105mm3,优化结果比初始值体积减少为：ω＝1－(6.1063×105/7.3247×105)×100%＝16.6%
所以优化后的体积比未优化前减少了16.6%，说明优化是很成功的。
学习机械优化设计课程体会
《机械优化设计》是将机械工程设计问题转化为最优化问题，然后选择适当的最优化方法，利用电子计算机从满足要求的可行设计方案中自动寻找实现预期目标的最优化设计方案。
最为机械专业的一名学生，本课程，掌握最优化问题的基本解决方法，从多个可能的方案中选出最合适的、能实现预定最优目标的最优方案有着很现实的意义，为今后的工程实际提供了良好的理论储备。
而在机械优化设计基本理论学习的基础上，再使用美国Math Works开发的Matlab软件，及其附带的优化工具箱作为最优化问题的运算工具，通过对行星齿轮减速器的优化求解，可以看出运用Mtalab优化工具箱来求解优化问题，计算方便、快捷，高效的处理了涵盖各种难度的最优化问题，着实丰富了我的本课程的学习。
不过在初期学习使用的过程中还是出现了很多问题，比如非线性约束如何书写，上、下限如何选择，函数格式如何书写，变量未定义等等或大或小的问题，每次运行都会提示出现错误的地方，比如函数编写的时候没写对、少了括号、分号等，导致无法计算结果。
在不断排除错误、重新编写程序的过程中，渐渐的对Mtalab熟悉起来，出现的Error不断消失，最终在显示器上通过自己从无到有的知识积累才得到的优化结果，感觉一切的郁闷、烦恼都不翼而飞。
所以虽然这门课对数学水平有一定的要求，原理的推导复杂、诡异，各种各样的优化过程更是看的头晕目眩、似懂非懂，不过通过计算机程序还是能便捷的实现各种优化方法。
这门课程让我加深了对前人的敬佩，懂得了一些优化方法的简单计算计算过程和原理，不过有了软件一切计算都交给了计算机，让我省去了优化计算过程，也加深了学好英语的决心，因为一切先进的软件大都是英文的，看不太懂让我在这次设计过程中多走了很多弯路。
总之学习机械优化设计课程，深有体会的就是2个字“头疼”。
参考文献
1、《机械优化设计》[M]，孙靖民、梁迎春主编，北京，机械出版社，2006.12。
2、《基于Matlab的行星齿轮减速器的结构优化设计》[J]，周
小灵，湖南工业大学，《机电产品开发与创新》第22卷第2期,2009.3。
3、《基于Matlab的行星轮减速器优化设计》[J],刘雷敏、李友兴、杨小娟，长江大学，《机械工程师》2009年第9期，2009.7。
4、《机械设计》[M]，濮良贵、纪名刚主编，北京，高等教育出版社，2006.5。