关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
目的是求取VFTO波形在50MHz范围内的幅频曲线
怎么感觉结果怪怪的？

clear all;
clc;
% ^6 v' s: W1 K# t6 a
/ z! _/ i7 d1 R+ F7 T& q3 E4 \%Reading the inout waveform
fid=fopen('wave.dat','r+');
( U  ?; A2 V  B, \5 Y1 ID=fscanf(fid,'%f');
; w; r0 z! }4 \& T: M
4 y3 I! l" r9 D' l0 Ffor i=1:numel(D)/2
    t(i)=D(2*i-1);
    x(i)=D(2*i);
end   

* a2 u* f  c2 u5 c. O9 r8 ?3 d8 O
1 o. ^+ q1 }' `/ v+ Zfs=500e+6;
p=log2(numel(D)/2);
number=2^p; %65536;

y=fft(x,number);
n=length(y);
$ q) N) c1 ^# b9 ?4 k; e* X7 d: E( p% Lm=abs(y)*2/n;
6 b8 A5 Y# N' C9 F! [1 ?f=(0:n/2-1)*fs/n;
6 H# s( x) ?% I. p7 P; C8 G  b6 a& \. Pplot(f, m(1:n/2),'k');

yanxl

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

monkeydl

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