Zjištění frekvenční charakteristiky pomocí „měření“ v MATLABu a SIMULINKu

1) nadefinujte v m-file proměnné
- sys pro soustavu - ze zadání
- frekv pro frekvenci - například 1
- perioda - vypočtěte z frekvence
- Tmax - z důvodu přechodného děje je nutné cca 20 period
- krok - cca 50 vzorků na periodu (u nízkých frekvencí má matlab problém s iteracemi)

2) model v simulinku
- generátor sinusu - frekvenci nastavit na frekv (v rad/s => je to w)
- soustava (LTI system) - s přenosem sys
- výstup lze na osciloskom, nebo do workspace - vstup i výstup - nastavte na array
- nastavit parametry simulace na fixed step s krokem krok a maximální čas nastavit na Tmax

3) v m-file spustit simulaci a vytisknout do grafu vstupní a výstupní signál

4) zjištění polohy a velikosti maxima výstupního signálu
- z výstupního signálu vezměte část jeho konce (ustálený stav) - cca 2 periody aby určitě bylo maximum
- příkazem max zjistěte hodnotu maximu a jeho polohu

5) zjištění polohy maxima vstupního signálu
- hledáme nejbližší maximum (k počátku) ve vstupním signálu (například tak, že sjedem do minima a potom vyjedeme do maxima)
- zjistíme velikost a polohu maxima

6) výpočet parametrů frekvenční charakteristiky
- z hodnot velikostí maxim zjistíme zesílení-amplitudu
- z hodnot poloh maxim zjistíme úhel-fázi posunu (posun o jednu periodu je 360 stupňů)
- z amplitudy a fáze vypočteme reálnou a imaginární složku

7) propočet celé charakteristiky
- zvolíme frekvence, při kterých budeme počítat (pomocí příkazu logspace)
- vytvoříme výstupní proměnné, do kterých budeme ukládat výstupní hodnoty
- přidáme cyklus tak, aby se simulace a výpočty prováděly pro všechny frekvence (je výhodné mít celočíselnou proměnnou index - vybere frekvenci a uloží se výstupní hodnoty na daný index)

8. zobrazíme do jednoho grafu výsledek pomocí nyquist získané hodnoty (re, im). Srovnáme výsledek

Poslední změna: 2009-03-04