exa2_10.m

来自「离散信号及其Matlab实现,包括单位冲激序列 ,单位阶跃序列,斜坡序列等」· M 代码 · 共 20 行

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