comp.lang.idl-pvwave archive: archive » Re: CHISQR

Home » Public Forums » archive » Re: CHISQR_CVF question.

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: CHISQR_CVF question. [message #67631]

Thu, 20 August 2009 09:57

David Fanning
Messages: 11724
Registered: August 2001

Senior Member

R.G. Stockwell writes:

> Craig, Sorry but I am a bit confused here.
>
> using the +1 direction is the "inverse" FFT here isn't it?
> and hence it lacks the 1/N normalization that occurs on the "forward" FFT.
> Is that right?
>
> Also, total(err^2) happens to be equal to the length here, so i looks like
> you are doing an inverse FFT ^2, and then dividing by len.
>
> BUT, that is the same as doing the forward FFT (with 1/N), squaring it, then
> multiplying
> by len.
>
> So, it almost looks like this just happens to be by coincidence the same as
> pow = fft(lc, /forward)*length
>
> And you have a factor of 2, which is coincidentally also the power of your
> spectrum. and it appears that again this may have just coincidentally
> cancelled out.
>
>
> basically, I am starting with a normalization of the spectrum as:
>
> d = 120*randomn(seed,len)
> spe = fft(d)
> pspe = abs(spe[0:len/2-1])^2
>
> ; normalize wrt length and variance, so we always get the same result
> pspe = pspe*(len)
> pspe = pspe/stddev(d)^2
>
>
> with this normalization, the mean of my spectrum is always the same.
> (as i vary the length of the time series, and as i vary the standard
> deviation,
> above i have a stdev of 120).
>
> Are you saying that there should be a factor of 2 in my power spectrum,
> i.e. I need a final line that states pspe = pspe*2?
> Because, when I do this, I do get the expected result. By expected I mean I
> calculate the number of points above the cutoff level (90%) and I find
> approximately
> 10% above, 90% below. ditto 95%, 99%.
>
> But, I want to justify that factor of 2.

You know what? I'm just going to stick with mastering
that down-the-line backhand, thank you very much!

Cheers,

David

--
David Fanning, Ph.D.
Coyote's Guide to IDL Programming (www.dfanning.com)
Sepore ma de ni thui. ("Perhaps thou speakest truth.")

Report message to a moderator

[Message index]

		Re: CHISQR_CVF question. By: David Fanning on Thu, 20 August 2009 09:57
		Re: CHISQR_CVF question. By: R.G. Stockwell on Thu, 20 August 2009 09:53
		Re: CHISQR_CVF question. By: Craig Markwardt on Thu, 20 August 2009 09:14
		Re: CHISQR_CVF question. By: R.G. Stockwell on Thu, 20 August 2009 08:21
		Re: CHISQR_CVF question. By: Craig Markwardt on Thu, 20 August 2009 01:03
		Re: CHISQR_CVF question. By: R.G. Stockwell on Wed, 19 August 2009 13:42
		Re: CHISQR_CVF question. By: pgrigis on Wed, 19 August 2009 12:29
		Re: CHISQR_CVF question. By: Craig Markwardt on Sat, 22 August 2009 10:51

Previous Topic:	Re: Loop Multiple Files and Rename
Next Topic:	Re: Voxels in IDL

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Nov 28 11:07:21 PST 2025

Total time taken to generate the page: 0.00270 seconds