关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
目的是求取VFTO波形在50MHz范围内的幅频曲线
怎么感觉结果怪怪的？
; G. d! Q  R; O! w1 ~% J
( J9 X3 L/ X9 l  R, i# o- cclear all;
1 \6 C+ C2 L' d1 n# F1 Q# z) _* Rclc;

5 Z' U9 p4 u2 y0 x- l4 S%Reading the inout waveform
fid=fopen('wave.dat','r+');
9 Q; {& a+ M* ZD=fscanf(fid,'%f');
! v, v2 l! S+ y0 f
for i=1:numel(D)/2
: j5 k3 p6 s$ _& ~" A3 U    t(i)=D(2*i-1);
$ K; y8 F) V. I+ ]1 P9 ?    x(i)=D(2*i);
7 p+ }5 @2 P# g" ]: ]# _1 Oend   
4 ^: F+ d6 D$ `0 [

fs=500e+6;
2 {" p( B) z( c+ F, @% C5 I" F8 hp=log2(numel(D)/2);
' q3 n  Q# S0 Ynumber=2^p; %65536;

8 E. S8 U2 B# p# G& yy=fft(x,number);
n=length(y);
0 o; Y: y# ?  H/ o$ k* B5 u* Y1 {m=abs(y)*2/n;
f=(0:n/2-1)*fs/n;
6 K" Y2 Z/ P' w% Gplot(f, m(1:n/2),'k');

yanxl

hehe，大家都忙，没人答理呢
欢迎讨论

monkeydl

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