|
楼主 |
发表于 2007-4-25 09:57
|
显示全部楼层
本帖最后由 牛小贱 于 2014-6-30 23:57 编辑
恩,谢谢 gghhjj 和 无水1324 两位前辈的指导,对于一些基本问题已经有所了解,按理论来将这个图对应的运动应该是周期运动而不是拟周期运动,而且poincare图也就是那样的,没有未表示出的点,从这点来说也应该是周期运动而不是拟周期运动,但是,现在问题的关键是频率幅值图,是否是有问题?因为我看过人家和我取的同样的参数的时候出来的频率幅值图,并没有这么多细小的连续的频谱,所以现在很纳闷,不知道该如何解释这幅频图?再次感谢两位的指导!附上频率幅值图的程序!大家帮分析下,也许是程序出错导致的,谢谢了!
幅值频率图程序:
- t0=0;tf=700;
- t=[0:4:tf];
- x0=[0.001 0 0.001 0];
- options = odeset('RelTol',1e-6,'AbsTol',[1e-6 1e-6 1e-6 1e-6]);
- [T,X]=ode45('111',t,x0,options);
- fs=200; %采样频率
- N=140001; %采样点数
- x=X(:,2); %生成信号(Y轴方向的位移)
- figure(1)%subplot(2,1,1);
- plot(T,x);
- xlabel('t/s');ylabel('Y');
- title('时间位移图');
- axis([0 700 -3 3]);
- xf=fft(x,N);
- mag=abs(xf);
- df=fs/N; %频率分辨率Hz
- %绘制双边幅值谱
- f=(0:N-1)*df; %频域序列
- figure(2)%subplot(2,1,2);
- plot(f,abs(xf));
- xlabel('Frequency(HZ)');
- ylabel('Magnitude');
- title('频率幅值图');
- axis([0 100 0 120]);
- %N=fs*t
- hold off
复制代码
|
|