clear all
Fp=1000; dp=0.05; e2=((1-dp)^-2 -1); op=j*2*pi*Fp; n=3;
s=e2*[16 0 -24 0 9 0 0].*op.^[0:2*n]; s=s+[zeros(1,2*n) op^(2*n)];
kk=1; N=1000; r=roots(s)
for k=1:2*n
if real(r(k)) <=0, rr(kk) = r(k); kk=kk+1; end
end
%mo=poly(rr); ja=(-1)^n* sqrt(1/(1+e2)) *prod(rr); tf(ja, mo)
mo=poly(rr); ja=(-1)^n* 1 *prod(rr); tf(ja, mo)
figure(2); subplot(121); zplane(ja,mo)
subplot(122); [H,w]=freqs(ja,mo,N); plot(w/(2*pi), abs(H)); xlim([0 2*Fp])
hold on; plot(Fp, 1-dp, 'ro'); hold off; grid on

이게 체비세프 3차 코드에요.
이걸 조금 바꿔서 4차, 5차로 만들고싶은데 좀 도와주세요~
멜주소 적어주시면 사례는 제대로 하겠습니다!