6028

Кореляційний аналіз сигналів

Лабораторная работа

Коммуникация, связь, радиоэлектроника и цифровые приборы

Кореляційний аналіз сигналів Метароботи: набути навичок кореляційного аналізу сигналів у середовищі MatLAB. Порядокроботи 1. Сформувати вектор відліків часу тривалістю 1 с для частоти дискретизації 256 Гц. Сформувати сигнал випадкового б...

Украинкский

2012-12-27

289 KB

17 чел.

Кореляційний аналіз сигналів

Мета роботи: набути навичок кореляційного аналізу сигналів у середовищі MatLAB.

Короткі теоретичні відомості

Див.: стор. 281 – 347,     Айфичер, Э. Цифровая обработка сигналов. Практический подход /

Э. Айфичер, Б. Джервис. – М. : Издательский дом «Вильямс», 2008. – 992 с. – ISBN 978-5-8459-

0710-3.

Порядок роботи

1. Сформувати вектор відліків часу тривалістю 1 с для частоти дискретизації 256 Гц.

Сформувати сигнал випадкового білого гаусівського шуму  (функція randn). Розрахувати та

побудувати графік автокореляційної функції за формулою (5.1).

2. Сформувати вектор відліків часу тривалістю 1 с для частоти дискретизації 256 Гц. Сформувати дискретний аналог сигналу X(t)=5cos(2pi50t)+2cos(2pi100t). Побудувати графік автокореляційної функції.

3. Сформувати вектор відліків часу тривалістю 1 с для частоти дискретизації 256 Гц. Сформувати послідовність прямокутних імпульсів частотою 10 Гц. Побудувати графік автокореляційної функції.

4. Розрахувати та побудувати графік коефіцієнту взаємної кореляції (формула 5.8) сигналів

п. 2 та 3.

5. Розрахувати та побудувати графіки взаємнокореляційних функцій для пар сигналів: ЕКГ

та плетизмограма, ЕКГ з різних каналів.  

6. Розрахувати та побудувати графіки автокореляційних функцій для оцифрованих сигналів

електрокардіограми, електроенцефалограми, прочитаної з файлу, а також ЕЕГ здорової та хворої людини, сигналів артеріального та внутрішньочерепного тиску та плетизмограми.

7*. Побудувати функцію, яка за допомогою кореляційного аналізу знаходить час затримки

відносно початку координат появи в шумовому сигналі зашумленного прямокутного імпульса з відомою шириною.


1.

x1=0:1/256:1;

x2=randn(1,257);

x3=zeros(1,10);

x2=[x2, x3];

k=10;

for n=1:k

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n-1);

end

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('З сигналом випадкового білого гаусівського шуму');

xlabel('j+1');

ylabel('r12');

2.

x1=0:1/256:1;

t=0:1/256:1;

for i=1:257

X(i)=5*cos(2*pi*50*t(i))+2*cos(2*pi*100*t(i));

end

x2=X;

x3=zeros(1,10);

x2=[x2, x3];

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('З сигналом X(t)=5cos(2pi50t)+2cos(2pi100t)');

xlabel('j+1');

ylabel('r12');

3.

x1=0:1/256:1;

t = 0:.004:1;

x2= SQUARE(2*pi*10*t);

x3=zeros(1,20);

x2=[x2, x3];

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('З послідовністю прямокутних імпульсів частотою 10 Гц ');

xlabel('j+1');

ylabel('r12');

4.

t1=0:1/256:1;

for i=1:257

X(i)=5*cos(2*pi*50*t1(i))+2*cos(2*pi*100*t1(i));

end

x1=X;

t2 = 0:.004:1;

x2= SQUARE(2*pi*10*t2);

x3=zeros(1,20);

x2=[x2, x3];

k=10;

m=0;

while (m<k)

summa=0;

for l=1:257

  summa=summa+x1(l)*x2(l+m);

end

m=m+1;

r12(m)=summa/257;

end

r12;

summa1=0;

summa2=0;

for n=1:257

 summa1=summa1+x1(n)*x1(n);

 summa2=summa2+x2(n)*x2(n);

end

znam=(1/257)*sqrt(summa1*summa2);

for j=1:10

ro12(j)=r12(j)/znam;

end

ro12

plot(ro12)

grid;

title(' графік коефіцієнту взаємної кореляції (формула 5.8) сигналів п. 2 та 3');

xlabel('j+1');

ylabel('ro12');

5.

load ('D:\flash 11-11-2009\ECG_rec\ecg_2.mat');

x1=d;

clear d;

load ('D:\flash 11-11-2009\ECG_rec\ecg_16.mat');

x2=d;

clear d;

x3=zeros(1,10);

x2=[x2, x3];

k=10;

m=0;

while (m<k)

summa=0;

for l=1:4096

  summa=summa+x1(l)*x2(l+m);

end

m=m+1;

r12(m)=summa/4096;

end

r12;

summa1=0;

summa2=0;

for n=1:4096

 summa1=summa1+x1(n)*x1(n);

 summa2=summa2+x2(n)*x2(n);

end

znam=(1/4096)*sqrt(summa1*summa2);

for j=1:10

ro12(j)=r12(j)/znam;

end

ro12

plot(ro12)

grid;

title('Графік взаємнокореляційної функції для пари сигналів ЕКГ з різних каналів');

xlabel('j+1');

ylabel('ro12');

6.

а)

x1=0:1/256:1;

load ('D:\flash 11-11-2009\ECG_rec\ecg_2.mat');

x2=d;

clear d;

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12;

plot(r12)

grid;

title('Графік автокореляційної функції для сигналу електрокардіограми ');

xlabel('j+1');

ylabel('r12');

б)

x1=0:1/256:1;

load ('D:\flash 11-11-2009\EEG_healthy\eeg_healthy_2.mat');

x2=sig;

clear sig;

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('Графік автокореляційної функції для ЕЕГ здорової людини ');

xlabel('j+1');

ylabel('r12');

в)

x1=0:1/256:1;

load ('D:\flash 11-11-2009\EEG_sick\eeg_sick_7.mat');

x2=sig;

clear sig;

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('Графік автокореляційної функції для ЕЕГ хворої людини ');

xlabel('j+1');

ylabel('r12');

г)

x1=0:1/256:1;

fid=fopen('D:\flash 11-11-2009\TBI_ICP.txt');

x2=fscanf(fid,'%f');

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('Графік автокореляційної функції для сигналу артеріального тиску');

xlabel('j+1');

ylabel('r12');

д)

x1=0:1/256:1;

x2 = textread('D:\flash 11-11-2009\TBI_ABP.txt');

k=10;

n=0;

while (n<k)

summa=0;

for i=1:257

  summa=summa+x1(i)*x2(i+n);

end

n=n+1;

r12(n)=summa/257;

end

r12

plot(r12)

grid;

title('Графік автокореляційної функції для сигналу внутрішньочерепного тиску');

xlabel('j+1');

ylabel('r12');


 

А также другие работы, которые могут Вас заинтересовать

56106. Словник-довідник труднощів української мови (вживання архаїзмів та історизмів) 1.31 MB
  Головне - дати дітям поштовх до роздуму, викликати бажання підносити рівень культури власного мовлення, дбайливо поводитись з найбільшим скарбом народу - його мовою, самовдосконалюватися, інтелектуально зростати.
56107. Проблема збагачення, уточнення й активізації словникового запасу молодших школярів 117.5 KB
  Пізнаючи значення слова діти пізнають світ розкриваючи смислові звязки слів у контексті встановлюють взаємозвязок між предметами і явищами навколишньої дійсності внаслідок чого в них розвивається логічне мислення.
56108. Інтерактивний підхід до словникової роботи на уроках української мови в початкових класах 56 KB
  Рекомендації щодо словникової роботи на уроках української мови в початкових класах ІІІ. Вступ Словникова робота на уроках української мови в початкових класах є однією зі складових різноманітної роботи з розвитку мовлення учнів.
56109. О, слово дзвінке озиваєшся в славі до нас ти крізь давні століття імляві. Година спілкування 223.5 KB
  Вступне слово вчителя. Подорож в слово. 1й учень: Слово школа є в усіх словянських мовах; було введено ще в давньоруській мові перший запис датується 1388р. У давньоруській мові було слово схоле від якого походить наше школа от тільки значення його у греків було іншим: спочатку затримка потім заняття у вільний час далі читання а вже пізніше лекція сучасне школа.
56110. Слухаймо Шевченкове слово 128 KB
  У центрі сцени на стіні заквітчаний і прикрашений рушником портрет Тараса Шевченка. Благословен той день і час Коли прослалась килимами Земля яку сходив Тарас Малими босими ногами Земля яку скропив Тарас Дрібними росамисльозами.
56111. Словограй. Логопедичне свято для молодших школярів 51 KB
  У змаганнях беруть участь 2 команди: Букварик і Читайлик.Про які букви ці загадочки В цім будиночку живе Весела як ведмедик буква З нею поруч Вінні-Пух Це її найкращий друг. Упізнали букву Е Нам знайома буква ця Ніби планка на стовпцях.
56112. Счастливый случай Математическая игра в 8 классе 89 KB
  Цель: формировать умение мыслить логически быстро принимать решения правильно и четко выражать свои мысли; стимулировать познавательную активность учащихся; воспитывать ответственность самостоятельность коллективизм.
56114. Інтерактивна дошка SMART Board на уроках математики 1.07 MB
  До таких форм можна віднести поперше мережеве навчання подруге навчання самостійне на основі якогонебудь комплекту навчально-методичних матеріалів УМК. Активізувати увагу учнів можна вже на етапі написання теми уроку.