t[GvZ\tg
Ԋu̐l񂩂A̐l̐g̎gA
UAʑZo邽߂̃\tgQłB

e
fftcyc.exe:t[GϊŁAf[^Ɋ܂܂ǧ񋓂܂B
dftcyc.exe:fftcyc.exeŎZoǧ₪ۂɂ邩mFA
           ̐UƈʑvZ܂B
subtract.txt:f[^␳t@C̃Tv

es2000.txt:Tvf[^
  iJPL-Horizons,zʒuAD2000-25004Ԋu,ωW)
 *) 1- 7s:2000N11̒ʓ(0h-UTC)
 0) 9-18s:zo|ϑzo(b)
 1)22-29s:zo|ϑzo|(359.9905x/N̐)|(719.981x/N̐)(b)
 2)33-39s:z(b)
 3)43-51s:zn(mAU)
 4)56-62s:zn|ϋ|(359.9905x/N̐)(mAU)
setting0.ini:es2000.txtpݒt@C[0)]
setting1.ini:es2000.txtpݒt@C[1)]
setting2.ini:es2000.txtpݒt@C[2)]
setting3.ini:es2000.txtpݒt@C[3)]
setting4.ini:es2000.txtpݒt@C[4)]

===============================================

fftcyc.exe(iK)
ʏ̍t[Gϊ(FFT)́Af[^(2̗ݏɌ)
ŜPƂgyт̐{̎g
vZo܂B

ŁAf[^ŜAJn_炵Ȃ琔FFTs܂B
Ƃ200,000Tṽf[^ɂāAJn_10,000_
炵ȂA2^17=131,072TvFFT7s܂B
eFFTł131,072TvPƂgyт̐{
ĝ݂Ńt[GvZ܂AV̌vZɂA
gł̐Uƈʑrĕω܂B
̔rɂAFFTł̎g\Q`R
g߂鎖o܂i^Ȃ΁jB

gp@
1)܂̓f[^t@CɊւf[^͂܂B
RgƂĖs̐擪A
vZɗpf[^es̊JnʒuƕAt@CłB
uSubtract Filev͂Ƃ肠󔒂ɂĂĂB

2)L͌Check DataƁAt@CPǂݍŁA
f[^ƍőlAŏl\܂B
KsKv͂܂񂪁Aȍ~Őݒ肷鐔l
肷ł̎QllƂȂ܂B

3)PFFTŗpf[^îQ̗ݏ̎wjw肵܂B
Rf[^ŜKv܂B
傫XFFŤvZx͗ǂȂ܂A
f[^Ŝɋ߂FFT̉񐔂ȂȂ܂B

4)Qڈȍ~FFTvZۂ́A擪ʒũVtg
ݒ肵܂B
 int{(Sf[^|[3)Œ߂f[^])[Vtg]}
FFTs܂B
VtǵA΁uɒ[ɑ傫ȐU̐ꍇ
듮vh܂AFFT̉񐔂鎖ŐUEʑω
΂₷Ȃ܂A
ԂɁAʑωȂďo̓t@C
LĂ܂̂ŁAx̑傫͕Kvł傤B
̂߁ASf[^3)Œ߂f[^߉߂
FFT񐔂mۏoȂȂ܂B

5)px\ۂ̒PʂƂāAuPvɑ
f[^͂܂B
Ƃ΂PP̃TvOłPNPPʂ
ꍇ́A̐l365i܂365.25jƂ܂B
̎APNpx360deg/cycleƂȂ܂B

6)f[^擪̌_ɑ_܂ł̎Pʐ
͂܂B
Ƃ΁APʂ365.25Af[^擪1980N11A
2000N11_Ƃꍇ20łB
fftcyc.exePڂ̌vZł͂قǏdvȍڂł͂܂񂪁A
dftcyc.exeł̈ʑvZAZofftcyc.exe
ꍇɕKvƂȂ܂B

7)FFTł̓f[^[3)Œ߂]̔Ɠ̎gɂāA
t[Gϊ܂Ƃ߂čs܂B131,072TvFFTs΁A
65,536̎gɂăt[Gϊs܂B
FFTs߁ASo͂ƂƂłȂʂɂȂ܂B
臒lw肵A傫ȐU̐
o܂BFFTsAw肵臒lPł
gt@Cɏo͂܂B2)ŕ\ꂽ
őlAŏlA̍ȂǂQlɌ߂ĂB

8)o͂t@Cw肵܂B

9)ȏ͍̓ڂ́Aݒt@CƂĕۑł܂B
܂Aqdftcyc.exeǂݎA
Kvڂ̂܂܎gp鎖ł܂B

10)uCalculatevƌvZăt@Co͂܂B
FFT̂͂̂łAFFT̉񐔎ł
ӊOƎԂ܂B
܂AFFTvẐ߂ɋȔzp̂ŁA
ɋNĂ\tg̓ɉe^ꍇ܂B

o̓f[^ɂ
o̓f[^͎̂悤ȐlA臒l𒴂px̐
т܂B

Number of data 45656 
Length of FFT 32768 
Step for FFT-block 3652 
--------------------------------
Avg-Swing(St-D)/D-Swing,Avg-Phase(St-D)/D-Phase:     -0.0014cycle----:     39.9932cycle----:     79.9877cycle----:    119.9822cycle----
>>>>   298(   298.9503deg)-> 299.267170518419 deg
   0.1358(   0.0023), 124.6561deg(  14.1136deg):   0.1324, 105.3560dg:   0.1368, 118.8698dg:   0.1386, 131.0289dg:   0.1355, 143.3695dg
   0.0010(   0.0031),  12.6712deg(   0.6004deg):   0.0000,   0.0000dg:   0.0044,  13.5138dg:   0.0019,  12.1591dg:  -0.0032,  12.3406dg
------------------------
>>>>   314(   315.0014deg)-> 315.557157933504 deg
   0.1172(   0.0036), 262.5140deg(  24.8259deg):   0.1232, 229.2900dg:   0.1165, 251.3999dg:   0.1139, 273.3911dg:   0.1152, 295.9750dg
  -0.0027(   0.0033),  22.2283deg(   0.2561deg):   0.0000,   0.0000dg:  -0.0067,  22.1099dg:  -0.0026,  21.9911dg:   0.0013,  22.5840dg
>>>>   315(   316.0046deg)-> 315.559010845009 deg
   0.1484(   0.0023), 112.9935deg( 141.0746deg):   0.1445,  49.8200dg:   0.1493,  31.6868dg:   0.1505,  14.1059dg:   0.1494, 356.3611dg
   0.0017(   0.0025), -17.8196deg(   0.2316deg):   0.0000,   0.0000dg:   0.0049, -18.1332dg:   0.0012, -17.5809dg:  -0.0011, -17.7448dg
------------------------
>>>>  4816(  4831.3586deg)-> 4832.01960102154 deg
   0.2445(   0.0001), 161.5849deg(  29.5561deg):   0.2444, 121.9321dg:   0.2446, 148.3630dg:   0.2444, 174.8083dg:   0.2444, 201.2363dg
  -0.0000(   0.0001),  26.4347deg(   0.0076deg):   0.0000,   0.0000dg:   0.0001,  26.4308dg:  -0.0001,  26.4453dg:  -0.0000,  26.4280dg
>>>>  4817(  4832.3618deg)-> 4832.01996209175 deg
   0.4727(   0.0001), 281.4432deg(  15.2872deg):   0.4728, 301.9534dg:   0.4726, 288.2808dg:   0.4727, 274.6041dg:   0.4726, 260.9346dg
  -0.0001(   0.0001), -13.6729deg(   0.0030deg):   0.0000,   0.0000dg:  -0.0002, -13.6725dg:   0.0000, -13.6767dg:  -0.0000, -13.6695dg
>>>>  4818(  4833.3650deg)-> 4832.02002352936 deg
   0.1202(   0.0000), 221.2386deg(  60.1421deg):   0.1202, 301.9365dg:   0.1201, 248.1269dg:   0.1202, 194.3322dg:   0.1202, 140.5589dg
   0.0000(   0.0001), -53.7925deg(   0.0149deg):   0.0000,   0.0000dg:  -0.0001, -53.8096dg:   0.0001, -53.7947dg:  -0.0000, -53.7733dg
------------------------

u>>>>vŎn܂śApxɊւf[^łB
ŏ̐FFTvZł̊px̒ʂԍA
̎()FFTvZɗppxł
(w肵PʂP360xƂꍇ̐l)B
̎́u->v̌̊px
uۂɊ܂܂Ǝvpx̌vłB
L̗̏Q`RԖځAS`UԖڂɂāA
ꂼꏬ_ȉQڂ܂łقڈvĂ̂ŁA
̒x̐x͂Ǝv܂B

AK̐ƕۏ؂Ăł܂B
ȂoĂ\ƂĂ

a)rI߂Q̎gx̐UĂꍇA
̒Ԃ̎g̐悤Ɍ鎖܂B
Q̎g߂قǁA܂FFT̃f[^Zق
o₷Ȃ܂BRȏ̎g߂ɏW܂Ăꍇ
ƕGȌۂN܂B

b)ɒ[ɑ傫ȐUꍇA◣ꂽpxɂĂ
̉eo܂BuU傫AꂽpxvƁA
QFFTvZł̈ʑω}180x𒴂Ă܂܂B
}180x𒴂ʑpω͌ołȂ̂ŁAʂƂ
uvZԈႢv\܂ĩGCAXۂłjB

tɁAۂɂ͂gȂP[XƂāA

c)ɒ[ɑ傫ȐU̐ꍇɁAɋ߂g
ȐU̐BĂ܂܂B

b)c)ɊւẮAdftcyc.exeő傫ȐU̎g
Uƈʑm肵Af[^獷ŁA
ēxsƁAuvZԈႢv̕A
uBĂȐv܂B
a)́AQ̎gFFT\ȏ㗣Ă
ooꍇ܂A߉߂ƎʏoȂȂ܂B
߉߂ꍇ́Adftcyc.exeōׂώ@Kv܂B


pxf[^̉QśAt[Gϊ̌ʁiUAʑj
\Ă܂B
PsڍA
U̕ϒliW΍jAʑ̕ϒliW΍jFPڂ̐UAʑFQڂ̐UAʑFc
QsڍA
Uω̕ϒliW΍jAʑω̕ϒliW΍jF0,0:PڂQڂ܂ł̐UωAʑωFc

Psڂ͐UȊO܂dvł͂܂BUڈxłB
Uƈʑ̐mȒldftcyc.exeŋ߂܂B
Qsڂ́AUω̑傫W΍傫ꍇA
uςȎNĂvƋ^܂Bʑω̕W΍
傫ꍇƂ₵łB
uςȎvȂ̂́A̐l͕Ȃ̂ŁA
dftcyc.exeł̎głŕ\Ă݂Kv܂B

=========================================================

dftcyc.exe(iK)
w肵pxɑ΂闣Ut[GϊAf[^Ƃ
vZA摜\܂̓t@Co͂܂B
܂Sf[^ɑ΂闣Ut[Gϊ̌ʂ\܂B
łUt[GΐAX̒`ɊÂ
fȃt[GvZłixƂjB

g
1)u~~~~~~~~v́Afftcyc.exeƓl̐ݒ荀ڂłB
fftcyc.exeŕۑݒt@Cǂݍ߂Ί܂B

2)t[GvẐ߂̊px͂܂B
fftcyc.exěvZʂŏoĂȂǂ܂B

3)P̃t[GϊŌvZf[^cyclePʂœ͂܂B
Disp(摜\p)File(t@Co͗p)܂A
t@Co͗p̎ŉ摜\ł̐Fς悤
Ă܂
(ĉ摜\݂̂ł0ȊO̒lKv܂)B
ZƂꂾ_ċȐ炵ȐɂȂ܂A
f[^ł͖\܂܂B
Ƃꂾuωv̂ŁA
Ȃ܂A_܂B
܂Af[^̎Ƃ̌ˍň̃GCAXۂ
N鎖܂Blςĉx悢ł傤B

4)摜\̏ꍇ͐Uɍ킹Ĕ{͂܂B
摜̈͏cƂ-200`+200ł̂ŁAU1x̏ꍇ
100{xǂł傤B

5)uDrawvΉ摜\n߂āAŌɑSf[^ɂ
t[GvŽʁiUƈʑj̕\gɕ\܂B

6)o̓t@Cw肵āuCalc.v΁AvZʂ
t@Cɏo͂܂B
܂ASf[^ɂt[GvŽʂ̕\g
\AquPole datavF̊ۂŐ}܂B


摜̉
a)_ȊÔقڂP_ɗ܂ĂAP_̎͂
xȂA
u̎gmɂBvZʂ̐Uƈʑ
Ƀf[^΁A̎g͂قڎ菜v
b)_̎ĂȂ
u̎g͂ȂvB

c)_ȊÔP_̎ŉ~`ꍇ́A
u̎g͂BāA̋߂̎g̐B
w߂̎gxƂ́A悻
wt[GvZɗpĂpx{~`̂PPʓ̉]pxxvB
d)_𒆐Sɉ~`ꍇ
uȂ΃t[GvZɗpĂg͒߂A
EȂ΍߂BCʂ́A~`̂PPʓ̉]pxvB

e)ɒPUꍇ́A
u̎ĝ͂ȂB̎g𒆐Sɋ
Q̎gB̂Q̊px
wt[GvZɗpĂpx}(PŮpx)xv
f)A܂͌_ʂ₩ȉ~ʏłقڈ葬ɓꍇA
f1)UP֐ŕωĂB
f2)e)̃P[Xɑ邪APU̎߂
  ^̂悤ɌB
f3)ĩf[^ߎ̃f[^Ăꍇj
  ߎ̊pxۂƂ͔ɂĂB̂
  w(ߎ̈ʑƒւ̂PPʓ̈ړ)(ߎ̐U)x(rad/cycle)
g)uT[W͗lv`ꍇ
uQi܂͂ȏj̎gAvZɗpĂg
  ΂ɂv

d),f)ȂǂŐUʑ̕ω𐳊mɒm肽ꍇ
t@Co͂sāA̓emF܂B


t@Co
t@Co̗͂łB

Angular Speed: 359.990502 
cycle-----)Swing------,Phase-------:d-Swing----,d-Phase-----:d-Vec(Lng)-,d-Vec(Ang)--
   50.0068)    0.17189,   9.2889deg:-----.-----,----.----deg:-----.-----,----.----deg
  150.0041)    0.49746, 295.0495deg:    0.32558, 279.4850deg:    0.48018, 274.8971deg
  250.0014)    0.08836, 179.2322deg:   -0.40911,-109.5371deg:    0.54182, 123.4911deg
  349.9986)    0.04557,  85.0143deg:   -0.04279, -87.9371deg:    0.10236,  25.5931deg
cycle-----)Swing------,Phase-------:d-Swing----,d-Phase-----:d-Vec(Lng)-,d-Vec(Ang)--
Average---)    0.07398, 175.4160deg:   -0.04211,  27.3369deg:    0.37479, 141.3271deg
          )---Standard Deviation---:    0.29994, 178.5135deg:    0.19427, 102.5564deg
Pole data-)    0.22634, 287.0868deg

[cycléǍ_OƂāAt[GvZ̒S_
Pʂŕ\̂łB
Averagesł̐Uƈʑ́ASf[^Ńt[GvZlłB

Q񂪐UωƈʑωłBۂ̃t[GvZ
vZԂɒ[t鎖o܂񂪁AUƈʑ̕ω
w肵vZԂɍ킹ĕ␳Ă܂B
̂ߒPȍƂ͔ɈقȂĂ܂B
L̗ႾƁA_t߂؂`ŒɕωĂ邽߁A
Uωƈʑω̕ϒlEW΍ςȒlɂȂĂ܂B

EQ񂪁AUEʑ̕fʕ\ɂωxNg
ɍWŕ\̂łBQ񂪌_S
UEʑ̕ω𑪂̂ł̂ɑ΂āA
͊e_̕ω݂̂łB
Û̂ωP[XȂǂł𗘗p܂B

͂exȊŐ߂exꍇɁA
摜`Ōƌʂ`ꍇ܂B
̌ʂ̒S_߂ʂuPole datavłB
{͍ŏ@pČvẐׂłA
蔲vZȂ̂łȂ덷܂B
QlxɎgĂB
̓IɐƁA
t@Co͂ŌvZ_(̗ł͂S_)dȂȂ悤
Q_̃yAQgIсAeyA̓񓙕̌_߂܂B
̌vZ\ȓ_̑gōsAdݕtςo͂܂B
(̗ł[1,22,3],[1,22,4],[1,23,4],
[1,32,4],[1,33,4],[2,3][3,4]̂Uʂ)


=======================================================
vZœꂽđ삷

fftcyc.exe  dftcyc.exe łx傫U̐ɂ
pxAʑAU܂B
Ȑ傫ȐɉBĂꍇ܂A
傫Ȑ̐UωĂꍇ̂ŁA
Zof[^ɂ̐Ȁł
ڂׂ܂B

܂Af[^f[^␳t@C쐬܂B
eLXgt@CɎ̂悤Ƀf[^ׂ܂B

"deg=3"
359.9905,267.52574deg,6892.3386,-0.1737,0,0
719.9810,265.0823deg,+72.006,0,0,0
#

܂"deg=n"ŐU̎߂܂B
ȍ~APsɂЂƂ̊pxłB擪
upx(deg/cyc),ʑ(deg),U(萔)A,U(1W)B,U(2W)C,U(3W)Dc(ȉw肵܂)v
łB̏ꍇɍĺA_T(cycle)̓_ɂ
(A+B*T+C*T^2+D*T^3+c)*cos(*T+)ɂȂ܂B
ۂɂ́Adftcyc.exepxƌvZʂRs[y[XgāA
U12WvZċ߂Ė߂邩AȂ
Ƃ肠u,0,0,0,cvŖ߂ĊesR{ɂ܂B
ԈႦƃG[Ŏ~܂܂B

ő40sA#܂(end of file)ŏIłB
#gΈȍ~Rgɏo܂B


̂悤ɂč쐬f[^␳t@CA
fftcyc.exe,dftcyc.exeɂāuSubtract FilevƂ
w肵ŁAꂼ̑sƁA
̃f[^elpČvZs܂B
fftcyc.exéuCheckv΁Aőlŏl͕␳̒l
\܂AuCalculatevł͊e
f[^ŌvZso͂܂B
čꂽpx̐͑ł؂덷x܂
Ȃ܂B
iUpxȂǂωꍇ́Aɉʂ
o͂܂Bj

ȍ~Åpx̐ł̌vZsAUƈʑ
vZ΁AɃf[^␳t@CɏA
őlŏlɏAɏȐ
vZ鎖o܂B

܂t߂ƒxȂ̂ŁA
Kȏŕ␳D荞ŌvZďo͂A
̃f[^t@Cő삵悢Ǝv܂B

=======================================================
Reference
A.V.Aho, J.E. Hopcroft and J.D. Ullman:The Design and Analysis of Computer Algorithms
        (a:ASY̐݌vƉII,菺O,쉺_,TCGX)
=======================================================
H20/07/20:Ƃ肠[X
H20/07/22:dftcyc.exẽt@Co͂ɕωxNg(ɍW\)ǉ
H20/07/26:f[^␳@\ǉAɍ킹ČvZʏo͌`ύXAoO
H21/04/19:dftcyc.exe̐FύX,dft̒[,ɍW\,
          t@C̎ݒ,
          t@C̋󔒍s,"Pole data"vZ
=======================================================
Apr.19th,2009

ftcenter@mth.biglobe.ne.jp
http://www2s.biglobe.ne.jp/~ftcenter/
