%《數(shù)字信號處理教程——MATLAB釋義與實(shí)現(xiàn)》第三章例3.2.6程序hc326
% DTFT的線性性質(zhì)檢驗(yàn)
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
x1 = randn(1,11); x2 = rand(1,11); n = 0:10;
alpha = 2; beta = 3;
k = 0:500; w = (pi/500)*k;
X1 = x1 * (exp(-j*pi/500)).^(n'*k);  % x1 的 DTFT
X2 = x2 * (exp(-j*pi/500)).^(n'*k);  % x2 的 DTFT
x = alpha*x1 + beta*x2;              % x1 和 x2 的線性組合
X = x * (exp(-j*pi/500)).^(n'*k);    % x 的 DTFT
% 校驗(yàn)
Xc = alpha*X1 + beta*X2;       % 把X1 和 X2 的線性組合稱為Xc
error = max(abs(X-Xc))          % 全頻段內(nèi)最大誤差絕對值