matlab(四连杆优化设计)

grad =
*
)
%函数在最优点处梯度信息
hessian =
…
>>
5）结果分析
① 采用 fmincon 求解的最优值： =[; ]; ；
采用算法：中型算法（mediun-scale）。
%函数在最优点处海塞矩阵
~
这与课本给出的最优解： =[;]， =相比，计算精度更高，最优解的数值更精确，故计 算准确度高。 ② 用 matlab 绘制输入——输出曲线关系图
<
上图中（单位为“度”）蓝色的线代表曲柄摇杆机构的实际输出角与输入角的关系，红色的 线代表理想输出角与输入角的关系。 可以看出：
}
实际输出和理论输出曲线之间存在线性误差，其最大线性误差为 ，误差在允许的范围 之内，故结果的可信度也较大，运用 matlab 优化工具箱计算所得结果正确。
小结
通过结合实际问题的分析，计算，求解，更加深入地了解和掌握机械优化设计的过 程和步骤，比较重要的步骤是数学模型的建立，以及设计变量的选取，以及数学模型的 尺度变换，根据机构实际工作需要，建立目标函数的约束条件等等，当数学模型建好以 后，剩下的工作可以再 matlab 里面完成，而 matlab 里面的优化工具箱，给用户提供了 多种优化函数，使用者只需要将数学模型按要求编写成子程序嵌入已有的优化程序即 可。
⑵精度的控制：为了得到更加精确地解，需要设置优化函数的控制精度，函数本身默认精 度为 1e-4，精度比较低，通过 options 的设置将精度提高到 1e-9,这样得到的结果更精确。
&
以上两点通过设置options参数即可： options=optimset('largescale','off','display','off','Algorithm','active-set','TolF un',1e-9);
在设计过程中也遇到一些困难，比如说在在用 matlab 计算时，计算机已知处于 busy 状态，得不到函数的最优解，最后反复的检查，终于找的了其原因，是由于初始点选择 不恰当引起的，如果初始点选择得好，可以节省计算时间和计算空间，故初始点的选取 比较重要。
:
！
、
&
、
附录
1. 编写目标函数 M 文件：
,'TolFun',1e-9); %采用中型算法，设计精度为1e-9
[x,fval,exitflag,output]=fmincon(@myfun,x0,A,b,[],[],lb,[],@constrain,o
ptions)
%调用fmincon函数
机械优化设计在 matlab 中的应用
东南大学机械工程学院**
一 优化设计目的：
在生活和工作中，人们对于同一个问题往往会提出多个解决方案，并通过各方面的论 证从中提取最佳方案。最优化方法就是专门研究如何从多个方案中科学合理地提取出最佳方 案的科学。由于优化问题无所不在，目前最优化方法的应用和研究已经深入到了生产和科研 的各个领域，如土木工程、机械工程、化学工程、运输调度、生产控制、经济规划、经济管 理等，并取得了显著的经济效益和社会效益。
二 优化设计步骤：
1.机械优化设计的全过程一般可以分为如下几个步骤： 1)建立优化设计的数学模型；
'
2)选择适当的优化方法； 3)编写计算机程序； 4)准备必要的初始数据并伤及计算； 5)对计算机求得的结果进行必要的分析。 其中建立优化设计数学模型是首要的和关键的一步，它是取得正确结果的前提。优化方 法的选取取决于数学模型的特点，例如优化问题规模的大小，目标函数和约束函数的性态以 及计算精度等。在比较各种可供选用的优化方法时，需要考虑的一个重要因素是计算机执行 这些程序所花费的时间和费用，也即计算效率。 2.建立数学模型的基本原则与步骤 ① 设计变量的确定； 设计变量是指在优化设计的过程中，不断进行修改，调整，一直处于变化的参数称为 设计变量。设计变量的全体实际上是一组变量，可用一个列向量表示：
3）约束条件的确定
① 曲柄摇杆机构应满足曲柄存在条件，可得
￥
②曲柄摇杆机构的传动角
可得
'
把约束条件简化( )
=5
=
=
，
【
其中 满足条件，故最后一共有两个设计变量（
4）优化计算
），7 个约束条件。
)
!
① 此问题的图解见上图，有 7 个约束条件构成了改优化模型的可行域，而最优解在可行域 内。
②优化方法选择： 该问题属于一般的约束非线性最优化类型，可以使用 matlab 优化工具箱里面的
式中
—期望输出角，
；
m — 输入角等分数；
—实际输出角，由下图得：
<
a) 0≤ <
b)
（3） 式中
{
将输入角分成 30 等分，并用近似公式计算，可得目标函数的表达式：
f(x)= 由题意知，传动角的变化范围是
，
，则上式中变量的最后形式可以写成：
将 带入（4） （5）得：
·
为当
时的理想输出角，其值在题目中已经给出：
所有的程序编好以后，在命令窗口输入：youhua 得到的 matlab 的运行结果如下： >> youhua
x=
>
%最优解
fval =
*
%目标函数最优点的值
exitflag = 5
%标志值，’5’表示重要方向导数小于规定的容许范 围并且约束违背小于
output =
`
iterations: 12
%迭代次数
message: [1x780 char]
%跳出信息
lambda =
lower: [2x1 double] upper: [2x1 double] eqlin: [0x1 double]
\
eqnonlin: [0x1 double] ineqlin: [5x1 double]
ineqnonlin: [2x1 double]
funcCount: 40
%函数的评价次数
lssteplength: 1
stepsize:
algorithm: 'medium-scale: SQP, Quasi-Newton, line-search'
%采用的
中型算法
firstorderopt:
%一阶最优性条件
constrviolation:
—
d(i)=acos((r^2+24)/(10*r));
if a(i)<=pi
e(i)=pi-c(i)-d(i);
%计算摇杆输出的实际值
else if a(i)<=2*pi
e(i)=pi-c(i)+d(i);
end
end
a(1)=a0;
f=f+((b(i)-e(i))^2)*(a(i)-a(i-1));
function f=myfun(x)
f=0;
%函数f赋初值
a0=acos(((1+x(1))^2-x(2)^2+25)/(10*(1+x(1)))); %初始计算点曲柄和摇杆的角度
b0=acos(((1+x(1))^2-x(2)^2-25)/(10*x(2)));
i=2;
while(i<=31)
%设置迭代次数为30次
a(i)=a0+(pi/2)*(i/30);
% 计算曲柄各分度的角度值
b(i)=b0+2*(a(i)-a0)^2/(3*pi);
% 计算摇杆各分度的角度值
r=sqrt(26-10*cos(a(i)));
c(i)=acos((r^2+x(2)^2-x(1)^2)/(2*x(2)*r));
三 实例分析 （机械优化设计 P241 页例 8-5）
设计一曲柄摇杆机构如图，要求：
曲柄 从
且
已
知
范围内变化。
分析：
1） ' 2） 设计变量的确定
决定机构尺寸的各杆长度，以及当摇杆按已知运动规律开始运行时，曲柄所载的位置角 应列为设计变量，即：
X=
=
考虑到机构的杆长按比例变化时，不会改变其运动规律，因此在计算时常取 ，而其他杆长则按比例取为 的倍数。若取曲柄的初始位置角为极位角，则 及相应摇杆 位置角 均为杆长的函数，几何图形关系如右图,其关系式为：
lb=[1;1];
%设计变量的下界
x0=[4;2];
%迭代初始点
A=[-1,0;0,-1;-1,-1;1,-1;-1,1];
%线性不等式约束
b=[-1;-1;-6;4;4];
options=optimset('largescale','off','display','off','Algohm','active-set'
‘fmincon’函数进行求解。 ⑴ fmincon 里面算法的选择：fmincon 里面一共提供了‘largescale’，‘ 'medium-scale’
两 种 算 法 ， 由 于 此 问 题 只 有 两 个 设 计 变 量 ， 维 数 较 低 ， 故 采 用 ‘ medium-scale ’ 算 法 。 ‘medium-scale’算法是采用 SQP，算法中 Hessian 阵可以通过 BFGS 迭代，初始 Hessian 阵任给。注意 BFGS 公式中 q 项是需要计算目标函数梯度得到的。所以 Hessian 矩阵的近似 计算是需要用到有限差分法。在采用‘'medium-scale’算法时，需提供其设计变量的初始 点 的信息，而初始点的选择也将影响计算得收敛性和收敛速度，如果初始点选择得不恰 当，可能最后函数不能收敛，得不到计算结果。