模擬有源濾波器設計的MATLAB實現(xiàn)

PAGEword文檔可自由編輯word文檔可自由編輯摘要本文闡述了濾波器的基本概念，介紹了模擬有源濾波器的設計原理和逼近理論。其中包括巴特沃思逼近、切比雪夫I型逼近、切比雪夫Ⅱ型逼近、橢圓函數(shù)逼近和貝塞爾逼近。并研究了模擬有源濾波器的設計流程及性能測試。綜合了傳統(tǒng)的硬件設計方法與軟件編程技術,由MATLAB仿真出了各種濾波器逼近技術的幅頻特性曲線并進行了實例分析。對巴特沃思濾波器實例的研究仿真，由程序快速的得到了最小階數(shù)和截止頻率，取代了傳統(tǒng)繁復的計算；方便的實現(xiàn)了由模擬低通濾波器向高通、低通和帶阻濾波器的轉換；對四運放復雜電路進行了設計仿真，通過求取其不同點的輸出傳遞函數(shù)，模擬了二階低通、高通、帶通和帶阻濾波器的幅頻特性曲線并得到了較好的仿真結果。關鍵詞：模擬有源濾波器；逼近理論；幅頻特性；MATLAB程序設計ABSTRACTThispaperdescribesthebasicconceptsoffiltersintroducedanalogactivefilterdesignprinciplesandapproximationtheory.IncludingButterworthapproximation,ChebyshevtypeIapproximation,ChebyshevⅡtypeapproximation,ellipticfunctionapproximationandBezierapproximation.Thefirstauthorstudiedtheanalogactivefilterdesignprocessandperformancetesting.Combinesthetraditionalhardwaredesignmethodsandsoftwareprogrammingtechniques，theMATLABsimulationofavarietyoffilterapproximationtechniqueofamplitude-frequencycharacteristiccurveandanillustrativeexample.Butterworthfilterinstanceonsimulationstudies,bytheprogramtoquicklygettheminimumorderandthecutofffrequency,replacingthetraditionalcomplexcalculations;convenientlyachievedbytheanaloglow-passfiltertothehigh-pass,low-passandband-stopfiltertheconversion;complexonquadopampcircuitdesignsimulation,throughitsdifferentpointsstrikeoutputtransferfunctiontosimulatethesecond-orderlow-pass,highpass,bandpassandband-stopfilterfrequencycharacteristiccurveandgetbettersimulationresults.Keywords:AnalogandActiveFilter;TheoryofApproximation,Amplitude-frequencycharacteristics，MATLABprogramdesignCompiling

and

organizing

data

After

you

have

established

the

purpose

of

the

report,

you

need

to

compile

and

organize

the

information

needed

to

support

it.

The

gathering

of

information

may

have

given

you

a

lot

of

materials,

but

you

need

to

be

a

firm

editor

and

retain

only

the

essential

data

and

throw

out

the

rest.

Consider

your

readers,

and

think

about

how

much

background

information

they

will

need.

Writing

the

report

A

report

consists

of

five

parts:

title,

introduction,

findings,

conclusion

and

recommendations.

1.

Title

This

should

run

no

longer

than

one

line.

2.

Introduction

This

tells

the

reader

what

the

purpose

and

objective

of

the

report

are.

It

might

also

give

the

reader

some

background

information

on

the

subject.

The

purpose/objective/aim

of

this

report

is

to...

This

report

aims

to/is

intended

to...

3.

Findings

This

is

the

main

part

of

the

report.

It

tells

the

reader

what

you

have

found

in

your

information

gathering.

It

is

important

to

differentiate

between

fact

and

opinion.

Which

of

the

following

phrases

report

facts

and

which

report

opinions?

We

found

that...

(could

be

both

fact

and

opinion)

It

clearly

shows

that...

(fact)

It

was

found

that...(fact)

We

discovered

that...(fact)

We

observed

that...(opinion)

There

is

evidence

that...(fact)

4.

Conclusion

This

part

tells

the

reader

about

the

results

of

the

report

based

on

the

findings.

It's

concluded/decided/agreed/felt

that...

It

can

be

seen

that...

In

conclusion...

No

conclusions

were

reached

regarding?...(關于一事未得出任何結論)

We

can

conclude

that...

5.

Recommendations

Finally,

recommendations

are

made

on

what

future

actions

need

to

be

taken.

Based

on

our

findings,

we

(would)

recommend

that...

It

is

recommend-ed/proposed/suggested

that?...

It

seems

to

me

that...title,

introduction,

findings,

conclusion

and

recommendations.PAGEword文檔可自由編輯1Introduction濾波是信號處理的一種最基本而重要的技術，利用濾波可從復雜的信號中提取所需要的信號，抑制不需要的部分。所謂濾波器是具有一定傳輸特性的信號處理裝置。根據(jù)濾波器所處理的信號不同，濾波器可分為模擬濾波器和數(shù)字濾波器兩類。模擬濾波器在測試系統(tǒng)或專用儀器儀表中是一種常用的變換裝置。例如：帶通濾波器用作頻譜分析儀中的選頻裝置；低通濾波器用作數(shù)字信號分析系統(tǒng)中的抗頻混濾波；高通濾波器被用于聲發(fā)射檢測儀中剔除低頻干擾噪聲；帶阻濾波器用作電渦流測振儀中的陷波器，等等。與模擬濾波器相對應，在離散系統(tǒng)中廣泛應用數(shù)字濾波器。它的作用是利用離散時間系統(tǒng)的特性對輸入信號波形或頻率進行加工處理?；蛘哒f，把輸入信號變成一定的輸出信號，從而達到改變信號頻譜的目的。數(shù)字濾波器一般可以用兩種方法來實現(xiàn)：一種方法是用數(shù)字硬件裝配成一臺專門的設備，這種設備稱為數(shù)字信號處理機；另一種方法就是直接利用通用計算機，將所需要的運算編成程序讓通用計算機來完成，即利用計算機軟件來實現(xiàn)。本論文首先介紹了濾波器的濾波原理以及模擬濾波器設計方法。重點介紹了模擬濾波器的設計和仿真。系統(tǒng)研究了模擬濾波器（包括巴特沃斯濾波器、切比雪夫濾波器、橢圓函數(shù)濾波器和貝塞爾濾波器）的設計原理和方法，并在此基礎上論述了模擬濾波器（包括低通、高通、帶通、帶阻）的設計。在此基礎上，用MATLAB虛擬實現(xiàn)模擬濾波器。此設計擴展性好，便于調節(jié)濾波器的性能，可以根據(jù)不同的要求在MATLAB上加以實現(xiàn)。通過MATLAB的仿真與實現(xiàn)，可以看出傳統(tǒng)的模擬濾波器設計方法繁瑣且不直觀，而MATLAB具有較嚴謹?shù)目茖W計算和圖形顯示這一優(yōu)點，使設計結果顯示的更加直觀，而且對濾波器的精度也有了很大的提高，能更好的達到預期效果。同時，又對模擬濾波器低通至高通、帶通、帶阻的轉換進行了理論上的闡述。2濾波器的基本概念2.1濾波原理濾波器是一種選頻裝置，可以使信號中特定的頻率成分通過，而極大地衰減其它頻率成分。在測試裝置中，利用濾波器的這種選頻作用，可以濾除干擾噪聲或進行頻譜分析。廣義地講，任何一種信息傳輸?shù)耐ǖ溃劫|）都可視為是一種濾波器。因為，任何裝置的響應特性都是激勵頻率的函數(shù)，都可用頻域函數(shù)描述其傳輸特性。因此，構成測試系統(tǒng)的任何一個環(huán)節(jié)，諸如機械系統(tǒng)、電氣網絡、儀器儀表甚至連接導線等等，都將在一定頻率范圍內，按其頻域特性，對所通過的信號進行變換與處理。按照濾波器處理信號的性質分為，模擬濾波器和數(shù)字濾波器。本文所述內容屬于模擬濾波范圍。主要介紹模擬濾波器（連續(xù)時不變系統(tǒng)）原理、種類、數(shù)學模型、主要參數(shù)、ＲＣ濾波器設計。盡管數(shù)字濾波技術已得到廣泛應用，但模擬濾波在自動檢測、自動控制以及電子測量儀器中仍被廣泛應用。2.2濾波器分類2.2.1根據(jù)濾波器的選頻作用分類⑴低通濾波器從頻率之間，幅頻特性平直，它可以使信號中低于的頻率成分幾乎不受衰減地通過，而高于f2的頻率成分受到極大地衰減。圖1低通濾波器圖2高通濾波器⑵高通濾波器與低通濾波相反，從頻率，其幅頻特性平直。它使信號中高于的頻率成分幾乎不受衰減地通過，而低于的頻率成分將受到極大地衰減。⑶帶通濾波器它的通頻帶在之間。它使信號中高于而低于的頻率成分可以不受衰減地通過，而其它成分受到衰減。⑷帶阻濾波器與帶通濾波相反，阻帶在頻率之間。它使信號中高于而低于的頻率成分受到衰減，其余頻率成分的信號幾乎不受衰減地通過。圖3帶通濾波器圖4帶阻濾波器低通濾波器和高通濾波器是濾波器的兩種最基本的形式，其它的濾波器都可以分解為這兩種類型的濾波器，例如：低通濾波器與高通濾波器的串聯(lián)為帶通濾波器，低通濾波器與高通濾波器的并聯(lián)為帶阻濾波器。圖5低通濾波器與高通濾波器的串聯(lián)圖6低通濾波器與高通濾波器的并聯(lián)2.2.2根據(jù)“最佳逼近特性”標準分類巴特沃斯濾波器從幅頻特性提出要求，而不考慮相頻特性。巴特沃斯濾波器具有最大平坦幅度特性，其幅頻響應表達式為：（1）n為濾波器的階數(shù)；wc為濾波器的截止角頻率，當w=wc時，|H(wc)|2=1/2,所以，wc對應的是濾波器的-3db點。巴特沃思低通濾波器是以巴特沃思函數(shù)作為濾波器的傳遞函數(shù)H(s)，以最高階泰勒級數(shù)的形式逼近濾波器的理想矩形特性。⑵切比雪夫濾波器切貝雪夫濾波器也是從幅頻特性方面提出逼近要求的，其幅頻響應表達式為：（2）ε是決定通帶波紋大小的波動系數(shù)，0<ε<1，波紋的產生是由于實際濾波網絡中含有電抗元件；wc是通帶截止頻率，Tn是n階切貝雪夫多項式。與巴特沃斯逼近特性相比較，這種特性雖然在通帶內有起伏，但對同樣的n值在進入阻帶以后衰減更陡峭，更接近理想情況。ε值越小，通帶起伏越小，截止頻率點衰減的分貝值也越小，但進入阻帶后衰減特性變化緩慢。切貝雪夫濾波器與巴特沃斯濾波器進行比較，切貝雪夫濾波器的通帶有波紋，過渡帶輕陡直，因此，在不允許通帶內有紋波的情況下，巴特沃斯型更可??；從相頻響應來看，巴特沃斯型要優(yōu)于切貝雪夫型，通過上面二圖比較可以看出，前者的相頻響應更接近于直線。⑶貝塞爾濾波器只滿足相頻特性而不關心幅頻特性。貝塞爾濾波器又稱最平時延或恒時延濾波器。其相移和頻率成正比，即為一線性關系。但是由于它的幅頻特性欠佳，而往往限制了它的應用。2.2.3按濾波器元件性質無源濾波器(R、L、C)有源濾波器(含運放)3模擬有源濾波器的設計原理3.1信號無失真?zhèn)鬏數(shù)臈l件信號通過一個系統(tǒng)，其響應若不失真，則這個系統(tǒng)稱為無失真?zhèn)鬏斚到y(tǒng)。實際中常常需要無失真?zhèn)鬏斚到y(tǒng)，例如，高保真音響系統(tǒng)、示波器等。無失真?zhèn)鬏斚到y(tǒng)對輸入信號的響應應是的準確的復制品，或者說兩者的大小和出現(xiàn)的時間可以不同(見圖7a、b)，但它們的波形變化規(guī)律應相同，即（3）式中、為常數(shù)，是系統(tǒng)增益，為延遲時間。上式可以改寫為（4）這說明無失真?zhèn)鬏斚到y(tǒng)的單位沖激響應為（5）對上式進行傅里葉變換，則有（6）由上式可知，無失真?zhèn)鬏斚到y(tǒng)應滿足下面兩個條件：(1)幅頻特性（7）即系統(tǒng)增益為一個常數(shù)，與信號頻率無關，如圖7(c)所示。(2)相頻特性（8）即系統(tǒng)的相位和與信號頻率成正比，如圖7(d)所示。(c) (d)圖7無失真?zhèn)鬏斚到y(tǒng)實際系統(tǒng)，幅頻特性和相頻特性均不滿足無失真?zhèn)鬏斚到y(tǒng)的兩個條件。工程上，只要信號在占有的頻率范圍內，系統(tǒng)的幅頻特性和相頻特性基本上滿足無失真?zhèn)鬏斚到y(tǒng)的兩個條件，就可以認為是無失真?zhèn)鬏斚到y(tǒng)了。3.2理想濾波器的特性濾波器是一個選頻裝置。理想濾波器應能無失真地傳輸有用信號，而又能完全抑制無用信號。有用信號和無用信號往往占有不同的頻帶。信號能通過濾波器的頻帶稱為通帶(Passband)。信號被抑制的頻帶稱為阻帶（Stopband）。理想濾波器頻率特性可寫為：但通過后面的講解我們會知道理想濾波器是物理不可實現(xiàn)系統(tǒng)。實際濾波器的頻率特性只能“逼近”理想濾波器。圖8為低通濾波器的幅頻特性示意圖。圖8低通濾波器的幅頻特性示意圖可見，濾波器的幅頻響應在通帶內不是完全平直的，而是呈波紋變化；在阻帶內，幅頻特性也不為零，而是衰減至某個值；在通帶和阻帶之間存在一個過渡帶，而不是突然下降。通常，實際設計要求濾波器的技術指標包括通帶波紋Rp(Passbandripple)(dB)、阻帶衰減(Stopbandattenuation)(dB)、通帶邊界頻率、阻帶邊界頻率、過渡帶寬。濾波器的通帶波紋Rp為相對于頻率響應最大點（一般為1）的下降，因此下降越少說明通帶越平直，濾波器的濾波效果越好（通常為1-5dB）。濾波器的阻帶衰減Rs也是相對于頻率響應最大點（一般為1）的下降，因此下降越多說明信號在阻帶內越不容易通過，因此濾波效果越好(通常要大于15dB)。過渡帶寬越窄，濾波器的頻率特性越接近于直角矩形特性，濾波效果也越好。3.3模擬濾波器傳遞函數(shù)設計原理模擬濾波器的設計理論通常在Laplace域內進行討論，本節(jié)內容我們只討論其應用。模擬濾波器的技術指標可由平方幅值響應函數(shù)形式給出，而和傳遞函數(shù)H(s)存在下面關系：（9）即（10）當給定模擬濾波器的技術指標后，由求出A(-s2)，再適當?shù)胤峙淞銟O點可求出H(s)。為了使濾波器穩(wěn)定，H(s)的極點必須落在s平面左半平面，這是因為對于一個濾波器的極點p，就可以寫成的形式，其逆Laplace變換（對應于時間域）為，若p>0，則隨著時間增大至無窮，該濾波器的輸出將出現(xiàn)不穩(wěn)定。濾波器的零點選擇可任取A(-s2)的一半零點，這是因為濾波器對Laplace域表示的傳遞函數(shù)并無特殊要求，但如果要求H(s)具有最小相位，零點也必須選擇在s左半平面。3.3.1模擬濾波器設計步驟用戶對設計的濾波器提出設計要求，我們可以針對濾波器的設計要求設計濾波器。通常用戶對模擬濾波器提出的要求有：（1）濾波器的性能指標，包括截止頻率(對于低通和高通)或上下邊界頻率、，通帶波紋、阻帶衰減等；（2）濾波器的類型，通常為Butterworth、ChebyshevI、ChebyshevII、Elliptic或Bessel濾波器。我們根據(jù)濾波器的類型通常按下列步驟設計濾波器。(1)給定模擬濾波器的性能指標，如截止頻率(對于低通和高通)或上下邊界頻率；通帶波紋、阻帶衰減以及濾波器類型等（用戶給定）。(2)確定濾波器階數(shù)(3)設計模擬低通原型濾波器。MATLAB信號處理工具箱的濾波器原型設計函數(shù)有butterap、cheb1ap、cheb2ap、ellipap、besselap。(4)按頻率變換設計模擬濾波器（低通、高通、帶通、帶阻）。MATLAB信號處理工具箱的頻率變換函數(shù)有l(wèi)p2lp、lp2hp、lp2bp、lp2bs。3.3.2模擬濾波器設計函數(shù)上面濾波器的設計步驟比較麻煩，根據(jù)設計要求求解濾波器的最小階數(shù)和邊界頻率之后需要設計模擬原型濾波器并進行頻率轉換。其實MATLAB將這一系列的過程組合成了更為方便的設計函數(shù)：butter,cheby1,cheby2,ellip,besself。這些函數(shù)稱為模擬濾波器完全設計函數(shù)。用戶在求得濾波器的最小階數(shù)和截止頻率之后只需調用一次完全設計函數(shù)就可以自動完成所有設計過程，編程十分簡單。這些工具函數(shù)適用于模擬濾波器的設計，但同樣也適用于數(shù)字濾波器。本節(jié)只討論這些函數(shù)在模擬濾波器設計中的應用。但要注意，MATLAB是將上述一系列的步驟復合而已，并不是一種新的設計模擬濾波器的方法。[b,a]=butter(n,wn[,'ftype'],'s')[z,p,k]=butter(n,wn[,'ftype'],'s')[b,a]=cheby1(n,Rp,wn[,'ftype'],'s')[z,p,k]=cheby1(n,Rp,wn[,'ftype'],'s')[b,a]=cheby2(n,Rs,wn[,'ftype'],'s')[z,p,k]=cheby2(n,Rs,wn[,'ftype'],'s')[b,a]=ellip(n,Rp,Rs,wn[,'ftype'],'s')[z,p,k]=cheby2(n,Rp,Rs,wn[,'ftype'],'s')[b,a]=besself(n,wn[,'ftype'],'s')[z,p,k]=besself(n,wn[,'ftype'],'s')在上面的調用方式中，n為濾波器的階數(shù)，wn為濾波器的截止頻率，單位rad/s（wn>0）；’s’為模擬濾波器，缺省時為數(shù)字濾波?！甪type’濾波器的類型可取為：‘high’高通濾波器，截止頻率為wn?！畇top’帶阻濾波器，截止頻率為wn=[w1,w2](w<w2)。‘ftype’缺省時為低通或帶通濾波器。當設計帶通濾波器時，截止頻率也為wn=[w1,w2](w1<w2)。a,b分別為濾波器的傳遞函數(shù)分子和分母多項式系數(shù)向量；z,p,k分別為濾波器的零點、極點和增益。Rp,Rs分別為所設計濾波器的通帶波紋和阻帶衰減,單位為dB。濾波器的傳遞函數(shù)具有下面的形式：（11）若濾波器為帶通或帶阻型，則濾波器的階數(shù)為2n，否則階數(shù)為n。3.4濾波器最小階數(shù)選擇前面所述的模擬濾波器設計中，濾波器階數(shù)是我們在編程時任意指定的。其實它是決定濾波器品質的主要參數(shù)之一。通常在滿足性能指標的前提下，階數(shù)應該盡可能小，以滿足易于實現(xiàn)、提高運算速度的要求。而在濾波器階數(shù)和濾波器性能之間存在一定的函數(shù)關系，我們通過這一函數(shù)關系可以求出滿足濾波性能指標的最低階數(shù)。下面介紹常用低通濾波器最小階數(shù)的確定原理及MATLAB實現(xiàn)；接著介紹MATLAB信號處理工具箱中用來計算最小階數(shù)和截止頻率的工具函數(shù)。3.4.1濾波器最小階數(shù)選擇原理以Butterworth低通模擬濾波器為例模擬低通濾波器的設計指標為：通帶邊界頻率，阻帶邊界頻率，通帶波紋Rp(dB)、阻帶衰減Rs(dB)。當時（12）以截止頻率（幅值下降3dB）為1，化為相對的相對頻率，則上式可寫成：(13)

同理，當時，(14)綜合上面兩式可得：式中，N應向上取整，則N為該濾波器的最小階數(shù)。表示為：或(15)3.4.2濾波器最小階數(shù)選擇函數(shù)上面給出了Butterworth濾波器的最小階數(shù)和截止頻率的選擇公式及程序。其實MATLAB工具箱中運用濾波器的最小階數(shù)選擇公式給出了濾波器最小階數(shù)選擇函數(shù)。幾種濾波器最小階數(shù)的選擇函數(shù)如下：[n,wc]=buttord(wp,ws,Rp,Rs,'s');Butterworth濾波器[n,wc]=cheb1ord(wp,ws,Rp,Rs,'s');ChebyshevI濾波器[n,wc]=cheb2ord(wp,ws,Rp,Rs,'s');ChebyshevII濾波器[n,wc]=ellipord(wp,ws,Rp,Rs,'s');Elliptic濾波器式中，wp為通帶邊界頻率，ws為阻帶邊界頻率，單位為rad/s。Rp,Rs分別為通帶波紋和阻帶衰減，單位為dB。二者分別表示通帶內的最大允許幅值損失和阻帶下降的分貝數(shù)。's'表示模擬濾波器（缺省時，該函數(shù)適用于數(shù)字濾波器）；函數(shù)返回值n為模擬濾波器的最小階數(shù)；wc為模擬濾波器的截止頻率，單位為rad/s。這四個函數(shù)適用于低通、高通、帶通、帶阻濾波器。若wp<ws，對應于低通模擬濾波器，當wp>ws時對應于高通模擬濾波器，對于帶通和帶阻濾波器存在兩個過渡帶，wp和ws均應該為含有兩個元素的向量，分別表示兩個過渡帶的邊界頻率。即：wp=[通帶下界頻率,通帶上界頻率]，ws=[阻帶下界頻率,阻帶上界頻率]。對于帶通濾波器，這四個頻帶界線的大小排列為:阻帶下界頻率<通帶下界頻率<通帶上界頻率<阻帶上界頻率；對于帶阻濾波器，這四個頻帶界線的大小排列為:通帶下界頻率<阻帶下界頻率<阻帶上界頻率<通帶上界頻率。這時返回值wc包括兩個元素（第一個元素小于第二個元素），分別為通帶和阻帶之間的界線頻率。函數(shù)自動判斷是帶通還是帶阻濾波器。3.5模擬濾波器的性能測試濾波器設計好之后，一般要進行各方面的測試。在正式設計濾波器之前，我們先介紹如何測試濾波器的性能。對于模擬濾波器，在本章第二節(jié)我們已采用函數(shù)freqs來求模擬濾波器的頻率響應，這里我們詳細介紹該函數(shù)，若其調用格式為：式中，b,a分別為模擬濾波器傳遞函數(shù)分子和分母多項式系數(shù)；h為對應頻率點的傳遞函數(shù)值。上面的表示中，[,w]和[,n]表示可有可無的參數(shù)，本書中的該類表示均為可有可無的參數(shù),由用戶根據(jù)需要確定。w為頻率點的值,n為頻率點數(shù)。若n=128，則用128個頻率點來給出此模擬濾波器的頻率特性（給定頻率點的傳遞函數(shù)值），缺省時為200。若該函數(shù)不寫輸出變量，則執(zhí)行后繪出該濾波器的幅頻響應和相位響應圖。此函數(shù)模擬濾波器的傳遞函數(shù)形式為：MATLAB工具箱還提供了兩個函數(shù)abs和angle，由頻率響應求幅頻響應和相頻響應。其中angle的輸出單位為rad?？刹捎胷ad2deg函數(shù)轉化為度。另外注意函數(shù)的幅頻響應經常用分貝(dB)來表示。求出的幅頻響應后，通過下式轉換為分貝:(dB)。我們知道，除了用傳遞函數(shù)描述濾波器特性外，還可用脈沖(沖激)響應來描述濾波器，因為在模擬濾波器中，脈沖響應與傳遞函數(shù)是Laplace變換對。此外還可以用階躍響應（輸入一個階躍時系統(tǒng)的輸出）來描述濾波器特性。下面我們介紹在MATLAB中如何得到模擬濾波器的脈沖響應和階躍響應。將濾波器的傳遞函數(shù)表示成分子和分母多項式系數(shù)的形式，如分子和分母多項式的系數(shù)為[b,a]，則在MATLAB中用H=[tf(b,a)]來表示此模擬濾波器，采用[[y,t]=]impulse(H)給出該系統(tǒng)的模擬脈沖響應。采用[[y,t]=]step(H)來得到該系統(tǒng)的階躍響應。這兩個函數(shù)與freqs一樣，若沒有輸出則程序自動繪圖模擬該濾波器的脈沖響應或階躍響應。輸出值y,t分別為該濾波器的脈沖響應或階躍響應及其對應的時間序列。另外，我們還可以運用一個輸入信號模擬該濾波器的輸出。若給定濾波器的輸入值序列和對應的時間序列為u,t，則系統(tǒng)的輸出可用y=lsim(H,u,t)來模擬,y為對應t的輸出。若該函數(shù)沒有輸出變量則程序自動繪圖顯示。3.6頻率變換模擬原型濾波器(下一章將講到)均是截止頻率為1的濾波器，在實際設計中是很難遇到的，然而它是設計其他各類濾波器的基礎。我們通常遇到的是截止頻率任意的低通濾波器、高通濾波器、帶通濾波器和帶阻濾波器。如何由低通原型模擬濾波器設計這些濾波器呢？這就要用到我們今天要講的頻率變換。所謂頻率變換是指各類濾波器（低通、高通、帶通、帶阻）和低通濾波器原型的傳遞函數(shù)中頻率自變量之間的變換關系。通過頻率變換，我們可以從模擬低通濾波器原型獲得模擬的低通濾波器、高通濾波器、帶通濾波器和帶阻濾波器，再借助于s域至z域的變換關系又可以設計各類后面所講的無限沖激響應數(shù)字濾波器，這是濾波器設計的重要方法之一。3.6.1頻率變換工具函數(shù)MATLAB信號處理工具箱有l(wèi)p2lp,lp2hp,lp2bp,lp2bs四個頻率變換函數(shù)。下面分別敘述其使用方法及各參量的意義。(1)函數(shù)lp2lp用于實現(xiàn)用低通模擬原型濾波器至低通濾波器的頻率變換，調用格式為：其中，a,b為模擬原型濾波器的分母和分子多項式的系數(shù)，為低通濾波器所期望的截止頻率（rad/s），若給定的單位為Hz，應乘以。bt,at為返回的低通濾波器的分母和分子多項式的系數(shù)。該函數(shù)將模擬原型濾波器傳遞函數(shù)模擬原型執(zhí)行下面變換：式中，H(p)為低通原型濾波器傳遞函數(shù)，H(s)為低通濾波器傳遞函數(shù)。該項操作可以執(zhí)行模擬原型濾波器由截止頻率為1到指定截止頻率的變換，其原理討論以超出本課程的范圍，可參看其他信號處理參考書。(2)函數(shù)lp2hp用于實現(xiàn)由低通模擬濾波器至高通濾波器的頻率變換。調用格式：式中，為高通模擬濾波器所期望的截止頻率（rad/s），若給定的頻率單位為Hz，應乘以。該函數(shù)將模擬原型濾波器傳遞函數(shù)執(zhí)行下面變換：(3)函數(shù)lp2bp用于實現(xiàn)由低通模擬原型濾波器至帶通濾波器的頻率變換。調用格式：式中，為帶通濾波器的中心頻率(rad/s)，Bw為帶通濾波器帶寬(rad/s)。而式中,為帶通濾波器的下邊界頻率，為帶通濾波器上邊界頻率。若給定的邊界頻率為Hz需乘以。該函數(shù)將模擬原型濾波器傳遞函數(shù)執(zhí)行下面變換運算：(4)函數(shù)lp2bs用于實現(xiàn)由低通模擬原型濾波器至帶阻濾波器的頻率變換。調用格式：式中，為帶阻濾波器的中心頻率(rad/s)，Bw為帶阻濾波器帶寬(rad/s)。而式中,為帶阻濾波器的下邊界頻率，為帶阻濾波器上邊界頻率。若給定的邊界頻率為Hz需乘以。該函數(shù)將模擬原型濾波器傳遞函數(shù)執(zhí)行下面變換運算：注意：輸出的帶阻濾波器和帶通濾波器是濾波器原型階數(shù)的2倍。4模擬濾波器的逼近理論由于實際的濾波器與理想的濾波器存在差距，所以我們每個濾波器都有一個原型與之對應。本章介紹常用的模擬原型濾波器的主要特點及其MATLAB實現(xiàn)，包括Butterworth、ChebyshevI,ChebyshevII,Elliptical、Bessel原型低通濾波器的設計。模擬原型濾波器指的是截止頻率為1的濾波器。后面所講的各類模擬濾波器和數(shù)字濾波器可通過這些低通原型濾波器變換得到。4.1Butterworth濾波器Butterworth模擬低通濾波器的平方幅頻響應函數(shù)為：（16）式中，為低通濾波器的截止頻率，N為濾波器的階數(shù)。Butterworth濾波器的特點：通帶內具有最大平坦的頻率特性，且隨著頻率增大平滑單調下降；階數(shù)愈高，特性愈接近矩形，過渡帶愈窄。傳遞函數(shù)無零點，極點等距離地分布在以為半徑的圓周上。所謂濾波器的零點就是將該點的值代入傳遞函數(shù)后，傳遞函數(shù)的值為零。所謂函數(shù)的極點就是將該點的值代入傳遞函數(shù)后，傳遞函數(shù)的值為無窮大。濾波器的增益是指傳遞函數(shù)表達式前的系數(shù)。若系統(tǒng)的傳遞函數(shù)表示為：(17)則濾波器零點為，極點為，濾波器的增益為K。這里所說的零點和極點分布在一個圓上為Laplace域中的形式，感興趣的同學可查看數(shù)學中的Laplace變換。MATLAB信號處理工具箱提供Butterworth模擬低通濾波器原型設計函數(shù)buttap,函數(shù)調用形式為：式中，n為butterworth濾波器階數(shù)；z,p,k分別為濾波器的零點、極點和增益。Butterworth濾波器傳遞函數(shù)具有下面的形式：（18）濾波器沒有零點，極點為，濾波器的增益為K。在有關模擬濾波器設計的MATLAB程序中，經常遇到一些特定函數(shù)。為將模擬原型濾波器函數(shù)（如buttap）設計出的零點z,極點p和增益k形式轉換為傳遞函數(shù)（transferfunction）形式；其中，b為濾波器傳遞函數(shù)分子多項式系數(shù)，a為濾波器傳遞函數(shù)分母多項式系數(shù)。求出傳遞函數(shù)形式（分子和分母多項式的系數(shù)為b,a）表示的濾波器的對應于頻率點w的復數(shù)頻率響應H（包括實部和虛部），這里w為一個矢量，表示對應的角頻率。若該函數(shù)不寫輸出變量，則執(zhí)行后繪出該濾波器的幅頻響應和相頻響應圖。圖9Butterworth濾波器原型平方幅頻圖由圖9知，Butterworth濾波器的幅頻平方特性隨著頻率單調下降。隨著濾波器階數(shù)的增大，其幅頻特性逐漸接近矩形。4.1.1Butterworth濾波器應用設計一個模擬帶通濾波器，我們將設計指標改為：通帶頻率：2000~3000Hz，兩側過渡帶寬500Hz，通帶波紋1dB，阻帶衰減100dB。檢驗信號，其中f1=100Hz，f2=2500Hz，f3=3900Hz。信號的采樣頻率為10000Hz。程序如下：wp=[20003000]*2*pi;ws=[15003500]*2*pi;Rp=1;Rs=100;%濾波器設計參數(shù)，對于給定HZ應乘以[N,Wn]=buttord(wp,ws,Rp,Rs,'s');%求得濾波器的最小階數(shù)和截止頻率w=linspace(1,4000,1000)*2*pi;%設置繪制頻率響應的頻率點[b,a]=butter(N,Wn,'s');%設計模擬Butterworth濾波器H=freqs(b,a,w);%計算給定頻率點的復數(shù)頻率響應magH=abs(H);phaH=unwrap(angle(H));%計算幅頻響應和相頻響應plot(w/(2*pi),20*log10(magH));%以頻率為橫坐標繪制幅頻響應xlabel('Frequency/Hz');ylabel('Magnitude/dB');title('Butterworthanalogbandpassfilter');holdon;plot([15001500],ylim);plot([35003500],ylim);%繪制阻帶邊界gridonfigure(2)dt=1/10000;%模擬信號采樣間隔f1=100;f2=2500;f3=3900;%輸入信號的三個頻率成分t=0:dt:0.04;%給定模擬時間段x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t)+0.5*sin(2*pi*f3*t);%輸入信號H=[tf(b,a)];%濾波器在MATLAB系統(tǒng)中的表示[y,t1]=lsim(H,x,t);%模擬輸出subplot(211),plot(t,x),title('Inputsignal')%繪出輸入信號subplot(212),plot(t1,y)%繪制輸出信號title('Outputsignal'),xlabel('Time/s')圖10模擬濾波器幅頻響應，阻帶邊界由圖10可見，濾波器的幅頻響應完全符合要求。在阻帶邊界確實下降到100dB。3900Hz和100Hz的頻率均在阻帶內，只有25000Hz的頻率在通帶內。圖11濾波器的輸入信號和輸出信號由圖11可見，通過模擬系統(tǒng)后，3900Hz的高頻成分和100Hz的低頻成分被濾除了，輸出只含有2500Hz頻率成分的振動。4.2切比雪夫（chebyshev）濾波器巴特沃茲濾波器在通帶內幅度特性是單調下降的，如果階次一定，則在靠近截止處，幅度下降很多，或者說，為了使通帶內的衰減足夠小，需要的階次N很高，為了克服這一缺點，采用切比雪夫多項式來逼近所希望的。切比雪夫濾波器的在通帶范圍內是等幅起伏的，所以在同樣的通常內衰減要求下，其階數(shù)較巴特沃茲濾波器要小。4.2.1ChebyshevI型ChebvshevI型模擬低通濾波器的平方幅值響應函數(shù)為：（19）式中，為小于1的正數(shù)，表示通帶內的幅值波紋情況；為截止頻率，N為Chebyshev多項式階數(shù)，為Chebyshev多項式，定義為：（20）ChebyshevI型濾波器特點是：通帶內具有等波紋起伏特性，而在阻帶內則單調下降，且具有更大衰減特性；階數(shù)愈高，特性愈接近矩形。傳遞函數(shù)沒有零點，極點分布在一個橢圓上。注意，這里所說的零點和極點分布在一個橢圓上為Laplace域中的形式。MATLAB信號處理工具箱函數(shù)cheb1ap設計N階ChebyshevI型模擬低通濾波器原型。因為ChebyshevI型模擬原型濾波器通帶內具有波紋特點，所以在調用時需增加通帶波紋所能下降的最低限度Rp（dB）。由于Rp為相對于增益1的下降分貝數(shù)，所以該值越小濾波器在通帶內越接近1，濾波器的性能越好，通常取值范圍為1-5（dB）。此函數(shù)的調用格式為：[z.p,k]=cheb1ap(N,Rp)式中，N為濾波器的階數(shù)，Rp為通帶波紋，單位為dB，z,p,k分別為濾波器的零點、極點和增益。圖12ChebyshevI型模擬原型濾波器平方幅頻圖與Butterworth濾波器相比（圖9），可以看到：在相同的階數(shù)下，ChebyshevI型模擬原型濾波器具有更窄（更陡）的過渡帶。但這種特性是以犧牲了通帶的單調平滑特性（而成為波紋狀）為代價的。如果我們不想犧牲通帶內的單調平滑特性有無辦法設計階數(shù)較小而過渡帶較窄的濾波器呢？這就涉及到ChebyshevII型濾波器。4.2.2ChebyshevII型濾波器ChebyshevII型低通模擬濾波器的平方幅值響應函數(shù)為：（21）式中各項參數(shù)的意義同上。ChebyshevII型模擬濾波器的特點是：阻帶內具有等波紋的起伏特性，而在通帶內是單調、平滑的，階數(shù)愈高，頻率特性曲線愈接近矩形。傳遞函數(shù)既有極點又有零點。ChebyshevII型模擬濾波器零點為[z(1),z(2),…,z(nz)]，極點為[p(1),p(2),…,p(np)]，濾波器的增益為K。MATLAB信號處理工具箱提供函數(shù)cheb2ap設計N階ChebyshevII型模擬濾波器的原型。由于ChebyshevII型濾波器阻帶內有波紋，通帶內單調平滑，則必須給定阻帶衰減Rs，前面已經講過，通帶波紋和阻帶衰減都是相對于增益1的下降，因此，Rp和Rs越大則與通帶增益1的差距越大。因此，阻帶衰減Rs越大對應的濾波器濾波效果越好，通常Rs的取值要大于16dB。該函數(shù)通常調用格式為：[z,p,k]=cheb2ap(N,Rs)式中，N為濾波器的階數(shù)；Rs為阻帶波紋，單位dB;z,p,k為濾波器的零極點和增益。圖13ChebyshevII型模擬原型濾波器平方幅頻圖可見ChebyshevII型濾波器在通帶內是單調平滑的，而阻帶內卻出現(xiàn)了波紋。隨著濾波器階數(shù)的增高，其幅頻特性越接近矩形。由前面的講解可知，Butterworth濾波器具有單調平滑的頻率特性，但在過渡帶寬相同的條件下，所需的濾波器階數(shù)較高（即設計的濾波器濾波速度較慢，成本較高）。而ChebyshevI和II型濾波器在通帶或阻帶內有波紋出現(xiàn)，但在過渡帶寬相同的條件下，所需的濾波器階數(shù)較低（即設計的濾波器濾波速度較快，成本較低）。4.2.3應用舉例為輸入信號，設計一個5階的chebyshevI型帶通濾波器，通帶波紋3dB，下邊界頻率500rad/s，上邊界頻率1000rad/s，繪制幅頻響應圖。給出該濾波器的脈沖響應、階躍響應。模擬濾波器的輸出并給出模擬輸入信號和模擬輸出信號的Fourier振幅譜。程序如下：N=5;Rp=3;w1=500*pi;w2=1000*pi;%濾波器階數(shù)、阻帶衰減[z,p,k]=cheb1ap(N,Rp);%設計ChebyshevI型原型低通濾波器[b,a]=zp2tf(z,p,k);%轉換為傳遞函數(shù)形式Wo=sqrt(w1*w2);%中心頻率Bw=w2-w1;%頻帶寬度[bt,at]=lp2bp(b,a,Wo,Bw);%頻率轉換[h,w]=freqs(bt,at);%計算復數(shù)頻率響應figure(1)subplot(2,2,1),semilogy(w/pi,abs(h));%繪制幅頻響應xlabel('Angularfrequency/pi');gridon;title('Magnitude');subplot(2,2,2),plot(w/pi,angle(h)*180/pi);%繪制相頻響應xlabel('Angularfrequency/pi');ylabel('Phase/^o');gridon;title('Phaseangle')H=[tf(bt,at)];%在MATLAB中表示此濾波器[h1,t1]=impulse(H);%繪出系統(tǒng)的脈沖響應圖subplot(2,2,3),plot(t1,h1);xlabel('Time/s');title('Impluseresponse')[h2,t2]=step(H);%繪出系統(tǒng)的階躍響應圖subplot(2,2,4),plot(t2,h2);xlabel('Time/s');title('Stepresponse')figure(2)dt=1/2000;t=0:dt:0.1;%給出模擬濾波器輸出的時間范圍u=sin(2*pi*30*t)+0.5*cos(2*pi*300*t)+2*sin(2*pi*800*t);subplot(221),plot(t,u)%繪制模擬輸入信號xlabel('Time/s');title('Inputsignal')[ys,ts]=lsim(H,u,t);%模擬系統(tǒng)的輸入u時的輸出subplot(222),plot(ts,ys)%繪制模擬輸出信號xlabel('Time/s'),title('Outputsignal');subplot(223),plot((0:length(u)-1)/length(u)*4000,abs(fft(u)));%繪輸入信號振幅譜title('Spectraofinputsignal'),xlabel('Angularfrequency/pi')Subplot(224)Y=fft(ys);plot((1:length(Y))/length(Y)*4000,abs(Y));%繪制輸出信號振幅譜title('Spectraofoutputsignal')xlabel('Angularfrequency/pi')圖14ChebyshevI型濾波器的幅頻響應、相頻響應、脈沖相應和階躍響應圖15濾波器模擬輸入和輸出信號的時間域和振幅譜程序輸出結果見圖14和圖15。圖14給出了該程序的幅頻圖、相頻圖、脈沖響應和階躍響應。幅頻圖清楚地給出了的通帶范圍和阻帶范圍。圖15給出了模擬輸入、輸出信號的時間域波形及頻率域振幅譜。可以看到通過該濾波器處理后，成功地濾除了通帶范圍之外的信號。但要注意，由該濾波器的相頻特性可知，該濾波器并不是線性相位。4.3橢圓濾波器橢圓(Elliptic)模擬低通原型濾波器的平方幅值響應函數(shù)為：（22）式中，為小于1的正數(shù)，表示波紋情況；為低通濾波器的截止頻率（Cutofffrequency），N為濾波器的階數(shù)，為橢圓函數(shù)，其定義已超出本課程的范圍，我們直接利用。橢圓濾波器的特點：在通帶和阻帶內均具有等波紋起伏特性，與以上濾波器原型相比，相同的性能指標所需的階數(shù)最小，但相頻響應具有明顯的非線性。MATLAB信號處理工具箱提供Elliptic模擬低通濾波器原型設計函數(shù)ellipap。由于橢圓濾波器通帶和阻帶均有波紋，因此函數(shù)的輸入?yún)?shù)中通帶波紋和阻帶衰減均需給出。該函數(shù)調用形式為：式中，n為butterworth濾波器階數(shù)；Rp,Rs分別為通帶波紋和阻帶衰減，單位dB,通常濾波器的通帶波紋的范圍為1-5dB，阻帶衰減的范圍大于15dB。z,p,k分別為濾波器的零點、極點和增益。圖16橢圓模擬原型濾波器平方幅頻圖可見階數(shù)為6的橢圓濾波器的過渡帶已相當窄（陡），但這種特性的獲得是以犧牲通帶和阻帶的單調平滑特性為代價的?？梢钥吹綖V波器的階數(shù)越高平方幅頻響應越接近于矩形。4.3.1應用舉例設計一個高通橢圓濾波器，設計性能指標為：通帶邊界頻率wp=1500Hz,阻帶邊界頻率ws=1000Hz,通帶波紋Rp=1dB,阻帶衰減Rs=100dB。假設一個信號，其中f1=400Hz,f2=1600Hz。信號的采樣頻率為10000Hz。試將原信號與通過該濾波器的模擬信號進行比較。wp=1500*2*pi;ws=1000*2*pi;Rp=1;Rs=100;%濾波器設計參數(shù)[N,Wn]=ellipord(wp,ws,Rp,Rs,'s');%求濾波器的最小階數(shù)和截止頻率w=linspace(1,3000,1000)*2*pi;%給出計算復數(shù)頻率響應的頻率點[b,a]=ellip(N,Rp,Rs,Wn,'high','s');%設計高通橢圓濾波器H=freqs(b,a,w);magH=abs(H);phaH=unwrap(angle(H));figure(1)plot(w/(2*pi),20*log10(magH));%以頻率為橫軸繪幅頻響應xlabel('Frequency/Hz');ylabel('Magnitude/dB');title('Ellipanaloghighpassfilter');holdon;plot([10001000],ylim);%阻帶邊界gridonfigure(2)%給出另一個圖形窗口dt=1/10000;f1=400;f2=1600;%信號中所含頻率成分t=0:dt:0.04;%時間序列x=3*sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t);%輸入信號H=[tf(b,a)];[y,t1]=lsim(H,x,t);%模擬輸出subplot(211),plot(t,x),title('Inputsignal')%繪輸入信號subplot(212),plot(t1,y)%繪制輸出信號title('Outputsignal'),xlabel('Time/s')圖17濾波器的幅頻響應圖圖18濾波器的輸入和輸出信號程序運行結果見圖17和圖18?？梢娝O計的高通濾波器完全符合要求，在阻帶邊界下降的分貝數(shù)達100dB。并且通帶和阻帶均有波紋，這正是橢圓濾波器的特點。將該濾波器模擬輸入400Hz和1600Hz的信號后，輸出結果完全濾除了400HZ低頻成分，只含有1600HZ的信號。4.4Bessel濾波器前面講過的各類原型濾波器均沒有繪出其相位隨頻率的變化特性（相頻特性）。在后面的數(shù)字信號處理學習中將會看到它們的相位特性是非線性的。本節(jié)所介紹的Bessel濾波器就能最大限度地減少相頻特性的非線性，使得通帶內通過的信號形狀不變(拷貝不走樣)。Bessel模擬低通濾波器的特點是在零頻時具有最平坦的群延遲，并在整個通帶內群延遲幾乎不變。在零頻時的群延遲為。由于這一特點，Bessel模擬濾波器通帶內保持信號形狀不變。但數(shù)字Bessel濾波器沒有平坦特性，因此MATLAB信號處理工具箱只有模擬Bessel濾波器設計函數(shù)。函數(shù)besselap用于設計Bessel模擬低通濾波器原型，調用格式為：式中，N為濾波器的階數(shù)，小于25。z,p,k為濾波器的零點、極點和增益。圖19Bessel模擬原型濾波器相頻圖可見，Bessel濾波器具有最優(yōu)的線性相頻特點，但這個特點的獲得是以犧牲窄過渡帶為代價的,即濾波器的幅頻平方特性與矩形特性相差甚遠。4.4.1應用舉例設計一個六階Bessel低通模擬濾波器，截止頻率為2000rad/s,繪制濾波器的頻率特性圖。假設一個信號，其中f1=400Hz,f2=1000Hz。信號的采樣頻率為10000Hz。試將原信號與通過該濾波器的模擬信號進行比較。N=6;%濾波器階數(shù)Wn=2000;%邊界頻率[b,a]=besself(N,Wn);%設計Bessel濾波器figure(1)freqs(b,a);%繪出濾波器的復數(shù)頻率特性figure(2)dt=1/10000;%信號采樣間隔f1=100;f2=1000;%輸入兩個頻率成分t=0:dt:0.1;x=7*sin(2*pi*f1*t)+6*cos(2*pi*f2*t);%輸入信號H=[tf(b,a)];%濾波器在NATLAB系統(tǒng)中的表示[y,t1]=lsim(H,x,t);%模擬濾波器的輸出subplot(211),plot(t,x),title('Inputsignal')%繪出輸入信號subplot(212),plot(t1,y)%繪出輸出信號title('Outputsignal'),xlabel('Time/s')圖20Bessel濾波器的頻率特性圖21Bessel濾波器的輸入和輸出程序運行結果為圖20和圖21。雖然模擬Bessel濾波器的幅頻特性不如其他濾波器好，但其相頻特性為線性的。因此圖21中的輸出信號與輸入的低頻信號形狀不變，只是有延遲。前面介紹了那么多個濾波器原型，我們可以得出以下結論：Butterworth濾波器在通帶和阻帶內均具有平滑單調的特點，但在相同過渡帶寬的條件下，該濾波器所需的階數(shù)最多。ChebyshevI和II型濾波器在通帶或阻帶內具有波紋，但在相同過渡帶寬的條件下，該濾波器所需的階數(shù)比Butterworth濾波器要少。橢圓濾波器在通帶和阻帶內均有波紋出現(xiàn)，但在相同過渡帶寬的條件下，該濾波器所需的階數(shù)最少。Bessel濾波器具有最寬的過渡帶，但具有最優(yōu)的線性相頻特性。因此沒有絕對“好”的濾波器，要根據(jù)解決問題的不同選擇不同的濾波器，因此，每一種濾波器的設計方法我們都要熟練掌握。5模擬有源濾波器設計舉例5.1有源一階濾波電路5.1.1一階濾波電路如果在一階RC低通電路的輸出端，再加上一個電壓跟隨器，使之與負載很好隔離開來，就構成一個簡單的一階有源RC低通濾波電路，由于電壓跟隨器的輸入阻抗很高，輸出阻抗很低，因此，其負載能力很強。

如果希望電路不僅有濾波功能，而且能起放大作用，則只要將電路中的電壓跟隨器改為同相比例放大電路即可。如圖22。圖22一階低通濾波電路傳遞函數(shù)

RC低通電路的傳遞函數(shù)為(23)

對于電壓跟隨器，其通帶電壓增益等于同相比例放大電路的電壓增益，即

(24)因此，可導出電路的傳遞函數(shù)為(25)式中=1/(RC)，稱為特征角頻率。

由于傳遞函數(shù)中分母為s的一次冪，故上述濾波電路稱為一階低通有源濾波路。（2）幅頻響應

對于實際的頻率來說，式25中的s可用代入，由此可得(26)

(27)圖23幅頻響應5.1.2應用實例若圖22中則由上述公式可得傳遞函數(shù)為,特征角頻率為。對傳遞函數(shù)用MATLAB仿真，求系統(tǒng)的幅頻響應和相頻響應其程序為：b=[10];%濾波器傳遞函數(shù)分子多項式系數(shù)a=[18];%濾波器傳遞函數(shù)分母多項式系數(shù)figure(1),freqs(b,a)%第一種輸出方法[h,w]=freqs(b,a);%計算濾波器的復數(shù)頻率響應mag=abs(h);pha=angle(h);%得到濾波器的幅頻和相頻響應figure(2),subplot(2,1,1),loglog(w,mag);%運用雙對數(shù)坐標繪制幅頻響應gridon;xlabel('Angularfrequency');ylabel('Magnitude');subplot(2,1,2),semilogx(w,pha*180/pi)%運用半對數(shù)坐標繪相頻響應gridon;xlabel('Angularfrequency');ylabel('Phase/degrees');圖24兩種方法對應的相頻響應和幅頻響應5.2雙二階濾波電路5.2.1雙二階電路所謂“雙二階電路”指的是能夠實現(xiàn)雙二階函數(shù)（即傳輸函數(shù)的分子、分母都是二階函數(shù)）的電路。雙二階電路可以實現(xiàn)各種類型的二階函數(shù)，例如二階低通、高通、帶通、具有傳輸零點以及全通等濾波器。5.2.2應用實例我們知道雙二階電路能構成多種濾波器，現(xiàn)有如下多運放二階電路，下面我們將通過以下電路進行討論仿真。圖25ken多運放二階電路我們分析其主電路是由一個相加器A1，和兩個積分器A2、A3組成的，由上圖可以看出，第一個放大器A1的輸出:（28）第二個放大器的輸出和第三個放大器的輸出都是反相積分器，它們電壓的關系為:由以上三式可求得（29）（30）（31）上述的H（s）都是一般的雙二階函數(shù)，我們將函數(shù)分母用極點頻率和等效品質因數(shù)值來表示，即將改寫為標準形式：（32）經變形即可得出每一個傳遞函數(shù)的等效品質因素Q和截止頻率：中：中： 中： 不同的參數(shù)將導致不同的頻率特性，為了簡單起見，不妨設圖25的電阻為，電容為。則有傳遞函數(shù)為：；=1；；；；；對傳遞函數(shù)用MATLAB仿真，其程序為：b=[001];a=[111];%濾波器傳遞函數(shù)分母多項式系數(shù)figure(1),freqs(b,a)；%第一種方法[h,w]=freqs(b,a);%計算濾波器的復數(shù)頻率響應mag=abs(h);pha=angle(h);%得到濾波器的幅頻和相頻響應figure(2),subplot(2,1,1),loglog(w,mag);%雙對數(shù)坐標繪制幅頻響應gridon;xlabel('Angularfrequency');ylabel('Magnitude');subplot(2,1,2),semilogx(w,pha*180/pi)%運用半對數(shù)坐標繪相頻響應gridon;xlabel('Angularfrequency');ylabel('Phase/degrees');圖形如圖26所示:圖26V3端輸出低通濾波器的相頻特性修改程序為：b=[0-10];%濾波器傳遞函數(shù)分子多項式系數(shù)其圖形如圖27圖27V2端輸出帶通濾波器的相頻特性修改程序為：b=[100];%濾波器傳遞函數(shù)分子多項式系數(shù)其圖形如圖28圖28V1端輸出高通濾波器的相頻特性由圖26-28我們看出，由不同點取輸出，可得到不同類型的傳遞函數(shù)，由第一級取輸出得高通函數(shù)；第二級取輸出得帶通函數(shù)；由第三級輸出則得低通函數(shù)。經上可知低通與高通串聯(lián)可的得到帶通濾波器，與第一章所講內容相同。5.3二階有源帶阻濾波電路帶通濾波電路是用來使某一頻段的信號通過，而讓該頻段以外的信號全部抑制或衰減。與帶通濾波電路相反，帶阻濾波電路是用來抑制或衰減某一頻段的信號，而讓該頻段以外的所有信號通過。在第一章中，我們提及過高通濾波器和低通濾波器并聯(lián)就可以組成帶阻濾波器。我們這里要討論的是另一種方案，即雙T帶阻濾波器。電路如圖29所示，其中。由節(jié)點導納方程不難導出電路的傳遞函數(shù)為