Kmitočtové a přenosové charakteristiky
> | restart; |
Definice přenosové funkce
> | H:=1/(s^2+sqrt(2)*s+1); |
Dosazení za Laplaceův operátor (toto značení je "anglické", u nás obvykle používáme označení ) komplexní kmitočet .
> | H_f:=subs(s=I*omega,H); |
Výpočet modulu.
> | AH:=evalc(abs(H_f)); |
Definice fáze jako funkce - pro procvičení.
> | PH:=unapply(evalc(argument(H_f)),omega); |
Modulová frekvenční charakteristika.
> | plot(AH,omega=0..3); |
Fázová frekvenční charakteristika ve stupních.
> | plot(PH(omega)*180/Pi,omega=0..3); |
Počítání s logaritmy
Přirozený logaritmus
> | ln(10.0); |
Obecný logaritmus - standartně přirozený
> | log(10.0); |
> | log[exp(1)](10.0); |
Volitelně desítkový
> | log[10](10.0); |
Desítkový logaritmus
> | log10(10.0); |
Modulová charakteristika v dB
> | plot(20*log10(AH),omega=0..3); |
Modulová charakteristika s logaritmickým měřítkem na ose kmitočtu.
> | plots[semilogplot](20*log10(AH),omega=0.1..10); |
Fázová chrakteristika. Logaritmické měřítko na ose x (kmitočtu) je uděláno jiným způsoem.
> | plot(PH(10^x),x=-1..1); |
Použití funkce semilogplot dává špatné výsledky - nerovnoměrně rozdělené body, což je patrné zejména pro velký rozsah hodnot kmitočtu.
> | plots[semilogplot](PH(omega),omega=0.01..100,style=point); |
S několika dalšími parametry můžeme výsledek opravit.
> | plots[semilogplot](PH(omega),omega=0.01..100,style=point,sample=[seq(10^(i/10),i=-20..20)],adaptive=false); |
Zpět na začátek sekce .