Fs=8000; %采样频率 f=[1200 1600]; a=[1 0]; dev=[0.01 0.1]; [n,fo,ao,w]=remezord(f,a,dev,Fs); %估计滤波器的参数 b=remez(n,fo,ao,w); %用 remez 函数设计滤波器 freqz(b,1,1024,Fs); %滤波器的频率响应, 如图 1 t=(0:100)/Fs; x=sin(2pit1000)+sin(2pit2000); %设置输入信号 x x1=filter(b,1,x); %对信号 x 进行滤波,得输出信号 x1 figure subplot(211);plot(t,x);grid; subplot(212);plot(t,x1);grid; 依次解释每句代码的含义

好的，下面是每句代码的含义：

1. Fs=8000; %采样频率
定义采样频率为 8000 Hz。

2. f=[1200 1600]; a=[1 0]; dev=[0.01 0.1]; [n,fo,ao,w]=remezord(f,a,dev,Fs);
定义截止频率为 1200 Hz 和 1600 Hz，通带增益为 1，阻带增益为 0，通带和阻带的最大衰减量分别为 0.01 和 0.1。使用 remezord 函数估计出滤波器的阶数、截止频率、增益和权重。

3. b=remez(n,fo,ao,w);
使用 remez 函数根据上一步得到的参数设计滤波器，并将其系数保存在 b 中。

4. freqz(b,1,1024,Fs);
使用 freqz 函数计算滤波器的频率响应，并绘制其幅度和相位响应。其中，1 表示滤波器的单位增益响应为 1，1024 表示计算频率响应的点数，Fs 表示采样频率。

5. t=(0:100)/Fs;
生成一个 101 个元素的时间向量 t，其取值范围为 0 到 100/Fs 秒。

6. x=sin(2*pit1000)+sin(2*pit2000); %设置输入信号 x
生成一个输入信号 x，其为 1000 Hz 和 2000 Hz 两个正弦波的叠加。

7. x1=filter(b,1,x);
使用 filter 函数对输入信号 x 进行滤波，得到输出信号 x1。

8. figure
创建一个新的图形窗口。

9. subplot(211);plot(t,x);grid;
在当前图形窗口中创建一个 2 行 1 列的子图，并将第一个子图设为当前子图。在该子图中绘制输入信号 x，并打开网格线。

10. subplot(212);plot(t,x1);grid;
在当前图形窗口中将第二个子图设为当前子图。在该子图中绘制输出信号 x1，并打开网格线。