邱少 发表于 2016-6-15 10:50

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

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

Pparis 发表于 2016-6-15 11:16

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

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

我不是小偷 发表于 2016-6-16 11:01

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

Edinburgh 发表于 2016-6-16 15:28

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
就类似于这样
页: 1 [2]
查看完整版本: 请教下这种相图怎么回事