cw_evaluate

内容

連続ウェーブレット変換の計算

これはOriginProのみの機能です。

1. cw_evaluate type:=dgauss param:=3 lb:=-8 ub:=8 npts:=200

2. cw_evaluate type:=morlet param:=8

3. cw_evaluate type:=maxhat npts:=50 oy:=(col(1), col(2))

表示名	変数名	I/O とデータ型	デフォルト値	説明
ウェーブレットの種類	type	入力 int	morlet	ウェーブレットの種類を指定します。オプションリスト morlet:Morlet Morletウェーブレット dgauss:DGauss Gaussianウェーブレットの導関数 mexhat:MexHat Mexican Hatウェーブレット
波数	param	入力 int	5	Morletウェーブレットの場合、この変数は、5より小さくなく、20より大きくないウェーブ数を指定します。DGaussウェーブレットの場合、この変数は、9より小さい正の正数であるような導関数の次数を指定します。
下側境界	lb	入力 double	10	評価するポイントの下側境界を指定します。
上側境界	ub	入力 double	10	評価するポイントの上側境界を指定します。
ウェーブレット関数のポイント数	npts	入力 int	100	評価するポイント数を指定します。
ウェーブレット関数	oy	出力 XYRange	<新規>	出力範囲を指定します。

このXファンクションは、与えられた評価ポイントにおいて選択したウェーブレットを計算します。評価ポイントの数はユーザが指定します。そして評価ポイントの範囲は、上側境界と下側境界で決まります。

この関数を使って、Morlet, Mexican Hat, Gaussian導関数の3種類のウェーブレットを評価できます。それぞれの定義は下記をご覧ください。

Morletウェーブレットは
のように定義されます。ここで k はウェーブ数です。

Mexican Hatウェーブレットは次のようになります。

そして、Gaussianウェーブレットはガウス関数の p階の微分で、これは
のように定義されます。ここで p は微分の階数です。

1. 範囲-8から7までの150の評価ポイントに対するMexican Hatウェーブレットの値を計算するには、次のスクリプトコマンドを実行します。

cw_evaluate type:=2 lb:=-8 ub:=7 npts:=150

2. ウェーブレットの種類を変更するには、ダイアログのType ドロップダウンリストを使うこともできます。これはコマンドウィンドウに次のように入力して開くことができます。

cw_evaluate -d;

または、スクリプト行にどのウェーブレットを使用するかを指定します。例えば次のように入力します。

cw_evalute type:=1

計算は3つのNAG関数nag_morlet(), nag_dgauss(), nag_mexhat()を使って組み込まれています。

Torrence, C and Compo, G... 1998.A practical guide to wavelet analysis, Bull.Amer.Meteor.Soc.

Ingrid Daubechies.1992.Ten Lectures On Wavelets.SIAM.