わかりやすい https://tech-blog.cerevo.com/archives/10159/

clear; close all; N = 61; A = 0; b = 0; syms w; M = (N-1)/2; c = [w];%dummy for n=1:M c(end+1) = 2*cos(w*n); end c(1) = 1; c = fliplr(c); c = transpose(c); omega = linspace(0, pi, N+1); %遮断周波数の定義 delta_w = 0.1*pi;%転移域の幅 wc = 0…