Delayed step/ impulse response

99 次查看（过去 30 天）
Daniel Stankowski
Daniel Stankowski 2020-1-14
回答： Paul 2020-9-20
Hi everyone,
I would like to plot an impulse and step response of some arbitrary system sys1.
Normally I would use just a command impulse(sys1,t) or step(sys1,t) however the signals in my case are delayed.
1) The impulse is represented as: diract(t-5)
2) step is represeneted as: 1(t)-2*1(t)(t-tsw)
I know some people use an exponential function for it, but I am unsure why and how,
I would greatly appreciate any help

请先登录，再回答此问题。

回答（2 个）

Raj
Raj 2020-1-14
You have your system model and input signals equation. You can just generate your input signals and get the system response for those signals using lsim command.
  2 个评论
Daniel Stankowski
Daniel Stankowski 2020-1-16
Thanks for help but impulse as well as step does not work well with lsim.
I have solved the problem by using pure delay element.
So to plot impulse represented as diract(t-20) i did the folowing
impulse(exp(-20*s)*Tf,t)
where Tf is a closed loop transfer function of mine.
Thanks anyway!
Bill Tubbs
Bill Tubbs 2020-9-20
It seems a shame that this wouldn't be possible as an option in stepDataOptions, similar to specifying the step amplitude. Something like this perhaps:
opt = stepDataOptions('StepTime',10,'StepAmplitude',2);
step(sys,opt)

请先登录，再进行评论。

Paul
Paul 2020-9-20
Use the time invariance and lineraity properties of an LTI system. One approach is to generate the the nominal impulse or step response and then use apppropriate time shifing and superposition of the nominal outputs.
sys1 = tf(1,[1 1]);
tsw = 4.3;
t = 0:.01:10; % Important to define tnom such that tsw is one element of t. verify this
sum(t == tsw);
% nominal impulse response
ynom = impulse(sys1,t);
% delayed impulse response
yimp = 0*ynom;
yimp(t >=tsw) = ynom(1:sum(t >= tsw));
% nominal step response
ynom = step(sys1,t);
% first term
y1 = ynom;
% delayed and scaled second term
y2 = 0*ynom;
y2(t >=tsw) = 2*ynom(1:sum(t >= tsw));
% total output
ystep = y1 + y2;
plot(t,[yimp ystep]),grid

类别

在 Help Center 和 File Exchange 中查找有关 Frequency-Domain Analysis 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by