对数频谱图的求解

sy0124

支路电流的大小如图所示： ，我想达到以下图的目的 ，我该怎么做，首先在网上找了个对数频谱图的程序
fs=100;%设定采样频率
N=128;
n=0:N-1;
t=n/fs;
f0=10;
%进行FFT变换并做频谱图
y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
figure(1);
subplot(232);
plot(f,mag);%做频谱图
axis([0,100,0,80]);
xlabel('频率(Hz)');
ylabel('幅值');
title('正弦信号y=2*pi*10t幅频谱图N=128');
grid;

求出的结果为 ，我修改了下程序，把电流的数值导出来之后a=[1 2 ];x=a(:,1);这样求不出结果，按上述程序求出的结果也得不到负值，求大神指点