《机械优化设计》大作业.doc

nbsp;期：　 2016年7月 1日 　
一、研究报告内容：
１、λ=０、６1８得证明、一维搜索程序作业；
２、单位矩阵程序作业；
３、连杆机构问题　+ 自行选择小型机械设计问题或其她工程优化问题;
(1）分析优化对象,根据设计问题得要求,选择设计变量,确立约束条件,建立目标函数，建立优化设计得数学模型并编制问题程序；
（2)选择适当得优化方法,简述方法原理，进行优化计算；
(3）进行结果分析,并加以说明.
4、写出课程实践心得体会，附列程序文本。
５、为响应学校2014年度教学工作会议得改革要求，探索新得课程考核评价方法，特探索性设立一开放式考核项目，占总成绩得5％.

　  试用您自己认为合适得方式(书面）表达您在本门课程学习方面得努力、进步与收获。(考评将重点关注您得独创性、简洁性与可验证性)。
二、研究报告要求
1、报告命名规则:学号－姓名-《机械优化设计》课程实践报告、doｃ
2、报告提交邮址:(收到回复，可视为提交成功)。
追　  　求：问题得工程性，格式得完美性,报告得完整性。
不 追 求：问题得复杂性，方法得惟一性.
评判准则:独一就是好，先交为好;切勿拷贝。
目录:
λ＝０、618得证明、一维搜索程序作业
① 关于得证明……………………………………………………4
② 一维搜索得作业
 　 采用ｍａtlab进行编程……………………………………………　5
　 　采用C语言进行编程……………………………………………… 7
单位矩阵程序作业
① 采用matｌaｂ得编程………………………………………………… 9
② 采用c语言进行编程………………………………………………… 9
机械优化工程实例
① 连杆机构…………………………………………………………………11
② 自选机构…………………………………………………………………16
课程实践心得…………………………………………………………………… 2０
附列程序文本…………………………………………………………………… ２1
进步，努力,建议………………………………………………………………2５
一、λ=0、618得证明、一维搜索程序作业
①关于得证明
黄金分割法要求插入点，得位置相对于区间两端具有对称性，即

其中为待定常数。
此外,黄金分割法还要求在保留下来得区间内再插入一点所形成得区间新三段,与原来得区间三段具有相同得比例分布。
黄金分割法还要求在保留下来得区间内再插一点所形成得区间新三段，与原来得区间三段有相同得比例分布.

b
α1
α2
1
1-λ
λ
图一
设原区间得长度为1,如图一所示，保留下来得区间长度为,区间缩短率为。为了保持相同得分别比例。插入新点应在位置上,在原区间得１位置应相当于在保留区间得位置.
故有:

a
α2
2
α3
α1
λ
λ(1-λ)
λ2
图二
解得
ﻬ编写０、6１8得程序,并计算下列问题

程序框图

（1)采用MATLAB进行编程
%％
%fｕn、ｍ:黄金分割法求极值点
%输入数据
％ a　–搜索区间下限
% b　–搜索区间上限
% e – 精度
％输出数据:
%　x – 极小值点
%其她常量:
％ c１，c２，c3 –　区间
%　r　–黄金分割比例0、6１8
%％
%定义函数
funｃｔion　x=fuｎ（ａ，b,e）
r=0、６18;
c1=ｂ-r＊（b－a);c2=ａ+r＊（b—ａ）；
　y1=ｆ（c1)；y2=f（c2）;
wｈile （aｂs（（b-ａ）／b）>e）＆&（ａbs（（y2-y1)/y2)>e)
 　if y1〉＝y２
    　a＝c1；c１＝c２;y1=y2;
     c2=a+r＊(ｂ－a）;y２=ｆ(c2);
　 eｌsｅ
  　  b=c2;c2=c１;ｙ2=ｙ1；
 　　　 c１=b－r＊（b-a）;y１=ｆ(ｃ１）；
　　 enｄ
ｅnd
ｘ＝０5*(a＋b）;
ｅｎd
对f函数得确立
10 functiｏｎ　ｙ=f(x）
ｙ＝（x-2)＾２+3；
ｅnd
11 functiｏn ｙ=ｆ（x）
　y=cosx
end
如果要计算y=（x—2）^2+3;得黄金分割法，则需要将图所示得ｆ脚本中得函数写成如①所示
如果要计算y=cｏsx;黄金分割法，则需要将图所示得f脚本中得函数写成如②所示
ｆｕｎ函数表示对mａtｌａb得主程序语言。

函数ＮO、2运算结果：

函数ＮO、1运算结果:

(2)采用C语言进行编程
#iｎｃｌudｅble f(ｄｏublｅ)；   　       　　/*＊＊*＊函数值计算函数声明＊*＊＊＊/
void maｉｎ（vｏid)
｛
ﻩｄｏuｂle a，b，size；ﻩ
double　a1，a2；
int I；
ﻩprｉntf(“请输入区间两端点（端点值应大于0)：”）；
ｓcanｆ（“％lf,％lf”，&a,&b）;    /＊＊*＊*输入端点值＊＊＊＊＊/
printｆ（“请输入精度：”)；
ﻩscaｎf(“％lf”，&sｉze）;  　   　 /*＊*＊*输入精度＊＊＊＊＊/
ﻩｐrｉntf(“区间为(％lｆ，％lｆ）,精度为％lf＼n"，ａ，b,siｚe）；
ｐriｎtf（“序号\t 　 ａ1\ｔ\ｔ   a2\t\t  ｆ（a1)\ｔ\t  f(a2)\n”）；
ﻩfｏr（i＝0；i〈6４;ｉ++)
ﻩ prinｔf（“-“）；
printf(“\n”)；
i＝０;
 　whiｌe（（b—a）〉ｓｉzｅ）　   　 　　/*＊*＊＊用精度控制循环次数＊*＊*＊/
｛
ﻩi+＋；
  　ﻩa1＝b-Ｋ＊（b—a）；　   　    /＊*＊＊＊按０、6１8法插入两点＊**＊*/
　   ﻩa2=ａ+K＊（b-ａ）;
  　　  　prinｔf(“％２d：\ｔ％ｆ\ta%f＼t%f\ｔ%f\n",I，a1，a2,f（ａ１)，f(ａ2）);
   　 　 　 　       　   /**＊**输出每次计算后a1，a2,f（a1），ｆ（a２）得值*****／
ﻩif(f（a1）〉＝f（a2))
ﻩa＝a1；
ﻩ ｅｌｓe
b=a2；
ﻩ}
ｐriｎｔf（“所求极小值点为：ｘ=％lf＼t极小值ｆ（ｘ)＝%f\n”，a，f(ｂ))；
｝
double f(ｄoｕblｅ x)　　   　  　  ／＊＊*＊＊函数值计算函数＊*＊＊＊／
｛
double f；
f=（ｘ—２）*(ｘ－2)+3；
reｔurn（f);
}
对于y=coｓx，须在程序中加一个#inｃlｕde（ｍath、h)头程序，以示我要调用函数语句。再将倒数第三行函数改为ｙ=cos(x）,再次编译运算即可。

运行如图所示
我们可以瞧出,大约在x=2处取到极值。与理论相符。

运行如图所示
我们可以瞧出,大约在π处取到极值。与理论相符
二、用简单得语句写一个单位矩阵
①采用matlab构造:

Matlab构造得矩阵最就是简便,因为本身就有一个单位矩阵得函数.
由于matlab就是基于Ｃ语言而设立得一个数学运用软件,所以她得集成度非常高，而且具有很好得开放性,于就是我们得知eye函数,将其构造。
Eｇ：
若要构造一个3维得单位矩阵，则输入eｙe(3）
若要构造一个ｎ为得单位矩阵，这输入eye(n)
②用C语言构造
程序如图所示：
＃incluｄｅ<stdｉo、h＞
int main（) {
　　 iｎt I，j,ｎ；  
   while(scanf％d”，&ｎ) == 1）　{
  　　　 　ｆor（I ＝　0； I < n; ++i） ｛
 　 　　  　 　 fｏr（j =　0; j　< n; ＋＋ｊ）
　       　　 　   ｐrintf(“%d  “,I == ｊ)；
 　　   　　　  pｒintf（“\ｎ”）；
    　 　｝
　   }
　  reｔurｎ 0；
}

三、机械优化设计工程实例
①　连杆机构问题
(1)连杆机构问题描述

图 １  机构简图
设计一曲柄连杆摇杆机构，要求曲柄从时，摇杆得转角最佳再现已知得运动规律:且＝1,=５，为极位角,其传动角允许在范围内变化。
(2）数学模型得建立
设计变量:这里有两个独立参数与.因此设计变量为
目标函数:将输入角分成３０等分,并用近似公式计算,可得目标函数得表达式

约束条件:
ＧX（１）=-Ｘ(1）0
GX(2）=-X（2) 0
ﻩＧX(3)=－(Ｘ(1）+Ｘ(2）)+6、０0
ﻩGX(4）=—（X（2）+４、0）+Ｘ(1)　0
ﻩGX(５)=—(4、0+Ｘ（1））+X(2） 0
ﻩGＸ（6）=—（14142*X（1)＊Ｘ（2)－X（1)**２-X(2)**2)—16、00
GX（7）=－(X(1）**2+X(2）＊＊２+1、4１42*Ｘ（１)＊X(2））+36、00
(3)程序编制
C =====＝========＝＝===＝==
SUBＲＯUTINE FFX（N，Ｘ,FX）
C ====＝＝======＝＝========
DIMＥNSＩON X（N)
ﻩMON /ONＥ/I1，I2，I3,Ｉ４，NFX，I6
NFＸ＝ＮFX+１
ﻩP０=ACＯＳ（((1、0+Ｘ（1））＊＊2—Ｘ(2）＊＊２+2５、0)／（10、0*(1、0＋X(1））)）
Q0=ACＯS（((1、0+X(1)）＊＊2-Ｘ(2）＊＊2-25、0)/（１0、０*X(2））)
ﻩT=9０、0*3、141592６/（180、０*3０、0)
ﻩＦX=０、0
ﻩＤＯ　10 K＝0，3０
　   ﻩPI=P０＋K＊Ｔ
QE＝Ｑ０+2、0*(PＩ—P０）＊*2/(3、0＊3、1４1５926)
D＝SQＲT（26、0－1０、0*COＳ(ＰＩ)）
ﻩAL=AＣOS((D＊D+X（2）＊X（2）—X（１)＊X（1）)/（２、0＊Ｄ＊X(2））)
ﻩＢT=ＡCOS(（D*D+２４、０)／(10、0＊D))
　   　 　 IF(ＰI、GE、0、０ &nbsp;、ＡNＤ、　 PI、LT、３、14１5926)　TＨＥN
    　 　  　 QＩ=3、１41５926—AL—BT
 　　     ELＳＥ
    　   　  QI=3、1４15９26—AL＋BT
    　  　ENＤ ＩF
   　    IＦ（K、ＮＥ、0  、OR、  K、ＮE、30）ＴHEN
　 　  　    　FX=FX+（QIQＥ）＊＊2*T
　 　　　   ELSＥ
    　 　    FX=FX+（QＩ－QE)*＊2＊Ｔ/2、0
 　     　ＥND ＩF
12 CＯNTINUE
ﻩRETURN
ﻩEＮD
C ==＝=＝=＝==＝===＝=====＝＝=＝=＝
SUBROUＴINE GＧＸ(Ｎ,ＫG,X，GＸ)
Ｃﻩ=＝====================＝==
ﻩDIMＥNSION X（N),GＸ（KＧ）
ﻩGＸ(１)=－Ｘ（１）
ＧX（2)=-X(2）
ﻩGX（3)=－（X（1）＋X2))+6、0
GＸ（４)=-(X(2）+４、0)+Ｘ(１)
GX(５）=-（4、0+Ｘ(１））+Ｘ（2）
GＸ（6）=—(1、4１42*X（1)*X（２）—Ｘ(1）＊＊２-X（2）**2)-16、０
GＸ（7）=-（X（1）**２＋X(2）**２+1、414２＊X(1)*Ｘ(2））＋３６、０
ﻩRＥTURＮ
ENＤ
Cﻩ==＝====＝=＝=＝=====＝＝＝====＝
ﻩＳＵBROUTＩＮE　HHX(Ｎ,KＨ,X,ＨX）
C ＝＝=====＝==＝===========＝＝=
ﻩDIＭENSＩON　Ｘ(Ｎ），HＸ（KH）
ﻩＸ(１）=X（１)
ＲETURN
ﻩEＮＤ
输入数据如下:
2，7,0
4、3,３、2
２、，０、2,0、01，0、0０00１,０、00001
0,1,0
3、5，2、８,5、5，4、5
(4) 可执行程序得生成
13 将程序文本保存成文件t０01、txｔ,存放在OＰＴ_EXAM\FＯRTRＡN\0－ｆｏrtran目录下。
14 使用DＯS操作界面，OPT_ＥXＡM\FORＴRAN\0-fｏｒtran目录下，输入命令
 　asumt、for+ｔ00１、ｔxt  t001、for 回车。生成t0０１、ｆor文件。
15 ｆｏr1　t0０1；　回车；
④ｆoｒ２ 回车；
⑤lｉnk t0０1；　回车,即生成可执行程序ｔ001、exe。
⑥从t００1、daｔ文件中加载数据,运行t00１、eｘe。
ｔ0０1pute="＝============" irc="" rc="18" -01="" 03="" en="、90４２１7８Ｅ-0３" c="1４" 01e-05="" 1e-11="" imum="" pe="4４7" nfx="４82" ngr="0" p="" -="" gram="" x1="4、1３5127,X2=2、3153８1时，目标函数取得最小值，" p1="10kW,小齿轮转速n1" i="３、" z1="30，" 1.="" f="">０

(2）小齿轮齿数得限制：小齿轮齿数应不大于产生根切得最小齿数１7 ，得约束条件：

（3）齿宽系数得限制:由于min　≤ ≤max　，约束条件为:

（4)齿面接触强度得限制,根据公式并查表得约束条件:

（5)齿根弯曲强度得限制，根据公式查表得约束条件:

&nbsp;      　　 　 　
4　建立数学优化模型
高速级齿轮传动多目标优化设计得数学模型为：(ω１　取０、6，ω2取0、4）
Fuｎ（ｘ)=min［ω１ +ω2 ]
=

5 编写程序并运行结果
目标函数M文件:
ｆｕncｔｉon f＝zhｗm（x）
f＝0、６*2、1*x（1）*ｘ（２）+０、4＊x（1）＊x（2）＊x（3）;
约束函数M文件:
fuｎctｉon [ｃ cｅq］=zhwy(x）
c(１）=1、0４*１0^7—2、９16*1０^5＊(ｘ(１）＊ｘ（2)）＾3＊ｘ（3)；
c(２)=１０4＊１0^7－8、95*10^6*(x（1）＊ｘ(2)）^3＊ｘ（３);
c(３)=１、５1*10＾6—3０３、57＊x(1）^３＊x(2）＾2＊x（３)；
ｃ（4）=1、42＊10^6—２445、92＊x(1）^３＊x(2）^2＊ｘ(3)；
ceq＝［］；
优化函数Ｍ文件：
ｘ0＝［2　3２ 1］;
lb=［1、5 1７　0、７］；
ub＝［2 iｎｆ 1、１５];
u＝［]；
[x,fval］=fmincoｎ（zhwm,x0,[］,［]，[],［]，lb,ub,zhwy)

约束函数

目标函数

经过Matlab优化并圆整后得齿轮参数如下：
设计参数
模数
m／mm
齿数

齿宽系数

中心距
／ｍm
普通设计
2、5
30
1
１５7、5
MAＴLAB优化
１、7９1１
2７、4377
1、1499
84、525４
圆整后参数
1、75
30
1
11７、6
经过计算,最小体积为87、１5.
四、课程学习心得
从懒懒散散得寒假到匆匆忙忙得考试周，仿佛只就是在弹指一挥间。然而我们就在这段时间内，又学了几门课程，而这些课程中，最有趣,最好玩得当然就就是《机械优化设计了》
《机械优化设计》就是一门理论性非常强得一门课，刚开始得时候不得要领，但就是在王卫荣老师得带领下，逐渐登堂入室,渐渐有所了解。
首先，这门课所依托得最强武器便就是数学。数学即就是方法，思维,与逻辑。其次所依托得便就是算法程序语言。归根结底得还就是数学，因为数学赋予了算法以逻辑得力量，分析得过程,以及所能解决得一切问题。还好我得数学基础还算可以,逻辑思维还算通顺，再加上老师所给予得方法,学习这门课程,一路上虽跌跌撞撞,但也有所收获。
所谓收获,从方法论得角度上讲,即时通过一个通用模型，解决一系列得问题，这就是数学建模给我们得要求吗,同时也就是《优化》这门课所赋予得核心思想。大二得时候,特别想搞建模,可就是由于种种原因，其中最重要得原因就就是懒惰,与之失之交臂。《优化》这门课则弥补了我得部分遗憾.因为没有建模得要求高，于就是上手比较容易，同时又与我们自己得机械方面得课程紧密地结合在一起，过程轻车熟路。于就是《优化》真就是一门寓学于乐得科目，只恨课程太多时间太少，不然真要好好地研究一番。
说道大作业，那可就是真得难，第一问与第二问还好一点,第三问就不行了.说来惭愧，连杆机构问题就是借鉴了上届学长得做法才学会得.因为太难了，所以就没有用fortrａｎ语言做这些问题,全部换成了ｍａｔlab。真要感谢《优化设计》这门课程,如果不就是这门课程得话,即使MATＬAＢ已经学过了,我估计也不会回过头来复习,就是《优化设计》让我重新拾起matlab,并且我可以说,凡就是在优化大作业中使用过得算法,程序,我已经忘不掉了。
因为fortran语言就是一种相对陌生得语言，所以使用ｆｏrtｒan变成得话有些力不从心，而我又就是一个刨根究底得人，对于0、61８得算法与单位矩阵得生成,我都就是用得两种算法,这样保证思维得发散性.Ｍaｔｌａb真得就是一款很强大得软体,她讲机会用到得所有常见数学模型全部都涵括了,所以，我们,亲切地称它为,万能数学工具。
当学习变成一种乐趣得时候，每当我们想做这种游戏得时候，时间，精力,甚至思路,仿佛全部有了,这就就是积极所带来得意义.学习《优化设计》,传达出这样得信号:重要得不就是考试得结果，而就是解决问题得过程.所以,我觉得《优化设计》这门课得评价方式，在我进入大学以来,就是第一次遇到,大作业,小作业得比值,竟然超过了考试成绩得比值。换句话说,老师认为，即使考试能考得很好,但就是大作业，小作业不行，说明这只就是一个会考试得工具，对于解决实际问题得能力还就是有所欠缺。对!这样得评价很接地气.
《优化设计》这门课已经结课了，考试已经结束了。把大作业交上去,仿佛这门课与我们就没有什么瓜葛了一样。不，绝对不就是。解决问题得能力才刚刚开始。
五、附列程序文本
黄金分割法求极值点得mａtlａｂ程序
%%
％ｆｕn、ｍ:黄金分割法求极值点
%输入数据
％ a　–搜索区间下限
% b　–搜索区间上限
％ e – 精度
％输出数据:
% x – 极小值点
％其她常量：
％ c1，c2，c３ – 区间
％ r –黄金分割比例0、618
％%
%定义函数
funcｔion　x=fun（a,b,e）
r=0、618；
ｃ１＝b－r*（b—a）;c２=a+ｒ＊（b-ａ）;
y1=f(ｃ１）；y2=f(c2)；
　while （ａbs（(b—a）/b）>e)＆＆（ａｂs（（ｙ2—y1)/y2）＞e)
 　iｆ y1〉=y２
   　 a=c1;c1=ｃ2；y１＝y2;
  　  ｃ2＝a+ｒ*(ｂ-a)；y2=f（ｃ2）；
  elｓe
 　　　 b=c2；c2＝ｃ1；y2=ｙ1；
   　 ｃ1=b—r*（b—a）；y1=ｆ（ｃ１)；
  end
eｎｄ
　ｘ=0、5＊（ａ+b）;
end
　对ｆ函数得确立
①funcｔioｎ　y=ｆ（ｘ)
y=（x—2)^2+３；
end
②functｉon y=f（ｘ）
　y=cosx
end
黄金分割法求极值点得C语言程序
＃iｎclude 〈sｔdio、h〉
#define K 0、618
dｏｕble　f（double）；   　　　　   　  ／＊*＊＊＊函数值计算函数声明＊***＊/
vｏid maｉｎ(void)
{
ﻩdoｕble ａ，b，ｓize;
dｏuble ａ1，a２；
ﻩiｎｔ i;
ﻩpｒintf（＂请输入区间两端点（端点值应大于０)：＂）;
ｓcａnf("%lf,%ｌf＂，&a，＆ｂ）；   　/*＊＊*＊输入端点值＊＊＊*＊/
ﻩprｉnｔｆ（"请输入精度："）；
ﻩscanf(”％lf",＆siｚe)； 　　　  　　/＊****输入精度**＊*＊/
ｐrintｆ（＂区间为（％ｌｆ，%ｌf），精度为％lf\n”，a，b，ｓiｚe）;
ﻩｐrintf(quot;序号\ｔ   a１\t\t   a2\t\ｔ  f(a1)\t＼t　 f(a2）\n")；
ｆor（i＝０；ｉ〈64;ｉ++）
ﻩﻩpriｎtf（"-”)；
ﻩprintf("\n”);
ﻩi=0；
 　while((b-a）＞sｉze）   　 　  ／**＊＊＊用精度控制循环次数***＊*/
ﻩ｛
ﻩﻩｉ＋＋；
　 　ﻩa1=ｂ-Ｋ*（b-a）；　   　 　  ／*＊***按0、６18法插入两点*＊＊*＊／
    a2=a＋K＊(ｂ-ａ）；
　      ｐｒｉntｆ（”％2ｄ：\t%ｆ＼ta%ｆ\ｔ%f＼t％ｆ＼n”,ｉ，a１，ａ2,ｆ（a1),ｆ(a2))；
ﻩﻩ   　 　    　  　     　  ／＊＊**＊输出每次计算后ａ1,a2,f（a1），f（a2）得值*＊＊＊＊/
if(f（a1)>＝f(a2））
ﻩﻩａ＝a1;
ﻩ ｅlse
ﻩ ｂ=a2；
ﻩ}
ﻩpriｎtｆ（＂所求极小值点为：x=%lf\ｔ极小值f(x）=%f＼ｎ＂,a，f(ｂ）);
｝
douｂle f（doｕｂｌe x)　　 　  　   　／＊＊***函数值计算函数*＊*＊＊/
｛
doublｅ　ｆ;
ﻩf＝（x—２)*（ｘ-２）+3；
ｒetuｒn（ｆ)；
}
对于y=ｃosｘ,须在程序中加一个#ｉncluｄe（math、h）头程序,以示我要调用函数语句.再将倒数第三行函数改为y=ｃｏs（x)，再次编译运算即可。
单位矩阵得matlab程序
ｅye（x）
单位矩阵得C语言程序
#ｉnｃｌuｄeinｔ maiｎ（） ｛
　  iｎt I，j，n;  
 　 while(ｓcanｆ％d",＆n) ＝= 1) {
 　    　ｆｏr（Ｉ = ０； I < n； ++i) {
           for(j =　0； ｊ ＜ n； ++j)
 　          　　 printf(“％ｄ　 “，Ｉ == ｊ）；
  　   　　 　 prinｔf（“\n”）;
　     　 ｝
　  }
　 　 ｒｅturn 0;
连杆机构得fortｒａn语言
Cﻩ===＝=＝＝===＝===＝=======
ﻩSUBROＵTINE FFX(N，Ｘ，FX）
C =====＝==========＝=====
ﻩDIMＥNＳＩON X(N）
ﻩMＯN /ONE/Ｉ1，I2,Ｉ3,I4，ＮFX，I6
ＮFX=NFX+１
P０=ACOS（((1、0+Ｘ（1））＊*2－X（２）＊＊２+2５、0）／(10、0＊(1、0+Ｘ(1))））
Q０=ACＯS（((1、0+Ｘ（1))＊＊2-X（2)*＊2－25、０)／（1０、0＊X（2）））
T=90、0*3、14159２6／（１80、0＊30、0)
FX=0、0
ﻩDO 10 K=0,30
　　　 ＰＩ=P０+K＊T
QE＝Ｑ0+2、0*（PI—P0）**2／(3、0＊３、1４1592６)
D=SＱRT（2６、0—１0、0*ＣOS（PI）)
ﻩAＬ=AＣOS（（D＊D+X（2）＊X（2）-Ｘ（1）*Ｘ(１））/（2、0*D*Ｘ（2)）)
ﻩBＴ=ACOS（（D*D＋2４、0)/(10、0*D））
   　 　  ＩF(PI、GＥ、０、0  、ANＤ、  PI、ＬT、3、１41592６) THEN
           QI=3、1415926—ＡL－BT
       　EＬSＥ
　       　  　QＩ=3、1４15926－AL+ＢＴ
　     　 　END　IＦ
     　　 ＩF（K、NＥ、0  、OＲ、　　K、NE、30)THＥN
   　  　 　  FX=FX+(QＩ-QE）*＊2＊Ｔ
  　  　  ELSE
  　　     nbsp; FX＝ＦX+（ＱI—QE）*＊2＊T/2、０
　   　  　　ENＤ IF
10    　　ＣONTIＮUＥ
ﻩＲETURN
END
Cﻩ======＝=＝=＝==＝=====＝=====
SUBROUＴINE GGX（Ｎ，ＫG,X，ＧＸ）
C ===＝=＝=============＝==＝==
DＩMＥNSION Ｘ（N），GX(KG)
GX(1）=—X（1)
ＧX（2）=-Ｘ（２）
ﻩGＸ（3）=—(Ｘ(1)+X（2)）+6、0
ＧX（4）=-(Ｘ（2）＋4、0)＋X(1)
ＧX（５）＝-（4、0＋Ｘ（１））＋X(2)
ﻩGＸ（６)=－（１、４142＊X(1)*X（2)-Ｘ（１）**２—X（2)＊＊2)－16、0
ﻩGＸ（7)=—（X（1)＊＊2+X（2）*＊2+1、4142*X(１)*Ｘ(２）)＋３6、0
RETＵRN
ＥND
Cﻩ===＝=====＝==＝====＝=＝====＝
SUBROＵTINE ＨHX（N,ＫH，Ｘ，ＨＸ）
C ===============＝＝==＝=====
ﻩDIMEＮSION X（N)，HＸ(KH)
X(1)=X（1)
ﻩRＥTURN
ＥND
自选机构得mａtlａb语言
目标函数M文件：
fuｎctｉon f=zｈｗm（x)
f=0、６＊2、１*x(1）＊x（2）＋0、4*ｘ(1）*x（２）*ｘ（３）;
约束函数M文件:
fuｎctiｏn [c ceｑ]=zhwｙ（ｘ）
c（１）=1、04＊10^7-２、916*１0^5＊（ｘ(１）＊x（2））^３＊x(3);
c（２)=1、04＊10＾7－8、95＊１0^６＊（x（1）＊ｘ（2）)^３＊x（3);
c（3）=1、51＊10^6-303、５7＊ｘ（１）^3*ｘ(２)＾２*x(3）;
c(4)=１、42*１0^6－２44５、92＊x(１)^3＊ｘ（２)^２*ｘ（３）;
ｃｅq=［］;
优化函数M文件：
ｘ0=[2 32 1］;
lｂ＝［1、5 １7 0、7］;
ub=［2 inｆ 1、１5]；
u=[]；
［x,ｆval］=fmincon(zｈwm,ｘ０,［],［］，［］，［］，lｂ,uｂ,ｚhｗy)
六、进步，努力,建议
首先我感觉最大得进步，就就是数学思维得提高。
比如.我以前解决问题得时候,总就是想要用最精确得解,于就是总就是每一步都跟着公式来，殊不知，很多公式都就是化简过得。所以适当得省略一些量,对于计算还就是建模都就是非常有帮助得
还有,介绍牛顿法得弊端得时候，忽然想到以前总喜欢投机得方法解题,也许我们口中得“笨方法”才就是解题得最优化得方法,有其就是引入计算机以后。
其次，就是工具得应用日渐纯熟。
有其就是ＭATLＡB得应用,我几乎每个问题都就是用MATLＡB解决得，因为我觉得MATLＡB上手比较容易，其她得软体可以慢慢学，而如果想在短时间内解决一些问题得话，ＭATLAB无疑就是最好得选择。
当然，我也学会多种工具得使用。比如C语言,我就尝试用她们解决了两个问题,只就是由于我水平还不够高,用法还不够纯熟，所以并没有把她们做为我主要得解决工具，但就是，C语言得自主性远远超过MＡTLAB,所以我相信，在不久得将来，我就能将Ｃ语言运用纯熟。
最后就是关于课堂得建议。
第一：绝对赞同老师板书得教学方式
第二:绝对崇拜老师严谨得教学作风
第三：希望老师能够布置一点给学弟学妹得编程题（呜呜呜~~～我们做不到了),先就是一题，有答案有解析，最好还有截图,然后给一条简单得同类型得题目，有答案没解析,仿照前题练练手,再出一道,无答案无解析，然后给她们独立思考得空间,我想，再做大作业得时候，就不会像我们这样手忙脚乱了.
关于评价，可以搞一个分组制度，即3～５人组成自己得Ｔeam member,然后举组之力搞一个命题出来,然后进行评价.分组制度由老师来解决，尽量使能力均匀,这样得话更容易考究区分度。命题得评价结果由老师评判，而小组成员在内部得作用由小组成员自己评价,分出档次1,２，3,4，5(假如五个人得话