育才学习网1. 怎样用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;
2. 怎样用MATLAB画正弦函数以及怎样看频谱图
sinx中角频率omega为1，周期为2*pi/omega=2*pi,2*pi/步长=628，你花的就是周期2*pi的函数
fft就是把时域上的曲线变成频域上的曲线，所以fft变换后横坐标是频率，比如你的这个函数，注意到频域图中第三个点出出现最大值，计算2/1257*100=0.159，其实它就是sinx的频率，即1/(2*pi)=0.159，不过你只用stem(abs(fft(y1)））；这条语句显然不能表示横坐标的情况，一般横坐标是0~采样频率，并且由于采样定理，进一步缩小为0~采样频率/2.，也就是说，横轴坐标应该是
f=(0:1256)/1257*100;
纵坐标上，你用了abs，所以应该是幅值。
另外，fft指快速傅立叶变换，其变换点是2^n个，你虽然没有指定，不过程序仍然是这么做的，所以你下次根据数据多少制定变换点数就行了，比如1024
最后，可以参见matlab中关于fft的帮助