关于MATLAB的fft变换

yanxl

ATP仿真得到的VFTO波形见附件，变换程序如下。
6 X  x. F" O/ Q  ^目的是求取VFTO波形在50MHz范围内的幅频曲线
4 M+ Z. I+ m0 ?3 {& S6 r4 y/ a怎么感觉结果怪怪的？

clear all;
  j& x- f4 D2 ?% wclc;
& n0 t; f2 R2 l- B( ]; Z; I. Y
%Reading the inout waveform
fid=fopen('wave.dat','r+');
( g& N( u/ `" W) i5 @1 `! \D=fscanf(fid,'%f');

9 Z, O1 f- `6 S) C# Mfor i=1:numel(D)/2
9 c: I' r" r5 R. ]5 B    t(i)=D(2*i-1);
    x(i)=D(2*i);
  k5 _! B; k6 d) k; y+ h) Z- [end   

1 ], o- r% X' h/ k+ [
7 _% x" i5 G, {fs=500e+6;
p=log2(numel(D)/2);
# c/ Q( }- `4 |, c( ~9 ?number=2^p; %65536;

y=fft(x,number);
n=length(y);
m=abs(y)*2/n;
f=(0:n/2-1)*fs/n;
$ H2 `$ {$ l/ f3 _0 i+ H( S* eplot(f, m(1:n/2),'k');

yanxl

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

monkeydl

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