社会实践课程设计

2022-07-14

第一篇：社会实践课程设计

课程设计实践心得

一、课程设计的目的

掌握电力系统潮流计算的基本原理。掌握并熟练运用MATLAB仿真软件。

采用MATLAB仿真软件，做出系统接线图的潮流计算仿真结果。

二、课程设计要求

1、在读懂程序的基础上画出潮流计算基本流程图

2、通过输入数据，进行潮流计算输出结果

3、对不同的负荷变化，分析潮流分布，写出分析说明。

4、对不同的负荷变化，进行潮流的调节控制，并说明调节控制的方法，并列表表示调节控制的参数变化。

5、打印利用MATLAB进行潮流分析绘制的系统图，以及潮流分布图。

三、课程设计基本内容

1. 对给定的网络查找潮流计算所需的各元件等值参数，画出等值电路图。

2. 输入各支路数据，各节点数据利用给定的程序进行在变电所在某一负荷情况下的潮流计算，并对计算结果进行分析。

3. 跟随变电所负荷按一定比例发生变化，进行潮流计算分析。 1) 4个变电所的负荷同时以2%的比例增大; 2) 4个变电所的负荷同时以2%的比例下降;

3) 1和4号变电所的负荷同时以2%的比例下降，而2和3号变电所的负荷同时以2%的比例上升;

4. 在不同的负荷情况下，分析潮流计算的结果，如果各母线电压不满足要求，进行电压的调整。(变电所低压母线电压10KV要求调整范围在9.5-10.5之间;电压35KV要求调整范围在35-36之间)。

5. 轮流断开环网一回线，分析潮流的分布。

6. 利用MATLAB软件，进行绘制系统图进行上述各种情况潮流的分析，并进行结果的比较。

7. 最终形成课程设计成品说明书。

四、问题分析

1、节点设置及分类

2设为节点1，10，根据系统图可知此系统为两端供电网路，将母线1，将变电所

1、

2、

3、4的高低压侧分别设为节点

2、

3、

4、

5、

6、

7、

8、9。并且，将节点1设为平衡节点，将节点10设为PV节点，其余节点设为PQ节点。

2、参数求取

设定基准值SB100MVA，UB220KV，所以ZB资料，计算发电厂、变压器及线路的参数。

(1)运用下列公式计算变压器参数：

SBUB2484根据题目原始

PUU%URX100S 1000S22kNKNT2TNNZ(2)计算线路参数

RTXTTZTZZTB

ZRjX(rjx)L

(3)变电所负荷分别为：

变电所1

SL=50+j30.987

变电所2

SL=40+j27.79 变电所3

SL=50+j30.987

变电所4

SL=60+j37.18

3、计算变压器分接头变比

变压器有5个抽头，电压调节范围为UN2*2.5%，UN对应的分接头开始时设变压器高压侧接主接头,降压变压器5个分接头时的非标准变比k*以备调压时选用

对变电所低压母线为35KV时，非标准变比与10KV时相同。

4、计算方法

利用牛顿拉夫逊法进行求解，用MATLAB软件编程，可以求解系统潮流分布根据题目的不同要求对参数进行调整，通过调节变压器变比和发电厂的电压，求解出合理的潮流分布，最后用PSAT进行潮流分析，将两者进行比较。

5、牛顿—拉夫逊法

(1)牛顿—拉夫逊法概要

首先对一般的牛顿—拉夫逊法作简单的说明。已知一个变量X函数为：

到此方程时，由适当的近似值X(n1)(0)220(15%)100.95511220 220(12.5%)10k*20.93211220 22010k*30.90911220 220(12.5%)10k*40.88611220 220(15%)10k*50.86411220 k*1f(X)0

出发，根据：

(n)X反复进行计算，当X的牛顿—拉夫逊法。

(n)Xf(X(n))(n1,2,......)f(X(n))

满足适当的收敛条件就是上面方程的根。这样的方法就是所谓

(n)这一方法还可以做下面的解释，设第n次迭代得到的解语真值之差，即X时，则：

把f(X(n))在X(n)的误差为附近对用泰勒级数展开

)f(X(n)f(X(n))0

)f(X(n)f(X(n))22!f(X(n))......0

上式省略去2以后部分

f(X(n))f(X(n))0

(n)X的误差可以近似由上式计算出来。

f(X(n)) (n)f(X)比较两式，可以看出牛顿—拉夫逊法的休整量和X用同样的方法考虑，给出n个变量的n个方程：

(n)的误差的一次项相等。

f1(X1,X2,,Xn)0f(X,X,,X)0212n fn(X1,X2,,Xn)0对其近似解X得修正量X可以通过解下边的方程来确定：

11f1f1f1,,Xn)f1(X1,X2xxx12nf(X,X,,X)X1f2n212f2f2X2x1x2 xnXfnfnfnnfn(X1,X2,,X)nxnx1x2fn,X2,,Xn的值。式中等号右边的矩阵都是对于X1这一矩阵称为雅可比(JACOBI)

xn,X2,,Xn后，得到如下关系矩阵。按上述得到的修正向量X1XnXn Xn,X2,,Xn更接近真实值。这一步在收敛到希望的值以前重复进行，一般要反这比X1复计算满足

maxX1n1X1n1,X2n1X2n1,,Xnn1Xnn1

为预先规定的小正数，Xnn1是第n次迭代Xn(2)牛顿法的框图及求解过程

用牛顿法计算潮流时，有以下的步骤： ①给这各节点电压初始值e(0),f(0);

的近似值。

②将以上电压初始值代入公式，求修正方程的常数项向量P(0),Q(0),(V2)(0); ③将电压初始值在带入上述公式，求出修正方程中系数矩阵的各元素。 ④解修正方程式e(0),f(0);

f(0)f(0);

⑥将e(1)，f(1)在带入方程式，求出P(1),Q(1),(V2)(1); ⑤修正各节点电压e(1)e(0)e(0)，f⑦检验是否收敛，即maxPi(1)(k),Qi(k)

如果收敛，迭代到此结束，进一步计算各线路潮流和平衡节点功率，并打印输出结果。如果不收敛，转回②进行下次迭代计算，直到收敛为止。

五、问题求解

电压是衡量电力系统电能质量的标准之一。电压过高或过低，都将对人身及其用电设备产生重大的影响。保证用户的电压接近额定值是电力系统调度的基本任务之一。当系统的电压偏离允许值时，电力系统必须应用电压调节技术调节系统电压的大小，使其维持在允许值范围内。本文经过手算形成了等值电路图，并编写好了程序得出节点电压标幺值，使其满足所要求的调整范围。

我们首先对给定的程序输入部分作了简要的分析，程序开始需要我们确定输入节点数、支路数、平衡母线号、支路参数矩阵、节点参数矩阵。

(1)为了保证整个系统潮流计算的完整性，我们把凡具有母线及发电机处均选作节点，这样，可以确定有10个节点，节点号见等值电路图。

(2)确定完节点及编号后，各条支路也相应确定了，网络中总计有13条支路，我们对各支路参数进行了计算。根据所给实际电路图和题中的已知条件，有以下公式计算各输电线路的阻抗和对地支路电容的标幺值和变压器的阻抗标幺值。选择电压基准值为UB=220KV和功率基准值SB=100MVA，所以ZB

六、误差分析

UB=484。 SB以系统在正常运行的情况下为例。利用PSAT软件，比对MATLAB中的B

1、B2矩阵，得到各个节点的电压状况。将PSAT软件运行的结果对MATLAB运行的结果进行校验。同时利用该软件对图形颜色的变化的描绘和数据对结果进行呈现。得到PSAT软件运行的结果各节点的电压标幺值，并与MATLAB程序运行得到的结果进行比较。

七、心得体会及总结

此次课程设计使我在潮流计算、MATLAB的使用方面均有所提高，但也暴露出了一些问题;理论知识的储备还是不足，对MATLAB的性能和特点还不能有一个全面的把握，对MATLAB中PSAT工具箱使用不够熟悉，相信通过以后的学习能弥补这些不足，从而达到一个新的层次。潮流计算是电力系统的最基本、最常用的分析计算。用以研究系统规划和运行中提出的各种问题。对规划中的电力系统，通过潮流计算可以检验所提出的电力系统规划方案能否满足各种运行方式的要求;对运行中的电力系统，通过潮流计算可以预知各种负荷变化和网络结构的改变会不会危及系统的安全，系统中所有母线的电压是否在允许的范围以内，系统中各种元件(线路、变压器等)是否会出现过负荷，以及可能出现过负荷时应事先采取哪些预防措施等。实际的情况远比我们计算的情况复杂，这让我深刻了解了潮流计算的重要性。精准的潮流计算不仅可以使电网处于稳定且平衡的状态，这对电网的安全运行起到关键性的作用，还可以实现运行的经济性。所以我认为学好电力系统的关键在于学好潮流计算，电力系统的安全运行不是儿戏，一个很小的失误就可能造成整个电网的崩溃，这给国家和社会带来的危害和损失之大是可想而知的。学好电力系统这门课、学好潮流计算是我们电专业学生的当务之急。没有坚实的基础就不可能建造出高楼大厦，所以，要想在以后的工作岗位上有所成就，就必须打好基础，努力提高自己的专业素养。

通过此次课程设使得我更加扎实的掌握了有关潮流计算和MATLAB方面相关的知识，在设计过程中虽然遇到了一些问题，但经过一次又一次的思考，一遍又一遍的检查终于找出了原因所在，也暴露出了前期我在这方面的知识欠缺和经验不足。实践出真知，通过亲自动手制作，使我们掌握的知识不再是纸上谈兵。过而能改，善莫大焉。在课程设计过程中，我们不断发现错误，不断改正不断领悟和不断获取。最终的检测调试环节本身就是在践行“过而能改，善莫大焉”的知行观。这次课程设计终于顺利完成了。在设计中遇到了很多问题，最后在老师的指导下。终于游逆而解。在今后社会的发展和学习实践过程中，一定要不懈努力不能遇到问题就想到要退缩，一定要不厌其烦的发现问题所在，然后一一进行解决。只有这样，才能成功的做成想做的事，才能在今后的道路上披荆斩棘，而不是知难而退，那样永远不可能收获成功，收获喜悦。也永远不可能得到社会及他人对你的认可。

课程设计诚然是一门专业课，给我很多专业知识以及专业技能上的提升，同时又是一门讲道课;一门辩思课，给了我许多道，给了我很多思，给了我莫大的空间。同时设计让我感触很深。使我对抽象的理论有了具体的认识。通过这次课程设计。我掌握了MATLAB中各种元器件的识别和使用，熟悉了和了解了潮流计算的基本原理。以及如何提高实验的求解精度和收敛性能。掌握了MATLAB软件的方法和技术。

我认为，在这学期的实验中，不仅培养了独立思考、动手操作的能力。在各种其它能力上也都有了提高。更重要的是，在实验课上，我们学会了很多学习的方法。而这是日后最实用的，真的是受益匪浅。要面对社会的挑战，只有不断的学习、实践、再学习、再实践。这对于我们的将来也有很大的帮助。以后不管有多苦，我想我们都能变苦为乐。找寻有趣的事情，发现其中珍贵的事情。就像中国提倡的艰苦奋斗一样，我们都可以在实验结束之后变的更加成熟，会面对需要面对的事情。

回顾起此课程设计，至今我仍感慨颇多。从理论到实践，在这段日子里，可以说得是苦多于甜，但是可以学到很多很多的东西，同时不仅可以巩固了以前所学过的知识。而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的。只有把所学的理论知识与实践相结合起来从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题，可以说得是困难重重，但可喜的是最终都得到了解决。

实验过程中，也对团队精神的进行了考察，让我们在合作起来更加默契，在成功后一起体会喜悦的心情。果然是团结就是力量。只有互相之间默契融洽的配合才能换来最终完美的结果。此次设计也让我明白了思路即出路，有什么不懂不明白的地方要及时请教或上网查询。只要认真钻研、动脑思考、动手实践，就没有弄不懂的知识，收获颇丰。课程设计是培养学生综合运用所学知识,发现,提出，分析和解决实际问题，锻炼实践能力的重要环节,是对学生实际工作能力的具体训练和考察过程。随着科学技术发展的日新日异，MATLB已经成为当前电力系统重要的仿真软件之一。因此作为二十一世纪的电气专业的大学来说掌握MATLB的开发技术是十分重要的。回顾起此次课程设计，至今我仍感慨颇多，的确，从选题到定稿，从理论到实践，在整整两星期的日子里，可以说得是苦多于甜，但是可以学到很多很多的东西。同时不仅可以巩固了以前所学过的知识，而且学到了很多在书本上所没有学到过的知识。通过这次课程设计使我懂得了理论。

与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。在设计的过程中遇到问题，可以说得是困难重重，这毕竟第一次做的，难免会遇到过各种各样的问题，同时在设计的过程中发现了自己的不足之处。对以前所学过的知识理解得不够深刻，掌握得不够牢固。比如说MATLAB中原理图的搭建以及潮流计算的基本原理掌握得不好。通过这次课程设计之后，一定把以前所学过的知识重新温故。

这次课程设计终于顺利完成了，在设计中遇到了很多编程问题。最后在陈老师的辛勤指导下，终于游逆而解。同时，在陈老师的身上我学得到很多实用的知识。

整个设计通过了软件调试和图形搭建以及程序设计。我想这对于自己以后的学习和工作都会有很大的帮助。在这次设计中遇到了很多实际性的问题，在实际设计中才发现书本上理论性的东西与在实际运用中的还是有一定的出入的。所以有些问题不但要深入地理解，而且要不断地更正以前的错误思维。一切问题必须要靠自己一点一滴的解决，而在解决的过程当中你会发现自己在飞速的提升。对于教材管理系统，其程序是比较简单的，主要是解决程序设计中的问题。而程序设计是一个很灵活的东西，它反映了你解决问题的逻辑思维和创新能力，它才是一个设计的灵魂所在。因此在整个设计过程中大部分时间是用在MATLAB图形的搭建和参数的设置。其中很多子程序是可以借鉴书本上的，但怎样衔接各个子程序才是关键的问题所在，这需要对系统的结构很熟悉。因此可以说系统的设计是软件和硬件的结合，二者是密不可分的。通过这次课程设计我也发现了自身存在的不足之处，虽然感觉理论上已经掌握，但在运用到实践的过程中仍有意想不到的困惑，经过一番努力才得以解决。这也激发了我今后努力学习的兴趣，我想这将对我以后的学习产生积极的影响。其次，这次课程设计让我充分认识到团队合作的重要性，只有分工协作才能保证整个项目的有条不紊。另外在课程设计的过程中，当我们碰到不明白的问题时，陈老师总是耐心的讲解，给我们的设计以极大的帮助，使我们获益匪浅。因此非常感谢老师的教导。通过这次设计，我懂得了学习的重要性，了解到理论知识与实践相结合的重要意义，学会了坚持、耐心和努力，这将为自己今后的学习和工作做出了最好的榜样。我觉得作为一名电气工程专业的学生，这次课程设计是很有意义的。更重要的是如何把自己平时所学的东西应用到实际中。虽然自己对于这门课程懂的并不多，很多基础的东西都还没有很好的掌握，觉得很难，也没有很有效的办法通过自身去理解。但是通过这次的课程设计的沟通和各个小组同学探讨，渐渐对这门课程逐渐产生了些许的兴趣，自己开始主动学习并逐步从基础慢慢开始弄懂它。

第二篇：《机械优化设计》课程实践

合肥工业大学研究报告

班级：机设10-5 学号：20100523 姓名：吴亮宏授课老师：王卫荣日期：

2012年5月7日

一维搜索程序作业(0.618法 )

a=2，b=6，d=0.00001 y=(x-4)*(x-4)-7 #include #include void main() { float a,b,c=0.618,aa[3],y[3],d; scanf("%f,%f,%f",&a,&b,&d); aa[1]=b-c*(b-a); aa[2]=a+c*(b-a);

y[1]=(aa[1]-4)*(aa[1]-4)-7;

y[2]=(aa[2]-4)*(aa[2]-4)-7;

do{ if(y[1]>y[2])

{ a=aa[1];aa[1]=aa[2];y[1]=y[2]; aa[2]=a+c*(b-a);

y[2]=(aa[2]-4)*(aa[2]-4)-7;

}

else

{ b=aa[2];aa[2]=aa[1];y[2]=y[1]; aa[1]=b-c*(b-a);

y[1]=(aa[2]-4)*(aa[2]-4)-7;

}

}while(fabs((b-a)/b)>d); aa[0]=(a+b)/2;

y[0]=(aa[2]-4)*(aa[2]-4)-7; printf("a*=%f ",aa[0]); printf("y=%f ",y[0]);

}

2、单位矩阵程序作业

#include void main () { inti,j,p ,n; printf("please input a number");

} scanf("%d",&n); for(i=0;i

{

if(i==j) p=1; else p=0; printf("%2d",p); if(j==n-1) printf(" ");

}

3曲柄摇杆机构优化设计(约束随机法)

课本P241例8-5 1程序文本

procedure ffx; //目标函数 var

p0,q0,T,PI,QE,D,AL,BT,QI:real;

K:integer; test:string; begin with form1.hfgd do

begin

NFX:=NFX+1;

p0:=arccos((sqr(1.0+X[1])-sqr(X[2])+25.0)/(10.0*(1.0+X[1])));

//φ0

q0:=arccos((sqr(1.0+X[1])-sqr(X[2])-25.0)/(10.0*X[2]));

//ψ0

T:=90.0/30.0*(3.1415926/180.0); //90度分30等分

FX:=0.0;

For K:=0 To 30 do;

//30次循环计算 begin

PI:=p0+K*T; //φ曲柄转角

QE:=Q0+2.0*sqr(PI-p0)*2/(3.0*3.1415926); //理想输出角

D:=SQRT(26.0-10.0*COS(pI));

//辅助线长度

AL:=ArcCos((D*D+X[2]*X[2]-X[1]*X[1])/(2.0*D*X[2])); //α角

BT:=arccos((D*D+24.0)/(10.0*D)); //β角 IF ((PI>=0.0) AND (PI<3.1415926)) //ψ(i)实际输出角 THEN

QI:=3.1415926-AL-BT ELSE

QI:=3.1415926-AL+BT; IF((K<>0) OR (K<>30)) THEN

FX:=FX+sqr(QI-QE)*T //理想输出角与实际输出角曲线包围面积 ELSE FX:=FX+sqr(QI-QE)*T/2.0; end; end; end; procedure ggx; //约束函数 begin with form1.hfgd do

begin GX[1]:=-X[1]; GX[2]:=-X[2]; GX[3]:=-(X[1]+X[2])+6.0; //3,4,5为曲柄存在条件 GX[4]:=-(X[2]+4.0)+X[1]; GX[5]:=-(X[1]+4.0)+X[2]; GX[6]:=-(1.4142*X[1]*X[2]-X[1]*X[1]-X[2]*X[2])-16.0; //6,7为传动角条件

GX[7]:=-(X[1]*X[1]+X[2]*X[2]+1.4142*X[1]*X[2])+36.0; end; end; 2输入截图

3结果

4工程实例

：制造一体积为200 m3长度不小于5m，不带上盖的箱盒，试确定箱盒的长x1，宽x2，高x3，使箱盒用料最省。 (1)、设计变量的选择与目标函数的建立设计变量为箱盒的长x1，宽x2，高x3 本设计要求箱盒用料最省，故以箱盒的表面积为目标函数 F(x)=x1x2+2(x2x3+x1x3) (2)、约束条件

根据体积要求---- x1x2x3=200 根据体积要求---- x1>5 x2>0

x3>0 g1(x)=- x1-5≤0 g2(x)=- x2≤0 g3(x)=- x3≤0 g4(x)=x1x2x3-200.001 g4(x)=199.999-x1x2x3

(3)、选择方法约束随机法 (4)、文本程序

procedure ffx; //目标函数 begin with form1.hfgd do

begin fx:= x[1]*x[2]+2*(x[2]*x[3]+x[1]*x[3]); end; end; procedure ggx; //约束函数 begin with form1.hfgd do

begin gx[1]:=-x[1]+5; gx[2]:=-x[2]; gx[3]:=-x[3]; gx[4]:=x[1]*x[2]*x[3]-200.0001; gx[5]:=199.9999-x[1]*x[2]*x[3]; end; end;

{procedurehhx; begin with form1.hfgd do

begin hx[1]:=x[1]*x[2]*x[3]-200; end; end;

} End .输入截图结果截图

5心得体会

这次对机械优化的学习，学生感受颇多。

通过上机实验，我深深的感觉到机械优化设计是一门实践性很强的课程，它考验了我们的各种能力和素质。首先我们要有一定的编程能力，这样才能把我们所要表达的以程序的形式表达出来。例如以c语言来编写单位矩阵;其次就是对机械优化设计的一些算法如惩罚函数法的了解，这样才能选择合适的算法去解决问题，得出最优解;另外这也是对我们的灵活运用优化设计方法解决工程设计问题的能力的考验。

在刚刚开始学习这门课程时，感觉这门课的数学味道好浓，感觉就是在讲数学课。因为没什么阅历，这不清楚这门课程的应用性。慢慢讲道后面，特别是听过一些例题后，对于其在工程实际中的应用才有些感受。传统的机械设计方法比较粗放，依赖于设计者的经验，往往考虑到安全性等一些原因，而不得不牺牲一些性能和经济性，这对设计产品是不利的。而优化设计将程序化的思想引入了设计工程，借助于计算机的计算，不断优化设计结果，不仅设计结果更优，效率也更高。优化设计是工程设计的发展方向，传统的依赖设计者经验的设计方法已不能适应社会发展了。努力学习这门课，对于我们以后的发展非常有用。

这是我第一次使用delphi语言来编写程序。由于自己在电脑方面的知识匮乏，所以用起来相当吃力。在调试程序的过程中，不断尝试，再用惩罚函数法处理小盒问题时，出现了四个错误，但就是不知道怎么改，于是就慢慢试，最后发现是自己的输入表达出了问题。

通过这次计算机实验，突然感觉其实程序是个挺好玩的东西。以前感觉语言啊，程序啊都是是一种很模糊的概念在我的脑海中，看到他们头皮都发麻。所以昨晚刚刚开始时感觉什么都不懂，无从下手，都不知道怎么打开程序。程序真是一种神奇的东西吧，我想这也是我第一次体会到用程序的乐趣吧。当你用排除一个个错误后，最终程序顺利运行，得出了结果。那在喜不自禁的感觉真是让人很舒服。我喜欢这种小小的成就感。

在这门课的学习中，我被老师深厚的学术功力和风趣的教学方法所折服。本来觉得这门课挺枯燥的，可是在您的课堂上，却欢声不断，这确实是老师的个人魅力。感学老师这些天来对我们的无私教诲，不仅教我们课本上的知识，还和我们交流您出国访问的心得，这扩大了我的视野。还记得您在访问日本的报告中还做了一首小诗，虽然不记得诗句，可还是被老师的人文素养所打动。在我们这些年轻人当中，已经很是见有人会作诗啦。作为一名工科学生，提高自己的人文素质也是同样重要的。这就是我在做完试验后的一些小小感想。谢谢王老师!

吴亮宏

20100523

2013

5 7

第三篇：《机械优化设计》课程实践报告

(课程实践报告封面模版)

合肥工业大学

《机械优化设计》课程实践

研究报告

班级：机设10 -04 学号： 20100495 姓名：李健授课老师：王卫荣日期： 2012年月日

一主要内容

1、一维搜索程序作业

A.λ = 0.618的证明..........................................1 B.编写用0.618法求函数极小值的程序..........................2

2、单位矩阵程序作业............................................4

3、其他工程优化问题..................................9 4连杆机构问题.....................................12

二实践心得体会...............................15

一: 主要内容

1. 一维搜索程序作业：

A.λ = 0.618的证明 (y2 > y1) 证明：0.618法要求插入点α

1、α2 的位置相对于区间 [a，b] 两端点具有对称性，即

已知 a1=a2 , 要求α1=α2 由于α1=b-λ(b-a) α2=a+λ(b-a) 12 若使α1=α2111222211则有：b-λ(b-a)=a+λ(b-a)= a1+λ(b-a) 121因此: b- a1=(λ+λ)( b- a1) 1

2( b- a1)( λ+λ-1)=0 1因为: b= a1 2所以: λ+λ-1=0 则有: 取方程正数解得

1若保留下来的区间为 [α1，b] ，根据插入点的对称性，也能推得同样的λ的值。其0.618法的程序框图如下:

B.编写用0.618法求函数极小值的程序例：(1)a=0 ，b=2π，f(x)=cox(x)

2 (2)a=0 ，b=10， f(x)=(x-2)+3(1)

#include #include void main(void) {

int i; float a1,a2,aa,y1,y2,ymin,e; float a=0,b=2*3.14159,n=0.618; a1=b-n*(b-a); a2=a+n*(b-a); print(“输入精度：”); scanf(“%f”,&e); for(i=0;i=10000;i=i++) {

y1=cos(a1); y2=cos(a2); if(y1

a=a1; a1=a2; a2=a+n*(b-a); } If(y1

b=a2; a2=a1; a1=b-n*(b-a); } if(fabs(b-a)/b

(2)

#include #include void main(void) {

int i; float a1,a2,aa,y1,y2,ymin,e; float a=0,b=10,n=0.618; a1=b-n*(b-a); a2=a+n*(b-a); print(“输入精度：”); scanf(“%f”,&e); for(i=0;i=10000;i=i++) {

y1=(a1-2)*(a1-2)+3; y2=(a2-2)*(a2-2)+3; if(y1>=y2) {

a=a1; a1=a2; a2=a+n*(b-a); } If(y1

b=a2; a2=a1; a1=b-n*(b-a); } if(fabs(b-a)/b

2.单位矩阵程序作业

编写生成单位矩阵的程序

程序文本

#include void main(void) { int a[100][100];

int N,i,j; printf("请输入所要输出矩阵的阶数(最多100阶):"); scanf("%d",&N); printf("输出的矩阵阶数为%d ",N); printf(" N "); /*****制作表头*****/

for(i=0;i

printf("%3d",i+1); printf(" ");

for(i=0;i

printf("---"); /*****分割线*****/

printf(" ");

for(i=0;i<100;i++) /*****数组赋值*****/

for(j=0;j<100;j++)

{

if(i==j)

a[i][j]=1;

else

a[i][j]=0;

}

for(i=0;i

{

printf("%2d:",i+1); /*****纵列序号*****/

for(j=0;j

{

printf("%3d",a[i][j]);

} printf(" "); } }

结果显示

从键盘输入9，显示9阶单位矩阵，结果如下

3. 其他工程优化问题

有一箱形盖板，已知长度

L=600mm，宽度b=60mm，厚度ts=0.5mm 承受最大单位载荷

4q=60N/cm，设箱形盖板的材料为铝合金，其弹性模量E710MPa，泊松比0.3，许用弯曲应力70MPa，许用剪应力45MPa，要求在满足强度、刚度和稳定性条件下，设计重量最轻的结构方案。数学模型的建立：

1)设计变量：取结构的翼板厚度tf 和高度h为设计变量，即

tfx1x

hx2 2)目标函数：取结构的总重量最轻为目标函数，计算公式为

f(x)260t20.5h120x1x2 3)约束条件：

g1(x)[]max[]7g2(x)1x1x210max457g3(x)c1x13x210max45g4(x)1.5121x1x210fmax32111x2104g5(x)x10g6(x)x204)、根据目标函数和约束条件在Delphi程序坏境下编制程序如下：

5)、利用Delphi程序解决工程优化问题

使用复合形发寻找最优点在Delphi程序下输入所需参数值：

6)、根据目标函数和约束条件以及输入的参数值使用Delphi程序进行计算找到工程优化问题的优化极值：

4 连杆机构问题描述

图 3-1 机构简图

设计一曲柄连杆摇杆机构，要求曲柄l1从l1从m090时，摇杆l3的转角最佳再现已知的运动规律：E0围内变化。

3.12 数学模型的建立 2(0)2且l1=1，l4=5，0为极位角，其传动角允许在45135范3lxx1x2l2l3 设计变量：这里有两个独立参数l2和3。因此设计变量为

Tt目标函数：将输入角分成30等分，并用近似公式计算，可得目标函数的表达式

fxiEiii12i130

约束条件：

GX(1)=-X(1)0 GX(2)=-X(2) 0 GX(3)=-(X(1)+X(2))+6.00 GX(4)=-(X(2)+4.0)+X(1) 0 GX(5)=-(4.0+X(1))+X(2) 0 GX(6)=-(1.4142*X(1)*X(2)-X(1)**2-X(2)**2)-16.00 GX(7)=-(X(1)**2+X(2)**2+1.4142*X(1)*X(2))+36.00

3.13 程序编制

procedure ffx;

var

p0,q0,T,PI,QE,D,AL,BT,QI:real;

K:integer;

test:string; begin

with form1.rand do

begin

NFX:=NFX+1;

p0:=arccos((sqr(1.0+X[1])-sqr(X[2])+25.0)/(10.0*(1.0+X[1])));

q0:=arccos((sqr(1.0+X[1])-sqr(X[2])-25.0)/(10.0*X[2]));

T:=90.0/30.0*(3.1415926/180.0);

FX:=0.0;

For K:=0 To 30 do;

begin

PI:=p0+K*T;

QE:=Q0+2.0*sqr(PI-p0)*2/(3.0*3.1415926);

D:=SQRT(26.0-10.0*COS(pI));

AL:=ArcCos((D*D+X[2]*X[2]-X[1]*X[1])/(2.0*D*X[2]));

BT:=arccos((D*D+24.0)/(10.0*D)); IF ((PI>=0.0) AND (PI<3.1415926)) THEN

QI:=3.1415926-AL-BT ELSE

QI:=3.1415926-AL+BT; IF((K<>0) OR (K<>30)) THEN

FX:=FX+sqr(QI-QE)*T ELSE FX:=FX+sqr(QI-QE)*T/2.0;

end;

end; end; procedure ggx; begin

with form1.rand do

begin

GX[1]:=-X[1];

GX[2]:=-X[2]; GX[3]:=-(X[1]+X[2])+6.0; GX[4]:=-(X[2]+4.0)+X[1]; GX[5]:=-(X[1]+4.0)+X[2]; GX[6]:=-(1.4142*X[1]*X[2]-X[1]*X[1]-X[2]*X[2])-16.0; GX[7]:=-(X[1]*X[1]+X[2]*X[2]+1.4142*X[1]*X[2])+36.0;

end; end; 3.14使用Delphi程序验证连杆机构问题方法：随机方向法。 Delphi程序的使用：

3.15验证结果显示

二:实践心得体会

总的看来,机械优化设计是适应生产现代化要求发展起来的,是一门崭新的学科。它是在现代机械设计理论的基础上提出的一种更科学的设计方法,它可使机械产品的设计质量达到更高的要求。因此,在加强现代机械设计理论研究的同时,还要进一步加强最优设计数学模型的研究,以便在近代数学、力学和物理学的新成就基础上,使其更能反映客观实际。同时机械优化设计的研究还必须与工程实践、数学力学理论、计算技术和电子计算机的应用等紧密联系起来,才能具有更广阔的发展前景。

作为21世纪的大学生，要使自己适应社会需求，首先在做任何事之前都应该有正确的态度看待问题，把这些想法作为促使自己进步的动力，再去学习课本知识，效果应该很不一样，有了想法就付诸行动，随着对课本内容的学习跟老师的讲解，发现并不是像自己在学期初想的那样困难，特别是在老师介绍了一些与机械优化设计相关的计算机语言和计算机软件后，真正体会到科学优化设计的强大跟简洁明了，与传统优化设计方法相比较，大大提高了设计效率和质量。作为机械专业的一名学生，本课程，掌握最优化问题的基本解决方法，从多个可能的方案中选出最合适的、能实现预定最优目标的最优方案有着很现实的意义，为今后的工程实际提供了良好的理论储备。在学完课程之后，反思自己在学习过程中的得失，深深体会到，不论在人生的哪个阶段，都要对自己负责，做任何事都要耐心，细致，“千里之行，始于足下”，学会在物欲横流的社会大潮中，坚持踏踏实实走好人生的每一步。