M=63 J=(M-1)/2 Wp=0.5*pi Ws=0.6*pi

Np=100 Ns=100

weightp=1 weights=100 通带和阻带的加权系数

WpN=linspace(0,Wp,Np)

WsN=linspace(Ws,pi-0.05,Ns)

Wm=[WpN WsN]

p=[ones(1,50),ones(1,50)/2]

d=[p,zeros(1,Ns)]

W=diag([weightp*ones(1,Np),weights*ones(1,Ns)])

C=cos(Wm*[0:J])

Q=diag(cos(Wm*0.5))

g=(W*Q*C)\(W*d) g=g

h=(g(J:-1:2)+g(J+1:-1:3))/4

h=[g(J+1)/4 h g(1)*0.5+g(2)/4]

h=[h fliplr(h)] W=linspace(0,pi-0.1,512)

H=freqz(h,1,W)

plot(W/pi,20*log10(abs(H)))

-filter design