風の周期性を調べてるっす

ぺんぎん

2009年08月15日 21:50


やってることは地味っす。
コレによってできるモノも・・・地味っす。
派手さは不要なんっすけどね。
どうもぺんぎんっす( ◎v◎ )


風って周期的に変化してるんっすかねぇ?
というわけで、調べてみてるっす。

llWindで返ってきたvectorのx要素,y要素をデータの実部、虚部として、
FFTに入れてみてるっす。

データ数512のFFTはヒープ・コリジョンが出なかったっす。
(Monoコンパイルで実部512要素、虚部512要素のlistを入力)
でも結果をメールで送信するとなると、「本文4K制限」にかかって、
データ数128が限界っす。
データ数512でやって、高周波成分をカットしたものを送信するっていう
手もあるっすね。


データ数128、10.0秒間隔で何回か実験してみたんっすけど、
低周波成分はもちろん残ってるんっすけど、高周波成分もかなり出てるっす。
もうちょっとやると面白い結果が得られそうっす。
「続きを読む」で1回分のFFT出力が見れるっす。
かなりの量っすから、覚悟して見てくださいっす。
10.0秒間隔
データ数 = 128

*注意*
逆変換するときにはスケーリングを考慮しないとダメっす。
自分の方法では128で割ればOKっす。

[実部]
-19.123089
-168.701233
66.628738
7.820580
-14.401188
-23.317215
-1.168656
-51.201679
-65.813614
-73.056374
-33.751213
-44.239685
-20.395493
-21.758551
-12.087668
-6.507188
-12.010941
-19.253149
-16.367975
-5.679879
-9.932274
-9.063625
-6.747901
-7.325851
-6.436904
-9.481175
-5.032416
-3.742121
-2.429868
-6.568584
-4.257585
-3.572488
-4.487659
-2.577536
-3.490913
-2.353604
-3.793070
-2.061213
-3.622206
-3.656405
-2.286545
-2.094496
-2.128628
-1.417548
-2.220318
-1.795951
-1.442102
-2.042731
-1.566714
-0.672022
-0.617702
-1.575737
-1.357228
-0.275834
-1.191645
-1.453022
-0.453846
0.449478
-0.626480
0.480705
0.340487
-0.335506
0.533697
0.255777
-0.275946
0.760658
0.322296
-0.252928
1.428037
1.322753
0.257809
0.019627
1.264328
1.279800
1.130836
1.980978
1.258453
1.888266
1.476460
2.592611
2.308464
2.207407
1.863877
2.655758
2.667846
2.534003
1.779533
3.164512
3.795562
4.046638
1.615408
2.710649
3.989767
4.068394
4.412798
5.075419
2.672479
4.277273
6.627113
7.310171
5.661512
3.533817
3.581910
7.448327
8.669477
7.427420
8.153700
6.150594
5.352184
5.218486
6.957374
17.198425
4.878880
16.297232
29.561916
-1.833489
-3.838338
1.526979
50.274208 ←ここから妙っす
89.566711
-76.201981
-28.835209
149.855774
142.880219
12.352022
-85.310120
-20.152678
16.657183
-77.544746
----------------
[虚部]
-77.544746
-111.998383
-204.805695
25.322128
-151.024506
18.810205
-29.359694
-82.801674
-57.430138
8.781224
19.766659
-22.603893
11.679066
12.978252
-1.455993
-6.672466
-4.962615
-7.294288
0.594550
6.182561
-2.777017
-1.604712
-0.179398
-2.961735
-6.376828
-1.275909
-0.045299
-1.049698
-5.193352
-4.222267
-2.583510
-4.145226
-2.973887
-2.904604
-3.659576
-4.407068
-3.753385
-2.937154
-4.026249
-3.536855
-3.353163
-3.432064
-4.575411
-4.208439
-2.845405
-4.354914
-3.673001
-3.806163
-3.944153
-4.292192
-4.286702
-4.538240
-3.327678
-3.052016
-4.542509
-3.252335
-3.552122
-3.403915
-5.388948
-3.772105
-3.325432
-3.636894
-4.224174
-3.221363
-3.974449
-4.061340
-3.928154
-3.630566
-4.288239
-3.996410
-3.944281
-3.924919
-2.654472
-3.785473
-3.654275
-2.393379
-3.700016
-4.482233
-2.714595
-3.484884
-4.467083
-3.883173
-4.161883
-3.572104
-4.367856
-3.039281
-3.346672
-2.400744
-2.477700
-5.152691
-4.959112
-1.609343
-1.907921
-2.739007
-2.537477
-4.693206
-2.412023
-1.107012
-1.908753
-3.992444
-6.186122
-4.403646
-0.612587
-0.531744
-2.524146
-5.544332
-1.421366
-6.128418
-2.087944
-1.606941
7.722182
-2.570473
-3.873180
6.023025
-12.607078
-22.348890
9.945086
25.591652 ←ここから妙っす
32.646225
-43.564278
-3.123683
122.679443
-31.000799
5.834145
-216.455750
-150.898468
-162.181030
-37.807247
スクリプト