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

syms t;
f = sin(t);
H = htrans(f)

H = $$-\mathrm{потому\; что}\left(x\right)$$

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

syms f(x) H(u);
f(x) = sinc(x);
H(u) = htrans(f,u)

H(u) = 
$$-\frac{\frac{\mathrm{потому\; что}\left(\pi \hspace{0.17em}u\right)}{u}-\frac{1}{u}}{\pi }$$

Постройте sinc функционируйте и его преобразование Гильберта.

fplot(f(x),[0 6])
hold on
fplot(H(u),[0 6])
legend('sinc(x)','H(u)')

Создайте синусоиду с положительной частотой на действительном пробеле.

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

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

Применяйтесь - сдвиг фазы с 90 степенями на положительную частотную составляющую с помощью преобразования Гильберта. Задайте независимую переменную как t и переменная преобразования как u.

H = htrans(y,t,u)

H = $$-A\hspace{0.17em}\mathrm{потому\; что}\left(5\hspace{0.17em}x+20\hspace{0.17em}\pi \hspace{0.17em}u\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}}$$

H = htrans(z)

H = $$A\hspace{0.17em}{\mathrm{e}}^{-10\hspace{0.17em}x\hspace{0.17em}\mathrm{i}}\hspace{0.17em}\mathrm{i}$$

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

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

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

Вычислите соответствующий аналитический сигнал $\mathit{F}\left(\mathit{s}\right)$ использование преобразования Гильберта.

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

Вычислите мгновенную частоту $\mathit{F}\left(\mathit{s}\right)$ использование

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

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

ans = $$75$$