clear; close all; global delta_p delta_s; delta_p = 1; delta_s = 1; global wc; wc = 0.5*pi; wd = 0.6*pi; N = 31; M = (N-1)/2; %STEP1: ωiを初期化 %ωiを初期化: 当分割 epsilon = 0.001; if mod(M,2)==1 %M=奇数 wi_left = linspace(0+epsilon,wc,(M…