Matlab中lsim函数使用

Matlab中lsim函数使用lsim函数:lsim函数是针对线性时不变模型,给定任意输入,得到任意输出。lsim函数表示任意输入函数的响应,连续系统对任意输入函数的响应可以利用lsim函数求取。语法(常用):1.分子分母形式lsim(num,den,u,t)2.传递函数形式lsim(sys,u,t)3.状态空间形式lsim(A,B,C,D,u,t)其中,u为由给定输入序列构成的矩阵,它的每列对应一个输入,每行对应一个新的时间点,其行数与时间t的长度相等,其它的用法与step函数相同。…

大家好,又见面了,我是你们的朋友全栈君。

目录

语法(常用)

说明

示例


        lsim函数:lsim函数是针对线性时不变模型,给定任意输入,得到任意输出。lsim函数表示任意输入函数的响应,连续系统对任意输入函数的响应可以利用lsim函数求取。

语法(常用)

分子分母形式

lsim(num,den,u,t)

传递函数形式

lsim(sys,u,t)

状态空间形式

lsim(A,B,C,D,u,t)

其它形式

lsim(sys,u,t,x0)

lsim(sys,u,t,x0,method)

y=lsim(sys,u,t,x0)

y=lsim(sys,u,t,x0,method)

        其中,u为由给定输入序列构成的矩阵,它的每列对应一个输入,每行对应一个新的时间点,其行数与时间t的长度相等,其它的用法与step函数相同。

说明

响应图

        lsim(sys,u,t)绘制动态系统模型sys对输入历史记录(t,u)的模拟时间响应。 向量t指定用于仿真的时间样本。 对于单输入系统,输入信号u是与t长度相同的向量。 对于多输入系统,u是一个数组,其行数与时间样本(length(t))一样多,而列数与sys的输入一样多。

        lsim(sys,u,t,x0)当sys是状态空间模型时,进一步指定初始状态值的向量x0。

        lsim(sys,u,t,x0,method)当sys是连续时间模型时,如何在样本之间插入输入值,method即插入输入值的方法。

响应数据

        y=lsim(sys,u,t)返回与输入y同时采样 的系统响应t。对于单输出系统,y是与长度相同的向量t。对于多输出系统, y是一个数组,其中的行与时间样本(length(t))一样多,列与的输出一样多sys。该语法不会生成图。

        y=lsim(sys,u,t,x0)当sys是状态空间模型时,进一步指定初始状态值的向量x0。

        y=lsim(sys,u,t,x0,method)当sys是连续时间模型时,如何在样本之间插入输入值,method即插入输入值的方法。

示例

示例1

%------------------------------------
        H(s)=[
                2s^2 + 5s + 1
               ---------------
                s^2 + 2s + 3
                     s - 1
                 -----------
                 s^2 + s + 5
              ]
%------------------------------------
%%
clc;
clear off;
H = [tf([2 5 1],[1 2 3]);tf([1 -1],[1 1 5])];
[u,t] = gensig('square',4,10,0.1);
lsim(H,u,t);

        如图所示:

Matlab中lsim函数使用

        当控制系统为状态空间表示时,则:

A = [-3 -1.5; 5 0];
B = [1; 0];
C = [0.5 1.5];
D = 0;
sys = ss(A,B,C,D);

[u,t] = gensig("square",10,20);
lsim(sys,u,t)
grid on

        采用方波输入,响应如图所示: 

Matlab中lsim函数使用

示例2

        如下所示的传递函数:

sys = tf(3,[1 2 3])

sys =
 
        3
  -------------
  s^2 + 2 s + 3

        它是一个连续的系统,则:

sys = tf(3,[1 2 3]);
t = 0:0.04:8;  % 201 points
u = max(0,min(t-1,1));
lsim(sys,u,t)
grid on

Matlab中lsim函数使用

        用lsim函数获取模拟的响应数据为:

y = lsim(sys,u,t);
size(y)
ans = 1×2

   201     1

        向量y包含t相应时间的模拟响应。

示例3

        当控制系统为离散型时,且该传递函数的采样时间为0.05 s。使用相同的采样时间来生成时间矢量t和阶跃信号u。如下所示:

sys = tf([0.06 0.05],[1 -1.56 0.67],0.05);
t = 0:0.05:4;  
u = max(0,min(t-1,1));
lsim(sys,u,t)

        如图所示:

Matlab中lsim函数使用

        当输入为正弦波时,即模拟系统对周期为1 s,持续时间为4 s的正弦波的响应。,则:

sys = tf([0.06 0.05],[1 -1.56 0.67],0.05);
t = 0:0.05:4;  
[u,t] = gensig("sine",1,4,0.05);
lsim(sys,u,t)

        其响应如图所示:

Matlab中lsim函数使用

示例4

        lsim允许在同一轴上绘制多个动态系统的模拟响应。 例如,比较带有PI控制器和PID控制器的系统的闭环响应。 创建系统的传递函数并调整控制器。如下所示:

H = tf(4,[1 10 25]);
C1 = pidtune(H,'PI');
C2 = pidtune(H,'PID');

%形成闭环系统
sys1 = feedback(H*C1,1);
sys2 = feedback(H*C2,1);

%绘制两个系统对周期为4 s的方波的响应。
[u,t] = gensig("square",4,12);
lsim(sys1,sys2,u,t)
grid on
legend("PI","PID")

        如图所示:

Matlab中lsim函数使用

        默认情况下,lsim绘制的每个系统选择不同的颜色。 可以使用LineSpec输入参数指定颜色和线条样式。

H = tf(4,[1 10 25]);
C1 = pidtune(H,'PI');
C2 = pidtune(H,'PID');

%形成闭环系统
sys1 = feedback(H*C1,1);
sys2 = feedback(H*C2,1);

%绘制两个系统对周期为4 s的方波的响应。
[u,t] = gensig("square",4,12);
lsim(sys1,"r--",sys2,"b",u,t)
 grid on
 legend("PI","PID")

        如图所示:

Matlab中lsim函数使用

示例5

        默认情况下,lsim假设在模拟开始时所有状态均为零,从而模拟模型。 模拟状态空间模型的响应时,当使用可选的x0输入参数指定非零初始状态值。 考虑以下两个状态的SISO状态空间模型。

A = [-1.5 -3;
      3   -1];
B = [1.3; 0];
C = [1.15 2.3];
D = 0;
          
sys = ss(A,B,C,D);

x0 = [-0.2 0.3];
t = 0:0.05:8;
u = zeros(length(t),1);
u(t>=2) = 1;
lsim(sys,u,t,x0)
grid on

        如图所示,图的前半部分显示了系统从初始状态值[-0.2 0.3]的自由演化。 在t = 2时,输入有一个阶跃变化,该图显示了系统从此时的状态值开始对这个新信号的响应:

Matlab中lsim函数使用

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/160814.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • maven找不到包但是确实引入了_idea写python好吗

    maven找不到包但是确实引入了_idea写python好吗直接上图,记住选中的位置要设置为Resouce,不然from导入时,会导致pycharm自动从cronjob根路径导入。放linux中运行时,涉及到调用的时候,就报错cronjob包找不到,因为只能找到当前执行的python文件下开始找,而不能找到父包设置位置在这里:…

  • PDF转换图片小工具(高清、无水印、支持随意页数)[通俗易懂]

    PDF转换图片小工具(高清、无水印、支持随意页数)[通俗易懂]疫情期间在家毕业需要,手写签名生成、成绩单的PDF文件需要加入到word中,经历了办理会员、限制5页等等的各种不方便,自己写了个小工具。平台:win764位vs15开发C#语言编写本文所编小工具,不限页数,支持高像素高清截图,随意页截取等截图如下所示:使用说明:选择起始页-终止页-像素级别(部分已默认,可以根据自己需要,一定得填写),点击按钮开始转换–选择需操作的PDF文件–选择转换完的图片的输出路径—程序提示转换成功即可。…

  • transactionscope mysql_TransactionScope事务操作

    transactionscope mysql_TransactionScope事务操作using(TransactionScopetrans=newTransactionScope()){try{InsertUserBase();//它插入不成功,自己回滚UserInfosuserInfo=newUserInfos{UserID=”1″,RealName=”zzl”,};db.UserInfos.InsertOnSubmit(userInfo);db.SubmitC…

  • 几何菱形包是什么牌子_菱形数学

    几何菱形包是什么牌子_菱形数学首页题库状态天梯课程比赛排行榜武则天 题目描述我的提交所有提交换下一题10049:凯的菱形时间限制:1000MS内存限制:131072KB提交总数:36AC总数:20问题描述打完一波农药,安琪拉感叹道:凯呀凯,你真是又帅又厉害!!!凯是真帅,他的菱形也是帅:最外层是字母A,然后里一层是B,再里一层是C。。。输入格式第一行输入整数n(1&lt…

  • 史上最全最详细的Anaconda安装教程[通俗易懂]

    史上最全最详细的Anaconda安装教程[通俗易懂]目录1.Anaconda简介2.Anaconda安装情况的选择2.1情况一2.1.1Anaconda的下载2.1.2测试安装2.1.3更改源2.1.4更新包2.1.5创建和管理虚拟环境2.2情况二2.2.1方法一:通过更改python.exe文件名2.2.2方法二:通过切换虚拟环境3.结束语1.Anaconda简介…

  • 测试19

    测试19一、Linux必备知识linux作为现在最流行的软件环境系统,一定需要掌握,目前的招聘要求都需要有linux能力。二、Shell脚本掌握shell脚本,包括shell基础与应用、shell逻辑控

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号