Міністерство освіти і науки України
Національний Університет “Львівська Політехніка”
Лабораторна робота №4
Із дисципліни:
«Методи та засоби обробки сигналів»
На тему:
АНАЛІЗ КОРЕЛЯЦІЙНИХ ХАРАКТЕРИСТИК СИГНАЛІВ
Львів – 2014
МЕТА РОБОТИ
Знайти автокореляційну та взаємокореляційну характеристику сигналів та дослідити їхні властивості.
ТЕОРЕТИЧНІ ВІДОМОСТІ
Кореляційний аналіз сигналів є однією з найважливіших операцій ЦОС. Він застосовується в таких методах обробки сигналів, як виявлення скритих періодичностей, визначення віддалі до об’єкту локації (може бути радіолокація і ехолокація), визначення швидкості рухомого об’єкту, згладжування шумів, тощо.
Коефіцієнт кореляції засвідчує міру подібності, міру взаємної залежності між сигналами і .
Для , маємо .
При цьому вважають: якщо велике, то між і є кореляція.
Для маємо .
При цьому вважають: якщо r велике по модулю, але від’ємне, то між і є кореляція, що розвивається в протилежному напрямку.
Якщо і , або хоча б одне з них є реалізаціями випадкового процесу з рівномірним розподілом на інтервалі : , при чому імовірність того, що чи набуде якогось значення є однакова, то маємо . При цьому вважають якщо r прямує до нуля, то процеси некорельовані, якщо ж , то ці два процеси корелюються і між ними існує зв’язок.
Автокореляційною функцією (АКФ) називають скалярний добуток сигналу та його копії (здебільшого затриманої в часі). Для неперіодичних сигналів АКФ обчислюють за виразами, (1) для аналогового і (2) для цифрового представлення
, (1) . (2)
де x - досліджуваний сигнал.
Автокореляційна функція застосовується для визначення часової затримки при проходженні сигналів крізь досліджувану систему.
Взаємокореляційною функцією (ВКФ) називають скалярний добуток двох сигналів. Для неперіодичних сигналів ВКФ обчислюють за виразами, (3) для аналогового і (4) для цифрового представлення. Взаємокореляційна функція застосовується для визначення подібності сигналів та розміщення їх на осі часу.
, (3) ,(4)
де x, y - досліджувані сигнали.
1. Завдання.
Сигнали задано наступним чином:
2. Аналітичне представлення сигналів.
Для сигналу X.
Виходячи із заданого графіка, сигнал описується фінітною функцією, тривалістю 5с. Перші 3с. сигнал має сталу амплітуду 10 од. Наступні 5с. – амплітуда змінюється до -9. Отже, аналітичний опис заданого графіком сигналу буде:
Для сигналу Y.
Виходячи із заданого графіка, сигнал описується фінітною функцією, тривалістю 14с. Перші 7с. сигнал має сталу амплітуду -3 од. Наступні 7с. – приймає значення 6.
.
3. Текст програми для виконання в SCILAB.
clc
//очистка екрану
clear
//очистка пам’яті
close;close;close;close
//закриття всіх графічних вікон
m=10;N=2^m;
//кількість точок
T1=5;
//інтервал дослідження першої функції
dt1=T1/N;
//крок дискретизації першої функції
t1=-1:dt1:T1-dt1;
//вектор часу для першої функції
t11=-1:dt1:2-dt1;
//вектор першої частини часу для першої функції
t12=2:dt1:T1-dt1;
//вектор другої частини часу для першої функції
x11=0*t11+10;
//обрахунок першої частини значень першої функції
x12=0*t12-9;
//обрахунок другої частини значень першої функції
x=[x11 x12];
//обчислення першої функції
//t1=resize_matrix(t1,1,1023);
T2=14;
//інтервал дослідження другої функції
dt2=T2/N;
//крок дискретизації другої функції
t21=0:dt2:7-dt2;
//вектор першої частини часу для другої функції
t22=7:dt2:14-dt2;
//вектор другої частини часу для другої функції
t2=[t21 t22]
//вектор часу для другої функції
y21=0*t21-3;
//обрахунок першої частини значень другої функції
y22=0*t22+6;
//обрахунок другої частини значень другої функції
y=[y21 y22];
//обчислення другої функції
figure(1,'BackgroundColor',[1,1,1]);
//вивід функцій на спільному графіку
plot2d2(t1,x,5);
plot2d2(t2,y,3)
legend('Вхідний сигнал - x(n)','Вхідний сигнал - y(n)')
title('Вхідні сигнали')
Bxx=(1/N)*xcorr(x,x);
//знаходження автокореляційної функції першого сигналу
figure(2,'BackgroundColor',[1,1,1]);
//вивід автокореляційної функції першого сигналу на графіку
plot(Bxx);
title('Автокореляційна функція сигналу X')
Byy=(1/N)*xcorr(y,y);
//знаходження автокореляційної функції другого сигналу
figure(3,'BackgroundColor',[1,1,1]);
//вивід автокореляційної функції другого сигналу на графіку
plot(Byy);
title('Автокореляційна функція сигналу Y')
Bxy=(1/(2*N))*xcorr(x,y);
//знаходження взаємокореляційної функції
figure(4,'BackgroundColor',[1,1,1]);
//вивід взаємокореляційної функції двох сигналів на графіку
plot(Bxy);
title('Взаємокореляційна функція сигналів X та Y')
4. Графіки вхідних сигналів на спільній часовій осі.
5. Графіки автокореляційних функцій для кожного із заданих сигналів.
6. Графік взаємокореліційної функції заданих сигналів.
7. Висновки.
З графіка взаємокореляційної функції видно, що вона досягає максимуму в момент часу n=716. Тому можна стверджувати, що вхідні функції є найбільш корельовані саме в цій точці. Також з форми взаємокореляційної функції видно, що коли один з вхідних сигналів рівний нулю, вона також рівна нулю.