关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
. i) `- }# {( f/ P- G3 j. c目的是求取VFTO波形在50MHz范围内的幅频曲线
怎么感觉结果怪怪的？
+ k! D, G% u+ e! g# ?- S9 u; [/ w# H
8 Z, z$ l1 p) K1 l- q: dclear all;
! N+ {9 u: ]4 W4 V/ yclc;
: f% Y' A. k" k& i9 L* O
! Q9 S! r4 ^3 l! j5 W. r2 N%Reading the inout waveform
fid=fopen('wave.dat','r+');
) a0 S: _2 J6 E! k* ]D=fscanf(fid,'%f');
4 h4 ^, g# `. a9 T, z# k
% W4 m/ T4 T( p" y, H! j8 efor i=1:numel(D)/2
! Q1 o, O+ r0 v1 K; r    t(i)=D(2*i-1);
- G' c. A* T5 z    x(i)=D(2*i);
end   

" F; y$ }7 J1 z+ ^& X
fs=500e+6;
2 T; e6 |5 @% ^* h+ D" up=log2(numel(D)/2);
0 }0 z! i' ]3 ?' b7 Z+ r' ^number=2^p; %65536;
3 y2 Q2 w6 V+ s- Z
y=fft(x,number);
n=length(y);
m=abs(y)*2/n;
7 ?/ G, ]( X; I4 q2 of=(0:n/2-1)*fs/n;
. G1 j5 T7 d: m4 Q; Rplot(f, m(1:n/2),'k');

yanxl

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

monkeydl

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