【文章內(nèi)容簡(jiǎn)介】 擊響應(yīng)長(zhǎng)度9，即濾波器階數(shù)為8b=fir1(N1,wc/pi,39。high39。,boxcar(N))。n=0:1:8。 %調(diào)用fir1函數(shù)subplot(1,2,1)。stem(n,b,39。.39。) %輸出單位脈沖沖擊響應(yīng)axis([0,8,1])。 grid。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。沖激響應(yīng)39。)。[h,w]=freqz(b,1,512) %調(diào)用函數(shù)求解頻率響應(yīng)subplot(1,2,2)。plot(w/pi,20*log10(abs(h)))。grid。axis([0,1,100,10])。xlabel(39。w/pi39。)。ylabel(39。幅度dB39。)。title(39。幅頻特性39。)。濾波器沖激響應(yīng)及幅頻特性輸出圖形如下： 圖31 基于矩形窗高通FIR濾波器特性由前面分析知基于矩形窗的FIR濾波器阻帶最小衰減是一定的，為21dB,是無(wú)法達(dá)到題設(shè)要求的100dB的。觀察幅頻特性圖形可知實(shí)際仿真結(jié)果驗(yàn)證了這一結(jié)論。通帶最大衰減1dB符合預(yù)期要求。 基于海明窗高通FIR濾波器同理利用表21求得濾波器階數(shù)為32.MATLAB源程序如下：%以下采用海明窗設(shè)計(jì)clear。close all。N=33。wc=*pi。 %沖擊響應(yīng)長(zhǎng)度33，即濾波器級(jí)數(shù)為32b=fir1(N1,wc/pi,39。high39。,hamming(N))。n=0:1:32。 %調(diào)用fir1函數(shù)subplot(1,2,1)。stem(n,b,39。.39。) %輸出單位脈沖沖擊響應(yīng)axis([0,32,1])。 grid。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。沖激響應(yīng)39。)。[h,w]=freqz(b,1,512) %調(diào)用函數(shù)求解頻率響應(yīng)subplot(1,2,2)。plot(w/pi,20*log10(abs(h)))。grid。axis([0,1,100,10])。xlabel(39。w/pi39。)。ylabel(39。幅度dB39。)。title(39。幅頻特性39。)。濾波器沖激響應(yīng)及幅頻特性圖形如下：圖32 基于海明窗高通FIR濾波器特性由前面分析知基于海明窗的FIR濾波器阻帶最小衰減是一定的，為53dB,是無(wú)法達(dá)到題設(shè)要求的100dB的。觀察幅頻特性圖形可知實(shí)際仿真結(jié)果驗(yàn)證了這一結(jié)論。通帶最大衰減1dB符合預(yù)期要求。 基于漢寧窗高通FIR濾波器同理利用圖21可知濾波器階次選為30。源程序如下：%以下采用漢寧窗設(shè)計(jì)clear。close all。N=31。wc=*pi。 %沖擊響應(yīng)長(zhǎng)度31，即濾波器級(jí)數(shù)為30b=fir1(N1,wc/pi,39。high39。,hanning(N))。n=0:1:30。 %調(diào)用fir1函數(shù)subplot(1,2,1)。stem(n,b,39。.39。) %輸出單位脈沖沖擊響應(yīng)axis([0,30,1])。 grid。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。沖激響應(yīng)39。)。[h,w]=freqz(b,1,512) %調(diào)用函數(shù)求解頻率響應(yīng)subplot(1,2,2)。plot(w/pi,20*log10(abs(h)))。grid。axis([0,1,100,10])。xlabel(39。w/pi39。)。ylabel(39。dB39。)。title(39。幅頻特性39。)。濾波器沖激響應(yīng)及幅頻特性圖形如下：圖33 基于漢寧窗高通FIR濾波器特性由前面分析知基于漢寧窗的FIR濾波器阻帶最小衰減是一定的，為44dB,是無(wú)法達(dá)到題設(shè)要求的100dB的。觀察幅頻特性圖形可知實(shí)際仿真結(jié)果驗(yàn)證了這一結(jié)論。通帶最大衰減1dB符合預(yù)期要求。 基于三角形窗高通FIR濾波器分析知濾波器階數(shù)為31，原理同上。MATLAB源程序如下：%以下采用三角形窗設(shè)計(jì)clear。close all。N=31。wc=*pi。 %沖擊響應(yīng)長(zhǎng)度31，即濾波器級(jí)數(shù)為30b=fir1(N1,wc/pi,39。high39。,triang(N))。n=0:1:30。 %調(diào)用fir1函數(shù)subplot(1,2,1)。stem(n,b,39。.39。) %輸出單位脈沖沖擊響應(yīng)axis([0,30,1])。 grid。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。沖激響應(yīng)39。)。[h,w]=freqz(b,1,512) %調(diào)用函數(shù)求解頻率響應(yīng)subplot(1,2,2)。plot(w/pi,20*log10(abs(h)))。grid。axis([0,1,100,10])。xlabel(39。w/pi39。)。ylabel(39。幅度dB39。)。title(39。幅頻特性39。)。濾波器沖激響應(yīng)及幅頻特性圖形如下： 圖34 基于三角形窗高通FIR濾波器特性由前面分析知基于三角形窗的FIR濾波器阻帶最小衰減是一定的，為25dB,是無(wú)法達(dá)到題設(shè)要求的100dB的。觀察幅頻特性圖形可知實(shí)際仿真結(jié)果驗(yàn)證了這一結(jié)論。通帶最大衰減1dB符合預(yù)期要求。 基于布拉克曼窗高通FIR濾波器分析知濾波器階數(shù)為54，原理同上。MATLAB源程序如下：%以下采用布拉克曼窗設(shè)計(jì)clear。close all。N=55。wc=*pi。 %沖擊響應(yīng)長(zhǎng)度55，即濾波器級(jí)數(shù)為54b=fir1(N1,wc/pi,39。high39。,blackman(N))。n=0:1:54。 %調(diào)用fir1函數(shù)subplot(1,2,1)。stem(n,b,39。.39。) %輸出單位脈沖沖擊響應(yīng)axis([0,54,1])。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。沖激響應(yīng)39。)。[h,w]=freqz(b,1,512) %調(diào)用函數(shù)求解頻率響應(yīng)subplot(1,2,2)。plot(w/pi,20*log10(abs(h)))。grid。axis([0,1,140,10])。xlabel(39。w/pi39。)。ylabel(39。幅度dB39。)。line([0,],[1,1])。line([,],[100,0])。line([,],[100,0])%劃線觀察通帶波紋與阻帶衰減是否符合要求title(39。幅頻特性39。)。濾波器沖激響應(yīng)及幅頻特性圖形如下： 圖35 基于布拉克曼窗高通FIR濾波器特性由前面分析知基于布拉克曼窗的FIR濾波器阻帶最小衰減是一定的，為74dB,是無(wú)法達(dá)到題設(shè)要求的100dB的。觀察幅頻特性圖形可知實(shí)際仿真結(jié)果驗(yàn)證了這一結(jié)論。通帶最大衰減1dB符合預(yù)期要求。 帶通FIR濾波器設(shè)計(jì)要求如下：數(shù)字帶通濾波器：，，通帶最大衰減為1dB，阻帶最小衰減為100dB。 基于矩形窗帶通FIR濾波器=,由圖21，窗寬度N=*=18 此處取濾波器階數(shù)為18.. %以下采用矩形窗設(shè)計(jì)clear。close all。N=19。wc=[*pi,*pi]。 %沖擊響應(yīng)長(zhǎng)度19，即濾波器級(jí)數(shù)為18b=fir1(N1,wc/pi,boxcar(N))。n=0:1:18。 %調(diào)用fir1函數(shù)subplot(1,2,1)。stem(n,b,39。.39。)。grid。 %輸出單位脈沖沖擊響應(yīng)axis([0,18,1])。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。沖激響應(yīng)39。)。[h,w]=freqz(b,1,512) %調(diào)用函數(shù)求解頻率響應(yīng)subplot(1,2,2)。plot(w/pi,20*log10(abs(h)))。grid。axis([0,1,100,10])。xlabel(39。w/pi39。)