ΠΠ°Π±ΠΎΡΠ°ΡΠΎΡΠ½Π°Ρ ΡΠ°Π±ΠΎΡΠ° № 1. ΠΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠ΅ ΠΌΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΠΎΠ±ΡΠ°Π±ΠΎΡΠΊΠΈ ΡΠΈΠ³Π½Π°Π»ΠΎΠ² Π»ΠΈΠ½Π΅ΠΉΠ½ΠΎΠΉ Π΄ΠΈΡΠΊΡΠ΅ΡΠ½ΠΎΠΉ ΡΠΈΡΡΠ΅ΠΌΠΎΠΉ Π² ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠΉ ΡΡΠ΅Π΄Π΅ MATLAB
1.1 Π’Π΅ΠΎΡΠ΅ΡΠΈΡΠ΅ΡΠΊΠΈΠ΅ ΡΠ²Π΅Π΄Π΅Π½ΡΡ
Π MATLAB ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠΉ ΠΌΠΎΠ΄Π΅Π»ΡΡ ΠΠΠ‘ Π½Π°Π·ΡΠ²Π°ΡΡ ΡΠΎΠΎΡΠ½ΠΎΡΠ΅Π½ΠΈΠ΅ Π²Ρ
ΠΎΠ΄/Π²ΡΡ
ΠΎΠ΄ Π² Π²ΠΈΠ΄Π΅ ΡΡΠ°Π²Π½Π΅Π½ΠΈΡ ΠΈΠ»ΠΈ ΡΠΈΡΡΠ΅ΠΌΡ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠΉ, ΠΊΠΎΡΠΎΡΡΠ΅ ΠΏΠΎΠ·Π²ΠΎΠ»ΡΡΡ Π²ΡΡΠΈΡΠ»ΠΈΡΡ ΡΠ΅Π°ΠΊΡΠΈΡ Π½Π° Π·Π°Π΄Π°Π½Π½ΠΎΠ΅ Π²ΠΎΠ·Π΄Π΅ΠΉΡΡΠ²ΠΈΠ΅.
ΠΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΡΠ°Π±ΠΎΡΡ ΠΠΠ‘ Π½Π° ΠΎΡΠ½ΠΎΠ²Π΅ ΡΠ°Π·Π½ΠΎΡΡΠ½ΠΎΠ³ΠΎ ΡΡΠ°Π²Π½Π΅Π½ΠΈΡ — Π²ΡΡΠΈΡΠ»Π΅Π½ΠΈΠ΅ ΡΠ΅Π°ΠΊΡΠΈΠΈ Π½Π° Π²Ρ
ΠΎΠ΄Π½ΠΎΠ΅ Π²ΠΎΠ·Π΄Π΅ΠΉΡΡΠ²ΠΈΠ΅ ΠΏΡΠΈ Π½ΡΠ»Π΅Π²ΡΡ
Π½Π°ΡΠ°Π»ΡΠ½ΡΡ
ΡΡΠ»ΠΎΠ²ΠΈΡΡ
— Π²ΡΠΏΠΎΠ»Π½ΡΠ΅ΡΡΡ Ρ ΠΏΠΎΠΌΠΎΡΡΡ ΡΡΠ½ΠΊΡΠΈΠΈ filter, ΡΠΎΡΠΌΠ°Ρ ΠΊΠΎΡΠΎΡΠΎΠΉ ΠΈΠΌΠ΅Π΅Ρ Π²ΠΈΠ΄: filter (b, a, x) Π³Π΄Π΅:
b — Π²Π΅ΠΊΡΠΎΡ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠΎΠ² Π² ΠΏΠΎΡΡΠ΄ΠΊΠ΅ ΠΈΡ
ΡΠ»Π΅Π΄ΠΎΠ²Π°Π½ΠΈΡ;
Π° — Π²Π΅ΠΊΡΠΎΡ ΠΊΠΎΡΡΡΠΈΡΠΈΠ΅Π½ΡΠΎΠ² Π² ΠΏΠΎΡΡΠ΄ΠΊΠ΅ ΠΈΡ
ΡΠ»Π΅Π΄ΠΎΠ²Π°Π½ΠΈΡ (ΠΏΠ΅ΡΠ²ΡΠΉ ΡΠ»Π΅ΠΌΠ΅Π½Ρ Π²ΡΠ΅Π³Π΄Π° ΡΠ°Π²Π΅Π½ 1);
Ρ
— Π²Π΅ΠΊΡΠΎΡ ΠΎΡΡΡΠ΅ΡΠΎΠ² Π²ΠΎΠ·Π΄Π΅ΠΉΡΡΠ²ΠΈΡ Ρ
(n).
1.2 ΠΠΎΡΡΠ΄ΠΎΠΊ Π²ΡΠΏΠΎΠ»Π½Π΅Π½ΠΈΡ ΡΠ°Π±ΠΎΡΡ
Π Π°ΡΡΡΠΈΡΠ°Π΅ΠΌ ΡΠ΅Π°ΠΊΡΠΈΡ ΠΠΠ₯-ΡΠΈΠ»ΡΡΡΠ° 2-Π³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ°, Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ ΡΠ°Π·Π½ΠΎΡΡΠ½ΡΠΌ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠ΅ΠΌ
y (n) = 0.1 x (n) + 0.7 x (n — 1) + 0.5 x (n — 2),
Π³Π΄Π΅: n = 32…64; wT = 0.7 ΡΠ°Π΄; x (n) = sin (wTn).
Π MATLAB Π² ΡΠ΅ΠΆΠΈΠΌΠ΅ ΠΏΡΡΠΌΡΡ
Π²ΡΡΠΈΡΠ»Π΅Π½ΠΈΠΉ ΡΡΠΎ Π²ΡΠ³Π»ΡΠ΄ΠΈΡ ΡΠ»Π΅Π΄ΡΡΡΠΈΠΌ ΠΎΠ±ΡΠ°Π·ΠΎΠΌ:
>> b= [0.1 0.7 0.5];
>> a= [1];
>> n=32: 64;
>> x=sin (0.7. *n);
>> y=filter (b, a, x);
>> plot (n, x, n, y,'—'), grid
>> hold on
>> stem (n, x)
>> stem (n, y)
>> gtext ('Vuhodnoy signal')
>> gtext ('Vhodnoy signal')
Π Π΅Π·ΡΠ»ΡΡΠ°ΡΡ ΡΠ°ΡΡΠ΅ΡΠ° ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Ρ Π½Π° ΡΠΈΡ. 1.1, Π³Π΄Π΅ ΠΊΡΠΎΠΌΠ΅ Π΄ΠΈΡΠΊΡΠ΅ΡΠ½ΡΡ
ΡΠΈΠ³Π½Π°Π»ΠΎΠ² ΠΈΠ·ΠΎΠ±ΡΠ°ΠΆΠ΅Π½Ρ ΠΈΡ
ΠΎΠ³ΠΈΠ±Π°ΡΡΠΈΠ΅.
Π ΠΈΡ. 1.1 ΠΡ
ΠΎΠ΄Π½ΠΎΠΉ ΠΈ Π²ΡΡ
ΠΎΠ΄Π½ΠΎΠΉ ΡΠΈΠ³Π½Π°Π»Ρ
Π Π°ΡΡΡΠΈΡΠ°Π΅ΠΌ ΡΠ΅Π°ΠΊΡΠΈΡ ΠΠΠ₯-ΡΠΈΠ»ΡΡΡΠ° 2-Π³ΠΎ ΠΏΠΎΡΡΠ΄ΠΊΠ°, Π·Π°Π΄Π°Π½Π½ΠΎΠ³ΠΎ ΡΠ°Π·Π½ΠΎΡΡΠ½ΡΠΌ ΡΡΠ°Π²Π½Π΅Π½ΠΈΠ΅ΠΌ
y (n) = 1.5 x (n) — 1,5 x (n — 1) + 0.5 x (n — 2) + 1.34 y (n — 1) — 0.9 y (n — 2),
Π³Π΄Π΅: n = 32…64; wT = 0.7 ΡΠ°Π΄; x (n) = sin (wTn).
Π MATLAB Π² ΡΠ΅ΠΆΠΈΠΌΠ΅ ΠΏΡΡΠΌΡΡ
Π²ΡΡΠΈΡΠ»Π΅Π½ΠΈΠΉ ΡΡΠΎ Π²ΡΠ³Π»ΡΠ΄ΠΈΡ ΡΠ»Π΅Π΄ΡΡΡΠΈΠΌ ΠΎΠ±ΡΠ°Π·ΠΎΠΌ:
>> b= [1.5 — 1.5 0];
>> a= [1 1.34 — 0.9];
>> n=32: 64;
>> x=sin (0.7. *n);
>> y=filter (b, a, x);
>> plot (n, x, n, y,'-. '), grid
>> hold on
>> stem (n, x)
>> stem (n, y)
>> gtext ('Vuhodnoy signal')
>> gtext ('Vhodnoy signal')
Π Π΅Π·ΡΠ»ΡΡΠ°ΡΡ ΡΠ°ΡΡΠ΅ΡΠ° ΠΏΡΠ΅Π΄ΡΡΠ°Π²Π»Π΅Π½Ρ Π½Π° ΡΠΈΡ. 1.2, Π³Π΄Π΅ ΠΊΡΠΎΠΌΠ΅ Π΄ΠΈΡΠΊΡΠ΅ΡΠ½ΡΡ
ΡΠΈΠ³Π½Π°Π»ΠΎΠ² ΠΈΠ·ΠΎΠ±ΡΠ°ΠΆΠ΅Π½Ρ ΠΈΡ
ΠΎΠ³ΠΈΠ±Π°ΡΡΠΈΠ΅.
Π»ΠΈΠ½Π΅ΠΉΠ½Π°Ρ Π΄ΠΈΡΠΊΡΠ΅ΡΠ½Π°Ρ ΡΠΈΡΡΠ΅ΠΌΠ° ΡΠΈΠ³Π½Π°Π» Π ΠΈΡ. 1.2 ΠΡ
ΠΎΠ΄Π½ΠΎΠΉ ΠΈ Π²ΡΡ
ΠΎΠ΄Π½ΠΎΠΉ ΡΠΈΠ³Π½Π°Π»Ρ
ΠΡΠ²ΠΎΠ΄: Π Π»Π°Π±ΠΎΡΠ°ΡΠΎΡΠ½ΠΎΠΉ ΡΠ°Π±ΠΎΡΠ΅ Π±ΡΠ»ΠΎ Π²ΡΠΏΠΎΠ»Π½Π΅Π½ΠΎ ΠΌΠ°ΡΠ΅ΠΌΠ°ΡΠΈΡΠ΅ΡΠΊΠΎΠ΅ ΠΌΠΎΠ΄Π΅Π»ΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΠΎΠ±ΡΠ°Π±ΠΎΡΠΊΠΈ ΡΠΈΠ³Π½Π°Π»ΠΎΠ² Π»ΠΈΠ½Π΅ΠΉΠ½ΠΎΠΉ Π΄ΠΈΡΠΊΡΠ΅ΡΠ½ΠΎΠΉ ΡΠΈΡΡΠ΅ΠΌΠΎΠΉ Π² ΠΏΡΠΎΠ³ΡΠ°ΠΌΠΌΠ½ΠΎΠΉ ΡΡΠ΅Π΄Π΅ MATLAB.