How to calculate a Fourier series in Numpy?

by

In the end, the most simple thing (calculating the coefficient with a riemann sum) was the most portable/efficient/robust way to solve my problem:

import numpy as np
def cn(n):
   c = y*np.exp(-1j*2*n*np.pi*time/period)
   return c.sum()/c.size

def f(x, Nh):
   f = np.array([2*cn(i)*np.exp(1j*2*i*np.pi*x/period) for i in range(1,Nh+1)])
   return f.sum()

y2 = np.array([f(t,50).real for t in time])

plot(time, y)
plot(time, y2)

gives me:
alt text

More Related Contents:

  1. Plotting a fast Fourier transform in Python
  2. how to extract frequency associated with fft values in python
  3. Improving FFT performance in Python
  4. Why are frequency values rounded in signal using FFT?
  5. FFT-based 2D convolution and correlation in Python
  6. Comparing two NumPy arrays for equality, element-wise
  7. How to take column-slices of dataframe in pandas
  8. Finding indices of matches of one array in another array
  9. Creating a Pandas DataFrame from a Numpy array: How do I specify the index column and column headers?
  10. How to copy a 2D array into a 3rd dimension, N times?
  11. How do I convert a numpy array to (and display) an image?
  12. A fast way to find the largest N elements in an numpy array
  13. How to add items into a numpy array
  14. multithreaded blas in python/numpy
  15. shuffling/permutating a DataFrame in pandas
  16. SQL-like window functions in PANDAS: Row Numbering in Python Pandas Dataframe
  17. numpy array assignment problem
  18. Fill zero values of 1d numpy array with last non-zero values
  19. Weighted percentile using numpy
  20. pandas.read_csv from string or package data
  21. How to multiply two 2D RFFT arrays (FFTPACK) to be compatible with NumPy’s FFT?
  22. Good ways to “expand” a numpy ndarray?
  23. Find large number of consecutive values fulfilling condition in a numpy array
  24. Multiple matrix multiplication
  25. A faster alternative to Pandas `isin` function
  26. Can i set float128 as the standard float-array in numpy
  27. Why is numpy.prod() incorrectly returning negative results, or 0, for my long lists of natural numbers?
  28. Faster way to transform group with mean value in Pandas
  29. How to add trendline to a scatter plot
  30. Bin values based on ranges with pandas [duplicate]

Leave a Comment