请教下这种相图怎么回事

a9303

在论坛上看到过类似的相图，但是都是没有解释清楚就不了了之了。还有没有研究非线性的 解释下

Frank

估计是没有收敛，或者拟用于绘图的数据中包含了收敛过程的数据

Vickyvictoria

同意楼上的说法

wzx1993

楼主是做这方面研究的吗？请问可以请教您一些问题吗？

a9303

Frank 发表于 2016-4-27 14:03
估计是没有收敛，或者拟用于绘图的数据中包含了收敛过程的数据

意思是运行的时间不够长 没达到 稳定吗。或者是没找到混沌区间？

a9303

Vickyvictoria 发表于 2016-4-28 10:17
同意楼上的说法

谢谢

a9303

wzx1993 发表于 2016-4-30 11:17
楼主是做这方面研究的吗？请问可以请教您一些问题吗？

。刚开始学习 可以相互交流经验

a9303

Vickyvictoria 发表于 2016-4-28 10:17
同意楼上的说法

谢谢回复  这个附图是积分步长为0.1的相图是这种   标题附的图积分步长是0.01的 这两种图哪一个合适

wzx1993

[url][/url]

wzx1993

wzx1993 发表于 2016-5-3 16:02
[/url]

寻求大侠帮助呐！时间-位移，相图这样子算正常吗？

Vickyvictoria

a9303 发表于 2016-5-3 10:22
谢谢回复  这个附图是积分步长为0.1的相图是这种   标题附的图积分步长是0.01的 这两种图哪一个 ...

你找论坛的相关程序看看就知道了，绘图时都不是从零时刻开始的，都是从迭代一定的部署之后才开始取点的，这主要是就为了去掉未收敛的数据

Vickyvictoria

wzx1993 发表于 2016-5-3 16:03
寻求大侠帮助呐！时间-位移，相图这样子算正常吗？

单纯看图无法确认，但是你这个图和常规的差距还是不小的，最好能贴个代码看看，贴别是绘图部分的

wzx1993

Vickyvictoria 发表于 2016-5-20 10:43
单纯看图无法确认，但是你这个图和常规的差距还是不小的，最好能贴个代码看看，贴别是绘图部分的

%FFT图%%%%%%%%%%%%%%%%%%%%5
    Fs = 1/(t(2)-t(1));
    N = length(x(:,1));
    n = 0:N-1;
   
    y = fft(x(:,1),N);
    mag = abs(y);
    Fre = n*Fs/N;
    figure
    plot(Fre(1:N/2),mag(1:N/2));%x的FFT
    xlabel('frequency');
    ylabel('x amplitude');
   
    y = fft(x(:,3),N);
    mag = abs(y);
    Fre = n*Fs/N;
    figure
    plot(Fre(1:N/2),mag(1:N/2));%y的FFT
    xlabel('frequency');
    ylabel('y amplitude');
   
    y = fft(x(:,5),N);
    mag = abs(y);
    Fre = n*Fs/N;
    figure
    plot(Fre(1:N/2),mag(1:N/2));%z的FFT
    xlabel('frequency');
    ylabel('z amplitude');
   
    y = fft(x(:,7),N);
    mag = abs(y);
    Fre = n*Fs/N;
    figure
    plot(Fre(1:N/2),mag(1:N/2));%thetaY的FFT
    xlabel('frequency');
    ylabel('thetaY amplitude');
   
    y = fft(x(:,9),N);
    mag = abs(y);
    Fre = n*Fs/N;
    figure
    plot(Fre(1:N/2),mag(1:N/2));%thetaZ的FFT
    xlabel('frequency');
    ylabel('thetaZ amplitude');
   

%相图%%%%%%%%%%%%%%%%%%%%%
    figure
    plot(x(:,1),x(:,2),'LineWidth',2);%x方向相图
    xlabel('x');
    ylabel('dx/dt');
   
    figure
    plot(x(:,3),x(:,4),'LineWidth',2);%y方向相图
    xlabel('y');
    ylabel('dy/dt');
   
    figure
    plot(x(:,5),x(:,6),'LineWidth',2);%z方向相图
    xlabel('z');
    ylabel('dz/dt');
   
    figure
    plot(x(:,7),x(:,8),'LineWidth',2);%thetaY相图
    xlabel('thetaY');
    ylabel('dthetaY/dt');
   
   
    figure
    plot(x(:,9),x(:,10),'LineWidth',2);%thetaZ相图
    xlabel('thetaZ');
    ylabel('dthetaZ/dt');
   
    现在感觉fft和相图还有问题，您能帮助看看吗？

Edinburgh

wzx1993 发表于 2016-5-27 10:35
%FFT图%%%%%%%%%%%%%%%%%%%%5
    Fs = 1/(t(2)-t(1));
    N = length(x(:,1));

这程序咋不用循环

wzx1993

Edinburgh 发表于 2016-5-27 13:13
这程序咋不用循环

怎么循环？求助大侠