example2_10.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

x1=rand(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);
X2=x2*(exp(-j*pi/500)).^(n'*k);
x=alpha*x1+beta*x2;
X=x*(exp(-j*pi/500)).^(n'*k);
error=max(abs(X-alpha*X1-beta*X2));
disp(error);