哈工大机电系统控制基础大作业一

《机电系统控制基础》大作业一

基于MATLAB的机电控制系统响应分析

哈尔滨工业大学

2013年11月8日

1

作业题目

1. 用MATLAB 绘制系统2

()25

()()425

C s s R s s s Φ=

=++的单位阶跃响应曲线、单位斜坡响应曲线。 2. 用MATLAB 求系统2()25

()()425

C s s R s s s Φ==++的单位阶跃响应性能指标：上升时间、峰值时间、调节时间和超调量。

3. 数控直线运动工作平台位置控制示意图如下：

X i

伺服电机原理图如下：

L R

(1)假定电动机转子轴上的转动惯量为J 1，减速器输出轴上的转动惯量为J 2，减速器减速比为i ，滚珠丝杠的螺距为P ，试计算折算到电机主轴上的总的转动惯量J ；

(2)假定工作台质量m ，给定环节的传递函数为K a ，放大环节的传递函数为K b ，包括检测装置在内的反馈环节传递函数为K c ，电动机的反电势常数为K d ，电动机的电磁力矩常数为K m ，试建立该数控直线工作平台的数学模型，画出其控制系统框图；

(3)忽略电感L 时，令参数K a =K c =K d =R=J=1，K m =10，P/i =4π，利用MATLAB 分析kb 的取值对于系统的性能的影响。

2

题目1 （1）操作过程

在Matlab 中进行编程，然后进行仿真,得到曲线。

（2）仿真结果及分析

00.5

1

1.5

22.533.5t(sec)

x (t )

00.2

0.4

0.6

0.8

1

1.2

1.4

t(sec)

x (t )

左图为单位阶跃响应曲线

右图为单位斜坡响应曲线

（3）源代码

t = [0:0.01:3.5];u = t nG = [25]; dG = [1 4 25]; G = tf(nG,dG); [y,T] = step(G,t); z = lsim(G,u,t);

subplot(121),plot(t,z,'-')

grid;xlabel('t(sec)'),ylabel('x(t)'); subplot(122),plot(T,y,'-') xlabel('t(sec)'),ylabel('x(t)'); grid on ;

3 题目2

（1）操作过程

按照要求编写程序，将写好的程序输入matlab中得到所要求的结果。

(2)仿真结果

（3）源代码

t=[0:0.001:10];

yss=1,dta=0.02;

num=25;

den=[1,4,25];

G1=tf(num,den);

y1=step(G1,t);

r=1;while y1(r)

tr1=(r-1)*0.001;

[ymax,tp]=max(y1);tp1=(tp-1)*0.001;

mp1=(ymax-yss)/yss;

s=5001;while y1(s)>1-dta & y1(s)<1+dta;s=s-1;end

ts1=(s-1)*0.001;

[tr1 tp1 ts1 mp1]

4 题目3

(一)公式推导

（1）机械部分

由公式：

各级轴折算转动惯量： 222311123224

(

)()()n n n

J J J J n n n ∑=+++

工作台折算转动惯量： 2()2L J m π

= 得：

2

122211()2P J J J m i i π

=+

+ （2）电路部分

①电枢回路电压平衡方程：

()

()()a a a di t u t L

Ri t E dt

=++ （1）

式中()a E V 是电枢反电势，其大小与激磁磁通及转速成正，d ()(//)a E k w t V rad s =。

②电磁转矩方程：

()()m m a M t k i t =

（2）

③电动机上的转矩平衡方程：

()

()()m c dw t J

M t M t dt

=- （3）

J 是电动机和负载折合到电动机轴上的转动惯量。

由式（2-1）、式（2-2）和式（2-3）中消去中间变量()a i t 、a E 及()m M t 便可得到以()w t 为输出量，以()u t 为输入量的微分方程：

()()

()()22

22m d m d w t dw t LJ RJ k k w t k u t d t dt

++= （4） 对式（2-4）进行拉普拉斯变换，得：

22()2()()()m d m LJs W s RJsW s k k W s k U s ++=

所以

02()

()()22m m d

k W s G s U s LJs RJs k k =

=++ 而由：

()()W s s s θ=

故，

12

1

()22m m d k G s LJs RJs k k si

=

++ 式中i 为减速比。

又:

2()2p G s π

=

控制系统框图：

（二） （1）操作过程

由已知条件电感L 时，令参数K a =K c =K d =R=J=1，K m =10，P/i =4π，所以

系统的开环传递函数为：

310()(5)

b

k G s s s =

+

系统的闭环传递函数为：

210()510b

b

k G s s s k =

++

（2）仿真结果及分析

Ⅰ.系统的时域分析: