关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
目的是求取VFTO波形在50MHz范围内的幅频曲线
: ?7 a7 X' e) v怎么感觉结果怪怪的？
! c, G) A) b: ]# q4 _
clear all;
; u6 L# n9 n" M* R" j: l: fclc;
4 {, e' e0 T  V& N/ ~3 x8 G
  h" @1 A0 {6 K( t7 c, a. j. D2 j%Reading the inout waveform
' K; W! o5 t( V* ~fid=fopen('wave.dat','r+');
1 ]- W. U7 ]* r' _% n9 RD=fscanf(fid,'%f');

4 P. s  i, O3 O, [+ m: |for i=1:numel(D)/2
; `8 }  W  s4 A' @1 l    t(i)=D(2*i-1);
    x(i)=D(2*i);
end   

4 N+ e9 q% f" `! {* |  a
fs=500e+6;
3 _8 I; }2 ?! F% W: `# Mp=log2(numel(D)/2);
number=2^p; %65536;

: {% Q+ x+ m4 G$ Uy=fft(x,number);
n=length(y);
m=abs(y)*2/n;
$ ~/ X' n* F$ T/ v5 @  h  @. \  jf=(0:n/2-1)*fs/n;
9 o1 [" W5 T# D: N% V' g* }& yplot(f, m(1:n/2),'k');

yanxl

hehe，大家都忙，没人答理呢
4 |; S1 H! h; ~  l$ p& b3 s2 @欢迎讨论

monkeydl

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