求助：碰摩转子频谱分析的问题

yiquncw

用ode45求解碰摩转子方程，得出的结果要进行频谱分析，可是得出频谱图总是在8hz和4096hz处有一峰值，
计算转速是90000rpm,对应频率1500hz,但此处没有振动幅值，请高手指点！
相关命令行如下：
w=3000*pi;  % 转速值
T=2*pi; ts=[0:T/512:128*T];
y=z(24577:end,1);  %取出的时域数据，用后40周期的数据
y=y-mean(y); fs=4096; N=length(y);
zz=fft(y); f=(0:N-1)*fs/N; mag=2*abs(zz(1:N))/N;
power=mag.^2;

[ 本帖最后由 ChaChing 于 2009-12-5 21:04 编辑 ]

ChaChing

fs=1/(T/512) ???

yiquncw

这是我的计算结果，强迫振动频率为1800hz
采样频率为4096hz，为什么在4096hz处总有一个小的峰值，2300hz处的频率又能说明什么？
计算工况为带有重力影响的强迫振动转子，没有碰摩情况发生？
望有经验的讨论下！

ChaChing

采样频率为4096hz, 仅能看到2048Hz, 后半是负频率部分
建议LZ去"信号处理方法"逛逛, 相信可以找到答案