# Estimating the Nyquist rate when simulating data

I want to find the power spectrum $\tilde{Q}(\omega)$ of a quantity $Q(t)$. The values of $Q$ are computed at discrete time steps so I have a series $Q(\Delta t)$, $Q(2 \Delta t)$, ... of values, but the calculation at each time step is independent rather than (for example) needing to know the value of $Q(\Delta t)$ in order to compute $Q(2 \Delta t)$. I then compute the power spectrum of $\tilde{Q}(\omega)$ by an FFT of the signal.

My main question is what is a good way to choose $\Delta t$ i.e. estimate the Nyquist rate of $Q(t)$ without wasting a lot of computing time?

I'm also wondering, given that the calculations at different times are independent, is there a smarter or more efficient way to find $\tilde{Q}(\omega)$ than generating a time-series and performing an FFT.

• Hi! What is the method (or simply the formula if there is any) that you use while computing the samples $Q(\Delta t)$ of $Q(t)$ ? Or stated in signal processing terms is the signal $Q(t)$ band limited? If so what's its bandwidth? Feb 28 '17 at 15:49

## 1 Answer

You must know the method of computing Q(t), and if Q(t) = Q[f(t)], meaning that Q is a function of another function, like the output of a filter. Then, knowing the autocorrelation or power spectrum of f(t), one can use various methods to find the power spectrum of Q. If the function Q is linear with respect to f(t), then the solutions are simple.

However, it seems like you should clear up your needs. An FFT of the signal is simply its spectrum, not is power spectrum or its power spectral density. They are slightly different.

To answer your last question, it seems like you are formulating this in a probability/statistics sense when you say independent? If they are truly independent (i.e. one value of Q gives absolutely no clue about the next value), then your power spectral density represents white noise and it is flat (a constant). Again, it would be best if you clarified your definition of "independent". For example, if Q is a relatively smooth function of time, then it is not statistically "independent" as you say, and to find the spectrum or power spectrum density, you should analytically look at other functions that Q depends on first.