Předmět Analogové elektronické obvody (AEO) je zařazen na FEKT VUT v Brně, v bakalářském studijním programu Elektrotechnika, elektronika, komunikační a řídicí technika,v zimním semestru 2. ročníku, jako povinný v oboru Elektronika a sdělovací technika. Do kurzu jsou zařazena i počítačová cvičení (1h. × 13t. ≡ 1 × 2h. každé 2t.), výuka tedy probíhá všesti čtrnáctidenních cyklech.
Vydal: FEKT VUT Brno
Autor: UREL - Lubomír Brančík, Tomáš Dostál, Zdeněk Kolka
Strana 13 z 59
Vámi hledaný text obsahuje tato stránku dokumentu který není autorem určen k veřejnému šíření.
Razítko
zdroje
Razítko
OZ
.
• Pro odpor modulová fázová charakteristika vykreslena Obr.m
function H=mat(s,R1)
R2=1000;
R=1000;
C=1e-6;
H= [(1/R1+1/R) -1/R1 -1/R 0
0 1/R2 -1/R2 1
-1/R1 -1/R2 (1/R1+1/R2) 0
-1/R (1/R+s*C) 0
1 0
0 0];
• Pro snazší zobrazení frekvenčních charakteristik pro různé hodnoty odporu bude i
vlastní řešení obvodu formulováno formě (hlavní) funkce: apf.m. 2. Tímto způsobem obvod
chová pouze při rovnosti odporů R2, jak lze snadno přesvědčit voláním funkce apf při
jiné hodnotě argumentu.
Vidíme, zisk (modul přenosu) konstantní obvod pouze natáčí fázi (0o
– 180o
) jedná
se tedy všepropustný fázovací dvojbran (all pass filter APF). apf(500); zde
tedy pro odpor 500 Ω.Analogove_elektronicke_obvody_P- Počítačová cvičení 13
• Matici vytvoříme jako funkci Matlabu, samostatném souboru: mat.
function f=apf(R1)
n=200;
f1=1;
f2=10e3;
% děleni kmitočtového intervalu
F=logspace(log10(f1),log10(f2),n);
% budicí vektor (pouze u1, MMUN u1, i1!)
b=[0;0;0;0;1;0]; nezapomeňte, definujeme sloupcový vektor!
%opakované řešeni soustavy rovnic cyklus:
for i=1:n,
H=mat(2*pi*j*F(i),R1); voláni matice pro s=jω hodnotu R1
x=H\b; vektor neznámých
% x=[u1 izdroj ioz]
K(i)=x(2); zajímá nás u2
% protože u1=1, tak přímo přenos napětí
end;
% kreslení modulové fázové charakteristiky
subplot(2,1,1)
semilogx(F,20*log10(abs(K)),'m','LineWidth',2);
xlabel('f');
ylabel('|K| [dB]');
grid on;
if R1==1000, axis([1,1e4,-1,1]), end;
subplot(2,1,2);
semilogx(F,angle(K)*180/pi,'b','LineWidth',2);
xlabel('f');
ylabel('phase [deg]');
grid on;
%*****************************************************************************************
• příkazovém řádku MATLABu tato funkce volána příkazem např.7
