gtstspct β

96/09/08 Ƿ
97/02/01 Ƿ (äȽ­)

 ǽ

(˷٤)- 2Υڥȥ롣

 ֶ֡׼ֳ֥åɤǤɬ٤ǤʤƤ褤ʲ
   ǤϴñΤ٤ꤷʡסʡפȤäɽ
   Ȥ

 (gtstspct_old)Ȥΰ㤤

¸ǽ ʬν(ץ steady=T)

ʤʤäǽ ֶ֡׼ꤹ륪ץ sdim ѻߤ
裱ּȲ

嵭εǽѻߤϡտޤΤȤꡢSDIM  1 ʳǤǽ
ĤޤޤǤΥץβ¤ݤäȤöŪͳ
ΤǤ롣ºݤλѤˤäƤ ּ ǤФȤäƤ
ΤǡSDIM=1 ɬפʤä

 

ϤϰۤʤˤϰۤʤإåĤ(gtseqƤʤ)
ꤹ롣ֶ֡׼ˤϥǡ켡Ѥ()

  ()  gtstspct_old Ǥϥץ SDIM ˤֶ֡׼٤롣

׻Ѥ FFT Թ塢ּּȤ⥰åɿ 2,3,5 ǰʬ
ǤǤʤФʤʤ

 

 

Ϥե

  ڥȥ (ǥե̾gtool.out)
         (ǥե̾gtoolph.out)
  ȿե (@WVNWE* : *ˤϺȿ)
  ưե (@FREQ*-? : *ˤϥǡθĿ?ˤϲʬλ
                    Ĥ줿)

Ǥ롣ڥȥȰϸҡڤӼեϽϤʤ
ȤǤ(ץ -noph, -noax)

ȿʤʤؤȤʤ (M=-MMAX,-MMAX+1,...,0,...,MMAX)ư
ȤȤؤ¤֡κݡʬʿưᣰˤά
ȤǽǤ롣

 ΥޥɤǽˤĤäȤˤʬϤʤ褦ˤ
   ˽ϤǤ褦ɤηаޤǥեȤǤϽ
   ʤ褦ˤƤ( ץ  steady ι򻲾ȤΤ)
   ʬϤ뤫ʤǤĤ뼴ե뤬㤦Τ̵
   ̤ȻפΤ줹褦᤿

Ϥǡ(λ)Ǥ硢ϤϣˤʤäƤޤ
gtool Թ塢ǡˤʤǥեȤǤϿư
ȤƽϤ뤬ץ seq ˤꡢ̤μФ餷ˤ
뤳ȤǤ롣ƥإå ETTL ˼μΰ֤񤭹ࡣ

 Х: SEQ=2ޤ3ξ硢κǽΥåɤϼեΰ
   ΰ֤ꤹΤ gtext եϤˤ ETTL 
   ʤ

 ץ 

out[C]  ϥե̾(ڥȥ. ǥեȤ $GTTMPDIR/gtool.out)
outp[C] ϥե̾(. ǥեȤ $GTTMPDIR/gtoolph.out)
noax[L] T-> եϤʤ (ǥեȤ F)
noph[L] T -> Ϥʤ (ǥեȤ F)
steady[L] T -> ʬ(ء0)Ϥ (ǥեȤ F)
seq[I] (֣ǡΥڥȥˤΤͭ)
       1-3 ꤹȡư3ܤμȤƽϤˡ
       ȤΥǡ seq ܤμФ餷ˤ(seq=1ʤȿ
       Ф餹).
       ǥեȤ 4. ʤưФ餷ˤ.
rstrd[L] T -> FFT˥˥ȥɤ(ޤԤäƤ뤫ɤ
         ǧƤʤäΤʤ.  ǥեȤ F.)
rttrd[L] T -> FFT˥˥ȥɤ (Ʊ)
mis00[L] T -> s==0ʬ»ͤˤ(ꤹȶŪ
         steady=Tˤ) ǥեȤ F. Υץ s==0 ʬ
         礭ƥդ񤯾ǤäȤѤȤ褤.
log[L] إå STYP  2 ˤ
help[L], item[C], unit[C], title[C], dset[C], greset[L] ¿Υޥ
       ȶ̡unitϲꤷʤƤŬʽ񤭴Ԥ롣

 ڥȥȰ

ǤϿư򣰰ʾȤ one-sided ڥȥѤ롣

֥֡å

  x=x_1,x_2,..,x_I, t=t_1,t_2,..,t_J  (ֳ)

Ϳ줿 f_{ij}f(x_i,t_j) 

  f_{ij}= \sum_{l=0}^{J/2} \sum_{k=-I/2}^{I/2}
             E_{kl} cos( ki - lj + _{kl} )

ɽ롣k>0ʡk<0ʬǤ롣ְפ _{kl} 
֥ڥȥפ

  S_{kl} = ( E_{kl} )^2 / 2 (l0ξ)
           ( E_{kl} )^2     (l0ξ)

롣

  (f_{ij})^2 ʿ  \sum_{l=0}^{J/2} \sum_{k=-I/2}^{I/2} S_{kl}  (*)

ȤʤʥѡХˡʤsteady=F ˤʬФʤ
硢դϤ f_{ij} ʿͤΤʿͤˤʤ롣

S_{kl} ñ̤ f_{ij} Τˤʤ롣


 ֥ѥڥȥ̩١פ᤿Ȥ  (ʲ 97/02/01 ˵)

ܤФ褤ʲǤɸŪʻˡǤ롢١֤˴ؤƥա
ꥨѴȤơפˤϲѤɤҤ٤롣

ȿ sư f ΩѿȤ硢ѥڥȥ̩ P_{kl}  
(*) Τ

  (f_{ij})^2 ʿ  \sum_{l=0}^{J/2} \sum_{k=-I/2}^{I/2} P_{kl}sf
                                                                     (**)
    s  sʬ  1
    f  fʬ  1/T (Tϥǡδ)

ΩĤ褦ˤɤäơ

   P_{kl}  S_{kl}/(sf)  TS_{kl}

Ĥޤ T ܤɤ

  -- ڥȥ

u.data  v.data ȤĤΥǡΥڥȥ롣

  % gtstspct u.data outp:u.phase                 #   u.phase
                                                   ڥȥ  gtool.out
  % gtsqrt   out:u.amp                           # ڥȥΥ롼ȡu.amp
  % gtstspct v.data outp:v.phase
  % gtsqrt   out:v.amp
  % gtsub u.phase v.phase
  % gtcos                                        # cos(꺹)  gtool.out
  % gtmlt u.amp                                  #  u.amp ݤ
  % gtmlt v.amp  out:uv.cospect \
      unit:'m|2"s|-2"' title:'cospectrum(U,W)'   # v.amp ݤƽ夬ꡣ
                                                   u,w ñ̡m/sꤷ

ʾǥڥȥ(0 or  ʬ)ޤ롣Ѥ EP ե
ȿưʬۤη׻ȤäȤ롣


