% exa2-11_hilbert.m, for example 2-11
% to test hilbert.m,and to find the Hilbert transform of signal x(n)

clear all;
N=16;
n=0:N-1;
xn=sin(0.3*pi*n);

% 求 x 的希爾伯特變換；
y=hilbert(xn);
yn=imag(y)
% imag(y) is the Hilbert transform of x(n);

% 驗證正交性
add=sum(xn.*yn)

stem(yn,'.');  
title('序列Hilbert變換結果')