写一段MATLAB用频率采样法设计线性相位低通滤波器的代码，注释越多越好

好的，下面是MATLAB用频率采样法设计线性相位低通滤波器的代码，注释已经尽可能详细。

%% 频率采样法设计线性相位低通滤波器
clear; clc;

%% 参数设置
fs = 8000;              % 采样频率，单位Hz
fc = 1000;              % 截止频率，单位Hz
N = 51;                 % 滤波器阶数（必须为奇数）
M = (N-1)/2;            % 中心点
w_c = 2*pi*fc/fs;       % 数字滤波器的截止角频率

%% 频率采样
n = 0:M;
h_d = (w_c/pi)*sinc(w_c*(n-M)/pi);      % 理想的数字低通滤波器的冲激响应
w = linspace(0,pi,M+1);                 % 频率采样点
H_d = h_d.*exp(-1j.*M.*w);               % 频率采样后的理想数字低通滤波器频率响应

%% IFFT变换
h = ifft(H_d);

%% 线性相位处理
h = real(h);
h = h.*((-1).^n);                       % 线性相位处理

%% 绘制图形
freqz(h,1,1024,fs);                     % 绘制滤波器频率响应
title('Frequency Response')

代码中，首先设置了采样频率、截止频率和滤波器阶数等参数。接着，通过频率采样的方法，求出了理想的数字低通滤波器的冲激响应，并进行了频率采样，得到了理想数字低通滤波器的频率响应。然后，通过IFFT变换将频率响应转换为时域冲激响应，再进行线性相位处理，最后绘制出了滤波器的频率响应。