МЕТА РОБОТИ
Отримати практичні навики використання програми спектрального аналізу, дослідити спектри дискретних сигналів різної форми та визначити їх особливості.
КОРОТКІ ТЕОРЕТИЧНІ ВІДОМОСІ
Визначення спектральних складових дискретних (дискретизованих) сигналів.
Обробка та дослiдження сигналiв з використанням персональних ЕОМ вимагає їх дискретного цифрового представлення. При цьому сигнали описуються сукупнiстю N вiдлiкiв (xk, k=0,N-1) на заданому iнтервалi часу (0,T). Ця сукупнiсть вiдлiкiв може описувати дискретний сигнал Xд(t), або представляти миттєвi значення неперервного сигналу X(t) у певнi моменти часу. В останньому випадку розглядається дискретизована неперервна функцiя, яка при виконаннi певних умов буде адекватно представляти неперервну функцію з необхідною точнiстю (питання дискретизацii неперервних функцiй розглядаються в iншiй лабораторній роботi).
Якщо задану сукупнiсть виборок подумки повторити безмежну кiлькiсть разiв, то дослiджуваний сигнал можна вважати перiодичним. Для визначення спектру можна ввести певну математичну модель дискретного перiодичного сигналу i використати розклад у ряд Фур'є. Якщо сигнал неперервний, то за допомогою послiдовностi дельта-iмпульсiв можна отримати його дискретне представлення на iнтервалi (0,T).
(1)
де: xk = X(k*d) - вiдлiки у k-й точцi; d – інтервал дискретизації; N=T/d.
Дискретну модель можна представити комплексним рядом Фур'є:
(2)
з коефiцiєнтами
(3)
Пiдставивши (1) в (3) пiсля нескладних математичних перетворень отримаємо
(4)
або у тригонометричнiй формi
(5)
(6)
(7)
(8)
Необхiдно зауважити, що при обчисленнi кута з використанням арктангенса потрiбно враховувати знаки Cns та Сnс для правильного визначення квадранта.
Вказанi формули визначають послiдовнiсть коефiцiєнтiв спектральних складових заданого вiдлiками сигналу i описують дискретне перетворення Фур'є (ДПФ).
Основнi властивостi ДПФ:
ДПФ є лінійним перетворенням, тобто ДПФ суми сигналiв є сума коефiцiентiв ДПФ кожного з них, а змiна амплiтуд сигналу в М-разiв викликає таку ж змiну вiдповiдних коефiцiєнтiв С(n).
Кiлькiсть рiзних коефiцiєнтiв С(0),...,С(N-1) визначається кiлькiстю вiдлікiв N (якщо n=N, то С(n)=C(0), тобто сигнали i спектри перiодично повторюються).
Коефiцiєнт С(0) (нульова гармонiка, яка визначає постiйну складову є середнiм значенням всiх вiдлiкiв.
(9)
Якщо кiлькicть вiдлiкiв N - парне число, то
(10)
Якщо значення вiдлiкiв xk- дiйснi числа, то коефiцiенти ДПФ, номери яких симетричнi вiдносно N/2 утворюють комплекснi спряженi пари
(11)
Тому можна вважати, що коефiцiенти С(N/2+1),...C(N-1) вiдповiдають вiд'ємним частотам.
Вiдновлення початкового сигналу по коефіцієнтах ДПФ.
Якщо на основi заданих вiдлiкiв знайденi коефiцiєнти ДПФ (С(0),...,С(N/2)), то по цих коефiцiєнтах завжди можна вiдновити початковий сигнал Хд(t), або дискретизований сигнал X(t). Для такого сигналу ряд Фур'є записується скiнченою сумою
(12)
де: ¦Сi¦ - модуль амплiтуди вiдповiдної гармонiки, а (i - її фаза.
Зворотнє перетворення Фур'є.
Нехай коефiцiенти Сn, що утворюють ДПФ, заданi. Якщо у формулi (2) t = k*d i сумується скiнченна кiлькiсть членiв ряду, якi вiдповiдають iснуючим гармонiкам у спектрi сигналу, то отримуємо таку формулу для обчислення значень вiдлiкiв
(13)
Ця формула є зворотнiм дискретним перетворенням Фур'е (ЗДПФ). Формула прямого (2) та зворотнього (13) дискретного перетворення Фур'є є дискретними аналогами пари перетворень Фур'є для неперервного сигналу
ВИКОНАННЯ ЛАБОРАТОРНИХ ЗАВДАНЬ
ЗАВДАННЯ 1,2,3
Напишіть програму (чи вручну за допомогою текстового редактора) згенеруйте N значень відліків для дискретного представлення таких функцій:
- симетричний прямокутний імпульс (меандр);
- несиметричний прямокутний імпульс тривалістю (перші n точок будуть, а решта нулі, де n ваш варіант);
- симетричний трикутний імпульс;
- пилоподібній імпульс;
- два пилоподібних імпульси;
- пачку з 10 прямокутних імпульсів (меандр);
Амплітуда імпульсів рівна 1.0, період Т=2.0, кількість відліків N=20.
Симетричний прямокутний сигнал
Таблиця 1
K
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

T
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9

X
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1


Файл L11F.txt
0 .000000000000000D+00 .100000000000000D+01
1 .100000000000000D-00 .100000000000000D+01
2 .200000000000000D-00 .100000000000000D+01
3 .300000000000000D-00 .100000000000000D+01
4 .400000000000000D-00 .100000000000000D+01
5 .500000000000000D-00 .100000000000000D+01
6 .600000000000000D-00 .100000000000000D+01
7 .700000000000000D-00 .100000000000000D+01
8 .800000000000000D-00 .100000000000000D+01
9 .900000000000000D-00 .100000000000000D+01
10 .100000000000000D+01 .100000000000000D+01
11 .110000000000000D+01 .100000000000000D+01
12 .120000000000000D+01 .100000000000000D+01
13 .130000000000000D+01 .100000000000000D+01
14 .140000000000000D+01 .100000000000000D+01
15 .150000000000000D+01 .100000000000000D+01
16 .160000000000000D+01 .100000000000000D+01
17 .170000000000000D+01 .100000000000000D+01
18 .180000000000000D+01 .100000000000000D+01
19 .190000000000000D+01 .100000000000000D+01
Файл L11S.txt
0 .100000000000000D+01 .000000000000000D+00
1 .444955726205437D-16 .307917384359384D+01
2 .632924504528419D-16 -.223183949564558D+01
3 .161364719961073D-15 -.150193983749385D+01
4 .659742500852807D-16 -.195911504496736D+01
5 .486114085725742D-15 -.309589928705444D+01
6 .297179143833852D-15 .572459813818051D+00
7 .223497022435708D-15 -.251002686643035D+00
8 .620633538311818D-16 -.175064982658738D+01
9 .311617374027836D-15 -.126312934822112D+01
10 .238932146631306D-15 .159403143648292D+01
11 .311617374027836D-15 .126312934822112D+01
12 .620633538311818D-16 .175064982658738D+01
13 .223497022435708D-15 .251002686643035D+00
14 .297179143833852D-15 -.572459813818051D+00
15 .486114085725742D-15 .309589928705444D+01
16 .659742500852807D-16 .195911504496736D+01
17 .161364719961073D-15 .150193983749385D+01
18 .632924504528419D-16 .223183949564558D+01
19 .444955726205437D-16 -.307917384359384D+01
20 .100000000000000D+01 .000000000000000D+00
Графік згенерованого симетричного прямокутного сигналу

Несиметричний прямокутний сигнал
Таблиця 2
K
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

T
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9

X
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
0
0
0
0


Файл L12F.txt
0 .000000000000000D+00 .100000000000000D+01
1 .100000000000000D-00 .100000000000000D+01
2 .200000000000000D-00 .100000000000000D+01
3 .300000000000000D-00 .100000000000000D+01
4 .400000000000000D-00 .100000000000000D+01
5 .500000000000000D-00 .100000000000000D+01
6 .600000000000000D-00 .100000000000000D+01
7 .700000000000000D-00 .100000000000000D+01
8 .800000000000000D-00 .100000000000000D+01
9 .900000000000000D-00 .100000000000000D+01
10 .100000000000000D+01 .100000000000000D+01
11 .110000000000000D+01 .100000000000000D+01
12 .120000000000000D+01 .100000000000000D+01
13 .130000000000000D+01 .100000000000000D+01
14 .140000000000000D+01 .100000000000000D+01
15 .150000000000000D+01 .000000000000000D+00
16 .160000000000000D+01 .000000000000000D+00
17 .170000000000000D+01 .000000000000000D+00
18 .180000000000000D+01 .000000000000000D+00
19 .190000000000000D+01 .000000000000000D+00
Файл L12S.txt
0 .750000000000000D+00 .000000000000000D+00
1 .226007351067010D+00 -.219911485751286D+01
2 .161803398874990D+00 -.125663706143592D+01
3 .778768257917526D-01 -.314159265358979D+00
4 .526625020286505D-16 -.189254688119154D+01
5 .500000000000000D-01 -.157079632679490D+01
6 .618033988749895D-01 -.628318530717958D+00
7 .396802246667421D-01 .314159265358977D+00
8 .502674853860431D-16 -.168145354796879D+01
9 .357960478079796D-01 -.942477796076933D+00
10 .500000000000000D-01 -.857259176026502D-15
11 .357960478079796D-01 .942477796076933D+00
12 .502674853860431D-16 .168145354796879D+01
13 .396802246667421D-01 -.314159265358977D+00
14 .618033988749895D-01 .628318530717958D+00
15 .500000000000000D-01 .157079632679490D+01
16 .526625020286505D-16 .189254688119154D+01
17 .778768257917526D-01 .314159265358979D+00
18 .161803398874990D+00 .125663706143592D+01
19 .226007351067010D+00 .219911485751286D+01
20 .750000000000000D+00 .000000000000000D+00
Графік згенерованого несиметричного прямокутного сигналу

Симетричний трикутний сигнал
Таблиця 3
K
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

T
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9

X
0
1
2
3
4
5
6
7
8
9
10
9
8
7
6
5
4
3
2
1


Файл L13F.txt
0 .000000000000000D+00 .000000000000000D+00
1 .100000000000000D-00 .100000000000000D+00
2 .200000000000000D-00 .200000000000000D+00
3 .300000000000000D-00 .300000000000000D+00
4 .400000000000000D-00 .400000000000000D+00
5 .500000000000000D-00 .500000000000000D+00
6 .600000000000000D-00 .600000000000000D+00
7 .700000000000000D-00 .700000000000000D+00
8 .800000000000000D-00 .800000000000000D+00
9 .900000000000000D-00 .900000000000000D+00
10 .100000000000000D+01 .100000000000000D+01
11 .110000000000000D+01 .900000000000000D+00
12 .120000000000000D+01 .800000000000000D+00
13 .130000000000000D+01 .700000000000000D+00
14 .140000000000000D+01 .600000000000000D+00
15 .150000000000000D+01 .500000000000000D+00
16 .160000000000000D+01 .400000000000000D+00
17 .170000000000000D+01 .300000000000000D+00
18 .180000000000000D+01 .200000000000000D+00
19 .190000000000000D+01 .100000000000000D+00
Файл L13S.txt
0 .500000000000000D+00 .000000000000000D+00
1 .204317290945307D+00 -.314159265358979D+01
2 .662836909475731D-17 -.167567326552513D+01
3 .242591999815959D-01 -.314159265358979D+01
4 .910192728294575D-17 -.143698534932031D+01
5 .100000000000002D-01 -.314159265358979D+01
6 .526422125426503D-16 .685033939490525D+00
7 .629808091841240D-02 -.314159265358979D+01
8 .628343567325538D-17 -.146013910562100D+01
9 .512542815468448D-02 -.314159265358978D+01
10 .888853460195529D-16 .153175357455006D+01
11 .512542815468448D-02 .314159265358978D+01
12 .628343567325538D-17 .146013910562100D+01
13 .629808091841240D-02 .314159265358979D+01
14 .526422125426503D-16 -.685033939490525D+00
15 .100000000000002D-01 .314159265358979D+01
16 .910192728294575D-17 .143698534932031D+01
17 .242591999815959D-01 .314159265358979D+01
18 .662836909475731D-17 .167567326552513D+01
19 .204317290945307D+00 .314159265358979D+01
20 .500000000000000D+00 .000000000000000D+00
Графік згенерованого симетричного прямокутного сигналу

Пилоподібний сигнал
Таблиця 4
K
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

T
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9

X
0
1
2
3
4
5
6
7
8
9
10
0
0
0
0
0
0
0
0
0


Файл L14F.txt
0 .000000000000000D+00 .000000000000000D+00
1 .100000000000000D-00 .100000000000000D+00
2 .200000000000000D-00 .200000000000000D+00
3 .300000000000000D-00 .300000000000000D+00
4 .400000000000000D-00 .400000000000000D+00
5 .500000000000000D-00 .500000000000000D+00
6 .600000000000000D-00 .600000000000000D+00
7 .700000000000000D-00 .700000000000000D+00
8 .800000000000000D-00 .800000000000000D+00
9 .900000000000000D-00 .900000000000000D+00
10 .100000000000000D+01 .100000000000000D+01
11 .110000000000000D+01 .000000000000000D+00
12 .120000000000000D+01 .000000000000000D+00
13 .130000000000000D+01 .000000000000000D+00
14 .140000000000000D+01 .000000000000000D+00
15 .150000000000000D+01 .000000000000000D+00
16 .160000000000000D+01 .000000000000000D+00
17 .170000000000000D+01 .000000000000000D+00
18 .180000000000000D+01 .000000000000000D+00
19 .190000000000000D+01 .000000000000000D+00
Файл L14S.txt
0 .275000000000000D+00 .000000000000000D+00
1 .202691841193974D+00 -.224894139649341D+01
2 .809016994374947D-01 .125663706143592D+01
3 .615305387037701D-01 -.221859649147892D+01
4 .425325404176020D-01 .942477796076939D+00
5 .390512483795333D-01 -.244685437739309D+01
6 .309016994374947D-01 .628318530717959D+00
7 .308970644815262D-01 -.271664116530481D+01
8 .262865556059567D-01 .314159265358980D+00
9 .278456769830596D-01 -.299891061148339D+01
10 .250000000000000D-01 -.134711154313764D-14
11 .278456769830596D-01 .299891061148339D+01
12 .262865556059567D-01 -.314159265358980D+00
13 .308970644815262D-01 .271664116530481D+01
14 .309016994374947D-01 -.628318530717959D+00
15 .390512483795333D-01 .244685437739309D+01
16 .425325404176020D-01 -.942477796076939D+00
17 .615305387037701D-01 .221859649147892D+01
18 .809016994374947D-01 -.125663706143592D+01
19 .202691841193974D+00 .224894139649341D+01
20 .275000000000000D+00 .000000000000000D+00
Графік згенерованого пилоподібного сигнал

Два пилоподібних сигнали
Таблиця 5
K
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

T
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9

X
0
1
2
3
4
5
6
7
8
9
10
1
2
3
4
5
6
7
8
9


Файл L15F.txt
0 .000000000000000D+00 .000000000000000D+00
1 .100000000000000D-00 .100000000000000D+00
2 .200000000000000D-00 .200000000000000D+00
3 .300000000000000D-00 .300000000000000D+00
4 .400000000000000D-00 .400000000000000D+00
5 .500000000000000D-00 .500000000000000D+00
6 .600000000000000D-00 .600000000000000D+00
7 .700000000000000D-00 .700000000000000D+00
8 .800000000000000D-00 .800000000000000D+00
9 .900000000000000D-00 .900000000000000D+00
10 .100000000000000D+01 .100000000000000D+01
11 .110000000000000D+01 .100000000000000D+00
12 .120000000000000D+01 .200000000000000D+00
13 .130000000000000D+01 .300000000000000D+00
14 .140000000000000D+01 .400000000000000D+00
15 .150000000000000D+01 .500000000000000D+00
16 .160000000000000D+01 .600000000000000D+00
17 .170000000000000D+01 .700000000000000D+00
18 .180000000000000D+01 .800000000000000D+00
19 .190000000000000D+01 .900000000000000D+00
Файл L14S.txt
0 .500000000000000D+00 .000000000000000D+00
1 .500000000000000D-01 .314159265358979D+01
2 .153884176858763D+00 .157079632679490D+01
3 .500000000000000D-01 -.314159265358979D+01
4 .688190960235587D-01 .157079632679490D+01
5 .500000000000002D-01 -.314159265358979D+01
6 .363271264002682D-01 .157079632679489D+01
7 .499999999999998D-01 .314159265358979D+01
8 .162459848116453D-01 .157079632679490D+01
9 .500000000000000D-01 -.314159265358979D+01
10 .113446578469914D-15 .152184524553693D+01
11 .500000000000000D-01 .314159265358979D+01
12 .162459848116453D-01 -.157079632679490D+01
13 .499999999999998D-01 -.314159265358979D+01
14 .363271264002682D-01 -.157079632679489D+01
15 .500000000000002D-01 .314159265358979D+01
16 .688190960235587D-01 -.157079632679490D+01
17 .500000000000000D-01 .314159265358979D+01
18 .153884176858763D+00 -.157079632679490D+01
19 .500000000000000D-01 -.314159265358979D+01
20 .500000000000000D+00 .000000000000000D+00
Графік згенерованих двох пилоподібних сигналів

Пачка з 10 прямокутних сигналів
Таблиця 6
K
0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

T
0
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
1
1,1
1,2
1,3
1,4
1,5
1,6
1,7
1,8
1,9

X
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0
1
0


Файл L16F.txt

0 .000000000000000D+00 .100000000000000D+01
1 .100000000000000D-00 .000000000000000D+00
2 .200000000000000D-00 .100000000000000D+01
3 .300000000000000D-00 .000000000000000D+00
4 .400000000000000D-00 .100000000000000D+01
5 .500000000000000D-00 .000000000000000D+00
6 .600000000000000D-00 .100000000000000D+01
7 .700000000000000D-00 .000000000000000D+00
8 .800000000000000D-00 .100000000000000D+01
9 .900000000000000D-00 .000000000000000D+00
10 .100000000000000D+01 .100000000000000D+01
11 .110000000000000D+01 .000000000000000D+00
12 .120000000000000D+01 .100000000000000D+01
13 .130000000000000D+01 .000000000000000D+00
14 .140000000000000D+01 .100000000000000D+01
15 .150000000000000D+01 .000000000000000D+00
16 .160000000000000D+01 .100000000000000D+01
17 .170000000000000D+01 .000000000000000D+00
18 .180000000000000D+01 .100000000000000D+01
19 .190000000000000D+01 .000000000000000D+00
Файл L16S.txt
0 .500000000000000D+00 .000000000000000D+00
1 .447545209131181D-16 -.208994244104142D+01
2 .298936698014091D-16 -.195130270390726D+01
3 .289776716758409D-16 -.186225312127276D+01
4 .388578058618805D-16 -.157079632679490D+01
5 .306161329334488D-16 -.157079632679490D+01
6 .400296604248672D-16 -.982793723247329D+00
7 .230889022824037D-15 -.474399882791907D+00
8 .819612325799379D-16 -.493941368919581D+00
9 .334221388864417D-15 -.844153986113171D+00
10 .500000000000000D+00 -.110218417373595D-14
11 .334221388864417D-15 .844153986113171D+00
12 .819612325799379D-16 .493941368919581D+00
13 .230889022824037D-15 .474399882791907D+00
14 .400296604248672D-16 .982793723247329D+00
15 .306161329334488D-16 .157079632679490D+01
16 .388578058618805D-16 .157079632679490D+01
17 .289776716758409D-16 .186225312127276D+01
18 .298936698014091D-16 .195130270390726D+01
19 .447545209131181D-16 .208994244104142D+01
20 .500000000000000D+00 .000000000000000D+00
Графік згенерованих пачки з 10 прямокутних сигналів

Обчислюємо: - косинусоїдальну складову спектру вручну.

Обчислюємо: - синусоїдальну складову спектру вручну.

Обчислюємо

Обчислюємо фазу
(в радіанах)
Обчислюємо відновлення значення функції у точці третього відліку

Обчислене відновлене значення функції у точці третього відліку зійшлось з програмним обчисленням.

ВИСНОВОК
На цій лабораторній роботі ми отримали практичні навики використання програми спектрального аналізу RKN, через яку пропустили 6 даних сигналів і отримали їх спектри. Дослідили спектри цих дискретних сигналів різної форми (симетричний прямокутний, несиметричний прямокутний, симетричний трикутний, пилоподібний, два пилоподібних, пачку з 10 прямокутних) та визначили їх особливості та властивості. Всі обчислення, графіки, записали в звіт лабораторної роботи.