关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
  w: X+ n. o! d" o; N- t$ A目的是求取VFTO波形在50MHz范围内的幅频曲线
怎么感觉结果怪怪的？
  ^. P9 Y8 r: z& J7 B) Z& t! ^
clear all;
, b' E, j. U1 C9 z# L- L1 Xclc;

7 H% t  h* P" B& t%Reading the inout waveform
6 G3 Y- c% `. V7 a  q+ e  Q7 Q8 F7 Jfid=fopen('wave.dat','r+');
8 b4 `+ S* U2 r. @# V3 i" y0 Y$ TD=fscanf(fid,'%f');

for i=1:numel(D)/2
( e2 V9 c- w, ?# z3 Y( s    t(i)=D(2*i-1);
. w) k* |1 [; A" B5 o7 {3 q    x(i)=D(2*i);
+ c- \. z/ n, H+ d' g$ w% Lend   


0 u& c* T9 D; A$ _9 c6 t2 b& nfs=500e+6;
  `4 r3 T: Q( `3 [1 n$ t$ T3 }! kp=log2(numel(D)/2);
number=2^p; %65536;
5 |, L" ^& h' O/ Z9 ~
! j; g- ^. v$ f& oy=fft(x,number);
n=length(y);
% [. f$ `0 i# e# y+ Qm=abs(y)*2/n;
' U1 j! n" A/ Z+ y# `6 _+ p+ Vf=(0:n/2-1)*fs/n;
plot(f, m(1:n/2),'k');

yanxl

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

monkeydl

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