spectrum.m

现代信号处理讲义,包含许多的例题,解答等Qinghua University publishing house modern signal processing course

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;