好心人帮注释一下这个MATLAB程序吧，最好每句都注释，急急急！万分感谢！

Question

Fs = 90;
n = 0:1/Fs:199*1/Fs;
x=sin(2*pi*20*n)+sin(2*pi*30*n);
N=length(x);
Ft=(1:N)*Fs/N;
X=abs(fft(x)); %傅里叶变换
figure(1);
subplot(211)
plot(Ft,X);
xlabel('频率');
ylabel('幅度');
title('合适的采样频谱');
Fs = 40;
n = 0:1/Fs:199*1/Fs;
x = sin(2*pi*20*n)+sin(2*pi*30*n);
N=length(x);
Ft=(1:N)*Fs/N;
X = abs(fft(x)); %傅里叶变换
subplot(212)
plot(Ft,X);
xlabel('频率');
ylabel('幅度');
title('欠采样频谱');

Answer 1

Fs = 90; %给Fs附初值90
n = 0:1/Fs:199*1/Fs; %给n附不同的值，范围是0--199*1/90,步长为1/90
x=sin(2*pi*20*n)+sin(2*pi*30*n); %对于给定的n求正弦函数，赋值给x存储这个矩阵中
N=length(x); %求x矩阵的长度
Ft=(1:N)*Fs/N; %求1--N之间的值，步长为1，根据公式求的一组数给Ft
X=abs(fft(x)); %傅里叶变换
figure(1); %作第一个图
subplot(2，1，1) %画第二行，第一列，画到第一幅图上
plot(Ft,X); %做曲线，Ft的值作为横坐标，X作为纵坐标
xlabel('频率'); %x轴表示‘频率’
ylabel('幅度'); %y轴表示‘幅度’
title('合适的采样频谱'); %这曲线标题为‘合适的采样频谱’
后半部分我就没写，是一样的，其中subplot(212) %%是表示在第二幅图上作图

Answer 2

居然有USTC的孩儿在在求...发给你了已经:
Digital Image Processing Using MATLAB.