ltisys和ltiss函数使用示例学习及频域分析(含实现程序)

本文介绍了MATLAB中用于描述线性时不变系统的ltisys函数和ltiss函数,详细阐述了时不变系统的一般形式,并通过实例展示了如何将二阶系统转化为状态空间方程。同时,讲解了如何使用splot函数绘制频率响应图,包括伯德图的绘制,并提供了具体的代码示例。最后,给出了离散系统的频率响应图绘制方法。
摘要由CSDN通过智能技术生成

目录

1.ltisys函数和ltiss函数

2.时不变系统一般形式

3.例如,给定一个二阶系统:

4.绘制频率响应图


1.ltisys函数和ltiss函数

ltisys函数的作用是描述系统,ltiss函数作用是从系统中提取状态空间矩阵,下面给出ltisys和ltiss函数的使用格式。

当被控对象为状态空间形式时:

sys=ltisys(A,B,C,D,E) %描述系统
[A,B,C,D] = ltiss(sys)%提取系统并以状态空间矩阵形式呈现

当被控对象为传递函数形式时:

sys = ltisys('tf',n,d)

2.时不变系统一般形式

时不变连续系统:

时不变离散系统:

  注:E为可逆矩阵,所以当E为单位矩阵时便是我们最熟悉的状态空间形式。

改写成系统矩阵的形式:

3.例如,给定一个二阶系统:

 其转化为上述标准形式的状态空间方程为:

 其中:

 所以使用ltisys的描述代码如下:

ltisys([0 1;-k -f],[0;1],[1 0],0,[1 0;0 m])

4.绘制频率响应图

使用函数及格式:

splot(sys,choice,xrange)         %连续系统
splot(sys,T,choice,xrange)       %离散系统
%{ 
①choice:
'bo'  :   Bode plot
'sv'  :   singular value plot
'ny'  :   Nyquist plot (SISO)
'li'  :   Lin-log Nyquist plot (SISO)
'ni'  :   Black/Nichols chart (SISO)
 
'st'  :   step response
'im'  :   impulse response
'sq'  :   response to a square wave
'si'  :   response to a sine wave
②xrange即坐标x(频率响应)的范围
%}

 取上例中的m = 2, f = 0.01, k = 0.5,以绘制伯德图为例,演示代码如下:

sys = ltisys([0 1;-0.5 -0.01],[0;1],[1 0],0,[1 0;0 2]);
splot(sys,'bo')

 出图结果:

 当限定频率范围在(10^-1~10^1)时,可改写为:

splot(sys,'bo',logspace(-1,1,50))

程序运行结果:

离散系统使用也是类似的,如给定以下离散系统:

取采样间隔为0.01,演示程序如下:

splot(ltisys(0.1,0.2,-1),0.01,'bo')

 出图结果:

注:仅为便利自己学习,错误在所难免,如有兴趣的学者可以参考交流,谢谢!

参考资料:

《LMI Control Toolbox》

MATLAB help文档

  • 2
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Mr. 邹

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值