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) > > > > > See http://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 -~----------~----~----~----~------~----~------~--~---
