关于MATLAB的fft变换

yanxl · 发表于 2009-6-24 17:52:25

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

您需要登录才可以下载或查看，没有账号？立即加入

×

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

clear all;
clc;

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

for i=1:numel(D)/2
t(i)=D(2*i-1);
x(i)=D(2*i);
end

fs=500e+6;
p=log2(numel(D)/2);
number=2^p; %65536;

y=fft(x,number);
n=length(y);
m=abs(y)*2/n;
f=(0:n/2-1)*fs/n;
plot(f, m(1:n/2),'k');

yanxl · 发表于 2009-6-25 08:57:17

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

monkeydl · 发表于 2009-6-28 19:53:13

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

账号		自动登录	找回密码
密码			立即加入

[讨论] 关于MATLAB的fft变换

马上加入，结交更多好友，共享更多资料，让你轻松玩转电力研学社区！

相关帖子