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

Mickael mickael.paris@gmail....
Tue Apr 21 01:36:28 CDT 2009


2009/4/20 Ivo Maljevic <ivo.maljevic@gmail.com>

> 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@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@scipy.org
>> http://mail.scipy.org/mailman/listinfo/scipy-user
>>
>>
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
> Hi Ivo,

thank you very much for your help.

i will try your code today and post the result

thanks,

Mickael.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/scipy-user/attachments/20090421/fb0ca4c8/attachment.html 


More information about the SciPy-user mailing list