ihtrans

逆ヒルベルト変換

構文

f = ihtrans(H)

f = ihtrans(H,transVar)

f = ihtrans(H,var,transVar)

説明

f = ihtrans(H) はシンボリック関数 H の逆ヒルベルト変換を返します。既定では、独立変数は x、変換変数は t です。

f = ihtrans(H,transVar) は、t の代わりに transVar を変換変数として使用します。

f = ihtrans(H,var,transVar) は、x と t の代わりに var と transVar をそれぞれ独立変数および変換変数として使用します。

• すべての入力引数が同じサイズの配列である場合、ihtrans は要素単位で動作します。

• 1 つの入力がスカラーで他が同じサイズの配列の場合、ihtrans はスカラーを同じサイズの配列に拡張します。

• f が異なる独立変数を持つシンボリック式の配列である場合、var は、独立変数に対応する要素を持つシンボリック配列でなければなりません。

例

シンボリック式の逆ヒルベルト変換

cos(x) の逆ヒルベルト変換を計算します。既定では、この逆変換は t の関数を返します。

syms x;
f = cos(x);
H = ihtrans(f)

H = $$-\sin \left(t\right)$$

sinc 関数の逆ヒルベルト変換

関数 sinc(t) の逆ヒルベルト変換を計算します。これは sin(pi*t)/(pi*t) と等しくなります。s の関数として結果を表現します。

syms H(t) f(s);
H(t) = sinc(t);
f(s) = ihtrans(H,s)

f(s) = 
$$\frac{\frac{\cos \left(\pi \hspace{0.17em}s\right)}{s}-\frac{1}{s}}{\pi }$$

関数 sinc とその逆ヒルベルト変換をプロットします。

fplot(H(t),[0 6],'b')
hold on
fplot(f(s),[0 6],'r')
legend('sinc(t)','f(s)')

位相シフトの適用

正の周波数を持つ正弦波を実空間に作成します。

syms A x t u;
assume([x t],'real')
H = A*sin(2*pi*10*t + 5*x)

H = $$A\hspace{0.17em}\sin \left(5\hspace{0.17em}x+20\hspace{0.17em}\pi \hspace{0.17em}t\right)$$

逆ヒルベルト変換を使用して、90 度の位相シフトを正の周波数成分に適用します。独立変数を x、変換変数を u としてそれぞれ指定します。

f = ihtrans(H,x,u)

f = $$A\hspace{0.17em}\cos \left(5\hspace{0.17em}u+20\hspace{0.17em}\pi \hspace{0.17em}t\right)$$

ここで負の周波数を持つ複素信号を作成します。逆ヒルベルト変換を使用して、-90 度の位相シフトを負の周波数成分に適用します。

Z = A*exp(-1i*10*t)

Z = $$A\hspace{0.17em}{\mathrm{e}}^{-10\hspace{0.17em}t\hspace{0.17em}\mathrm{i}}$$

f = ihtrans(Z)

f = $$-A\hspace{0.17em}{\mathrm{e}}^{-10\hspace{0.17em}u\hspace{0.17em}\mathrm{i}}\hspace{0.17em}\mathrm{i}$$

瞬時周波数の計算

60 Hz と 90 Hz の 2 つの周波数成分を持つ実数値信号 $$f(s)$$ を作成します。

syms s f(x) F(t)
f(s) = sin(2*pi*60*s) + sin(2*pi*90*s)

f(s) = $$\sin \left(120\hspace{0.17em}\pi \hspace{0.17em}s\right)+\sin \left(180\hspace{0.17em}\pi \hspace{0.17em}s\right)$$

逆ヒルベルト変換を使用して、対応する解析信号 $$F(t)$$ を計算します。

F(t) = ihtrans(f(s),t) + 1i*f(t)

F(t) = $$\cos \left(120\hspace{0.17em}\pi \hspace{0.17em}t\right)+\cos \left(180\hspace{0.17em}\pi \hspace{0.17em}t\right)+\sin \left(120\hspace{0.17em}\pi \hspace{0.17em}t\right)\hspace{0.17em}\mathrm{i}+\sin \left(180\hspace{0.17em}\pi \hspace{0.17em}t\right)\hspace{0.17em}\mathrm{i}$$

$$F(t)$$ の瞬時周波数を、以下を使用して計算します。

$$f_{instant}(t)=\frac{1}{2\pi }\frac{d\varphi (t)}{dt},$$

ここで、$$\varphi \left(t\right)=\arg [F(t)]$$ は解析信号の瞬間位相です。

InstantFreq(t) = diff(angle(F(t)),t)/(2*pi);
assume(t,'real')
simplify(InstantFreq(t))

ans = $$75$$

入力引数

H — 入力

シンボリック式 | シンボリック関数 | シンボリック ベクトル | シンボリック行列

入力。シンボリック式、シンボリック関数、シンボリック ベクトル、またはシンボリック行列として指定します。

var — 独立変数

x (既定値) | シンボリック変数 | シンボリック ベクトル | シンボリック行列

独立変数。シンボリック変数、シンボリック ベクトルまたはシンボリック行列として指定します。この変数は通常、時間領域にあります。変数を指定しない場合、ihtrans は既定で x を使用します。H が x を含まない場合、ihtrans は関数 symvar を使用して独立変数を決定します。

transVar — 変換変数

t (既定値) | u | シンボリック変数 | シンボリック ベクトル | シンボリック行列

変換変数。シンボリック変数、シンボリック ベクトルまたはシンボリック行列として指定します。この変数は、var と同じ領域にあります。変数を指定しない場合、ihtrans は既定で t を使用します。t が H の独立変数である場合、ihtrans は変換変数 u を使用します。

出力引数

f — H の逆ヒルベルト変換

シンボリック式 | シンボリック関数 | シンボリック ベクトル | シンボリック行列

入力関数 H の逆ヒルベルト変換です。出力 f は transVar で指定される変数の関数です。

ihtrans が入力関数を変換できない場合は、未評価の呼び出しが返されます。元の式を返すには、htrans を使用して出力にヒルベルト変換を適用します。

詳細

逆ヒルベルト変換

変数 x の式 H = H(x) の点 t における逆ヒルベルト変換 f = f(t) は次のようになります。

$$f\left(t\right)=\frac{1}{\pi }\text{p}\text{.v}\text{.}{\displaystyle \underset{-\infty }{\overset{\infty }{\int }}\frac{H\left(x\right)}{x-t}dx.}$$

ここで、p.v. はコーシー主値積分を表します。関数 H(x) は複素数でもかまいませんが、x および t は実数でなければなりません。

ヒント

• ヒルベルト変換を計算するには、htrans を使用します。関数の逆ヒルベルト変換は、そのヒルベルト変換の負の値と等しくなります。

• 時間領域内の信号に対しては、逆ヒルベルト変換は 90 度の位相シフトを対応するフーリエ成分の負の周波数に適用します。正の周波数に対する -90 度の位相シフトも適用します。

• 実数値の信号 b は、その逆ヒルベルト変換 a = ihtrans(b) の調和共役です。逆ヒルベルト変換 a = real(z) と信号 b = imag(z) は、解析信号 z = a + 1i*b を形成します。

バージョン履歴

R2019a で導入

参考

htrans | fourier | ifourier | laplace | ilaplace