%《數字信號處理教程——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))          % 校驗差值