怎样用MATLAB分析 Sa（50*pi*t）的频谱（麻烦写下程序，急求。。。谢谢啦）

fs=1;%设定采样频率
N=100;%采样点数
n=0:N-1
t=n/fs;
f0=10;%设定正弦信号频率
x= Sa（50*pi*t）%生成正弦信号
figure(1);
subplot(2，3，1);
plot(t，x);%作正弦信号的时域波形
xlabel('t');
ylabel('y');
title('正弦信号y=0.5*pi*t时域波形');
grid;
%进行FFT变换并做频谱图
y=fft(x，N);%进行fft变换
mag=abs(y);%求幅值
f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换
figure(1);
subplot(2，2，2);
plot(f，mag);%做频谱图
xlabel('频率(Hz)');
ylabel('幅值');
title('正弦信号y=0.5*pi*t幅频谱图N=100');
grid;