? ? ? ?今天清明放假的第二天，早晨出去吃飯時天氣有些陰，十點多開始“清明時節(jié)雨紛紛”了。

母親遠(yuǎn)在他鄉(xiāng)看孫子，挺勞累的。父親照顧生病的爺爺……

? ? ? ?我打算今天把《DSP using MATLAB》第7章結(jié)束，剩下的幾個例子看不懂了，先跳過去。

ch7最后一個：

?

?

?

代碼：

ws = 0.6*pi; wp = 0.75*pi; Rp = 0.5; As = 50; [delta1, delta2] = db2delta(Rp, As); f = [ws, wp]/pi; m = [0, 1]; delta = [delta2, delta1]; delta[N, f, m, weights] = firpmord(f, m, delta); Nh = firpm(N, f, m, weights); [db, mag, pha, grd, w] = freqz_m(h, [1]); delta_w = 2*pi/1000; wsi = ws/delta_w; wpi = wp/delta_w;Asd = -max(db(1:1:wsi))N = N + 2 h = firpm(N, f, m, weights); [db, mag, pha, grd, w] = freqz_m(h, [1]); [Hr, ww, a, L] = Hr_Type1(h); Asd = -max(db(1:1:wsi))M = N + 1

　　運行結(jié)果：

? ? ? ??注意：代碼中讓N+2，是為了保證濾波器階數(shù)為偶數(shù)。最優(yōu)M=29。

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

總結(jié)

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

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