求MATLAB学的好的注释一下，注释有%号的地方就行，越详细越好！

Question

clc;close all;clear;t=(1:100)/100;w0=pi;x=sin(w0*t)+sin(3*w0*t)/3;figure(1);subplot(211)plot(x);xlabel('t');ylabel('x(t)');title('原时域连续信号f(t)');grid;N=512;n=0:511;Xk=abs(fft(x,N)); %subplot(212)figure(1);plot(n,Xk);axis([0 N min(Xk) max(Xk)]); %xlabel('频谱图');ylabel('|Xk|');%下面开始采样X=fft(x,512);w=(0:255)/256*500 %T=4*tx=sin(w0*t)+sin(3*w0*t)/3;figure(2);stem(x) %xlabel('n');ylabel('x(n)');title('采样信号波形图');grid;%设计采样信号频谱图figure(3);w==(0:255)/256*500 %plot(w,abs([X(1:256)]));xlabel('Hz');ylabel('频率响应幅度');grid;%设计低通滤波器[B,A]=butter(8,350/500); [H,w]=freqz(B,A,512,2000);figure(4);subplot(211);plot(w*2000/(2*pi),abs(H));xlabel('Hz');ylabel('频率响应幅度');title('低通滤波器波形图');grid;%恢复原信号y=filter(B,A,x); %figure(5);subplot(2,1,1);plot(y);xlabel('t');ylabel('x(t)');title('还原信号');grid;Y=fft(y,512);w=(0:255)/256*500; %subplot(2,1,2);plot(w,abs([Y(1:256)]));xlabel('Hz');ylabel('频率响应幅度');title('恢复后的频谱波形图');grid;

Answer 1

clc;
close all;
clear;
t=(1:100)/100;
w0=pi;
x=sin(w0*t)+sin(3*w0*t)/3;%x和t的关系式
figure(1);%图1
subplot(211)%分割图1
plot(x);%画图
xlabel('t');%x轴坐标名
ylabel('x(t)');%y轴坐标名
title('原时域连续信号f(t)');%图标题
grid;%加网格
N=512;
n=0:511;%产生0到511整数序列
Xk=abs(fft(x,N)); %傅里叶变换，求绝对值
subplot(212)%分割图2
figure(1);
plot(n,Xk);
axis([0 N min(Xk) max(Xk)]); %坐标轴范围为x轴0:N，y轴min(Xk):max(Xk)
xlabel('频谱图');
ylabel('|Xk|');%x轴、y轴的坐标名称
%下面开始采样
X=fft(x,512);%傅里叶变换
w=(0:255)/256*500；%产生一组数
T=4*t;
x=sin(w0*t)+sin(3*w0*t)/3;
figure(2);
stem(x);
xlabel('n');
ylabel('x(n)');
title('采样信号波形图');%图形标题
grid;
%设计采样信号频谱图
figure(3);
w==(0:255)/256*500
plot(w,abs([X(1:256)]));
xlabel('Hz');
ylabel('频率响应幅度');
grid;
%设计低通滤波器
[B,A]=butter(8,350/500);
[H,w]=freqz(B,A,512,2000);figure(4);
subplot(211);
plot(w*2000/(2*pi),abs(H));
xlabel('Hz');
ylabel('频率响应幅度');
title('低通滤波器波形图');
grid;
%恢复原信号
y=filter(B,A,x);
figure(5);
subplot(2,1,1);
plot(y);
xlabel('t');
ylabel('x(t)');
title('还原信号');
grid;
Y=fft(y,512);
w=(0:255)/256*500;
subplot(2,1,2);
plot(w,abs([Y(1:256)]));
xlabel('Hz');
ylabel('频率响应幅度');
title('恢复后的频谱波形图');
grid;