我的振动加速度频域积分的结果怎么是这样？请大家帮忙看看，谢谢
频域积分的程序来源于王济教授的书本和论坛，数据为发动机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=;
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）');
1900rpm对应的转频为31.67Hz，可以将高通频率设的高些，例如5Hz，再试一下。 westrongmc 发表于 2015-2-6 13:18
1900rpm对应的转频为31.67Hz，可以将高通频率设的高些，例如5Hz，再试一下。
哦，谢谢！
低频的值不准确么？
页:
[1]