我的振动加速度频域积分的结果怎么是这样？请大家帮忙看看，谢谢

卡卡卡帅

频域积分的程序来源于王济教授的书本和论坛，数据为发动机1900rpm振动数据。见附件
发动机振动理想中应该是在0值附近的一个强迫性周期振动，但是计算出来的结果飘的很远，这是什么原因？
fmin=1的设定是根据传感器设定的。
程序如下：
clear;
load top6
C=C(1:10240);
fni=C;
i=(-1)^0.5;
sf=10240;                          
fmin=1;                        
fmax=5120;                        
c=1;                             
it=2;                             
x=fni';                             
n=length(x);

t=0:1/sf:(n-1)/sf;   
subplot(211);
plot(t,x);
nfft=2^nextpow2(n);                 
y=fft(x,nfft);                     
df=sf/nfft;                        
ni=round(fmin/df+1);               
na=round(fmax/df+1);               
dw=2*pi*df;                        
w1=0:dw:2*pi*0.5*sf;           
w2=-2*pi*(0.5*sf-df):dw:-dw;         
w=[w1,w2];
w=w.^it;                           
a=zeros(1,nfft);
a(2:nfft-1)=y(2:nfft-1)./w(2:nfft-1);        
if it == 2
   y=-a;                              
else
   a1=imag(a); a2=real(a); y=a1-a2*i;
end
a=zeros(1,nfft);
a(ni:na)=y(ni:na);                     

a(nfft-na+1:nfft-ni+1)=y(nfft-na+1:nfft-ni+1);            

y=ifft(a,nfft);
y=real(y(1:n))*c;
a=y;
subplot(2,1,2);
plot(t,a);
xlabel('时间（s）');
ylabel('位移（m）');

westrongmc

1900rpm对应的转频为31.67Hz，可以将高通频率设的高些，例如5Hz，再试一下。

卡卡卡帅

westrongmc 发表于 2015-2-6 13:18
1900rpm对应的转频为31.67Hz，可以将高通频率设的高些，例如5Hz，再试一下。

哦，谢谢！
低频的值不准确么？