关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
9 C; ]1 X! g/ r1 |& v目的是求取VFTO波形在50MHz范围内的幅频曲线
3 e. P: t( `4 g7 n4 _6 W5 X* j# u9 x$ t怎么感觉结果怪怪的？
; f: T! F1 V& I$ P' [$ b& C
  f6 F# D8 O3 R' @8 {+ Tclear all;
: c$ l. `  I9 Q+ X: J( v0 vclc;
: A' G/ J4 b, i6 {9 h! J6 g6 J
& [/ x& {1 g# }. |: o+ B%Reading the inout waveform
fid=fopen('wave.dat','r+');
D=fscanf(fid,'%f');
( e5 z2 T, O( Z; |: A5 M2 i
8 e* f- a6 X  W% f- kfor i=1:numel(D)/2
9 t2 L  X& E4 I9 T7 R8 [    t(i)=D(2*i-1);
    x(i)=D(2*i);
: m3 V1 b4 I' C9 O, s$ w: v, i# eend   
# l6 |: ^- f; Y" D: U

fs=500e+6;
4 r( b# Z$ G& X# Lp=log2(numel(D)/2);
/ z3 S) `% y: z& D: i' cnumber=2^p; %65536;

y=fft(x,number);
) _% C" u* ~2 H( `n=length(y);
m=abs(y)*2/n;
  W6 Y5 A; b, ?" W& P, Of=(0:n/2-1)*fs/n;
plot(f, m(1:n/2),'k');

yanxl

hehe，大家都忙，没人答理呢
$ d4 E# y6 O9 [: F+ }9 Y- O2 I% p$ j- Z欢迎讨论

monkeydl

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