Cvičení 2_3 SNAP pro studenty

dou

abs(equalx(j*2*pi*1000))

ans =

9.583471831005770e-001

angle(equalx(j*2*pi*1000))*180/pi

ans =

2.180574342864527e+001
(je uveden přepočet úhlu na stupně z radiánů).

Koeficienty čitatele a jmenovatele přenosové funkce (tj. výsledky semisymbolické analýzy) získáme takto:

num=equalx('numer')

num =

4.804575000000000e+006 6.035950000000000e+010 9.350000000000000e+013

den=equalx('denom')

den =

4.804575000000000e+006 6.440950000000000e+010 2.435000000000000e+014

Který je který? Viz předchozí hlášku

polynomials are in the Matlab style, ie. C(1)*s^N + ... + C(N)*s + C(N+1)

Kmitočtovou charakteristiku získáme nejsnadněji pomocí klasické matlabovské funkce:

freqs(num,den)



Abychom mohli pracovat s vnitřními proměnnými funkce equalx, tj. s proměnnými R1, R2, …, krokovat je
atd., musíme je nejprve exportovat do prostoru globálních proměnných MATLABu, abychom na ně „viděli“
z příkazové řádky:

equalx(‘export’)

Budeme pracovat s proměnnými R1, R2, R3 a R4 (ovládání potenciometrů). Tyto proměnné zvolíme za
globální:

global R1 R2 R3 R4
V pracovním adresáři MATLABu založíme soubor KRESLENI.M:

flog=(1:0.1:5); tvorba log. kmitočtové osy
f=10.^flog; vektor kmitočtů od 10Hz do 100kHz, odstupňovaných logaritmicky
for R1=0:10000:50000, krokování R1od 0 do 50kΩ po 10kΩ
R2=50000-R1; dopočty R2, spolu s R1 tvoří potenciometr
K=equalx(j*2*pi.*f); výpočet kmitočtové char. ekvalizéru v daném kroku
Kdb=20*log10(abs(K)); přepočet na decibely
semilogx(f,Kdb); vykreslení kmitočtové charakteristiky pro daný R1 s log. stupnicí f
hold on; obrázky budou drženy na obrazovce, nebudou se přepisovat
end;

Soubor uložíme a v příkazové řádce MATLABu napíšeme jeho jméno
kresleni
a potvrdíme. Otevře se okno výsledků:

Pokud chceme zviditelnit pomocnou mřížku, zapíšeme
Grid

Vytváření třírozměrného grafu: Vytvoříme soubor KRESLENIx.M:

logf=(1:0.1:5); viz výše; celkem 41 členů
f=10.^logf; viz výše; celkem 41 členů
z=zeros(50,41); matice o 50 řádcích a 41 sloupcích vyplněná nulami
y=zeros(1,50); matice o 1 řádku a 50 sloupcích vyplněná nulami
for i=1:50, proběhne 50 cyklů
R1=1000*i; krokování a výpočet odporů = ovládání potenciometrů
R2=50000-R1;
R3=25000;
R4=50000-R3;
m=abs(equalx(j*2*pi.*f)); výpočet amplitudové kmitočtové charakteristiky
z(i,1:41)=20.*log10(m); zapsání 41 hodnot decibelové kmit. char. do řádku i matice z
y(1,i)=R1; zapsání odporu R1 do sloupce i vektoru y
end
surfc(logf,y,z) vykreslení funkce z – funkce proměnných logf a y; osa x je logaritmická
Před jeho spuštěním vymažeme předchozí obrázek. Výsledkem provedených instrukcí v souboru
KRESLENIx je třírozměrný graf:



N ose “x” je vynesen logaritmus kmitočtu (1 až 5, tedy kmitočty od 10 do 100000Hz), na ose “y” krokovaný
odpor R1 (od 0 do 50kohmů), na ose “z” přenos v decibelech.
Nyní je vhodné uspořádat okna MATLABu, obrázku a editační okno souboru KRESLENIx tak, aby se
nepřekrývala. V souboru KRESLENIx vždy změníme odpor R3, klikneme na ikonu diskety pro uložení,
přejdeme do okna MATLABu, zmáčkneme kurzorovou klávesu ↑, tím si vyvoláme poslední použitý příkaz
(kreslenix), potvrdíme a sledujeme změnu kmitočtové charakteristiky.