对位移做频谱分析，如何让它的工频为1

dj15175771160

转速n=20000,工频f=n/60=20000/60,可是f为什么是5、10有输出

程序如下：
% 2 频谱图
% fs=1200;                          %采样频率
% figure
% N1=length(Y1);                    %采样点数
% Fn1=(1:N1/2)*fs/N1;               %频率序列
% YY1=abs(fft(Y1))/N1*2;
% YYY1=YY1(1:N1/2);                  %单边普
% mag1=abs(YYY1);
% plot(Fn1,mag1,'r')
% xlabel('f');
% ylabel('幅值/Xr');
% axis([0,25,0,1])
% %
% figure
% N2=length(Y5);                    %采样点数
% Fn2=(1:N2/2)*fs/N2;               %频率序列
% YY2=abs(fft(Y5))/N2*2;
% YYY2=YY2(1:N2/2);                  %单边普
% mag2=abs(YYY2);
% plot(Fn2,mag2,'r')
% xlabel('f');
% ylabel('幅值/Xj');
% axis([0,25,0,1])

hcharlie

频域如此干净，将你的时域信号打出来看看。

dj15175771160

hcharlie

感谢能打出时域图。图上500毫秒39.5个波，79Hz，而20000/60=333Hz，时域图能对上吗？

dj15175771160

对不上，我也不知道怎么解决了

dj15175771160

hcharlie 发表于 2015-1-13 20:08
感谢能打出时域图。图上500毫秒39.5个波，79Hz，而20000/60=333Hz，时域图能对上吗？

不好意思，我的问题该怎么解决，请指教

hcharlie

你的位移振幅近似=1，单位是什么？比如1mm，频率79Hz时加速度可以算出，如果频率是333，频率高4倍多，加速度应该大16倍多。用加速度验证也是方法。
科学问题带不得一点含糊的，你的含糊处太明显了。

dj15175771160

没有单位，这是无量纲化之后的位移