求解系统的状态方程教学文案Word文档下载推荐.docx
《求解系统的状态方程教学文案Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《求解系统的状态方程教学文案Word文档下载推荐.docx（21页珍藏版）》请在冰豆网上搜索。
（1）求下列系统矩阵A对应的状态转移矩阵
（a）
（b）
代码：
symslambda
A=[lambda00;
0lambda0;
00lambda];
symst;
f=expm（A*t）
（c）
symst;
symslambda;
A=[lambda000;
0lambda10;
00lambda1;
000lambda];
（2）已知系统
a）用MATLAB求状态方程的解析解。
选择时间向量t，绘制系统的状态响应曲线。
观察并记录这些曲线。
（1）
A=[01;
-2-3];
B=[3;
0];
C=[11];
D=[0];
u=1;
f=expm（A*t）;
%状态转移矩阵
x0=0;
s1=f*B*u;
s2=int（s1,t,0,t）%状态方程解析解
状态曲线：
（2）A=[01;
-2-3];
f=expm（A*t）;
X0=[1;
t=[0:
0.5:
10];
fori=1:
length（t）;
g（i）=double（subs（f
（1）,t（i）））;
end
plot（t,g）
（3）状态转移矩阵
symslambda
A=[lambda00;
symst
b）计算系统在初始状态作用下状态响应和输出响应的数值解（用函数initial（））,绘制系统的状态响应曲线和输出响应曲线。
观察并记录这些响应曲线，然后将这一状态响应曲线与a）中状态响应曲线进行比较。
G=ss（A,B,C,D）;
t=[0:
x0=[1;
0]
[y0,t,x0]=initial（G,x0,t）;
plot（t,x0,'
-'
t,y0,'
）
c）根据b）中所得的状态响应的数值解，绘制系统的状态轨迹（用命令plot（x（:
1）,x（:
2）））。
记录系统状态转移的过程，结合a）和b）中的状态响应曲线分析这一过程。
0.01:
G=ss（A,B,C,D）
[y,t,x]=initial（G,x0,t）;
plot（x（:
1）,x（:
2））
2）令初始状态为零，输入为u（t）=1（t）.
[y,t,x]=step（G）;
plot（t,x）
b）计算系统在初始状态作用下状态响应和输出响应的数值解,绘制系统的状态响应曲线和输出响应曲线。
观察并记录这些响应曲线，然后将这一状态响应曲线与a）.中状态响应曲线进行比较。
-1];
c）根据b）中所得的状态响应的数值解，绘制系统的状态轨迹。
x0=[00];
绘制系统的状态响应曲线、输出响应曲线和状态轨迹。
观察和分析这些响应曲线和状态轨迹是否是
（1）和
（2）中的响应曲线和状态轨迹的叠加。
x0=[1-1];
4）令初始状态为零，输入为u（t）=3sin（5t）。
计算状态响应和输出响应的数值解（用函数
lsim（）），并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。
u=3*sin（5*t）;
[y,t,x]=lsim（G,u,t）;
（3）已知系统
1）当输入为u（t）=
（t）时，用函数initial（）和impulse（）求解系统的状态响应和输出响应的
数值解，并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。
状态响应：
A=[0,1,0;
0,0,1;
-6,-11,-6];
B=[0;
0;
1];
C=[6,0,0];
D=0;
x0=[1,0,-1];
u=ones（size（t））;
plot（t,x,t,y）
输出响应：
G=ss（A,B,C,D）;
[y,t,x]=lsim（G,u,t）;
plot（t,x）
2）当输入为u（t）=1（t）时，用函数initial（）和step（）求解系统的状态响应和输出响应的数
值解，并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。
状态响应
u=step（G）;
3）当输入为u（t）=t时，用函数initial（）和lsim（）求解系统的状态响应和输出响应的数值
解，并绘制系统的状态响应曲线、输出响应曲线和状态轨迹。
u=t;
输出响应
4）当输入为
时，用函数initial（）和lsim（）求解系统的状态响应和输出响应的数值解，并绘制系统的状态响应曲线、输出响应曲线和状态轨迹.
u=sin（t）;
（4）已知一个连续系统的状态方程是
若取采样周期秒0.05T=
1）试求相应的离散化状态空间模型；
代码:
symsT;
-25-4];
[Gz,Hz]=c2d（A,B,T）
G=
[exp（-2*T-21^（1/2）*T*i）/2+exp（-2*T+21^（1/2）*T*i）/2+（21^（1/2）*exp（-2*T-21^（1/2）*T*i）*i）/21-（21^（1/2）*exp（-2*T+21^（1/2）*T*i）*i）/21,（21^（1/2）*exp（-2*T-21^（1/2）*T*i）*i）/42-（21^（1/2）*exp（-2*T+21^（1/2）*T*i）*i）/42]
[-（21^（1/2）*exp（-2*T-21^（1/2）*T*i）*25*i）/42+（21^（1/2）*exp（-2*T+21^（1/2）*T*i）*25*i）/42,exp（-2*T-21^（1/2）*T*i）/2+exp（-2*T+21^（1/2）*T*i）/2-（21^（1/2）*exp（-2*T-21^（1/2）*T*i）*i）/21+（21^（1/2）*exp（-2*T+21^（1/2）*T*i）*i）/21]
H=
1/25-exp（-2*T+21^（1/2）*T*i）/50-（21^（1/2）*exp（-2*T-21^（1/2）*T*i）*i）/525+（21^（1/2）*exp（-2*T+21^（1/2）*T*i）*i）/525-exp（-2*T-21^（1/2）*T*i）/50
（21^（1/2）*exp（-2*T-21^（1/2）*T*i）*i）/42-（21^（1/2）*exp（-2*T+21^（1/2）*T*i）*i）/42
2）分析不同采样周期下，离散化状态空间模型的结果。
[Gz,Hz]=c2d（A,B,0.05）
Gz=
0.97090.0448
-1.12120.7915
Hz=
0.0012
0.0448
5、实验心得：
本次实验是求解系统的状态方程。
并且要通过编入代码、调试，掌握传递函数与状态空间表达式之间相互转换的方法，系统状态空间表达式与传递函数相互转换方法。
同时也掌握了状态空间表达式的相似变换和如何将状态空间表达式转换为对角标准型、约当标准型、能控标准型和能观测标准型的方法。
通过本次试验我初步学会了使用matlab解决基本的线带控制问题，对书本上的知识有了根深的理解。
升级会员 