请教下这种相图怎么回事
在论坛上看到过类似的相图，但是都是没有解释清楚就不了了之了。还有没有研究非线性的 解释下估计是没有收敛，或者拟用于绘图的数据中包含了收敛过程的数据 同意楼上的说法 楼主是做这方面研究的吗？请问可以请教您一些问题吗？ Frank 发表于 2016-4-27 14:03
估计是没有收敛，或者拟用于绘图的数据中包含了收敛过程的数据
意思是运行的时间不够长 没达到 稳定吗。或者是没找到混沌区间？ Vickyvictoria 发表于 2016-4-28 10:17
同意楼上的说法
{:{23}:} 谢谢 wzx1993 发表于 2016-4-30 11:17
楼主是做这方面研究的吗？请问可以请教您一些问题吗？
。刚开始学习 可以相互交流经验{:{05}:} Vickyvictoria 发表于 2016-4-28 10:17
同意楼上的说法
{:{05}:} 谢谢回复这个附图是积分步长为0.1的相图是这种 标题附的图积分步长是0.01的 这两种图哪一个合适{:{14}:} wzx1993 发表于 2016-5-3 16:02
寻求大侠帮助呐！时间-位移，相图这样子算正常吗？ a9303 发表于 2016-5-3 10:22
谢谢回复这个附图是积分步长为0.1的相图是这种 标题附的图积分步长是0.01的 这两种图哪一个 ...
你找论坛的相关程序看看就知道了，绘图时都不是从零时刻开始的，都是从迭代一定的部署之后才开始取点的，这主要是就为了去掉未收敛的数据 wzx1993 发表于 2016-5-3 16:03
寻求大侠帮助呐！时间-位移，相图这样子算正常吗？
单纯看图无法确认，但是你这个图和常规的差距还是不小的，最好能贴个代码看看，贴别是绘图部分的 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和相图还有问题，您能帮助看看吗？ wzx1993 发表于 2016-5-27 10:35
%FFT图%%%%%%%%%%%%%%%%%%%%5
Fs = 1/(t(2)-t(1));
N = length(x(:,1));
这程序咋不用循环 Edinburgh 发表于 2016-5-27 13:13
这程序咋不用循环
怎么循环？求助大侠
页:
[1]
2