最近在用matlab处理振动数据，加速度频域积分为位移，但是我不明白为什么最后是用逆变换的实部作为结果？

以上程序是对加速度二次积分为位移，将时域积分装换为频域积分， 实现方法为将时域信号进行FFT到频域，然后毁知对频谱做简单运算实现了频域上的积分 最后将频域经过IFFT到时域 实现裂培了频域上的积分。

时域积分对应在频域上的频谱函数为：∫f(x)dx →F(jw)/jw

所以时域上的一次积分对应在频域上的处理为：逆完FFT之后取纤源消虚部,因为加速度与速度之间的关系为jw

二次积分时因为是-w^2，所以有 y=-a 即为取实部

假设以矩形或梯形面积来拟合吧带慎键；

function yanshi(fname,a,b,n)

%定积分蠢巧孝卜演示程序

xi(1)=a

xi(n+1)=b

for i=1:n-1

xi(i+1)=a+(i+rand(1))*(b-a)/n

end

I=0

for i=1:length(xi)-1

sxi=xi(i)+rand(1)*(xi(i+1)-xi(i))

syi=feval(fname,sxi)

I=I+syi*(xi(i+1)-xi(i))

xx=[xi(i) xi(i) xi(i+1) xi(i+1) xi(i)]

yy=[0 syi syi 0 0]

fill(xx,yy,'c')

hold on

end

x=a:(b-a)/100:b

y=feval(fname,x)

plot(x,y,'r','markersize',20)

fprintf('n=%6d I=%12.5f\n',n,I)

hold off

function y=lirui(x)

y=1+x.^2

t = (1/16:1/8:1)'*2*pi

x = sin(t)

y = cos(t)

fill(x,y,'r')

axis square

---