FAQ
FAMOS
シーケンス / 関数top
└シーケンス / 関数
ID.s029
Q. ヒルベルト変換の関数はありますか？
A.
Ver4.0の時点では存在しません。メーカーに問い合わせたところ、バージョンアップで対応するという回答をもらいました。
バージョンアップでサポートされるまで以下のシーケンスで代用してください。
;---------------------------------------------
; ヒルベルト変換
;---------------------------------------------
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; サンプルデータの読み込み
Famos
Load C:\imc\dat\slope.dat __data
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 演算
FftOption 0 0
_fan = fft( __data )
_fan.p = _fan.p - 90
_t1 = ifft ( _fan )
hilbert = Sqrt(Sqr(__data)+Sqr(_t1)) < 包落線
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; 瞬間周波数
angle = atan2(_t1, __data)
angle1 = cutIndex(angle, 1, leng?(angle)-1)
angle2 = cutIndex(angle, 2, leng?(angle))
angle2=xoff(angle2,0)
angle = angle2-angle1
i=1
while i<=leng?(angle)
if angle[i]<-pi
angle[i]=angle[i]+pi2
end
if angle[i]>pi
angle[i]=angle[i]-pi2
end
i=i+1
end
Freq = angle/pi2
delete _t1
delete _fan
