对位移做频谱分析，如何让它的工频为1
程序如下：
% 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()
% %
% 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()
频域如此干净，将你的时域信号打出来看看。
本帖最后由 hcharlie 于 2015-1-13 20:10 编辑
感谢能打出时域图。图上500毫秒39.5个波，79Hz，而20000/60=333Hz，时域图能对上吗？ 对不上，我也不知道怎么解决了 hcharlie 发表于 2015-1-13 20:08
感谢能打出时域图。图上500毫秒39.5个波，79Hz，而20000/60=333Hz，时域图能对上吗？
不好意思，我的问题该怎么解决，请指教 本帖最后由 hcharlie 于 2015-1-13 22:12 编辑
你的位移振幅近似=1，单位是什么？比如1mm，频率79Hz时加速度可以算出，如果频率是333，频率高4倍多，加速度应该大16倍多。用加速度验证也是方法。
科学问题带不得一点含糊的，你的含糊处太明显了。 没有单位，这是无量纲化之后的位移
页:
[1]