关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
目的是求取VFTO波形在50MHz范围内的幅频曲线
& V  l  F# h/ _7 y7 J+ B0 o( P怎么感觉结果怪怪的？
) V: T, |. Q, k( s7 ~7 \+ J# h
' y- J4 S$ v7 k% \2 @) l+ Hclear all;
; E  y& y9 k" w' \clc;

%Reading the inout waveform
fid=fopen('wave.dat','r+');
D=fscanf(fid,'%f');

8 T+ k8 H: S) b3 z( @; U! u. P" t# Qfor i=1:numel(D)/2
9 V- o3 |6 R6 J    t(i)=D(2*i-1);
4 ?$ g8 j. E5 S! D) `- H    x(i)=D(2*i);
; O: e$ n+ g' a' Y6 J0 l3 W' xend   

4 }3 O0 s$ e3 I* ~8 z
; W- `: @' h! A2 `$ Z& f+ ]5 {+ wfs=500e+6;
) P( Q0 d0 S: g0 L& L( g$ @p=log2(numel(D)/2);
  Q6 q% N2 h+ b9 P6 |number=2^p; %65536;

y=fft(x,number);
3 i, a+ z1 A6 c  nn=length(y);
2 }2 V( l: k4 H1 V  i/ T" @m=abs(y)*2/n;
3 P* e$ z6 M" @; N# X3 ~( R( Cf=(0:n/2-1)*fs/n;
plot(f, m(1:n/2),'k');

yanxl

hehe，大家都忙，没人答理呢
: v7 j2 `8 {- c# i. ]# J+ Y6 G欢迎讨论

monkeydl

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