clear; close all; syms w; chebyshev = @(c1,c0) 2*w*c1 - c0; c0 = 1; c1 = w; cs = [c0 c1]; for n=2:10 cn_prev = cs(n-1); cn = cs(n); cn_next = chebyshev(cn, cn_prev); cs(end+1) = cn_next; end n = 5; xis = [0.1 0.5 1.0]; for idx=1:length(xis…

clear; close all; s = tf('s'); syms w; xis = [0.1 0.5 1.0]; n = 5; for idx=1:length(xis) xi = xis(idx); e = sqrt(10.^(xi/10)-1); tmp = sqrt(1+e^(-2))+e^(-1); sinhv = 1/2*(tmp^(1/n)-tmp^(-1/n)); coshv = 1/2*(tmp^(1/n)+tmp^(-1/n)); H_gain = …

clear; close all; s = tf('s'); Ns = [5,8,10]; syms w; for idx = 1:length(Ns) n = Ns(idx); D = 1; H_gain = 1; H_phi = 0; for jdx=1:n j = jdx; sigma_j = -sin((2*j-1)/(2*n)*pi); omega_j = cos((2*j-1)/(2*n)*pi); sj = sigma_j + i*omega_j; D = D…

clear; close all; s = tf('s'); Ns = 1:8; for idx = 1:length(Ns) n = Ns(idx); D = 1; for jdx=1:n j = jdx; sigma_j = -sin((2*j-1)/(2*n)*pi); omega_j = cos((2*j-1)/(2*n)*pi); sj = sigma_j + i*omega_j; D = D*(s-sj); end H = 1/D; subplot(4,2,n)…