hc328.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第三章例3.2.8程序hc328
% 验证DTFT的共轭性质
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
n = -5:10; x = rand(1,length(n)) + j*randn(1,length(n));
k = -100:100; w = (pi/100)*k;        % -pi 和 +pi 之间的频率
X = x * (exp(-j*pi/100)).^(n'*k);    %  x的 DTFT
y = conj(x);                         % 求信号共轭检验共轭特性
Y = y * (exp(-j*pi/100)).^(n'*k);    % y的DTFT 
Y_check = conj(fliplr(X));           % conj(X(-w))
error = max(abs(Y-Y_check))          % 校验差值