谁能给我解释一下这个MATLAB源程序的每一行是啥意思

Question

要求设计一个10阶的带通巴特沃斯数字滤波器，带通频率为100Hz到200Hz，采样频率为1000Hz，绘出该滤波器的幅频于相频特性
其源程序如下：clear all;
N=10;
Wn=[100 200]/500;
[b,a]=butter(N,Wn,’bandpass’);
freqz(b,a,128,1000)
figure(2)
[y,t]=impz(b,a,101);
stem(t,y
源程序错了 是这个
clear all;
N=10;
Wn=[100 200]/500;
[b,a]=butter(N,Wn,'bandpass');
freqz(b,a,128,1000)
figure(2)
[y,t]=impz(b,a,101);
stem(t,y)
figure(3)
plot(t,y);xlim([0,0.1]);title('滤波后的信号');
hold on;
plot(t,x1,'--r');xlim([0,0.1]);

Answer 1

程序很清晰啊，貌似没什么可解释的！
如果不很理解函数内容，用help自己看是最有收获的。

clear all; % 清除工作空间中的所有变量和函数
N=10;% 滤波器阶数
Wn=[100 200]/500;% 滤波器带通频率范围
[b,a]=butter(N,Wn,'bandpass');% 生成巴特沃斯滤波器的系数
freqz(b,a,128,1000)% 滤波器响应
figure(2)% 生成绘图窗口2
[y,t]=impz(b,a,101);%
stem(t,y)% 绘制离散点图
figure(3)% 生成绘图窗口3
plot(t,y);% 绘图
xlim([0,0.1]);% 给X轴限制范围
title('滤波后的信号');% 给当前窗口的图像加上标题
hold on;% 维持当前绘图窗口
plot(t,x1,'--r');% 绘图
xlim([0,0.1]);% 给X轴限制范围