代碼：

wp = 0.2*pi; ws = 0.3*pi; Rp = 0.25; As = 50;
[delta1, delta2] = db2delta(Rp, As);[N, f, m, weights] = firpmord([wp, ws]/pi, [1, 0], [delta1, delta2]);N
f
m
weightsh = firpm(N, f, m, weights);
[db, mag, pha, grd, w] = freqz_m(h, [1]);delta_w = 2*pi/1000; wsi = ws/delta_w+1; wpi = wp/delta_w;
Asd = -max(db(wsi:1:501))N = 46;
h = firpm(N, f, m, weights);
[db, mag, pha, grd, w] = freqz_m(h, [1]);
[Hr, ww, b, L] = Hr_Type2(h);delta_w = 2*pi/1000; wsi = ws/delta_w+1; wpi = wp/delta_w;
Asd = -max(db(wsi:1:501))M = N + 1

　　運行結(jié)果：

? ? ? ??注意，我們是在最大阻帶衰減As的值超過給定的50dB后停止迭代的，最優(yōu)的M=47。這個值要比采用窗函數(shù)（Kaiser窗時M=61）和頻率采樣方法（M=60）

要小得多。

轉(zhuǎn)載于:https://www.cnblogs.com/ky027wh-sx/p/6661508.html

總結(jié)

以上是生活随笔為你收集整理的《DSP using MATLAB》示例Example7.23的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。