这样的频谱图能说明什么问题？

胡晓宇

附件为我做的频谱图，从图上看不出主要的振动频率，感觉很乱，不知道是做法的问题还是振动的性质本身造成的？

007007boy

可能是混沌了吧。

hsfy919

个人感觉应该是混沌了，不知道你的系统是不是非自治的，看看庞加莱截面就知道了

ChaChing

LZ这图在实物上量测信号时好像时常遇到! 一般都是取样多次再平均处理!
非线性还在与各位学习!? 混沌? why?

[ 本帖最后由 无水1324 于 2010-3-22 17:48 编辑 ]

hsfy919

出现连续谱，说明响应中含有相当多周期分量存在，这时庞加莱截面由应该是许多无规则的点构成，也即微分方程的解失去了周期性，根据庞加莱理论，该系统可能是混沌，不知理解的是否正确，望指正

胡晓宇

看到楼上各位的回复，很感动，也很感谢。我做的是两个转子的动力学分析，的确有很多不同周期的激励作用，实在是不知道如何分析，找不出运动的特点，一直在迷茫

无水1324

理论来说，混沌和含噪音都会出现连续谱。同时你这里有多个谱线，说明里面还有多周期的振动信号。而且三个信号的频率之间也没有出现明显的倍数关系。
如果是理论分析的话，应该就是混沌的了。
如果是实际的检测信号，可能还要考虑噪音的影响。有可能是多周期运动。

ydiok

转子的转频是多少

胡晓宇

转速分别为12000rpm和9600rpm，好像和图也对应不上

胡晓宇

function myfrequencyplot
% all the files are save as t and s
clc;
clear;
load f4_n1_12000_n2_9600_5.00.mat;
% 读取文件，里面保存有时间序列t和积分结果s，s有18维
disp('three frequencies');

N=1024*4;%选择点数
xtemp1=s(end-N*4:4:end,1);%转子1在x方向振幅

dt = 6/100000;%积分时，积分区间为【0，6】，点数为100000
xa=fft(xtemp1,N);%转子1在x方向位移的傅里叶变化
xa(1)=[];

f = 0.25/dt*(0:N/2)/N;
figure;
set (gcf,'Position',[400,100,400,300], 'color','w')
plot(f(1:N/2),abs(xa(1:N/2)));
h = gca;
set(h,'FontSize',14);
xlabel('Frequency( Hz)')
ylabel('Vertical Amplitude(m)')