Dolní propust

Zadání aproximační úlohy - definice mezních kmitočtů.

>    fp1:=1000:fs1:=2000:ap1:=3:as1:=50:

Transformace tolerančního schématu daného filtrtu dolní propusti  (LP) na nornmovanou dolní propust (NLP).

>    infolevel[syntfil]:=3:   #zajistí vykreslení tolerančního schématu

>    x1:=LP2NLP(fp1,fs1,ap1,as1);

>    infolevel[syntfil]:=1:

Filter specification:

fp = 1000.000000 Hz

fs = 2000.000000 Hz

ap = 3.000000 dB

as = 50.000000 dB

[Maple Plot]

Specification of NLP:

Os = 2.000000 1/s

ap = 3.000000 dB

as = 50.000000 dB

x1 := 2., 3, 50

Výpočet stupňů a nových útlumů v nepropustném pásmu pro všechny typy standardních aproximací.

>    Nb1:=ButterworthNLPOrder(x1);

>    Butterworth_asnew(Nb1);

>    Nch1:=ChebyshevNLPOrder(x1);

>    Chebyshev_asnew(Nch1);

>    Nc1:=CauerNLPOrder(x1);

>    Cauer_asnew(Nc1);

Nb1 := 9, 2., 3

54.164791466084385603

Nch1 := 5, 2., 3

51.153580310131411741

Nc1 := 4, 2., 3

57.753956206760133983

Výpočet pólů a případně nul (pro eliptické filtry) přenosu NLP pro vybranou aproximaci. Toto je potřebné pro následnou kaskádní syntézu ARC (aktivní RC) struktury filtru.

>    Gc1,poles1:=ButterworthPoles(Nb1);

Gc1, poles1 := .99762834511098350277, vector([-.93994057277650027102+.34211039046977123336*I, -.93994057277650027102-.34211039046977123336*I, -.76624657543368573891+.64295721883725957044*I, -.766246575...
Gc1, poles1 := .99762834511098350277, vector([-.93994057277650027102+.34211039046977123336*I, -.93994057277650027102-.34211039046977123336*I, -.76624657543368573891+.64295721883725957044*I, -.766246575...
Gc1, poles1 := .99762834511098350277, vector([-.93994057277650027102+.34211039046977123336*I, -.93994057277650027102-.34211039046977123336*I, -.76624657543368573891+.64295721883725957044*I, -.766246575...

Výpočet provozního činitele přenosu a charakteristické funkce pro vybranou aproximaci. Toto je potřebné pro následnou syntézu příčkové LC struktury filtru.

>    g1,chf1:=Butterworth(Nb1,p);

g1, chf1 := .99999999999999999968+.99762834511098350277*p^9+5.7466286011093642717*p^8+16.551123692969231713*p^7+31.114145368128408135*p^6+41.931035867794812144*p^5+41.942100000449492500*p^4+31.13878166...
g1, chf1 := .99999999999999999968+.99762834511098350277*p^9+5.7466286011093642717*p^8+16.551123692969231713*p^7+31.114145368128408135*p^6+41.931035867794812144*p^5+41.942100000449492500*p^4+31.13878166...

>   

Výpočet některých ostatních aproximací

Výpočet modulu přenosové funkce dané NLP z provozního činitele přenosu a jeho vykreslení v závislosti na kmitočtu.

>    plot(MagnitudeH(1/g1)(omega),omega=0..5);

Výpočet modulu přenosu pro mez propustného pásma (měl by samozřejmě vyjít jako zadané ap1  - kontrola).

>    evalf(MagnitudeHdB(1/g1)(1));

Výpočet modulu přenosu pro mez nepropustného pásma (měl by samozřejmě vyjít jako zadané asnew  - kontrola).

>    evalf(MagnitudeHdB(1/g1)(x1[1]));

[Maple Plot]

-2.9999999999999999139

-54.164791466084385596

>   

Vykreslení modulu v dB, fáze a skupinového zpoždění NLP

>