that is wonderful!! by the way, I am looking forward to see those LTI and Fourier Transform features applied to our powerful symbolic system :D
Maurizio On Jul 26, 9:23 pm, Rafael Cardoso Dias Costa <rfc...@gmail.com> wrote: > The problem is solved! > > http://www.sagenb.org/home/pub/677/ > > Thanks, Jason! > Thanks, everyone! > > On Sat, Jul 25, 2009 at 10:40 PM, Jason Grout > <jason-s...@creativetrax.com>wrote: > > > > > > > Rafael Costa wrote: > > > Why Sage not execute the python code? > > > > from scipy import * > > > from pylab import * > > > > sample_rate=1000.00 > > > t=r_[0:0.6:1/sample_rate] > > > N=len(t) > > > s=sin(2*pi*50*t)+sin(2*pi*70*t+pi/4) > > > S=fft(s) > > > f=sample_rate*r_[0:(N/2)]/N > > > n=len(f) > > > plot(f,abs(S[0:n])/N) > > > Seehttp://sagenb.org/home/pub/676/for a working version: > > > from scipy import * > > from pylab import * > > sample_rate = 1000.0 > > t = r_[0:0.6:1/sample_rate].astype(float) > > N = len(t) > > s = sin(2*pi*50*t)+sin(2*pi*70*t+pi/4) > > S = fft(s) > > f=sample_rate*r_[0:(N/2)]/N > > n=len(f) > > line(zip(f,abs(S[0:n])/N)) > > > As Minh points out, the need for the "astype" above will go away once > > the three tickets (that are up for review!) are reviewed and merged; > > hopefully in the next few days. > > > Thanks, > > > Jason > > -- > Rafael Cardoso Dias Costa > Engenheiro de Controle e Automação - UFMG > Mestrando em Engenharia Química - Controle de Processos - UFRGS > +55 51 8541-0926 :: Porto Alegre - RS > +55 51 8430-1421 :: Porto Alegre - RS > +55 31 9165-5488 :: Belo Horizonte - MG > +55 31 3466-0146 :: Residência em Belo Horizonte - MG --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
