Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-09 Thread Anne Archibald
2008/8/8 Matthias Hillenbrand <[EMAIL PROTECTED]>: > My Gaussian beam and the lenses have a diameter of approximately 2^16 > array elements while the total array has a size of 2^18. The absolute > value of the Gaussian beam multiplied by the lenses converges to zero. > Out of this reason I assume

Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-07 Thread Matthias Hillenbrand
Hello, today I found time to do concentrate on the proposed zero padding. My Gaussian beam and the lenses have a diameter of approximately 2^16 array elements while the total array has a size of 2^18. The absolute value of the Gaussian beam multiplied by the lenses converges to zero. Out of this

Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-06 Thread Matthias Hillenbrand
Hello, > When you convolve two signals, of lengths N and M, you need to pad the > FFTs to length (N+M-1) before multiplication. > You can take a look at my linear position-invariant filtering code at: > http://mentat.za.net/hg/filter I understand your comments that I need zero padding when doin

Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-06 Thread Stéfan van der Walt
2008/8/6 Matthieu Brucher <[EMAIL PROTECTED]>: > Exactly. Using FFT to do a convolution should be done after some > signal processing readings ;) When you convolve two signals, of lengths N and M, you need to pad the FFTs to length (N+M-1) before multiplication. You can take a look at my linear p

Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-06 Thread Matthieu Brucher
ias Hillenbrand > נשלח: ד 06-אוגוסט-08 05:26 > אל: numpy-discussion@scipy.org > נושא: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT) > > Hello, > > I am trying to calculate the propagation of a Gaussian beam with an > angular spectrum propagation algo

Re: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-06 Thread Nadav Horesh
- מאת: [EMAIL PROTECTED] בשם Matthias Hillenbrand נשלח: ד 06-אוגוסט-08 05:26 אל: numpy-discussion@scipy.org נושא: [Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT) Hello, I am trying to calculate the propagation of a Gaussian beam with an angular spectrum propagation algorithm

[Numpy-discussion] Horizontal lines in diffraction image (NumPy FFT)

2008-08-05 Thread Matthias Hillenbrand
Hello, I am trying to calculate the propagation of a Gaussian beam with an angular spectrum propagation algorithm. My program does the following steps: 1. Generation and multiplication of the Gaussian beam, an aperture, and a lens -> u 2. FFT(u) 3. Multiplication of FFT(u) with the angular spectru