[求助]频谱分析

[ct8613]

本人菜鸟一个，刚开始学习MATLAB。编了一个有关频谱分析的程序，但是实现不了，还望大家帮忙。小弟在此谢过了。在线等。
Ts=0.0005;
t=0:Ts:0.1;
%采样信号频率是fs=2000；
fs=1/Ts;
%原信号频率是50；
x=60*sin(2*pi*50*t);
%N是总共的采样点数；
N=length(x);
Ix=fft(x);
w=2*pi*50,c=5.72e-5;
Ux=Ix/(i*w*c);
u=ifft(Ux);
Fx=fft(u^2);
%N个点有N-1个间隔；
df=fs/(N);
n=0:N-1;
f=n*df;
subplot(211);
plot(t,x);grid;
subplot(212);
plot(f,abs(Fx));grid;