北航计算机控制大作业,北航计算机控制系统大作业.docx
计算机控制系统
大作业
姓名:陈启航
学号:13031144
教师:周锐
日期:2016年6月1日
22
综合习题1
已知: ,
1) 试用 Z 变换、一阶向后差分、向前差分、零极点匹配、Tustin变换和预修正的Tustin(设关键频率=4)变换等方法将D(s)离散化,采样周期分别取为0.1s 和 0.4s;
2) 将 D(z)的零极点标在Z 平面图上
3) 计算D(jω)和各个D(ejωT )的幅频和相频特性并绘图,w由0~ 20rad,计算40 个点,应包括=4 点,每个T绘一张图(Z变换方法单画)
4) 计算 D(s)及T=0.1,T=0.4 时D(z)的单位脉冲响应,运行时间为4 秒
5) 结合所得的结果讨论分析各种离散化方法的特点
6) 写出报告,附上结果。
解:
(1) Z变换法:
a.离散化:
T=0.1s 时,
Dz=4zz-0.6703;
T=0.4s 时,
Dz=4zz-0.2019;
b.Dz的零极点
c. Djω和DejωT幅频相频特性曲线
连续系统:
T=0.1s时
T=0.4s时
d. Ds和Dz单位脉冲响应
Ds单位脉冲响应:
Dz单位脉冲响应:
T=0.1s时
T=0.4s时
(2) 各种离散化方法:
a.离散化后的D(z)
1、 一阶向后差分:
T=0.1s时
D(z)=0.2857zz-0.7143
T=0.4s时
D(z)=0.6154zz-0.3846
2、 一阶向前差分:
T=0.1s时
Dz=0.4z-0.6
T=0.4s时
Dz=1.6z+0.6
3、 零极点匹配
T=0.1s时
Dz=0.1648z+1z-0.6703
T=0.4s时
Dz=0.3991(z+1)z-0.2019
4、 Tustin变换
T=0.1s时
Dz=0.1667(z+1)z-0.6667
T=0.4s时
Dz=0.4444(z+1)z-0.1111
5、 预修正的Tustin变换(设关键频率=4)
T=0.1s时
Dz=0.1685(z+1)z-0.6629
T=0.4s时
Dz=0.5073(z+1)z+0.0146
b.Dz的零极点
1、一阶向后差分
2、一阶向前差分
3、零极点匹配
4、 Tustin变换
5、预修正的Tustin变换(设关键频率=4)
c. Djω和DejωT幅频相频特性曲线
1、一阶向后差分
T=0.1s时
T=0.4s时
2、一阶向前差分
T=0.1s时
T=0.4s时
3、零极点匹配
T=0.1s时
T=0.4s时
4、Tustin变换
T=0.1s时
T=0.4s时
5、 预修正的Tustin变换(设关键频率=4)
T=0.1s时
T=0.4s时
d. Ds和Dz单位脉冲响应
1、一阶向后差分
T=0.1s时
T=0.4s时
2、一阶向前差分
T=0.1s时
T=0.4s时
3、零极点匹配
T=0.1s时
T=0.4s时
4、Tustin变换
T=0.1s时
T=0.4s时
5、预修正的Tustin变换(设关键频率=4)
T=0.1s时
T=0.4s时
二、实验结果分析和总结:
在本题中,当采样周期T=0.4s时所有离散方法的都会出现频率混叠现象,使得采样信号失真。因为此采样周期不满足采样定理导致采样信号失真现象。
当满足采样定理时,各种离散化方法的特点如下:
① Z变换法
由Z变换的脉冲响应可看出,连续系统与离散后的系统的脉冲响应相同,故其可以应用在要求脉冲响应不变的场合。但是Z变换容易产生频率混叠,而且变换本身也十分麻烦,在多个环节进行串联时无法单独改变某一串联环节(即无串联性)。
②一阶向后差分法
稳定性保持不变,无混叠但畸变严重。在采样周期很小的情况下,所得的频率特性与原连续系统的频率特性比较接近但在采样周期较大时相差较大。
③一阶向前差分法
由于其映射关系,将原系统离散后可能不稳定。但在本题中,系统离散后依然稳定,所以用此方法也可以。在采样周期很小的情况下,所得的频率特性与原连续系统的频率特性比较接近但在采样周期较大时相差较大。
④零极点匹配法
当D(s)分子阶次低于分母时,D(z)中将带有(z+1)项,可以防止频率混叠,频率保持良好。但在使用单零点匹配时将带来相位的较大延迟,同时,零极点匹配法在使用中要求将稳态增益进行匹配,因此在使用中不是很方便。
⑤Tustin变换法
采样周期较小的情况下,在低频段范围内,s域与z域的频率近似保持线性,所以在此范围内Tustin变换频率失真较小。但随着频率的升高,Tustin变换频率范围变得很窄,高频失真严重。所以主要应用于低频环节离散化。
⑥修正Tustin变换法
由于此方法只能保证在指定频率点附近时的频率特性与原系统相近,但在其他频率点附近无法保证不失真,故其多用在要求在某些特征频率处离散前后频率特性保持不变的场合。
综合习题2
计算机伺服控制系统设计
1. 已知:
1)被控对象为一个带有均质圆盘负载的直流力矩电机,其伺服系统方框图如下:
2)D/A输出5v,电机转速ω=26rad/s
其中,电机传递函数为角速率ω/Δu;模拟控制其由KG1,KG2,KG3组成,数字控制器由采样、CPU(控制率)和D/A组成。给定参数如下:
电机传函Gs=θ(s)u(s)=kms(Tms+1),km=2 rad/s/V,Tm=0.1s
电机启动电压ΔuA=1.7v
测速机传递系数kω=1 v/rad/s
电位计最大转角为345,输出5v
功放KA=2=KG3
采样周期T=0.010s
2. 设计要求:
1) D/A输出120mv,电机启动:∆uA=1.7v
2) D/A输出5v,电机转速ω=26rad/s
3) 设计状态反馈增益,使系统闭环极点ξ≥0.9,ωn≥20rad/s
4) 设θ可测,设计降维观测器,观测器衰减速率是系统闭环衰减速率的4倍。
5) 求离散控制率D(z)
6) 将D(z)进行实现,配置适当的比例因子。
7) 编制程序框图,给出控制律的差分方程。
仿真验证调节器的控制效果。假设系统受到扰动,初试状态为:
初速
8)仿真验证调节器的控制效果。假设系统受到扰动,初试状态为:
初速度w=0,初始角度θ=10。看看是否经过一定时间后,系统状态
回到平衡的零态。
9) (选作)引进指令信号,设计相应的指令跟踪控制器,仿真给出闭
环系统的阶跃响应曲线。
解:
一)速度回路设计:
根据要求1),有
KG1∙KG3≥1.70.12=14.2≈15
取
KG1KG3=15
根据要求2),有
ωuDA=lims→0sKG1KG3kmTs+1+KG1KG2KG3Kωkm1s=KG1KG3km1+KG1KG2KG3Kωkm=265 rad/s /V
所以
KG2=0.159
二)状态反馈增益设计:
取ξ=0.9,ωn=20rad/s
则要求的闭环系统极点为
s1,2=-18j8.7178
映射到z平面中,离散域的极点为
z1,2=0.83210.0727i
离散域的特征方程为
z2-1.6642z+0.6977
对该系统,开环传递函数为:
Gs=300s2+105.4s
取转角θ=x1,角速度ω=x2为状态变量x=(x1,x2)T,系统状态方程为
x=Ax+Bu
y=Cx
其中
A=010-57.7
B=0300
C=[0 1]
则离散状态方程为
xn+1=Fxn+Gu(n)
y(n)=Cx(n)
F=eAT=L-1[(sI-A)-1]|t=T=10.007600.5616
G=0TeAtBdt=0.012492.2795
设状态反馈增益为:
k=k1k2
闭环系统特征方程为
zI-F+GK=z-1+0.01249k10.01249k2-0.00762.2795k1z-0.5616+2.2795k2=z2+0.01249k1+2.2795k2-1.5616z+0.0103k1-2.2795k2+0.5616
对应系数相等,得
k1=1.4699k2=-0.0531
k=[1.4699-0.0531]
闭环系统的状态方程为:
xk+1=0.98160.0083-3.35060.6826xk+0.012492.2795u(k)
yk=10x(k)
三)设计降维观测器:
降维观测器方程为
x2k+1=F22-LF12x2k+F21-LF11yk+G2-LG1uk+Ly(k+1)
其中,
F22=0.6826,F12=0.0083,F21=-3.3506,F11=0.9816,G2=2.2795,G1=0.01249
闭环系统极点:
s1,2=-18j8.7178
观测器极点:
σL=-184=-72
闭环极点:
z=eσLT=0.4868
则有
F22-LF12=0.4868
L=23.59
观测器方程:
x2k+1=0.4868x2k-26.5065yk+1.9839uk+23.59y(k+1)
设计速度回路:rk=0,所以uk=-kx(k)则
zθz-0.4868θz-0.1053θz=-26.5065θz-2.9161θz+23.59zθz
整理得
θzθz=23.59z-1.2472z-0.5921
四)离散控制律设计:
状态反馈输入:
uk=-kxk=-1.4699θk+0.0531θk
则有
Uz=-1.4699θz+0.0531θzθzθz=23.59z-1.2472z-0.5921
计算控制律D(z)得
Dz=U(z)θz=-0.2173z+3.1841z-0.5921
-0.6919
-0.5921
D/A
A/D
-0.2173
z-1
θt
ut
五)控制律的实现:
采用直接零极型编排:
配置比例因子:
稳态增益:
Dz|z=1=2.229
高频增益:
Dz|z=-1=0.093
A/D,D/A增益补偿取1,比例因子取22。
配置比例因子后结构图如下:
14
-14
-0.17298
-0.5921
D/A
A/D
-0.054325
z-1
θt
ut
4
六)控制算法实现:
算法1:
u1k=-0.054325ek+x(k-1)
u2k=14u1(k)-1 u1k14u1(k)14u1k-14
算法2:
xk=-0.17298ek-0.5921u1k
算法流程图如下:
A/D
X=0
u1=-0.054325e+x
u2=1
u2=-1
u2=4u1
D/A
x=-0.17298e+0.5921u1
u20.25?
u2-0.25?
是
否
是
否
展开阅读全文
北航计算机控制大作业,北航计算机控制系统大作业.docx相关推荐
- 201609计算机控制技术作业三,计算机控制技术大作业2015..doc
计算机控制技术大作业2015. 深圳大学考试答题纸 (以论文.报告等形式考核专用)二○一四 -二○一五 学年度第 2 学期 课程编号1700470001课程名称计算机控制技术主讲教师评分学 号姓名专业 ...
- 计算机控制系统第2次作业,201609计算机控制系统作业2
201609计算机控制系统作业2 单项选择题 第1题 在样条函数中,当加速度是常数直线段的时候,对应的位移是___. A.常数 B.直线 C.抛物线 D.三次样条 答案:C 第2题 步进电机细分的作用 ...
- 用计算机控制食品生产论文,计算机控制系统论文范文2篇
计算机控制是指用数字计算机对动态系统进行的控制,在生产中得到了广泛的应用.下面是秋天网小编为大家整理的计算机控制系统论文,供大家参考. 计算机控制系统论文范文一:计算机控制技术课程教学改革论文 1课程 ...
- 计算机控制系统设计题例题,计算机控制系统练习题..doc
<计算机控制系统> 学习中心: 专 业: 学 号: 姓 名: 作业三 第五章 5.1 已知连续系统的传递函数: (1) ; (2) 试用冲击不变法﹑阶跃不变法﹑零极点匹配法﹑双线性变换法. ...
- 计算机控制z变换例题,计算机控制系统总习题(及答案).ppt
<计算机控制系统总习题(及答案).ppt>由会员分享,可在线阅读,更多相关<计算机控制系统总习题(及答案).ppt(36页珍藏版)>请在人人文库网上搜索. 1.1,8255,微 ...
- 计算机控制系统第2次作业,《计算机控制系统》第二次作业答案
- <计算机控制系统>第二次作业答案 你的得分: 100.0 完成日期:2018年09月10日 01点51分 说明: 每道小题选项旁的标识是标准答案. 一.单项选择题.本大题共16个小题, ...
- 计算机控制实时的概念,计算机控制系统复习资料..doc
计算机控制系统复习资料. 计算机控制系统复习 第一章·绪论 一.计算机控制系统的基本结构,控制过程,实时概念? 基本结构: 控制过程:1数据采集及处理 2.实时控制和信息处理 实施概念:信号的输入计算 ...
- 计算机仿真作业题目,计算机仿真技术作业一
<计算机仿真技术作业一>由会员分享,可在线阅读,更多相关<计算机仿真技术作业一(14页珍藏版)>请在人人文库网上搜索. 1.11292058 杨苑 电气1110计算机仿真技术作 ...
- 计算机仿真作业三,计算机仿真技术作业三.doc
计算机仿真技术作业三.doc 计算机仿真技术作业三 题目:三相桥式全控整流电路仿真 姓名: 班级: 学号: 计算机仿真技术作业三 题目:三相桥式全控整流电路仿真 利用simpowersystems建立 ...
最新文章
- Visual C++ Windows 用来定位 DLL 的搜索路径
- JS代码的window.location属性详解
- linux写入二进制文件内容,linux – 从管道读取数据并写入标准输出,中间延迟.必须处理二进制文件...
- 1008 数组元素循环右移问题 (20分)
- 排序算法之快速排序(Java)
- linux net 伪指令码,1.2.2. ldr pWTCON
- 批量画同心不同半径圆lisp_【微课视频】青岛版数学六年级上册5.1圆的认识
- adb shell 运行时报错“adb server version (26) doesn‘t match this client (39); killing...“的解决方案
- IPTV软件提取后无法安装:显示安装失败,需求的共享用户签名错误!请大神帮助,谢谢
- python身份证号处理代码_Python实现身份证号码解析
- Cascading介绍
- maven中无法下载ojdbc
- css自适应图片样式,css怎么让图片自适应?css图片自适应大小的方法介绍
- 19 个接私活平台,有技术就有钱
- 大数据学习路线,芜湖 起飞~
- 麒麟系统开启关闭防火墙服务
- 咖说 | 揭秘佳士得首次拍卖的区块链艺术品Portraits of a Mind
- Java -mail 退信代码说明
- 使用openssl 生成RSA、SM2、ECC的P12证书的方法
- 数学建模 趣味数模问题(数模入门)