跪求：matlab频谱分析，如何得到功率谱图？ 离散数据信号为p，长度为n。请高手指教！！！！！！！！！！！

跪求：matlab频谱分析，如何得到功率谱图？ 离散数据信号为p，长度为n。请高手指教！！！！！！！！！！！

logistic混沌映射序列的迭代式为：x(n+1)=4*x(n)*(1-x(n)),初值x(0)为0.1。要用matlab写出完整的代码画出此序列的功率谱图。
x(1)=0.1;
for i = 2 : 1000
x(i)=4*x(i-1)*(1-x(i-1));
end

L=1000;
NFFT = 1024;
Y = fft(x,NFFT)/L;
f = 1/2*linspace(0,1,NFFT/2+1);

plot(f,2*abs(Y(1:NFFT/2+1)))
title('功率谱')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

搜一下：跪求：matlab频谱分析，如何得到功率谱图？ 离散数据信号为p，长度为n。请高手指教！！！！！！！！！！！