如何采用matlab对单片机采集的数据进行频谱分析？

N=2500; %2500个采样点
Fs=1000; %采样频率
t=(0:N-1)/Fs; %时域
X=[]; %要快速傅里叶变换的原始信号
Y=fft(X，N); %快速傅里叶变换
f=(0:N-1)*Fs/N; %频域的点数
subplot(2，1，1);
plot(f，abs(Y)) %频谱图
subplot(2，1，2);
plot(f，angle(Y)) %相位谱

将这些数据通过通信协议和电脑通信，然后电脑就可以处理这些数据了

采用函数subplot(N，M，N*M)
N为窗口水平分割的图像数目，M为窗口竖直分割的图像数目，后面的N*M为第几个图像
例如：
x=[1:10];
subplot(211)
plot(x，y，'r')
subplot(212)
plot(x，2*y，'m')