exa2_11.m

来自「新版Matlab 7辅助信号处理技术与应用」· M 代码 · 共 19 行

% 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变换结果')