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

製品詳細はこちら

＜＜ シーケンス / 関数に関するFAQ一覧へ戻る