Transforming a spatial average into a spectral average

  • Thread starter Thread starter nickthequick
  • Start date Start date
  • Tags Tags
    Average
nickthequick
Messages
39
Reaction score
0
Hi,

First off, apologies if this is in the wrong place - any directions on where this is more appropriate are appreciated.

I am trying to figure out how to convert a smoothing operation in physical space into one in Fourier space.

In physical space, with equally spaced points x_j, (j =1,2,3...), we define an averaging of a function f(x) as

\bar{f}_j =\frac{1}{16}(-f_{j-2} +4 f_{j-1} + 10 f_j + 4f_{j+1} -f_{j+2}),

where f_j=f(x_j). The operation here is a type of moving average that smooths a function that has an oscillation, of period 2, about a smooth function.


Now, I do not have information in physical space, rather I have the Fourier coefficients of f, i.e. I know the set \hat{f}_k, (k=1,2,...), and I would like to apply the averaging to these coefficients.


I'm wondering if there's an analog for this spatial averaging in Fourier space, but am very much out of my area of expertise.

Any tips, or even relevant references, would be much appreciated.

Nick
 
Physics news on Phys.org
It might be easier to understand this if you think about the inverse operation. Working backwards, your averaging process converts a single value ##\bar f_j = a## (and all the other ##\bar f## values = 0)into a "smeared out" set of 5 values ## f_{j-2} = -a/16##, ##f_{j-1} = a/4##, etc. and all the ##f## values = 0 except those 5.

The process is linear, because if you had for example ##\bar f_j = a## and ##\bar f_{j+1} = b## (and all the other ##\bar f## values = 0), working backwards would give you 6 non-zero values ## f_{j-2} = -a/16##, ##f_{j-1} = a/4 - b/16##, ... , ##f_{j+3} = -b/16##.

So, to see what this does in the frequency domain, you can do Fourier transforms of the output data ## ... 0, 0, 0, 0, a, 0, 0, 0, 0, ...## and the corresponding input data ## ... 0, 0, -a/16, a/4, 5a/8, a/4, -a/16, 0, 0, ...## and compare the FFTs term by term in the frequency domain.

And since the "amplitudes" of the two data sets are both proportional to a, you might as well take a = 1.

Look for information about digital signal processing for more details - this corresponds to a finite impulse response (FIR) filter and its frequency response. (The "impulse" is the single value of ##\bar f##, and "finite response" means there are a finite number (5) of non-zero values of ##f##.

This will also be covered in books and websites on image processing, but they will be more about processing 2D or 3D arrays of data, not your 1D case.
 
AlephZero,

Thanks, this is very helpful!Nick
 
##\textbf{Exercise 10}:## I came across the following solution online: Questions: 1. When the author states in "that ring (not sure if he is referring to ##R## or ##R/\mathfrak{p}##, but I am guessing the later) ##x_n x_{n+1}=0## for all odd $n$ and ##x_{n+1}## is invertible, so that ##x_n=0##" 2. How does ##x_nx_{n+1}=0## implies that ##x_{n+1}## is invertible and ##x_n=0##. I mean if the quotient ring ##R/\mathfrak{p}## is an integral domain, and ##x_{n+1}## is invertible then...
The following are taken from the two sources, 1) from this online page and the book An Introduction to Module Theory by: Ibrahim Assem, Flavio U. Coelho. In the Abelian Categories chapter in the module theory text on page 157, right after presenting IV.2.21 Definition, the authors states "Image and coimage may or may not exist, but if they do, then they are unique up to isomorphism (because so are kernels and cokernels). Also in the reference url page above, the authors present two...
When decomposing a representation ##\rho## of a finite group ##G## into irreducible representations, we can find the number of times the representation contains a particular irrep ##\rho_0## through the character inner product $$ \langle \chi, \chi_0\rangle = \frac{1}{|G|} \sum_{g\in G} \chi(g) \chi_0(g)^*$$ where ##\chi## and ##\chi_0## are the characters of ##\rho## and ##\rho_0##, respectively. Since all group elements in the same conjugacy class have the same characters, this may be...
Back
Top