[SciPy-user] Convolve an input with a seri of impulses

Ivo Maljevic ivo.maljevic at gmail.com
Mon Apr 20 16:05:50 EDT 2009


I never use LTI systems defined in this way, but out of curiosity I tried to
make some changes into my original code in a hope I will learn something and
help you in a way.

Try the attached script, it uses your code.

Ivo

2009/4/20 Mickael <mickael.paris at gmail.com>

> Hi,
>
> I'm trying to convolve a reference signal, like a step signal, with a
> sequence of impulses.  Impulses have different amplitudes and timing :
>
>
> ----------------------------------------------------------------------------
> zeta = 0.0193
> wn = 35037.0
> wd = wn*math.sqrt(1.0-zeta**2)
> K = math.exp(((-zeta*pi)/math.sqrt(1.0-zeta**2)))
>
> num = array([9.871])
> den = array([1,1354,1.228*10**9])
> H = signal.lti(num, den)
>
> # Impulses : amplitudes
> A1 = (1.0/(1.0+K))
> A2 = (K/(1.0+K))
> # Impulses : timing
> t1 = 0
> t2 = (pi/wd)
>
> -----------------------------------------------------------------------------
>
> I would like to apply to "H" the result of a step signal convolve with
> these two pulses but I'don't know to define the vector of the impulses.
> --->T,yout,xout = signal.lsim(H, entree, t)  <=== entree = scipy.signal.
> convolve(*A_step*, My_Impulses?, mode='full')
>
> Thanks,
>
> Mickael.****
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20090420/b9979d7d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: conv2.py
Type: text/x-python
Size: 858 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20090420/b9979d7d/attachment.py>


More information about the SciPy-User mailing list