[SciPy-user] Convolve an input with a seri of impulses
Mickael
mickael.paris at gmail.com
Tue Apr 21 02:36:28 EDT 2009
2009/4/20 Ivo Maljevic <ivo.maljevic at 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 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
>>
>>
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user at 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