hc326.m

来自「数字信号处理教程：MATLAB释义与实现 书中实例码源」· M 代码 · 共 15 行

%《数字信号处理教程——MATLAB释义与实现》第三章例3.2.6程序hc326
% DTFT的线性性质检验
% 电子工业出版社出版  陈怀琛编著 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
% 校验
Xc = alpha*X1 + beta*X2;       % 把X1 和 X2 的线性组合称为Xc
error = max(abs(X-Xc))          % 全频段内最大误差绝对值