- 征文：我，祝福祖国（精选合集）一号文库
- 征文：我，祝福祖国（精选合集）一号文库
- 最新国家开放大学电大《自然科学基础》形考任务2试题及答案一号文库
- XX社区软弱涣散党组织整顿方案一号文库
- 最新烟草工作心得体会 烟草店工作心得体会(汇总13篇)一号文库
- 2025年拒绝燃放烟花爆竹倡议书(模板10篇)一号文库
- “3+1”机制破解人居环境整治难题一号文库
- 党办：在市委党建领导小组会议上的讲话一号文库
- 集团转型发展和深化改革各项工作大会讲话五篇范文一号文库
- 县委书记在2025年全县领导干部会议上的发言（5篇范文）一号文库
- 部编版五年级下语文期末分类----16 文本积累测试题含参考答案一号文库
- 《消防整改报告》一号文库
- 部编教一材年级上册知识点复习一号文库
- 2025+ 2025年农村人居环境整治专项行动实施方案一号文库
- 儿童电台策划案一号文库
- 2025年县委常委书记到村支部讲党课材料（推荐5篇）一号文库
- 建筑工程收费标准一号文库
- 项目建设整改工作报告一号文库
- 苏教版语文三年级下册第三单元第9课《少年王勃》同步练习(I)卷一号文库
- 人教小学美术二上《第14课做一道拿手菜》[黄秀芳]【市一等奖】优质课一号文库
- XX市教育局XX市2025年度教育法治建设工作要点一号文库
- 梅学堂：苏教版数学5年级下册复学摸底测试卷2（含答案）一号文库
- 2025年XX市高考期间餐饮服务食品安全保障工作方案一号文库
- 消防安全整治百日攻坚行动工作方案一号文库
- 乡村人才振兴存在问题及建议一号文库
- 工商联组政协委员座谈会讲话提纲一号文库
- 在全市青年干部培训班开班仪式上的讲话：关于学习一号文库
- 学党史有感——星星之火可以燎原一号文库
- “五动”教学法——高效的课堂一号文库
桥式吊车运动控制系统的建模及MATLAB仿真(附程序)
简介
桥式起重机是横架于车间、仓库及露天堆场的上方，用来吊运各种物体的机械设备，通常称为“天车”或
“吊车”。它是机械工业、冶金工业和化学工业中应用最广泛的一种起重机械。实际生产中的桥式吊车（天车）类似，是一个MIMO复杂控制系统。
桥式吊车系统由三部分组成：桥架驱动系统，车体驱动系统和重物装吊系统。其工作流程为：先将重物起吊至预先设定好的高度，然后吊车运动将重物运到想要放置的位置上方，最后把重物下放到想要放置的位置上。
确定要研究的系统为桥式吊车运动控制系统
桥式吊车系统工作示意图见下图1：
p
M
m
x
z
F
θ
mg
图1
桥式吊车工作示意图
对于如上桥式吊车控制系统，首先做如下假设：
1)
吊车的行走运动仅限于吊车一个自由度，即假设桥架不运动，只有吊车在桥架上行走。
2)
吊车行走时吊装重物的绳索长度不变。
图中，x坐标为水平方向，z坐标为垂直方向。重物的摆动是由吊车与重物的运动产生的，可以根据动力学有关规律建立吊车及重物的运动方程式。
1)
在水平方向，吊车和重物整体受力为F(t)，由牛顿第二定律得
（1）
2)
在垂直于绳索方向，重物受力为，由牛顿第二定律得
（2）
由吊车在行走时吊装重物的绳索长度不变的假设可得出下面两个关系式：
（3）
（4）
式中，为绳索长度。
由（3）可得
（5）
（5）代入（1）得：
（6）
同样由式（4）可得：
（7）
将（5）（7）代入（2）得
（8）
又尽量小，所以有如下近似式：，将（6），（8）线性化可得：
（9）
（10）
由（9）和（10）计算得
（11）
和
（12）
3）
吊车驱动装置的方程式。吊车由电动机驱动，简化的认为电动机是一个时间常数为的一阶惯性环节，即它产生的驱动力F(t)与其控制电压v(t)之间满足方程式：
（13）
其中K为放大系数。
选择系统的输入、输出变量和状态变量
选择5个状态变量分别为：，，；
输入变量为：；
两个输出变量为：。
建立状态空间描述
根据（11）（12）（13）式可得出描述吊车运动系统的状态空间表达式为：
选取适当参数:对一个实际的桥式吊车吊车运动系统，假定具有如下各具体参数：M=1000kg,m=4000kg,l=10m,K=100N/V。将它们代入上面的状态空间表达式得：
分析系统的稳定性
用特征值法。在MATLAB中输入以下程序：
eig(A)
ans
=
0
0
0
+
2.2136i
0
2.2136i
系统的5个开环特征值不全位于S左平面上，有4个位于虚轴上，所以系统为临界不稳定。
系统输出仿真波形如下图所示：
判断系统的能控性
使用MATLAB判断系统的能控性，输入以下程序：
A=[0
0
0
0;0
0
-39.2
0
0.001;0
0
0
0;0
0
-4.9
0
0.0001;0
0
0
0
-1];
B=[0;0;0;0;100];
C=[1
0
0
0
0;0
0
0
0];
rct=rank(ctrb(A,B))
rct
=
根据判别系统能控性的定理，该系统的能控性矩阵满秩，所以该系统是能控的。
采用状态反馈进行系统综合因为系统是能控的，所以，可以通过状态反馈来任意配置极点。例如将极点配置在：s1=-0.16-j0.16
s2=-0.16+j0.16
s3,s4,s5=-1。
在MATLAB中输入：
P=[-0.16+0.16i,-0.16-0.16i,-1,-1,-1];
K=acker(A,B,P)
求出状态反馈矩阵K：
K
=
0.52245
4.8327
-1420.7
-137.21
0.0232
在MATLAB中输入
A-B*K
ans
=
0
0
0
0
0
0
-39.2
0
0.001
0
0
0
0
0
0
-4.9
0
0.0001
-52.245
-483.27
1.4207e+005
13721
-3.32
因此综合后系统的状态空间描述为：
采用MATLAB编写m文件进行仿真。
运行仿真程序，得到仿真曲线如下图：
将极点配置在：s1=-0.2
s2=-0.2
s3,s4,s5=-1。计算出K=[0.40816
5.3061
-1438.1
-119.06
0.024]此时输出y的仿真曲线如下：
将极点配置在：s1=-0.16+0.16i
s2=-0.16-0.16i
s3=-1
s4=-2
s5=-3。计算出K=[3.1347
25.339
-2145.4
553.73
0.0532]此时输出y的仿真曲线如下：
实验结论
通过比较3组不同的极点配置下状态反馈系统的输出响应曲线和原系统的输出响应曲线可以看出，不同的极点配置对系统性能有一定的影响，但只要极点都配置在S左平面，就可以保证系统具有一定的动态和稳态性能。
附录：程序
%*******************************************
%桥是吊车运动控制系统极点配置设计及仿真
2025.12.30
lwd
%*******************************************
%建立状态空间表达式
A
=
[0,1,0,0,0;0,0,-39.2,0,10^(-3);0,0,0,1,0;0,0,-4.9,0,10^(-4);0,0,0,0,-1]
B
=
[0;0;0;0;100]
C
=
[1,0,0,0,0;0,0,1,0,0]
D
=
[0;0]
%分析系统稳定性
eig(A)
%求A的特征值，通过特征值在S平面分布判断系统稳定性
%系统仿真输出波形
[num,den]
=
ss2tf(A,B,C,D);
sys1
=
tf(num(1,:),den);
sys2
=
tf(num(2,:),den);
figure
step(sys1)
figure
step(sys2)
%判断系统的能控性
rct
=
rank(ctrb(A,B))
obs
=
rank(obsv(A,C))
%如果能空则进行极点配置设计
if
==
rct
%第一次极点配置设计
P
=
[-0.16+0.16*i,-0.16-0.16*i,-1,-1,-1];
K
=
acker(A,B,P)
A1
=
A-B*K
[num,den]
=
ss2tf(A1,B,C,D);
sys1
=
tf(num(1,:),den);
sys2
=
tf(num(2,:),den);
figure
step(sys1)
figure
step(sys2)
%第二次极点配置设计
P
=
[-0.2,-0.2,-1,-1,-1];
K
=
acker(A,B,P)
A1
=
A-B*K
[num,den]
=
ss2tf(A1,B,C,D);
sys1
=
tf(num(1,:),den);
sys2
=
tf(num(2,:),den);
figure
step(sys1)
figure
step(sys2)
%第三次极点配置设计
P
=
[-0.16+0.16*i,-0.16-0.16*i,-1,-2,-3];
K
=
acker(A,B,P)
A1
=
A-B*K
[num,den]
=
ss2tf(A1,B,C,D);
sys1
=
tf(num(1,:),den);
sys2
=
tf(num(2,:),den);
figure
step(sys1)
figure
step(sys2)
end