example2_2.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 11 行

%例2-2程序:example2_2
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);  %x1的FT
X2=x2*(exp(-j*pi/500)).^(n'*k);  %x2的FT
x=alpha*x1+beta*x2;              %x1与x2线性组合
X=x*(exp(-j*pi/500)).^(n'*k);    %x的FT
%校验
Xcheck=alpha*X1+beta*X2;         %X1与X2的线性组合
error=max(abs(X-Xcheck))         %两种方式所得结果的误差