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


 

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

21844. Монтаж каркасных многоэтажных зданий 78 KB
  Все размеры несущих и ограждающих конструкций кратны номинальным размерам с восприятием горизонтальных усилий жёсткими узлами рам; связевой – при которой колонны работают только на вертикальные нагрузки а горизонтальные воспринимаются системой вертикальных дисков и ядер жёсткости; и рамнопланировочного модуля 05м и высотного 06м. Сетка колонн кратна укрупнённому планировочному модулю 15м. Промышленные здания выполняются из железобетонных элементов с сеткой колонн 6х6 или 6х9м высотой этажей 3672м количество этажей от 2 до 12 с...
21845. Управление риском. Основные понятия и определения 97.5 KB
  Все участники проекта заинтересованы в том чтобы исключить возможность провала проекта или хотя бы убытка для себя. Вместе с тем ни у кого из них нет и не может быть уверенности в благополучном исходе проекта риск в любом реальном деле общепризнан. Под неопределенностью понимается неполнота или неточность информации об условиях реализации проекта в том числе о связанных с ними затратах и результатах. Неопределенность связанная с возможностью возникновения в ходе реализации проекта неблагоприятных ситуаций и последствий характеризуется...
21846. Реализация проекта и обратная связь 76.5 KB
  Реализация проекта и обратная связь. Организационные формы управления проектами Структуры управления проектами Это понятие включает в себя вопервых организационные формы и вовторых организационные структуры управления проектом Под организационной структурой управления проектом понимают совокупность взаимосвязанных органов управления расположенных на различных ступенях системы. Использование системы управления проектами предполагает создание специальной группы которая становится самостоятельным участником проекта или структурно входит...
21847. Управление качеством проекта 88 KB
  В чем же здесь роль управления проектами и чем новая система может помочь Оказывается весьма многим прежде всего потому что обеспечение требуемого уровня качества есть задача организационноэкономическая и проектменеджер оказался именно той фигурой которая в состоянии вооружившись методологией управления проектами навести должный порядок в сложном и потому зачастую хаотичном мире среды и окружения проекта. Задача обеспечения качества на должном уровне пронизывает весь жизненный цикл проекта. А потому задача этой главы показать не...
21848. Человеческие аспекты в управлении проектами 49 KB
  Создание проектной команды. В практике проектменеджмента используются два основных вида структуры проектной команды. Матричная форма структуры команды Эта форма применяется как правило для малых и средних проектов продолжительность жизненного цикла которых до двух лет этот критерий в разных странах варьируется от 05 до 2 лет. Преимущества: гибкость в организации и развитии команды; полномочия функциональных отделов фирмы не дублируются; наличие у членов команды уверенности в завтрашнем дне .
21849. Эффективность проекта 139.5 KB
  Эффективность проекта. Разработка концепции проекта. Любой инвестор прежде чем вложить деньги задается вопросом: в какой проект стоит вложить деньги средства сколько хотя бы примерно этих средств будет нужно когда вложенные средства начнут приносить доход сколько прибыли на вложенные средства можно получить и наконец если средства ограничены а они зачастую ограничены то где взять деньги для проекта Разработка концепции состоит из двух этапов: Формирование инвестиционного замысла идеи проекта; Анализ инвестиционных...
21850. Мир управления проектами 121 KB
  Мир управления проектами. Определение проекта. Проект является целенаправленной ориентированной во времени последовательностью как правило однократных комплексных и нерегулярно повторяющихся действий мероприятий или работ со следующими специфическими признаками: однократность и комплексность структуры проекта; сложность структуры проекта; специфичность содержательных и финансовых результатов; заданность сроков начала и окончания и тем самым заданность временной цели; нерегулярность осуществления. Проект это одноразовая совокупность...
21851. Разработка и планирование проекта 130 KB
  Разработка и планирование проекта. Планирование проекта. Сущность планирования состоит в обосновании целей и способов их удовлетворения на основе выявления детального комплекса работ определения эффективных методов и способов ресурсов всех видов необходимых для их выполнения и установления взаимодействия между организациямиучастниками проекта. Основная цель планирования интеграция всех участников проекта для выполнения комплекса работ обеспечивающих достижение конечных результатов проекта.
21852. Управление геомеханическими процессами при системах с искусственным поддержанием выработанного пространства: с магазинированием руды и креплением очистного пространства 433 KB
  Магазинирование полезного ископаемого накопление отбитого полезного ископаемого очистной выработке. Различают полное магазинирование полезного ископаемого если оно ведётся на всю высоту этажа блока или частичное слоевое если оно ведётся в пределах отдельных частей блока. Магазинирование полезного ископаемого составляет технологическую основу специального класса систем разработки. Отличительной особенностью этого класса систем: выемка полезного ископаемого в восходящем порядке; выпуск 3040 отбитой руды; поддержание боков...