对大量的离散数据进行fft变换

clear
clc
fs=500000;%采样频率
n=100000;%采样的数据
t=0:1/fs:0.2;
base=importdata('C:\Users\72910\Desktop\f100000.TXT');
y1=fft(base(:,2));
y2=fftshift(y1);
f=(0:n-1)*fs/n-fs/2;
figure(1)
plot(t,base(:,2),'r');%原始采样图
figure(2)
plot(f,abs(y2),'b');%频谱图
axis([0,10000,0,500000]);






错误是：
错误使用 plot
向量长度必须相同。
出错 shiyan2 (line 13)
plot(f,abs(y2),'b');%频谱图

长度问题已经解决了，但是做出来的第二个图是空的

为什么y1的计算结果是nan

你的数据里面出现了inf~

>>find(isinf(base(:,2)))

ans =

       45308
       45309
       55364
       55365
       60306
       65195
       75382
       75383
       75384
       80338
       80340
       85715
       85716
       85717
       85718
       85719
       90346
       90347
       95848
       95849
       95850
       95851
       95852

修改后频谱图如图所示：

是不是频率成分太多还是程序有问题。

godblessmee 发表于 2018-7-23 18:21
你的数据里面出现了inf~

>>find(isinf(base(:,2)))

谢谢。修改后显示的频谱是不是频率成分太多

jiangkai 发表于 2018-7-24 18:26
谢谢。修改后显示的频谱是不是频率成分太多

实际信号里面频率成分多是正常的，你的情况很正常

godblessmee 发表于 2018-7-24 18:31
实际信号里面频率成分多是正常的，你的情况很正常

好的，谢谢

能问一下修改后的程序嘛