matlab中的X=x*(exp(-j*pi/500)).^(n'*k);是什么意思?取自：n=-1:3;x=1:5;k=0:500;w=(pi/500)*k;X=x*(exp(-j*pi/500)).^(n'*k);magX=abs(X);angX=angle(X);subplot(2,1,1);plot(w/pi,magX);title('幅度响应');grid;ylabel('幅度');xlabel('以\pi为单

Question

matlab中的X=x*(exp(-j*pi/500)).^(n'*k);是什么意思?
取自：
n=-1:3;x=1:5;
k=0:500;w=(pi/500)*k;
X=x*(exp(-j*pi/500)).^(n'*k);
magX=abs(X);angX=angle(X);
subplot(2,1,1);plot(w/pi,magX);title('幅度响应');grid;
ylabel('幅度');xlabel('以\pi为单位的频率');
subplot(2,1,2);plot(w/pi,angX);title('相位响应');grid;
ylabel('相位/\pi');xlabel('以\pi为单位的频率');

Answer 1

X是x序列的离散傅里叶变换.
原本应该这么编：
for k=0:500
for n=-1:3
X(k+1)=X(k+1)+x(n+2)*(exp(-1j*pi/500)).^(n*k);
end
end
但原程序比较牛,直接用矩阵相乘代替了向量求和.