【正文】 offtitle(39。High Pass Filter39。)。subplot(2,2,3)plot(w,abs(h),39。r39。)。xlabel(39。f39。)。ylabel(39。H39。)。box offsubplot(224)plot(w,angle(h),39。m39。)。gridxlabel(39。f39。)。ylabel(39。phase39。)。grid on。box offclear all。Fs=2000。Ts=1/Fs。fs1=200。f1=300。f2=400。fsh=500。wp=2*pi*fs1*Ts。ws=2*pi*fsh*Ts。w1=2*pi*f1*Ts。w2=2*pi*f2*Ts。rp=3。rs=40。[n wn]=buttord(wp/pi,ws/pi,rp,rs)。[numz denz]=butter(n,[w1/pi w2/pi])。Gz_BP=tf(numz,denz,Ts)[h w]=freqz(numz,denz,128,Fs)。figure(2)subplot(2,2,[1 2])。plot(w,20*log10(abs(h)),39。r39。)。xlabel(39。f39。)。ylabel(39。20log(H)39。)。grid on。box offtitle(39。Band Pass Filter39。)。subplot(2,2,3)。plot(w,abs(h),39。r39。)。gridxlabel(39。f39。)。ylabel(39。H39。)。grid on。box offsubplot(2,2,4)。plot(w,angle(h),39。m39。)。xlabel(39。f39。)。ylabel(39。phase39。)。grid on。box offclear all實(shí)驗(yàn)四 FIR數(shù)字濾波器的設(shè)計(jì)實(shí)驗(yàn)?zāi)康?． 掌握FIR數(shù)字濾波器的設(shè)計(jì)方法與步驟；2． 理解系統(tǒng)頻率響應(yīng)的概念，學(xué)習(xí)編寫(xiě)計(jì)算系統(tǒng)頻率響應(yīng)的方法。實(shí)驗(yàn)原理 1. FIR數(shù)字濾波器設(shè)計(jì)的詳細(xì)內(nèi)容見(jiàn)教材第4章。(1) 窗口法：窗口法設(shè)計(jì)FIR數(shù)字濾波器的步驟： A 給出希望的濾波器頻率響應(yīng)函數(shù)； B 根據(jù)允許的過(guò)渡帶寬度及阻帶衰減確定所采用的窗函數(shù)和N值； C 做的逆傅里葉變換得hd(n)； D對(duì)hd(n)加窗處理得到有限長(zhǎng)序列 E對(duì)h(n)做傅里葉變換得到頻率響應(yīng)，用作為的逼近，并用給定的技術(shù)指標(biāo)來(lái)檢驗(yàn)。 (2) 切比雪夫一致逼近計(jì)算機(jī)輔助設(shè)計(jì)方法：用Remez算法實(shí)現(xiàn)FIR濾波器的等波紋逼近。 2. （與FIR數(shù)字濾波器設(shè)計(jì)相關(guān)MATLAB函數(shù)）。實(shí)驗(yàn)內(nèi)容1. 請(qǐng)選擇合適的窗函數(shù)及N來(lái)設(shè)計(jì)一個(gè)線性相位低通濾波器 要求其最小阻帶衰減為－45dB，過(guò)渡帶寬為，(1) 已知，求出并畫(huà)出曲線。(2) 保留原有軌跡，畫(huà)出用滿(mǎn)足所給條件的其他幾種窗函數(shù)設(shè)計(jì)出的曲線?！境绦虼a】clear all。close all。clcrs=45。w0=8*pi/51。w0_ham=。 N=2*pi*w0_ham/w0。M=ceil(N1)1。wc=*pi。hn=fir1(M,wc/pi,hamming(M+1))subplot(211)stem(0:M,hn,39。filled39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。Hamming Window39。)。[h w]=freqz(hn,1,256)。figure(1)subplot(212)plot(w,20*log10(abs(h)))。xlabel(39。w39。)。ylabel(39。20lg|H(e^jw)|39。)。title(39。Freqency Magnitude(dB)39。)。grid onfigure(2)。plot(w,20*log10(abs(h)),39。.39。)。clear all。w0=8*pi/51。 w0_kai=。 N=2*pi*w0_kai/w0。M=ceil(N1)1。wc=*pi。hn=fir1(M,wc/pi,kaiser(M+1))。[h w]=freqz(hn,1,256)。hold onplot(w,20*log10(abs(h)),39。g39。)。w0_han=。 N=2*pi*w0_han/w0。M=ceil(N1)1。hn=fir1(M,wc/pi,hanning(M+1))。[h w]=freqz(hn,1,256)。plot(w,20*log10(abs(h)),39。m:39。)。w0_box=。N=2*pi*w0_box/w0。M=ceil(N1)1。hn=fir1(M,wc/pi,boxcar(M+1))。[h w]=freqz(hn,1,256)。plot(w,20*log10(abs(h)),39。r39。)。xlabel(39。w39。)。ylabel(39。20lg|H(e^jw)|39。)。title(39。Freqency Magnitude(dB)39。)。legend(39。hamming39。,39。kaiser39。,39。hanning39。,39。boxcar39。)grid onclear all2 給定一理想低通FIR濾波器的頻率特性 希望用窗函數(shù)（矩形窗和哈明窗）法設(shè)計(jì)該濾波器，要求具有線性相位。假定濾波器系數(shù)的長(zhǎng)度為29，即M/2＝14?！境绦虼a】clear all。close all。clcN=29。M=N1。wc=*pi。hn_1=fir1(M,wc/pi,hamming(M+1))figure(1)subplot(211)stem(0:M,hn_1,39。filled39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。Hamming Window39。)。[h_1 w_1]=freqz(hn_1,1,256)。subplot(212)plot(w_1,abs(h_1))。xlabel(39。w39。)。ylabel(39。H39。)。grid onhn_2=fir1(M,wc/pi,boxcar(M+1))figure(2)subplot(211)stem(0:M,hn_2,39。filled39。,39。m39。)。xlabel(39。n39。)。ylabel(39。h(n)39。)。title(39。Boxcar Window39。)。[h_2 w_2]=freqz(hn_2,1,256)。subplot(212)plot(w_2,abs(h_2),39。m39。)。xlabel(39。w39。)。ylabel(39。H39。)。grid onfigure(3)plot(w_1,abs(h_1),w_2,abs(h_2),39。r.39。)。xlabel(39。w39。)。ylabel(39。H39。)。grid onlegend(39。hamming39。,39。boxcar39。)。title(39。contrast39。)1. 若不給自己設(shè)限，則人生中就沒(méi)有限制你發(fā)揮的藩籬。2. 若不是心寬似海，哪有人生風(fēng)平浪靜。在紛雜的塵世里，為自己留下一片純靜的心靈空間，不管是潮起潮落，也不管是陰晴圓缺，你都可以免去浮躁，義無(wú)反顧，勇往直前，輕松自如地走好人生路上的每一步3. 花一些時(shí)間，總會(huì)看清一些事。用一些事情，總會(huì)看清一些人。有時(shí)候覺(jué)得自己像個(gè)神經(jīng)病。既糾結(jié)了自己，又打擾了別人。努力過(guò)后，才知道許多事情，堅(jiān)持堅(jiān)持，就過(guò)來(lái)了。4. 歲月是無(wú)情的，假如你丟給它的是一片空白，它還給你的也是一片空白。歲月是有情的，假如你奉獻(xiàn)給她的是一些色彩，它奉獻(xiàn)給你的也是一些色彩。你必須努力，當(dāng)有一天驀然回首時(shí)，你的回憶里才會(huì)多一些色彩斑斕，少一些蒼白無(wú)力。只有你自己才能把歲月描畫(huà)成一幅難以忘懷的人生畫(huà)卷。