matlab对非正弦信号频谱分析,急,

Fs=2048;%采样频率
dx=1/Fs;
N=2048;
df=Fs/N;
x=(0:N-1)*dx;
%从y的表达式可看出信号最高频率为500/2=250Hz,选采样频率至少为2.56倍
y=(40/pi)*sin(100*pi*x)+(40/(3*pi))*sin(300*pi*x)+(40/(5*pi))*sin(500*pi*x);
figure(1)
plot(x(1:200),y(1:200));
figure(2)
Y=abs(fft(y));
f=(0:N/2)*df;
Y(1)=Y(1)/N;
Y(2:N/2)=Y(2:N/2)/(N/2);
Y(N/2+1)=Y(N/2+1)/N;
plot(f,Y(1:N/2+1))