微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何在 scilab 程序不是 xcos中模拟干扰

如何解决如何在 scilab 程序不是 xcos中模拟干扰

你好吗?。

我需要使用 scilab 模拟控制系统中的扰动,即 csim 函数用于通过使用阶跃、脉冲、斜坡或任何其他输入来模拟系统的响应,但是,我需要输入例如在 t = 0.5 秒内的干扰以查看系统行为。

这给我带来了另一个问题,因为我不知道如何让 csim 或 syslin 来确认两个不同的输入,或者就像定义两个系统一样简单,一个带有参考输入,另一个带有干扰入口并求和?.

预先感谢您的帮助。

解决方法

假设您有以下线性时不变系统 (A,B,C)

x'=A*x+B1*v+B2*d
y=C*x

B=[B1,B2],其中 v 是控制/输入,d 是干扰。如果你想要例如模拟阶跃响应和干扰,您必须定义自己的整体输入 [v;d] 并决定何时应用干扰。下面是一个例子:

function ud = step(t)
    ud = [1;0];
endfunction

function ud = input(t)
    ud = zeros(t);
    ud(1,:) = 1;
    s = abs(t-0.5);
    // 0.2 is the half-width of disturbance
    ud(2,:) = -0.1*(1-s/0.2).*(s<0.2)
endfunction

A = [-2 1;
      1 -2];
B1 = [1;0];
B2 = [0;4];
C = [0 1];

sl = syslin('c',A,[B1 B2],C);
t = linspace(0,5,1000);

x = csim(step,t,sl)
xd = csim(input,sl)

clf
plot(t,x,xd,input(t)(2,:))
legend('step','step and disturbance','disturbance',2)

我在这里进行了两次 csim 调用,一次用于通常的阶跃响应,第二次用于扰动阶跃响应。但是,我警告您 csim 使用的 ode 求解器:很容易错过不连续的输入,这就是我在这里应用帽子形干扰的原因。扰动输入的代码旨在允许向量时间输入,以便轻松绘制扰动图。 enter image description here

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