给出我的频谱图程序和结果，出现问题。

胡晓宇

结果出现在0点和最大点处，振幅无穷大，其他地方都为零。我的频谱图程序如下：
function myfrequencyplot
% all the files are save as t and s
clc;
clear;

load result_6000_30.mat;
%读取文件，里面保存有时间序列t和积分结果s，s有18维
fs = 20000;%采样频率
N = 1024*8;%采样点数
n = 0:N-1;
vt = n/fs;
vvt = t(1)+vt;
x = interp1(t,s(:,1),vvt,'linear');%插值获的等间距的数据点
y = fft(x,N);
mag = abs(y);
f = (0:length(y)-1)'*fs/length(y);
plot(f,mag);
xlabel('Frequency(Hz)');
ylabel('Magnitude');
title('N = 128');
grid;

ChaChing

直觉DC没去除! 同待高人路过!