spectrum.m

来自「现代信号处理讲义,包含许多的例题,解答等Qinghua University p」· M 代码 · 共 8 行

function y = spectrum(k);
temp = 0;
load mysave;
n = length(x);
for i = -(n-1):n-1
   temp = temp +autocorrelation(i)*exp(-2*pi*j*i*k/n);
end
y = temp;