错误使用 buttord 截止频率必须在 (0,1) 区间内。 出错 code (第 24 行) [n, Wn] = buttord(Wp, Ws, 3, 60); % 计算滤波器的阶数和归一化截止频率

非常抱歉，我的回答还是有误。在使用buttord函数计算滤波器阶数和归一化截止频率时，确实应该将通带边界频率和阻带边界频率都除以采样频率的一半，并且将第四个参数40改为60。但是，计算归一化截止频率时应该使用通带边界频率和阻带边界频率的平均值，而不是只使用通带边界频率。因此，下面这段代码应该是正确的：

Wp = 100; % 通带边界频率
Ws = 120; % 阻带边界频率
Rp = 3; % 通带最大衰减量
Rs = 60; % 阻带最小衰减量
Fs = 1000; % 采样频率
Wn = (Wp/(Fs/2) + Ws/(Fs/2))/2; % 归一化截止频率
[n, ~] = buttord(Wp/(Fs/2), Ws/(Fs/2), Rp, Rs); % 计算滤波器的阶数
[b, a] = butter(n, Wn, 'stop'); % 根据阶数和归一化截止频率设计滤波器

这样就可以得到一个截止频率为100Hz的巴特沃斯带阻滤波器了。再次对之前的回答表示深深地歉意。