关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
目的是求取VFTO波形在50MHz范围内的幅频曲线
怎么感觉结果怪怪的？

clear all;
) I4 R' @1 z/ L5 h# V& T5 Uclc;

  B- l: K: h$ x$ s/ Y* C%Reading the inout waveform
fid=fopen('wave.dat','r+');
5 d) K( V! l5 k+ ED=fscanf(fid,'%f');
  j( k' Y6 v) q- k! X) [; \( u
- \* G$ o/ }# y( J3 _3 Xfor i=1:numel(D)/2
! Y3 T4 o7 x& c- w7 g7 U    t(i)=D(2*i-1);
! y& t+ u0 S' l8 ^3 `  n    x(i)=D(2*i);
end   

" T. @. s, C/ T
0 e$ l* h% T2 `; Y4 B1 T* Hfs=500e+6;
p=log2(numel(D)/2);
1 f1 n3 q' S9 _) d- mnumber=2^p; %65536;
/ o  A& Y1 B# \$ u5 ^0 ~7 Y4 \
1 Z2 X% W! F. Hy=fft(x,number);
n=length(y);
7 u! Y' Z/ m5 ], im=abs(y)*2/n;
f=(0:n/2-1)*fs/n;
0 \7 j* A$ N' g. h* F) v: nplot(f, m(1:n/2),'k');

yanxl

hehe，大家都忙，没人答理呢
" m1 c. I! t+ Y% d. R. U7 {欢迎讨论

monkeydl

我感觉程序没有错，只是直流分量应该除以n而不是n/2