高頻課程設(shè)計(jì)報(bào)告--基于matlab的帶通濾波器設(shè)計(jì)

1、<p><b>　　課程設(shè)計(jì)</b></p><p>　　設(shè)計(jì)題目： 小信號(hào)調(diào)諧放大器電路設(shè)計(jì)、 </p><p>　　基于matlab的帶通濾波器設(shè)計(jì) </p><p>　　學(xué)生姓名： </p><p>　　學(xué) 號(hào)：

2、 </p><p>　　專業(yè)班級(jí)： </p><p>　　指導(dǎo)老師： </p><p>　　實(shí)驗(yàn)組員： </p&

3、gt;<p>　　20 13 年 7 月 6 日</p><p><b>　　目錄</b></p><p><b>　　一、硬件部分</b></p><p><b>　　1、原理1</b></p><p>　　1.1小信號(hào)調(diào)諧放大器的主要特點(diǎn)1</p&

4、gt;<p>　　1.2小信號(hào)調(diào)諧放大器的主要質(zhì)量指標(biāo)1</p><p>　　1.2.1諧振頻率1</p><p>　　1.2.2諧振增益（Av）1</p><p>　　1.2.3通頻帶1</p><p>　　1.2.4增益帶寬積2</p><p>　　1.2.5選擇性2</p>

5、<p>　　1.2.6噪聲系數(shù)3</p><p>　　1.3 單級(jí)單調(diào)諧回路諧振放大器電路原理3</p><p><b>　　電路圖4</b></p><p><b>　　實(shí)驗(yàn)過(guò)程4</b></p><p><b>　　實(shí)驗(yàn)結(jié)果6</b></p

6、><p><b>　　軟件部分</b></p><p><b>　　凱瑟窗7</b></p><p><b>　　1.1 原理7</b></p><p>　　1.2 參數(shù)設(shè)置7</p><p>　　1.3 matlab程序8</p>&

7、lt;p>　　1.4 實(shí)驗(yàn)圖像9</p><p><b>　　頻率采樣法11</b></p><p><b>　　2.1 原理11</b></p><p>　　2.2 參數(shù)設(shè)置11</p><p>　　2.3 matlab程序12</p><p>　　2.4

8、實(shí)驗(yàn)圖像13</p><p><b>　　三、實(shí)驗(yàn)心得15</b></p><p><b>　　四、參考文獻(xiàn)16</b></p><p><b>　　一、硬件部分</b></p><p>　　我們組硬件部分設(shè)計(jì)題目是：小信號(hào)調(diào)諧放大器的設(shè)計(jì)。</p><

9、;p><b>　　1、原理：</b></p><p>　　1.1 小信號(hào)調(diào)諧放大器的主要特點(diǎn)</p><p>　　晶體管集電極負(fù)載通常是一個(gè)由 LC組成的并聯(lián)諧振電路。由于 LC 并聯(lián)諧振回路的阻抗是隨著頻率變化而變化，理論上可以分析，并聯(lián)諧振在諧振頻率處呈現(xiàn)純阻，并達(dá)到最大值。即放大器在回路諧振頻率上將具有最大的電壓增益。若偏離諧振頻率，輸出增益減小?？傊?，

10、調(diào)諧放大器不僅具有對(duì)特定頻率信號(hào)的放大作用，同時(shí)也起著濾波和選頻的作用。</p><p>　　1.2小信號(hào)調(diào)諧放大器的主要質(zhì)量指標(biāo)</p><p>　　衡量小信號(hào)調(diào)諧放大器的主要質(zhì)量主要包括以下幾個(gè)方面：</p><p><b>　　1.2.1諧振頻率</b></p><p>　　放大器調(diào)諧回路諧振時(shí)所對(duì)應(yīng)的頻率稱為放大

11、器的諧振頻率，理論上，對(duì)于 LC組成的并聯(lián)諧振電路，諧振頻率 的表達(dá)式為：</p><p>　　式中，L 為調(diào)諧回路電感線圈的電感量；C 為調(diào)諧回路的總電容。</p><p>　　1.2.2諧振增益（Av）</p><p>　　放大器的諧振電壓增益放大倍數(shù)指：放大器處在在諧振頻率f0下，輸出電壓與輸入電壓之比。 </p><p>　　Av的測(cè)

12、量方法：當(dāng)諧振回路處于諧振狀態(tài)時(shí)，用高頻毫伏表測(cè)量輸入信號(hào)Vi和輸出信號(hào)Vo大小，利用下式計(jì)算：</p><p><b>　　dB</b></p><p><b>　　1.2.3通頻帶</b></p><p>　　由于諧振回路的選頻作用，當(dāng)工作頻率偏離諧振頻率時(shí)，放大器的電壓放大倍數(shù)下降，習(xí)慣上稱電壓放大倍數(shù)Av=Vo/V

13、i下降到諧振電壓放大倍數(shù)Avo的 0.707 倍時(shí)所對(duì)應(yīng)的頻率偏移稱為放大器的通頻帶帶寬BW，通常用2Δf0.7 表示。有時(shí)也稱2Δf0.1為 3dB 帶寬。通頻帶帶寬：</p><p>　　式中，Q為諧振回路的有載品質(zhì)因數(shù)。當(dāng)晶體管選定后，回路總電容為定值時(shí)，諧振電壓放大倍數(shù)fo與通頻帶BW的乘積為一常數(shù)。</p><p>　　圖1-1 放大器的通頻帶和諧振曲線</p>&

14、lt;p>　　1.2.4增益帶寬積</p><p>　　增益帶寬積BW?G也是通信電子電路的一個(gè)重要指標(biāo)，通常，增益帶寬積可以認(rèn)為是一個(gè)常數(shù)。放大器的總通頻帶寬度隨著放大級(jí)數(shù)的增加而變窄，BW越大，增益越小。二者是一對(duì)矛盾。 </p><p>　　不同電路中，放大器的通頻帶差異可能比較大。如：在設(shè)計(jì)電視機(jī)和收音機(jī)的中頻放大器時(shí)，對(duì)帶寬的考慮是不同的，普通的調(diào)幅無(wú)線電廣播所占帶寬是9k

15、Hz，而電視信號(hào)的帶寬需要6.5MHz，顯然，要獲得同樣的增益，中頻放大器的帶寬設(shè)計(jì)是完全不同的。</p><p><b>　　1.2.5選擇性</b></p><p>　　放大器從含有各種不同頻率的信號(hào)總和中選出有用信號(hào)，排除干擾信號(hào)的能力，稱為放大器的選擇性。選擇性的基本指標(biāo)是矩形系數(shù)。其中，定義矩形系數(shù)是電壓放大倍數(shù)下降到諧振時(shí)放大倍數(shù)的10％所對(duì)應(yīng)的頻率偏移和

16、電壓放大倍數(shù)下降為0.707時(shí)所對(duì)應(yīng)的頻率偏移2Δf0.1之比，即：</p><p>　　同樣還可以定義矩形系數(shù)，即：</p><p>　　顯然，矩形系數(shù)越接近1，曲線就越接近矩形，濾除鄰近波道干擾信號(hào)的能力愈強(qiáng)。</p><p><b>　　1.2.6噪聲系數(shù)</b></p><p>　　，NF越接近 1 越好。<

17、;/p><p>　　1.3單級(jí)單調(diào)諧回路諧振放大器電路原理</p><p>　　圖1-2 單級(jí)單調(diào)諧回路諧振放大器</p><p>　　上圖是一個(gè)單級(jí)單調(diào)諧回路諧振放大器的原理圖，理論上分析，諧振時(shí)電壓增益：</p><p>　　放大器的增益可用帶寬表示為：</p><p><b>　　，</b>&l

18、t;/p><p>　　其中單調(diào)諧放大器的選擇性用矩形系數(shù)來(lái)表示為：</p><p>　　所以單調(diào)諧放大器的矩形系數(shù)比 1大得多，選擇性比較差。</p><p><b>　　2、電路圖：</b></p><p><b>　　圖1-3 電路圖</b></p><p><b>

19、;　　圖1-4 仿真圖</b></p><p><b>　　3、實(shí)驗(yàn)過(guò)程：</b></p><p>　　我參與了電路搭建的工作，按照電路圖將元件連接好后，輸入信號(hào)465KHz，在示波器示波器顯示出有抖動(dòng)的波形，剛開始我們只是調(diào)節(jié)中周，但是波形依舊在抖動(dòng)，無(wú)法顯示正確波形，經(jīng)過(guò)陳麗娟老師的指點(diǎn)，我們才知道應(yīng)該是信號(hào)源和中周一起調(diào)節(jié)，在調(diào)節(jié)中周和信號(hào)源輸入頻率

20、的后，波形仍無(wú)好轉(zhuǎn)，于是更換了從別人那里借來(lái)的中周，情況還是沒(méi)有轉(zhuǎn)變，于是我們認(rèn)為是電路搭建出現(xiàn)了問(wèn)題，然后就重新搭建，但是波形和前面的沒(méi)有任何變化，后來(lái)問(wèn)了陳麗娟老師后，我們換了信號(hào)源和信號(hào)輸入的線，調(diào)節(jié)后，出現(xiàn)了正確的波形。</p><p>　　圖1-5實(shí) 物 電 路 圖</p><p><b>　　圖1-6工 作 臺(tái)</b></p><p&g

21、t;　　圖1-7輸入信號(hào)波形</p><p>　　圖1-8調(diào)節(jié)后的波形輸出</p><p><b>　　4、實(shí)驗(yàn)結(jié)果： </b></p><p>　　我們選擇的是輸入A=40mV，諧振頻率f=3.38Mhz，經(jīng)過(guò)放大器后輸出1.18V，f1=3.82MHz，f2=3.07Mhz，通頻帶BW=f1-f2=0.75MHz。</p>

22、<p><b>　　二、軟件部分 </b></p><p>　　軟件部分由xx和我負(fù)責(zé)，我們組選擇的是題目1-B，設(shè)計(jì)帶通濾波器，要求使用兩種方法。</p><p><b>　　凱瑟窗：</b></p><p><b>　　1.1原理：</b></p><p>　　

23、凱瑟窗函數(shù)是一種可調(diào)整的窗函數(shù)，是最有用且最優(yōu)的窗函數(shù)之一。通過(guò)調(diào)整控制參數(shù)可以達(dá)到不同的阻帶最小衰減,并提供最小的主瓣寬度，也就是最窄的過(guò)渡帶。反之，對(duì)于給定的指標(biāo)，凱瑟窗函數(shù)可以使濾波器階數(shù)最小。凱瑟窗函數(shù)由下式給出： </p><p>　　式中，β是調(diào)整參數(shù)，（β）表示零階第一類修正貝塞爾函數(shù)，可以用下式計(jì)算：</p><p>　　實(shí)際中取前2

24、0項(xiàng)就可以滿足精度要求。</p><p>　　參數(shù)β控制加窗設(shè)計(jì)濾波器的阻帶最小衰減。凱瑟給出的估算β和濾波器階數(shù)N的公式如下：</p><p>　　式中，=是數(shù)字濾波器過(guò)渡帶寬度。應(yīng)當(dāng)注意，因?yàn)樯鲜綖殡A數(shù)估算，所以必須對(duì)設(shè)計(jì)結(jié)果進(jìn)行檢驗(yàn)。</p><p><b>　　1.2參數(shù)設(shè)置：</b></p><p>　　根據(jù)題

25、目要求，需要濾除頻率為0.13和0.64的噪聲，通過(guò)頻率為0.3和0.46。因此，我設(shè)置的阻帶和通帶截止頻率分別為wsl=0.14,wpl=0.28,wph=0.48,wsh=0.6,Rs=60dB,故過(guò)渡帶寬為DB=wpl-wsl=0.14,其估算β的式子為0.112（-8.7）。由奈奎斯特抽樣定理得，fs>=2*320=640(Hz)，這里為了得到更好的抽樣效果，選取fs=1000Hz。</p><p>

26、;　　1.3 matlab程序：</p><p>　　%kaiser濾波器%</p><p>　　wsl= 0.14*pi; %臨界頻率采用模擬角頻率表示 </p><p>　　wsh= 0.6*pi; %臨界頻率采用模擬角頻率表示</p><p>　　wpl= 0.28*pi;

27、%臨界頻率采用模擬角頻率表示 </p><p>　　wph=0.48*pi; %臨界頻率采用模擬角頻率表示</p><p>　　rs=60; %阻帶最小衰減</p><p>　　DB=wpl-wsl; %過(guò)渡帶寬</p><p>　　B=0.112*(rs-8.7

28、); %計(jì)算凱瑟窗控制參數(shù)</p><p>　　M=ceil((rs-8)/2.285/DB);%計(jì)算所需h(n)長(zhǎng)度M，ceil(x)取大于等于x的最小整數(shù)</p><p>　　wc=[(wpl+wsl)/2/pi,(wph+wsh)/2/pi];%計(jì)算理想帶通濾波器通帶截止頻率（關(guān)于π歸一化）</p><p>　　h=fir1(M-1,wc,kaise

29、r(M,B)); %調(diào)用fir1計(jì)算帶通FIRDF的h(n)</p><p>　　Fs=1000; %抽樣頻率</p><p>　　[H,w1]=freqz(h,1,M,Fs); %求濾波器幅度響應(yīng)，設(shè)置最大幅度為1</p><p>　　plot(w1,abs(H)); %畫出濾波器幅度響應(yīng)</p><

30、;p>　　grid;xlabel('頻率/Hz');ylabel('幅值');</p><p>　　title('濾波器幅頻響應(yīng)'); %設(shè)置圖像窗口標(biāo)題</p><p>　　figure(2); %創(chuàng)建圖像窗口（2）</p><p>　　freqz(h);

31、%畫圖，濾波器幅度響應(yīng)(db)和相位響應(yīng)</p><p><b>　　%信號(hào)的處理%</b></p><p>　　t=0:0.001:1.999; %設(shè)置t變量范圍，和步長(zhǎng)</p><p>　　n=2000; %抽樣點(diǎn)數(shù)</p><p>　　Fs=1000;

32、 %抽樣頻率</p><p>　　y=cos(2*pi*65*t)+1.2*cos(2*pi*150*t)+2*sin(2*pi*230*t)+1.5*sin(2*pi*320*t);%輸入信號(hào)</p><p>　　y1=fft(y); %輸入信號(hào)的傅里葉變換</p><p>　　y2=fftshift(y1);

33、%調(diào)整頻譜圖</p><p>　　f=(0:1999)*Fs/n-Fs/2; %計(jì)算頻率f</p><p>　　figure(3); %創(chuàng)建圖像窗口（3）</p><p>　　plot(f,abs(y2)/Fs); %畫出輸入信號(hào)頻譜圖</p><p>　　title('輸入信號(hào)頻譜圖

34、9;); %設(shè)置圖像窗口標(biāo)題</p><p><b>　　%輸出的信號(hào)&</b></p><p>　　G=fftfilt(h,y); %輸入信號(hào)y通過(guò)濾波器</p><p>　　G1=fft(G); %濾波后，輸出信號(hào)傅里葉變換</p><p>　　G2=fftshift

35、(G1); %輸出信號(hào)傅里葉變換重新排布，使數(shù)據(jù)與頻率對(duì)應(yīng)</p><p>　　figure(4); %創(chuàng)建圖像窗口（4）</p><p>　　plot(f,abs(G2)/Fs); %畫出輸出信號(hào)頻譜圖</p><p>　　title('輸出信號(hào)頻譜圖'); %設(shè)置圖像窗口標(biāo)題</p>&l

36、t;p><b>　　1.4實(shí)驗(yàn)圖像：</b></p><p>　　圖1-1濾波器幅頻響應(yīng)</p><p>　　圖1-2濾波器幅度響應(yīng)(db)和相位響應(yīng)</p><p>　　圖1-3 輸入信號(hào)頻譜圖</p><p>　　圖1-4 輸出信號(hào)頻譜圖</p><p><b>　　2、頻率采

37、樣法：</b></p><p><b>　　2.1原理：</b></p><p>　　頻率采樣設(shè)計(jì)法屬FIR濾波器設(shè)計(jì)法中的一種, 該設(shè)計(jì)法可以設(shè)計(jì)任意形狀頻率響應(yīng)特性的FIR濾波器。</p><p>　　首先根據(jù)濾波器阻帶最小衰減值確定過(guò)渡帶采樣點(diǎn)m的值，然后求出采樣點(diǎn)N，根據(jù)濾波器通帶和阻帶截止頻率構(gòu)造出一個(gè)符幅序列。<

38、/p><p>　　Ak=[ ones(1,Np),T1,T2 ,zeros(1,Ns-2),T2 ,T1,ones(1,Np-1)];</p><p>　　與窗函數(shù)設(shè)計(jì)方法相同，首先構(gòu)造一個(gè)希望逼近的頻率響應(yīng)函數(shù)：</p><p>　　然后對(duì)在[0,2]上采樣N點(diǎn)得到</p><p>　　，k=0,1,2，...，N-1</p>&

39、lt;p>　　，n=0,1,2, ...,N-1</p><p>　　將h(n)作為設(shè)計(jì)的FIRDF的單位脈沖響應(yīng)，其系統(tǒng)函數(shù)為</p><p>　　以上即為頻率采樣法的基本思想。</p><p><b>　　2.2參數(shù)設(shè)置：</b></p><p>　　根據(jù)題目要求，需要濾除頻率為0.13和0.64的噪聲，通過(guò)

40、頻率為0.3和0.46。為了方便計(jì)算，我們將阻帶和通帶的截止頻率分別設(shè)置為wsl=0.15,wpl=0.25,wph=0.5,wsh=0.6, Rs=50dB。</p><p>　　表一：過(guò)渡帶采樣點(diǎn)個(gè)數(shù)m與濾波器阻帶最小衰減的經(jīng)驗(yàn)數(shù)據(jù)</p><p>　　所以，m=1，根據(jù)公式可以得出N可以取40，從文獻(xiàn)一中有</p><p>　　T1=0.1095，T2=0.5

41、9417456;由奈奎斯特抽樣定理得，fs>=2*320=640(Hz)，這</p><p>　　里為了得到更好的抽樣效果，選取fs=1000Hz。</p><p>　　2.3 Matlab程序：</p><p><b>　　%濾波器%</b></p><p>　　N=40; %詳見頻率采樣

42、法</p><p>　　alfa=(N-1)/2;</p><p><b>　　k=0:N-1;</b></p><p>　　w1=(2*pi/N)*k;</p><p>　　T1=0.1095;T2=0.59417456;%過(guò)渡帶頻率采樣值</p><p>　　hrs=[zeros(1,3),T

43、1,T2,ones(1,6),T2,T1,zeros(1,15),T1,T2,ones(1,6),T2,T1,zeros(1,2)];%構(gòu)造相應(yīng)符幅序列</p><p>　　hdr=[0,0,1,1,0,0]; wd1=[0,0.15,0.25,0.5,0.6,1];</p><p>　　k1=0:floor((N-1)/2); k2=floor((N-1)/2)+1:N-1;&l

44、t;/p><p>　　thetak=[-alfa*(2*pi)/N*k1,alfa*(2*pi/N*(N-k2))];%相位采樣向量</p><p>　　H=hrs.*exp(1j*thetak);%構(gòu)造頻域采樣向量</p><p>　　h=real(ifft(H)); %濾波器函數(shù)</p><p>　　n=0:1:N-1;</p>

45、<p>　　[db,mag,pha,grd,w] = freqz_m(h,1);</p><p>　　[Hr,ww,a,L] =hr_type2(h);</p><p>　　plot(w1/pi,hrs,'o',wd1,hdr);</p><p>　　axis([0,1,-0.1,1.1]);</p><p>　

46、　title('帶通：N=40，T1=0.109021， T2=0.59417456')</p><p>　　xlabel('頻率 （單位：pi）');ylabel('Hr（k）');</p><p>　　figure(2); plot(ww/pi,Hr,w1/pi,hrs,'o');</p><p>

47、　　axis([0,1,-0.1,1.1]);</p><p>　　title('振幅響應(yīng)');xlabel('頻率 （單位：pi）'); ylabel('Hr(w)');</p><p>　　figure(3); plot(w/pi,db); axis([0,1,-100,10]);</p><p>　　title(

48、'幅度響應(yīng)');xlabel('頻率 （單位：pi）'); ylabel('分貝');</p><p><b>　　%信號(hào)的處理%</b></p><p>　　t=0:0.001:1.999; %設(shè)置t變量范圍，和步長(zhǎng)</p><p>　　n=2000;

49、 %抽樣點(diǎn)數(shù)</p><p>　　Fs=1000; %抽樣頻率</p><p>　　y=cos(2*pi*65*t)+1.2*cos(2*pi*150*t)+2*sin(2*pi*230*t)+1.5*sin(2*pi*320*t);%輸入信號(hào)</p><p>　　y1=fft(y); %輸入信號(hào)的傅里葉變換

50、</p><p>　　y2=fftshift(y1); %調(diào)整頻譜圖</p><p>　　f=(0:1999)*Fs/n-Fs/2; %計(jì)算頻率f</p><p>　　figure(4); %創(chuàng)建圖像窗口（3）</p><p>　　plot(f,abs(y2)/Fs); %畫出輸入信號(hào)

51、頻譜圖</p><p>　　title('輸入信號(hào)頻譜圖'); %設(shè)置圖像窗口標(biāo)題</p><p><b>　　%輸出的信號(hào)&</b></p><p>　　G=fftfilt(h,y); %輸入信號(hào)y通過(guò)濾波器</p><p>　　G1=fft(G); %濾波后

52、，輸出信號(hào)傅里葉變換</p><p>　　G2=fftshift(G1); %輸出信號(hào)傅里葉變換重新排布，使數(shù)據(jù)與頻率對(duì)應(yīng)</p><p>　　figure(5); %創(chuàng)建圖像窗口（4）</p><p>　　plot(f,abs(G2)/Fs); %畫出輸出信號(hào)頻譜圖</p><p>　　title(&

53、#39;輸出信號(hào)頻譜圖'); %設(shè)置圖像窗口標(biāo)題 </p><p><b>　　2.4實(shí)驗(yàn)圖像：</b></p><p>　　圖2-1 構(gòu)造符幅序列 </p><p><b>　　圖2-2振幅響應(yīng)</b></p><p>　　圖2-3 幅度響應(yīng) </p><p>

54、　　圖2-4 輸入信號(hào)頻譜圖</p><p>　　圖2-5 輸出信號(hào)頻譜圖</p><p><b>　　實(shí)驗(yàn)心得：</b></p><p>　　這次課程設(shè)計(jì)有著不少的體會(huì)，在硬件方面，我們一調(diào)不出波形就重新搭電路、換中周，這種做法太沖動(dòng)，我們應(yīng)當(dāng)逐點(diǎn)檢測(cè)，在檢測(cè)電路沒(méi)有問(wèn)題的情況下再考慮實(shí)驗(yàn)儀器的問(wèn)題。軟件方面心得是較多的，之前學(xué)習(xí)信號(hào)處理

55、這門課的時(shí)候，注重的是理論知識(shí)，這次課設(shè)是對(duì)理論的實(shí)踐，讓我對(duì)matlab有了更多地了解。掌握了matlab在處理信號(hào)過(guò)程中的部分語(yǔ)法，但是在頻率采樣法中還存在個(gè)問(wèn)題未能解決，在構(gòu)造符幅序列時(shí)，對(duì)于zeros()和ones()函數(shù)括號(hào)里的數(shù)值怎么取還存在疑問(wèn)，一開始我是按照課本上低通濾波器中的公式來(lái)計(jì)算的 但是發(fā)現(xiàn)這樣做出來(lái)的濾波器有較大的變形，后來(lái)在參考資料上看到了另一種取法，試了之后還是不行，這個(gè)括號(hào)里的值如何去取，是一個(gè)后面需要

56、思考的問(wèn)題。</p><p><b>　　參考資料：</b></p><p>　　[1] 《數(shù)字信號(hào)處理——原理、實(shí)現(xiàn)及應(yīng)用》；高西全，丁玉美，闊永紅；電子工業(yè)出版社；2008. </p><p>　　[2]《基于MATLAB的頻率采樣法設(shè)計(jì)FIR濾波器》；龍安；廣西科學(xué)院學(xué)報(bào)；2007.2</p><p>　　[3]《