System impulse response and Convolution by matlab

2011 2 20
3 个回答
更新时间：2021 2 25
4 次查看（30 天）

1 个投票

Hi everyone, i am very new to matlab, and would like to konw how to obtain y(t) from below
Generate a unit step function as the input function, x(t), and an exponentially decay function as the impulse response function, h(t), such as h(t)=exp(-t/2) (note: 2 is the time constant of the system dynamic response). Using MATLAB to calculate the output of the system, y(t).
Thank you so much in advance

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

回答（3 个）

psyx21
psyx21 2011-4-29

1 个投票

hey friends what will be the peak acceleration response for 100g 6ms half sine pulse with zeta=0. I need a matlab script with natural frequency on x axis..thanks
Paulo Silva
Paulo Silva 2011-2-20

0 个投票

t=0:0.1:10;
u=0*t;
u(t>=0)=1;
h=exp(-t/2);
y=u.*h;
plot(t,y)

8 个评论
显示 6更早的评论 隐藏 6更早的评论

gedaa
gedaa 2011-2-21
thank you sooo much, but I have one more question if it's not 2 much, what is u=0*t and how did u come up with that?
Paulo Silva
Paulo Silva 2011-2-21
it's the same as
u=zeros(size(t))
It creates a variable with the same size as t but full of zeros. I learned that trick here at matlabcentral :)
David Young
David Young 2011-2-21
Paulo - I just wonder if this is right. In the time domain, don't you need to convolve the input and the impulse response rather than multiplying them?
Paulo Silva
Paulo Silva 2011-2-21
Yes you are correct, try
t=0:0.1:10;
u=0*t;
u(t>=0)=1;
h=exp(-t/2);
y=conv(u,h);
plot(y)
Paulo Silva
Paulo Silva 2011-2-21
Correction for the amplitude and time scale
T=0.1;
t=0:T:10;
u=0*t;
u(t>=0)=1;
h=exp(-t/2);
y=conv(u,h);
plot(t,T*y(1:numel(t)))
gedaa
gedaa 2011-2-21
ok, so if the problem asks for a pulse function instead of a unit step function, it would be the exact same codes except with the addition of PulseWidth=1, am i correct?
Paulo Silva
Paulo Silva 2011-2-21
I'm not sure about the amplitude at
plot(t,T*y(1:numel(t)))
with the step function is good compared to
step(tf([1],[1 1/2])) but with the impulse it's plot(t,y(1:numel(t))) without the T, I can't figure out why that happens.
gedaa
gedaa 2011-2-23
hey, paulo, i actually checked with my professor and he said that T doesn't need to be multiplied in the last line. Thank you soo much

请先登录，再进行评论。

Arpan Patel
Arpan Patel 2021-2-25

0 个投票

t=0:0.1:10;
u=0*t;
u(t>=0)=1;
h=exp(-t/2);
y=u.*h;
plot(t,y)

类别

在 帮助中心 和 File Exchange 中查找有关 Programming 的更多信息

标签

提问:

gedaa
gedaa
2011-2-20

回答：

Arpan Patel
Arpan Patel
2021-2-25

Community Treasure Hunt

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

Start Hunting!

Translated by