谁能帮忙用matlab画出logistic混沌映射序列的功率谱图?
来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/07/11 13:28:24
谁能帮忙用matlab画出logistic混沌映射序列的功率谱图?
logistic混沌映射序列的迭代式为:x(n+1)=4*x(n)*(1-x(n)),初值x(0)为0.1.要用matlab写出完整的代码画出此序列的功率谱图.悬赏不够的话可以加.
logistic混沌映射序列的迭代式为:x(n+1)=4*x(n)*(1-x(n)),初值x(0)为0.1.要用matlab写出完整的代码画出此序列的功率谱图.悬赏不够的话可以加.
![谁能帮忙用matlab画出logistic混沌映射序列的功率谱图?](/uploads/image/z/16546302-54-2.jpg?t=%E8%B0%81%E8%83%BD%E5%B8%AE%E5%BF%99%E7%94%A8matlab%E7%94%BB%E5%87%BAlogistic%E6%B7%B7%E6%B2%8C%E6%98%A0%E5%B0%84%E5%BA%8F%E5%88%97%E7%9A%84%E5%8A%9F%E7%8E%87%E8%B0%B1%E5%9B%BE%3F)
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)|')
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)|')