Вычислите обратное Гильбертово преобразование cos(x). По умолчанию обратное преобразование возвращает функцию t.

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

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

Вычислите обратное Гильбертово преобразование функции sinc, которая равна sin(πt)/πt. Выразите результат как функцию s.

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

f(s) = 
$$\frac{\frac{\mathrm{потому\; что}\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}\mathrm{потому\; что}\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}{e}^{-10\hspace{0.17em}t\hspace{0.17em}i}$$

f = ihtrans(Z)

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

Создайте сигнал с действительным знаком $$f(s)$$ с двумя частотными составляющими, 60 Гц и 90 Гц.

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) = $$\mathrm{потому\; что}\left(120\hspace{0.17em}\pi \hspace{0.17em}t\right)+\mathrm{потому\; что}\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}i+\sin \left(180\hspace{0.17em}\pi \hspace{0.17em}t\right)\hspace{0.17em}i$$

Вычислите мгновенную частоту $$F(t)$$ использование

где $$\varphi \left(t\right)=\mathrm{аргумент}[F(t)]$$ мгновенная фаза аналитического сигнала.

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

ans = $$75$$