I try to apply a periodical pulse signal to a pT1 filter and to convert the result to the time domain.
When the code below is executed, I receive a time domain function, which returns NaN after a certain point in time.
This time point is independent of the chosen pulse period - here 0.2s. Up to 0.2s all looks good.
A sinusodial signal works as expected....
A = 20e3;
B = 10e3;
H_filter = B/(C*A*B*s + (A + B));
T = 100e-3;
v = 10;
Hper_pulse = (1-exp(-T/v*s)) / (s*(1-exp(-T*s)));
in_t = ilaplace(Hper_pulse);
fplot(in_t, [0 0.25]);
out_s = Hper_pulse * H_filter;
out_t = ilaplace(out_s);
fplot(out_t, [0 0.25]);
Thanks in advance