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');


 

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

33195. Медицинская классификация Л.В.Неймана 13.88 KB
  Важное значение для правильного понимания особенностей психического развития детей с нарушениями слуха для своевременной диагностики и организации их обучения и воспитания в частности для определения типа учреждения в котором должен учиться ребенок имеет классификация таких детей. Проблема дифференциации лиц имеющих нарушения слуха интересовала как врачей так и сурдопедагогов. В нашей стране наибольшее распространение получила медицинская классификация нарушений слуха у детей предложенная Л. Если нарушение слуха распространяется на...
33196. ПСИХОЛОГИЧЕСКАЯ ХАРАКТЕРИСТИКА ГЛУХИХ ДЕТЕЙ 14.42 KB
  Отсутствие речи влечет за собой различные отклонения психического развития. Словесноречевую систему глухого в основном формируют экспрессивные и импрессивные виды письменной речи. Потребность в общении реализуется с помощью предметов действий рисунков и жестовой речи. Наиболее трудной задачей учебновоспитательного процесса является формирование словесной речи и речевого мышления.
33197. Характеристика слабослышащих детей 14.71 KB
  Возникновение тугоухости после сформирования речи не ограничивает дальнейшего достаточно интенсивного ее развития не смотря на её недостатки: относительная бедность словарного запаса сочетается с неправильным усвоением и употреблением понятий. Характерные отклонения в устной речи отражаются и на письменной. На степень сохранности речи влияют время наступления глухоты условия развития и роста ребенка уровень остаточного слуха и качество работы по развитию речи ребенка. Без специального воспитания дети потерявшие слух в 45 лет к поступлению...
33198. ПРОБЛЕМА КОМПЕНСАЦИИ ГЛУХОТЫ 14.85 KB
  Компенсация это возмещение недоразвитых или нарушенных психических функций путем использования сохранных или перестройки частично нарушенных функций. При компенсации психических функций возможно вовлечение в ее реализацию новых структур которые раньше не участвовали в осуществлении данных функций или выполняли при этом другую роль. Второй тип межсистемная компенсация которая осуществляется путем перестройки функциональных систем и включения в работу новых элементов из других структур выполнения ими несвойственных ранее функций....
33199. Закономерности психического развития детей с нарушением слуха 14.53 KB
  Общие закономерности: Закономерность соотношения биологических и социальных факторов в процессе психического развития ребенка. Процесс перехода от одной стадии психического развития к другой предполагает глубокое преобразование всех структурных компонентов психики. Неравномерность психического развития.
33200. МЕТОДЫ ИССЛЕДОВАНИЯ СЛУХА 19.41 KB
  Исследование нарушения слуха затруднено тем что: Ребенок не жалуется на отсутствие слуха. Родители не замечают отсутствие слуха и не принимают эффективных мер. Отсутствие унифицированных методов позволяющих получить достоверные сведения о снижении слуха в детском возрасте.
33201. Особенности зрительного восприятия детей с нарушением слуха 14.16 KB
  Особенности зрительного восприятия детей с нарушением слуха. Это связано с менее подробным анализом и синтезом предметов в прошлом опыте с замедленным формированием произвольности процесса восприятия. Для точного восприятия формы предмета важно выделить его контур. Особенности развития осмысленности восприятия отчетливо проявляются при анализе изображений и картин.
33202. Особенности двигательных ощущений. Осязание глухих и слабослышащих 15.25 KB
  Многие дети имеющие нарушения слуха отстают от нормально слышащих по развитию движений. На протяжении всего дошкольного возраста сохраняется некоторая неустойчивость трудность сохранения статичного и динамичного равновесия недостаточно точная координация неуверенность движений и относительно низкий уровень развития пространственной ориентировки. У большинства имеется отставание в развитии мелких движений пальцев рук артикуляционного аппарата. Замедленная по сравнению со слышащими скорость выполнения отдельных движений влияет на темп...
33203. Кожные ощущения и восприятия 13.83 KB
  Возникают при непосредственном контакте предмета с кожей подразделяются на 4 вида: тактильные вибрационные температурные и болевые. Наибольшее значение для компенсации слуха имеют вибрационные т. Вибрационные ощущения возникают при воздействии меньшей силы чем слуховые. Для того чтобы вибрационные ощущения смогли использовать как средство познания для детей необходимо проводить специальную работу.