请教下这种相图怎么回事

邱少

wzx1993 发表于 2016-5-27 21:39
怎么循环？求助大侠

            请先查看一下你的原始数据，时间序列是否真是等差序列，如果你的时间序列是某个仿真软件生成的，那么，很有可能不是等间距的时间。
        另外，你采用的是哪种离散数值迭代方法?   ode45？ Euler？改进型的Euler？
         后两者计算简单，但是可能会不收敛。

Pparis

邱少 发表于 2016-6-15 10:50
请先查看一下你的原始数据，时间序列是否真是等差序列，如果你的时间序列是某个仿真软件生成 ...

最好还是采用ode45函数，毕竟是解决数值问题的首选方法，精度也比较高！

我不是小偷

楼主是想要单周期的效果么？是不是某些参数的设置有问题？》

Edinburgh

wzx1993 发表于 2016-5-27 21:39
怎么循环？求助大侠

你看你程序中的前一部分  1 3 5 7 9  
可以用
for i=1:2:9
y = fft(x(:,i),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');
end
就类似于这样