求助：加速度计算位移问题

ldhltw

n=512;
n1=0:n-1;
n2=1:n;
t=n1*0.2;
%生成正弦信号
load d:\a512.txt;
A=a512;
figure(1);
subplot(231);
plot(t,A);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
grid;
%进行FFT变换并做频谱图
y=fft(A,n);%进行fft变换
mfz=abs(y);%求幅值,即是求模
mpl=2*pi*n1'*5/512;%圆频率
mxw=atan(imag(y)./real(y));%初相位
mazh=(mfz'.*(cos((mpl'.*n1*0.2)'+mxw))')';%变换后加速度值
figure(1);
subplot(232);
plot(t,mazh);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
grid;
dds(1)=0
ddx(1)=0
for j=2:n
    %计算位移幅值
    dds(j)=cos(mxw(j)-pi)*mfz(j)/(mpl(j)^2);
    %计算位移初相
    ddx(j)=sin(mxw(j)-pi)*mfz(j)/(mpl(j)^2);
end
dds=dds;
ddx=ddx;
for j=1:n
    %位移表达式
    dwy(j)=dds(j)+ddx(j)*i;
end
dwy=dwy';
d=ifft(dwy);%逆变换求位移
figure(1);
subplot(233);
plot(d);
最后给出的位移图不正确，请高手帮忙看看问题所在。

ChaChing

可否说说!
请勿一帖多发! 谢谢