Re: How to interpolation in frequency domain will not affect time domain signal?



On Jul 1, 9:20 am, "coly" <wolon...@xxxxxxxxx> wrote:
But what puzzled me is that my interpolation method in the spectrum leads
to a distorted time data. I want to know why.

There are many different interpolation methods. You can zero pad
things, fit parabolas to the data, use a zoom technique, etc. You
could use the DFT or inverse DFT to compute as many millions of time
or frequency domain points as you want. One of the nice things about
the DFT or IDFT is that you can use fractional values of 'n' (time
index) and 'k' (frequency index).

But no matter how you go about it, you're going to have the same
fundamental problem. Let me see if I can explain it in the following
way. Astronomers routinely capture millions of data points and do
FFTs on them. They also understand that they can express what the FFT
is doing as an N equation, N known and N unknown type of problem,
where N is in the millions.

Now suppose you said to them: "I can take 100 data points in time, FFT
it to get 100 frequency domain points, then zoom interpolate my 100
frequency domain points to create a million points. Then I can inverse
transform to obtain a million points in the time domain. So you
wouldn't have to get millions of samples, just 100 of them."

They’ll probably think many things about that kind of statement, none
of them good, and all of them reflecting badly on the person making
the statement. But, being polite, they might just say: "But that
makes no sense at all. We have a million knowns (our million data
samples), and we need to solve for a million unknowns (a million
frequency domain points). It's mathematically impossible to solve for
a million unknowns when all we have is 100 knowns, unless we impose
some ridiculous assumptions on our signal processing model to reduce
the million unknowns down to 100."

So just imagine that you actually had the million point astronomy
input data. Pick any 100 of those points you want, and do whatever
interpolating you wish and try to ‘create’ or 'estimate' all those
other points from your interpolation scheme. Can you begin to
appreciate now why your interpolated results will have errors when
compared to the actual data? Unless you can impose some extraordinary
assumptions to simplify the problem, it’s an impossible task no matter
what interpolation technique you use.

As Rune pointed out, you can use other methods, such as frequency
estimation. Basically, you’d be using a different signal processing
model to obtain the needed information from the available data.

Kevin McGee
.



Relevant Pages

  • Re: measuruement of peak amplitude of each cycle of sinusoid having varying frequency
    ... 3 point frequency domain interpolation to determine the frequency ... of a sinusoid only works for sinusoids. ... There is time interpolation in the time domain and frequency ...
    (comp.dsp)
  • Re: interpolation accuracy, oversampling and fractional interpolation
    ... This is implemented via and forward FFT into the ... This results in no change in the frequency domain characteristics though ... interpolation and the kind used in DSP ... audio work for up sampling. ...
    (comp.dsp)
  • Re: interpolation fft
    ... Interpolation between those is wrong, unless you smooth in the ... frequency domain to approximate a continuous function. ... You need to recalculate the 2nd set of data using the same FFT ...
    (comp.soft-sys.matlab)
  • Pilot based interpolation
    ... i'd like to ask about pilot based interpolation in channel estimation ... I have a block of reference cells which i need to ... interpolate in time and frequency domain. ... the help of pilot cells. ...
    (comp.soft-sys.matlab)
  • Pilot based interpolation
    ... i'd like to ask about pilot based interpolation in channel estimation ... I have a block of reference cells which i need to ... interpolate in time and frequency domain. ... the help of pilot cells. ...
    (comp.dsp)