现控实验报告.docx
- 文档编号:26142668
- 上传时间:2023-06-17
- 格式:DOCX
- 页数:8
- 大小:92.28KB
现控实验报告.docx
《现控实验报告.docx》由会员分享,可在线阅读,更多相关《现控实验报告.docx(8页珍藏版)》请在冰豆网上搜索。
现控实验报告
实验报告
实验项目名称:
线性系统的数学模型转换同组人
实验时间实验室指导教师
一、实验目的
1、学习系统状态空间模型的建立方法;
2、通过编程、上机调试,掌握系统状态空间模型之间转换的方法。
二、实验主要仪器与设备
1、PC计算机1台;
2、MATLAB6.X或MATLAB7.X软件1套。
三、实验原理
1、假设系统是单输入单输出系统(简称SISO),其输入、输出分别用u(t)、y(t)来表示,则线性系统的传递函数模型:
在MATLAB语言中,可以利用传递函数分子、分母多项式的系数向量进行描述,分子num、分母den多项式的系数向量分别为:
num=[bm,bm-1,…,b0],den=[1,an-1,…,a0]
这里分子、分母多项式系数按s的降幂排列。
2、假设系统是单输入单输出系统(简称SISO),其零极点模型为:
其中zi(i=1,2,…,m)和pi(i=1,2,…,n)分别为系统的零点和极点,K为系统增益。
[z]、[p]、[k]分别为系统的零极点和增益向量。
3、设系统的状态空间模型为:
在MATLAB中,系统状态空间用(A,B,C,D)矩阵组表示,MATLAB提供了建立状态空间模型的函数ss(),其常见的调用格式为:
sys=ss(A,B,C,D)。
4、将传递函数模型转换为状态空间模型的函数为tf2ss(),如[A,B,C,D]=tf2ss(num,den)。
5、将零极点模型转换为状态空间模型的函数为zp2ss(),如[A,B,C,D]=zp2ss(z,p,k)。
6、利用MATLAB提供的函数canon()可以将系统直接转化为对角型,其常用调用格式为:
[As,Bs,Cs,Ds,Ts]=canon(A,B,C,D,’mod’)
其中,A、B、C、D是变换前系统的状态空间实现,参烽’mod’表示转化成对角型,As、Bs、Cs、Ds是变换后的对角型,Ts表示所作的线性变换。
7、利用MATLAB提供的函数ss2ss()可以进行状态空间表达式的线性变换,其常用调用格式为:
[A1,B1,C1,D1]=ss2ss(A,B,C,D,inv(T))
四、预习要求
做实验前必须认真复习MATLAB仿真课程中相关的知识以及现代控制理论课本中关于状态空间表达式的建立的内容。
五、实验内容及步骤
1、利用[A,B,C,D]=tf2ss(num,den)将传递函数模型转换为状态空间模型。
2、利用[A,B,C,D]=zp2ss(z,p,k)将零极点模型转换为状态空间模型。
3、利用sys=ss(A,B,C,D)直接得出系统的状态空间模型。
4、利用[As,Bs,Cs,Ds,Ts]=canon(A,B,C,D,’mod’)将系统状态空间模型转化为对角型。
5、利用[A1,B1,C1,D1]=ss2ss(A,B,C,D,inv(T))实现系统状态空间表达式的线性变换。
六、思考题
1、变换矩阵T如何计算。
2、如何将系统的状态空间模型转换为传递函数模型?
3、如何将系统的状态空间模型转换为零极点模型?
4、系统的传递函数模型与零极点模型又是如何转换的?
七、实验报告
1、根椐已知的系统传递函数(自拟或选择现代控制理论课本中例题或习题),将传递函数模型转换为状态空间模型。
写出设计的MATLAB程序和运行结果。
2、根椐已知的系统零极点模型(自拟或选择现代控制理论课本中例题或习题),将零极点模型转换为状态空间模型。
写出设计的MATLAB程序和运行结果。
3、根据已知的系统矩阵A、控制矩阵B、输出矩阵C和直接转移矩阵D(自拟或选择现代控制理论课本中例题或习题)直接得出系统的状态空间模型。
写出设计的MATLAB程序和运行结果。
4、根据给定的系统状态空间表达式(自拟或选择现代控制理论课本中例题或习题),写出直接转化对角型和线性变换的程序代码和运行结果。
指导教师批阅及成绩
指导教师签名:
年月日
实验报告
实验项目名称:
能控能观性及其标准同组人
实验时间实验室指导教师
一、实验目的
1、掌握系统能控能观性判别方法及能控能观性标准型;
2、掌握利用MATLAB实现系统能控能观性的判别及能控、能观性标准型的转化。
二、实验主要仪器与设备
1、PC计算机1台;
2、MATLAB6.X或MATLAB7.X软件1套。
三、实验原理
1、能控性定义及其判别
设单输入n阶线性定常连续系统为:
。
若存在一个分段连续的控制函数u(t),能在有限的时间段[t0,tf]内把将系统从t0时刻的初始状态X(t0)转移到任意指定的终态X(tf),则称系统在t0时刻的状态X(t0)是能控的;如果系统的每一个状态都是能控的,则称系统是状态完全能控的;反之,只要有一个状态不能控,就称系统是不能控的。
n阶线性定常连续系统能控的充分必要条件为能控判别阵M=[BAB…An-1B]的秩是满秩的,即rank[M]=n。
2、能观性定义及其判别
若对任意给定的输入u(t),总能在有限的时间段[t0,tf]内,根据系统的输入u(t)及系统观测的y(t),能唯一地确定t0时刻的每一状态X(t0),则称系统在t0时刻是状态能观的。
n阶线性定常连续系统能观的充分必要条件为能观判别阵N=[CCA…CAn-1]T的秩是满秩的,即rank[N]=n。
3、能控标准型
控制系统的能控标准型有两种形式,分别称之为能控I型和能控II型。
对于能控I型
,其各矩阵形式为:
对于能控II型
,其各矩阵形式为:
若线性定常单输入系统:
是能控的,则可以通过非奇异矩阵
的线性变换,将系统
变换成能控I型
。
其中ai(i=0,1,2,…,n-1)是系统特征多项式的系数。
或通过非奇异矩阵
将系统
变换成能控II型
。
4、能观标准型
控制系统的能观标准型有两种形式,分别称之为能观I型和能观II型。
对于能观I型
,其各矩阵形式为:
对于能观II型
,其各矩阵形式为:
若线性定常单输入系统:
是能观的,则可以通过非奇异矩阵
的线性变换,将系统
变换成能观I型
。
或通过非奇异矩阵
将系统
变换成能观II型
。
其中ai(i=0,1,2,…,n-1)是系统特征多项式的系数。
四、预习要求
做实验前必须认真复习MATLAB仿真课程中相关的知识以及现代控制理论课本中关于系统的能控性、能观性判别、能控能观标准型以及能控性分解和能观性分解等的内容。
五、实验内容及步骤
1、利用M=ctrb(A,B),求系统的能控判别矩阵M=[B,AB,A2B,…],结合求M秩的函数rank(M),从而判别系统的能控性。
2、利用N=obsv(A,C)求系统的能观判别矩阵N=[C,CA,CA2,…]T,结合求N秩的函数rank(N),从而判别系统的能观性。
3、当系统能控性矩阵的秩小于系统的维数n时,利用[Ac,Bc,Cc]=ctrbf(A,B,C)对线性系统进行能控性分解。
4、当系统能控性矩阵的秩小于系统的维数n时,利用[Ao,Bo,Co]=obsvf(A,B,C)对线性系统进行能观性分解。
5、利用[A1,B1,C1,D1]=ss2ss(A,B,C,D,inv(T))实现系统状态空间表达式的线性变换。
六、思考题
1、Cc1中的βi与Cc2中的βi是不是同数值?
2、Bo1中的βi与Bo2中的βi是不是同数值?
七、实验报告
1、已知系统
的相应系统矩阵为
,
,
,试判断系统是否能控?
是否能观?
写出设计的MATLAB程序和运行结果。
2、已知系统
的相应系统矩阵为
,
,
,D=0,试判断它的能控性,如完全能控,将其转化为能控II型。
写出设计的MATLAB程序和运行结果。
3、已知系统
的相应系统矩阵为
,
,
,D=0,试判断它的能观性,如完全能观,将其转化为能观I型。
写出设计的MATLAB程序和运行结果。
4、已知系统
的相应系统矩阵为
,
,
,试对系统进行能控性结构分解和能观性结构分解写出设计的MATLAB程序和运行结果。
指导教师批阅及成绩
指导教师签名:
年月日
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 报告