试验二相平面作图文档格式.docx
- 文档编号:6500249
- 上传时间:2023-05-06
- 格式:DOCX
- 页数:11
- 大小:191.54KB
试验二相平面作图文档格式.docx
《试验二相平面作图文档格式.docx》由会员分享,可在线阅读,更多相关《试验二相平面作图文档格式.docx(11页珍藏版)》请在冰点文库上搜索。
连续系统命令:
(大部分s—>
z仍可用)
Gs=tf(num,den);
zpk(z,p,k);
zp2tftf2zpprintsys(Gs)
parallel,series,feedback,cloop
step(sysc),impulse,lsim
pzmap,rlocus,sgridbode,nyguistmargin
模型转换:
连续系统—离散系统
c2d(sysc,ts,’zoh/imp’),d2c
离散系统特有命令:
t=0:
ts:
tend;
step(sysd,t)
dstep(numd,dend,[LongNumber]),dimpulse,dlsim
pzmap,rlocus,dgriddbode,dnyguistmargin
其他有关命令:
helpdoc
roots(dend),
residueresiduez部分分式展开
conv(num1,num2)多项式乘法
deconv(numz,denz)多项式除法
zplane离散系统零极点图
figureplot(t,yc,’g’,t,yd,’r’)stem(yd)
实验二相平面作图
1、利用MATLAB完成控制系统的相平面图;
2、了解二阶系统相平面图的一般规律;
3、利用相平面图进行系统分析。
二、相平面作图实验方法
1、相关MATLAB命令
plot(t,x),plot(t,x,’:
g’,t,y,’-r’),plot(x,y)
subplot(n,m,i)
2、相平面作图方法
方法一:
逐步求解并作图
例如:
对二阶系统
(1)先求系统的
;
(2)再根据
判断系统是欠阻尼、临界阻尼还是过阻尼;
(3)然后求出系统的时间响应函数;
(4)再求出其导数。
num=[10];
den=[1210];
0.01:
6;
x=1-sqrt(9/10).*exp(-t).*sin(sqrt(11).*t)%时间响应函数
y=sqrt(9/10).*exp(-t).*sin(sqrt(11).*t)-sqrt(9/10).*exp(-t).*sqrt(11).*cos(sqrt(11).*t)%时间响应函数的导数
plot(x,y);
grid;
方法二:
利用状态空间模型的STEP求各阶导数
[a,b,c,d]=tf2ss(num,den);
sys=ss(a,b,c,d);
[y,t,x]=step(sys);
当sys是状态空间模型时,step可以计算得到输出向量y、时间向量t及状态向量x(n个状态,x(:
1)…x(:
n)分别与y、y一阶导数…y的n-1阶导数有关)
例:
二阶系统G(s)=10/(s2+2s+10)输入信号为r(t)=1(t),作该系统的相平面图。
[a,b,c,d]=tf2ss(num,den);
plot(y);
holdon;
plot(x(:
1));
2));
grid
subplot(2,2,1);
plot(t,x(:
subplot(2,2,2);
subplot(2,2,3);
2),x(:
subplot(2,2,4);
plot(-x(:
2)+0.1,-x(:
方法三:
利用Simulink建立系统模型
参见ver4p360.mdl及教材p360二阶系统分析
>
holdoff
plot(xx,yy);
axis([-1010-1010]);
%调整坐标
grid;
%画格
title('
a=-0.2,b=2,x0=0.1,y0=-0.1'
);
%加标题
holdon;
%在Simulink文件中改参数运行
%holdon后再plot,在原图上增加曲线
三、实验内容
对开环传递函数为G(s)=K/[(s+a)(s+b)]的单位负反馈系统,给出不同的K、a和b参数,使系统处于下列不同平衡点,用MATLAB画出系统误差的相平面图和单位阶跃响应曲线。
(1)稳定节点;
(2)稳定焦点;
(3)不稳定节点;
(4)不稳定焦点;
(5)中心点;
(6)鞍点。
四、实验报告要求
1、记录给定系统在不同参数下的相平面图;
2、分析不同参数下平衡点的性质及起系统动态响应的影响。
实验三用Matlab进行状态空间分析及设计
一、实验目的:
掌握使用MATLAB进行及状态空间分析及状态反馈控制系统的设计
二、主要函数及命令
(一)、状态空间描述及其转换
1.状态空间表达及显示
Ø
ss(A,B,C,D)%显示ABCD构成的状态空间模型
sys1=ss(A,B,C,D)%将ABCD状态空间模型赋给结构变量sys1
[A2,B2,C2,D2]=ssdata(sys1)%模型sys1中的矩阵赋给矩阵变量
printsys(A,B,C,D)%显示ABCD构成的状态空间模型
sys1%显示结构变量sys1即状态空间模型
2.模型转换
[num,den]=ss2tf(A,B,C,D)Ø
sys_ss=ss(A,D,C,D)
[A,B,C,D]=tf2ss(num,den)Ø
sys_tf=tf(num,den)
sys_ss2=ss2ss(sys_ss1,T)Ø
[A2,B2,C2,D2]=ssdata(ss2)
3.标准形式
sys_ss2=canon(sys_ss1,,modal,)%对角线,复数
sys_ss3=canon(sys_ss1,,companion,)%A为伴随矩阵
4.最小实观
[Am,Bm,Cm,Dm]=minreal(A,B,C,D)
sys_ssm=mineral(sys_ss)Ø
sys_tfm=minreal(sys_tf)
(二)、求A的特征值,特征向量
[W,Ad]=eig(A)Ø
eig(A)Ø
eig(sys_ss)Ø
det(A)
(三)、状态方程的解及状态转移矩阵
expm(A*t)%t必须已知如:
t=[0:
0.1:
5]
[y,x]=lsim(A,B,C,D,u,t,x0)%例u=0*t,oru=u+1,x0=[0;
0;
0]
plot(x)%另plot(x
(1)),plot(y)等
step(sys_ss)
(四)、可控可观性分析
Qc=ctrb(sys_ss)Ø
Qo=obsv(sys_ss)
Qc=ctrb(A,B)Ø
Qo=obsv(A,C)
n=length(A)Ø
nc=rank(Qc)
no=rank(Qo)Ø
no=rank(obsv(ss(a,b,c,d)))
[A1,B1,C1,T,K]=ctrbf(A,B,C)Ø
[A1,B1,C1,T,K]=obsv(A,B,C)
(五)状态反馈控制器的设计(极点配置法举例)
a=[-2–2.5–0.5;
100;
010]>
b=[1;
P=[-1,-2,-3]%希望配置的闭环极点
Ø
k=place(a,b,p)%求状态反馈矩阵
aclose=a-b*k%求状态反馈控制系统闭环状态矩阵
eig(a)%求开环状态矩阵特征值(开环极点)
eig(aclose)%求闭环状态矩阵特征值(闭环环极点)
三、实验及报告内容
1,系统状态空间模型如下:
;
(1)求其传递函数,由传递函数求系统的极点;
(2)由上述状态空间模型,求系统的特征值;
(3)求上述系统状态转移矩阵;
(4)求其在x0=[2;
1;
2],u为单位阶跃输入时x及y的响应;
(5)分析上述系统的可控性、可观性;
(6)将上述状态空间模型转换为其他标准形式;
(7)取T=[124;
010;
001]对上述状态空间模型进行变换,分析变换后的系统。
2,分析下列系统的可控性、可观性
(1)
(2)
(3)
3,系统状态空间模型如下,
(1)判别系统的可控性;
(2)设计状态反馈控制器使闭环极点为p=[-1,-10,-12];
(3)求出闭环系统的传递函数和动态方程;
(4)比较反馈前后系统的阶跃响应。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 试验 平面 作图