现代控制理论实验1资料.docx
- 文档编号:902596
- 上传时间:2022-10-13
- 格式:DOCX
- 页数:25
- 大小:320.82KB
现代控制理论实验1资料.docx
《现代控制理论实验1资料.docx》由会员分享,可在线阅读,更多相关《现代控制理论实验1资料.docx(25页珍藏版)》请在冰豆网上搜索。
现代控制理论实验1资料
河南工业大学《现代控制理论》实验报告
一、实验题目
线性系统状态空间表达式的建立以及线性变换
二、实验目的
1.掌握线性定常系统的状态空间表达式。
学会在MATLAB中建立状态空间模型的方法。
2.掌握传递函数与状态空间表达式之间相互转换的方法。
学会用MATLAB实现不同模型之间的相互转换。
3.熟悉系统的连接。
学会用MATLAB确定整个系统的状态空间表达式和传递函数。
4.掌握状态空间表达式的相似变换。
掌握将状态空间表达式转换为对角标准型、约当标准型、能控标准型和能观测标准型的方法。
学会用MATLAB进行线性变换。
三、实验过程及结果
1.已知系统的传递函数
(a)
传递函数:
状态空间模型:
传递函数:
零极点模型:
(b)
传递函数:
num=[168];den=[143];Gtf=tf(num,den)
零极点模型:
z=[-2-4];p=[-1-3];k=1;Gzpk=zpk(z,p,k)
状态空间模型:
Gss=ss(Gtf)
传递函数:
Gtf=tf(Gss)
零极点模型:
Gzpk=zpk(Gtf)
(c)
(1)建立系统的TF或ZPK模型。
TF模型:
num=[111];den=[16116];Gtf=tf(num,den)
ZPK模型:
Gzpk=zpk(Gtf)
(2)将给定传递函数用函数ss()转换为状态空间表达式。
再将得到的状态空间表达式用函数tf()转换为传递函数,并与原传递函数进行比较。
状态空间表达式:
Gss=ss(Gtf)
传递函数:
Gtf=tf(Gss)
该传递函数和原传递函数相同,没有变化。
(3)将给定传递函数用函数jordants()转换为对角标准型或约当标准型。
再将得到的对角标准型或约当标准型用函数tf()转换为传递函数,并与原传递函数进行比较。
对角标准型:
num=[111];den=[16116];Gj=jordants(num,den)
传递函数:
Gtf=tf(Gj)
(4)将给定传递函数用函数ctrlts()转换为能控标准型和能观测标准型。
再将得到的能控标准型和能观测标准型用函数tf()转换为传递函数,并与原传递函数进行比较。
能控标准型:
Gc=ctrlts(num,den)
传递函数:
Gtf=tf(Gc)
能观标准型:
Ao=(Gc.a)';Bo=(Gc.c)';Co=(Gc.b)';Do=Gc.d;Go=ss(Ao,Bo,Co,Do)
传递函数:
能控标准型与能观标准型转化为传递函数后与原函数相同,不会改变。
2.
已知系统的状态空间表达式
(a)
状态空间表达式:
A=[01;-5-6];B=[0;1];C=[11];D=0;Gss=ss(A,B,C,D)
对角标准型:
Gc1=canon(Gss,'modal')
(b)
状态空间表达式:
A=[010;302;-12-7-6];B=[2;1;7];C=[111];Gss=ss(A,B,C,0)
特征值:
eig(A)
对角标准型:
Gc=canon(Gss,'modal')
(d)
(1)建立给定系统的状态空间模型。
用函数eig()求出系统特征值。
用函数tf()和zpk()将这些状态空间表达式转换为传递函数,记录得到的传递函数和它的零极点。
比较系统的特征值和极点是否一致,为什么?
状态空间模型
A=[00-1;10-3;01-3];B=[1;1;0];C=[01-2];Gss=ss(A,B,C,0)
特征值:
eig(Gss)
传递函数:
Gtf=tf(Gss)
零极点模型:
Gzpk=zpk(Gss)
可知系统的极点是:
-1,-1,-1
系统的特征值和极点一致
(2)用函数canon()将给定状态空间表达式转换为对角标准型。
用函数eig()求出系统特征值。
比较这些特征值和
(1)中的特征值是否一致,为什么?
再用函数tf()和zpk()将对角标准型或约当标准型转换为传递函数。
比较这些传递函数和
(1)中的传递函数是否一致,为什么?
对角标准型:
Gc=canon(Gss,'modal')
特征值:
eig(Gc)
特征值和
(1)中的特征值相同
传递函数:
Gtf=tf(Gc)
零极点模型:
Gzpk=zpk(Gss)
可知,传递函数和
(1)中的传递函数一致,因为状态空间表达式不影响系统的输入输出特性,故传递函数不会因此发生改变。
(3)用函数ctrlts()将给定的状态空间表达式转换为能控标准型和能观测标准型。
用函数eig()求系统的特征值。
比较这些特征值和
(1)中的特征值是否一致,为什么?
再用函数tf()将它们转换为传递函数。
比较这些传递函数和
(1)中的传递函数是否一致,为什么?
能控标准型:
num=[121];den=[1331];Gc=ctrlts(num,den)
系统特征值:
eig(Gc)
能观标准型:
Ao=(Gc.a)';Bo=(Gc.c)';Co=(Gc.b)';Do=Gc.d;Go=ss(Ao,Bo,Co,Do)
系统特征值:
eig(Go)
这些特征值和
(1)中的特征值一致,因为非奇异变换不改变系统的原有性质,是等价变换。
3.已知两个子系统
(1)建立两个子系统的传递函数模型。
求它们串联、并联、反馈连接时,整个系统的传递函数模型。
然后将所得传递函数模型转换为状态空间模型。
第一个系统的传递函数:
num1=[11];den1=[144];Gtf1=tf(num1,den1)
第二个系统的传递函数:
num2=[35];den2=[16116];Gtf2=tf(num2,den2)
串联:
传递函数Gc=Gtf1*Gtf2
状态空间模型:
Gssc=ss(Gc)
并联:
传递函数Gb=Gtf1+Gtf2
状态空间模型:
Gssb=ss(Gb)
反馈传递函数:
Gf=feedback(Gtf1,Gtf2)
状态空间模型:
Gssf=ss(Gf)
(2)将两个子系统的传递函数模型转换为状态空间模型。
求它们串联、并联、反馈连接时,整个系统的状态空间模型。
然后将所得状态空间模型转换为传递函数模型。
比较
(1)和
(2)所得的相应的结果。
第一个子系统的状态空间模型:
Gss1=ss(Gtf1)
第一个子系统的状态空间模型:
Gss2=ss(Gtf2)
串联状态空间模型:
Gssc=Gss1*Gss2
传递函数:
Gtfc=tf(Gssc)
并联状态空间模型:
Gssb=Gss1+Gss2
传递函数:
Gtfb=tf(Gssb)
反馈状态空间模型:
Gssf=feedback(Gss1,Gss2)
传递函数:
Gtff=tf(Gssf)
(2)和
(1)的状态空间模型不同,但是传递函数相同
(3)
(2)中所得的整个系统的状态空间模型的系数矩阵与教材中推导出的整个系统的状态空间表达式的系数矩阵是符合的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 现代 控制 理论 实验 资料