exa2_10.m

新版Matlab 7辅助信号处理技术与应用

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

clear all;
N=25;
f=1/16;
x=sin(2*pi*f*[0:N-1]);
y=hilbert(x);
% 求 x 的希尔伯特变换；

subplot(211)
stem(x,'.');
hold on;
plot(zeros(size(x)));
title('原始序列')
subplot(212)
stem(imag(y),'.');        
% imag(y) is the Hilbert transform of x(n);
title('序列Hilbert变换结果')