怎么循环？求助大侠
请先查看一下你的原始数据，时间序列是否真是等差序列，如果你的时间序列是某个仿真软件生成的，那么，很有可能不是等间距的时间。
另外，你采用的是哪种离散数值迭代方法? ode45？ Euler？改进型的Euler？
后两者计算简单，但是可能会不收敛。 邱少 发表于 2016-6-15 10:50
请先查看一下你的原始数据，时间序列是否真是等差序列，如果你的时间序列是某个仿真软件生成 ...
最好还是采用ode45函数，毕竟是解决数值问题的首选方法，精度也比较高！ 楼主是想要单周期的效果么？是不是某些参数的设置有问题？》 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]