清空筛选条件
清空筛选条件

Plotting Power spectral density of Simple Rossler Attractor

1 次查看（过去 30 天）
ABHIJIT MAJUMDER
ABHIJIT MAJUMDER 2017-11-5
I have been trying to plot a figure of the power spectrum for Rossler system which I found in a paper entitled 'Power spectra and mixing properties of strange attractors' by Doyne Farmer. I wrote a MATLAB code to generate the figure but not getting the plot exactly what is in that paper. Here is my code-
close all;clc;
global a b c
a=0.2;
b=0.4;
c=5.7;
ts=0:.01:4000;
z0=[1 0 1];
opt=odeset('RelTol',10e-12);
[t,z]= ode45('System', ts, z0,opt);
Fs = 1000;
x=z(:,1);
nfft = 2^nextpow2(length(x));
Pxx = abs(fft(x,nfft)).^2/length(x)/Fs;
Hpsd = dspdata.psd(Pxx(1:length(Pxx)/2),'Fs',Fs/10);
figure()
plot(Hpsd)
subject to the function
function zdot=System(t,z)
global a b c
zdot=[ -(z(2)+z(3));
z(1)+a*z(2);
b+z(1)*z(3)-c*z(3)];
I, hereby, attach the paper also for your perusal. Please solve the error in my code.

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Parametric Spectral Estimation 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by