51132

Кореляційний аналіз сигналів. Властивості сигналів з використанням кореляційного аналізу

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

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

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

Украинкский

2015-01-19

199.85 KB

9 чел.

Національний технічний університет України

«Київський політехнічний інститут»

Факультет електроніки

Лабораторна робота № 7

з дисципліни «Теорія сигналів»

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

Виконав:  студент 3-го курсу

групи ДП-92

 Лонтковський С.А.

Київ – 2011

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

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

1. Сформувати вектор відліків часу тривалістю 1 с для частоти дискретизації 256 Гц. Сформувати сигнал випадкового білого гаусівського шуму (функція randn). Розрахувати та побудувати графік автокореляційної функції за формулою 5.1.

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

3. Сформувати вектор відліків часу тривалістю 10 с для частоти дискретизації 256 Гц. Сформувати сигнали прямокутний імпульс амплітуди 1 В. тривалості 1 сек. в момент часу 3 сек. (сигнал 1), та прямокутного імпульсу амплітуди 1 В тривалістю 1 сек. в момент часу 0 сек. (сигнал 2). Побудувати графік взаємнокореляційної функції (формула 5.8), зробити висновки. Повторити для сигналу 2 амплітуди 10 В.

4. Повторити п. 3 попередньо додавши до сигналу 1 шум з нульовим середнім значенням. Зробити висновки.

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

висновки щодо зв’язку між цими сигналами.

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

close all;

clear all;

clc;

 

%завдання 1-------------------------------------------

Fs=256;

T=1;

t=0:1/Fs:T;

s=randn(size(t));

for i=1:length(s)

   N=length(s);

   r(i)=(1/N)*sum(s(i:N).*s(1:(N-i+1)));

end;

figure

plot(1:length(t),r);xlabel('n');ylabel('r(n)');grid on

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

 

%завдання 2-------------------------------------------

Fs=256;

T=1;

t=0:1/Fs:T;

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

for i=1:length(s)

   N=length(s);

   r(i)=(1/N)*sum(s(i:N).*s(1:(N-i+1)));

end;

figure

plot(1:length(t),r);xlabel('n');ylabel('r(n)');grid on

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

 

% %завдання 3-------------------------------------------

T=10;

j=1;

Fs=256;

t=0:1/Fs:T-1/Fs;

am=[1 2 10];

s=rectpuls(t-3,1);

while j~=4

for i=1:length(s)

s1=am(j)*rectpuls(t,1);

   N=length(s);

   r(i)=(1/N)*sum(s(i:N).*s1(1:(N-i+1)));  

end;

subplot(3,1,j),plot(1:length(t),r);xlabel('n');ylabel('r(n)');grid on

title('Графік взаємнокореляційної функції сигналiв s i s1 ');j=j+1;

end;

 

%завдання 4--------------------------------------------

T=10;

j=1;

Fs=256;

t=0:1/Fs:T-1/Fs;

am=[1 2 10];

s=rectpuls(t-3,1);

re=rand(1,length(t));

while j~=4

for i=1:length(s)

   s1=am(j)*rectpuls(t,1)+re;

   N=length(s);

   r(i)=(1/N)*sum(s(i:N).*s1(1:(N-i+1)));  

end;

subplot(3,1,j),plot(1:length(t),r);xlabel('n');ylabel('r(n)');grid on

title('Графік взаємнокореляційної функції сигналiв s i s1 ');j=j+1;

end

 

%завдання 5-------------------------------------------

% load eeg_healthy_10

EKG=fopen('vavreschuk','r');

sig=fread(EKG,'int16');

for i=1:length(sig)

       N=length(sig);

   r(i)=(1/N)*sum(sig(i:N).*sig(1:(N-i+1)));  

end;

plot(1:length(sig),r);xlabel('n');ylabel('r(n)');grid on

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

Графіки:

1)        2)

3)       4)

5)


 

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

39564. Разработка модели анализа структуры затрат на услуги ООО «Противо Пожарная защита» 2.68 MB
  В экономической литературе можно найти множество вариантов снижения затрат на предприятии, однако комплексного подхода для разработки модели анализа структуры затрат на услуги предприятия нет. В связи с этим разработка такой модели является важным и актуальным шагом в изучении экономики предприятия.
39565. Загрузка синхронного транспортногой модуля STM – RR или STM – 1 1.62 MB
  Объединяются 1300 каналов тональной частоты и 6 каналов звукового вещания первого класса в системе с временным разделением каналов и 8ми разрядной импульснокодовой модуляцией. Рассчитать временные и частотные характеристики нарисовать структурную схему объединения и разделения каналов с учетом плезиохронной цифровой иерархии. Рассчитаем частоты дискретизации для каналов ТЧ и звукового вещания. Определим общее число каналов.
39566. Организационная схема предприятия «Вольво» ОАО «Белмагистральавтотранс» 377 KB
  Проверка давления включения и отключения регулятора давления 2. Проверка стояночного тормоза и запорного клапана 3. Проверка педали и тормозного клапана Проверки снаружи 4. Проверка крепления кабины гидросистемы запорных устройств решетки радиатора и механизма подъема кабины 5.
39567. Проект компрессора высокого давления с гладкой проточной частью для ТРДД с тягой 140 кН 1.34 MB
  Описание конструкции КВД 7 3.Описание конструкции КВД. КВД осевой левого вращения восьмиступенчатый приводится во вращение турбиной газогенератора. КВД состоит из следующих узлов: входного направляющего аппарата ВНА; корпуса с направляющим аппаратом 1 ст.
39568. БЛОК УПРАВЛЕНИЯ ВАКУУМНОЙ СИСТЕМЫ УСКОРИТЕЛЯ ЭЛЕКТРОНОВ УРТ-1 12.63 MB
  БЛОК УПРАВЛЕНИЯ ВАКУУМНОЙ СИСТЕМЫ УСКОРИТЕЛЯ ЭЛЕКТРОНОВ УРТ1 ВЫПУСКНАЯ КВАЛИФИКАЦИОННАЯ РАБОТА БАКАЛАВРА Пояснительная записка 210100. ЭКСПЕРИМЕНТАЛЬНАЯ УСТАНОВКА БЛОК УПРАВЛЕНИЯ ВАКУУМНОЙ СИСТЕМОЙ. Объект исследования: автоматизация контроля и управления вакуумной системы наносекундного ускорителя электронов. Цель работы: разработка блока управления вакуумной системой ускорителя электронов создание программного обеспечения для персонального компьютера предназначенного для управления и отладки блока управления.
39569. Проект мероприятий по увеличению объема реализации услуг предприятия туристического бизнеса ООО «Флагман Груп» г. Москва 238.77 KB
  Основным видом деятельности ООО «Флагман Груп» является осуществление туристической деятельности как внутри Российской Федерации, так и за рубежом.
39570. Проектирование гидроэлектростанции приплотинного типа на расчетный напор 158,75 м 439.28 KB
  В соответствии с номенклатурой осевые гидротурбины (Пр и ПЛ) применяют на напоры H=2…80м. Применение Пр гидротурбин выгодно при незначительных колебаниях напора и работе ГЭС в энергосистеме. Кроме того Пр экономически предпочтительнее чем Пл гидротурбины
39571. Совершенствование системы мотивации персонала (на примере РИА «Волекс») 127.25 KB
  В РИА «Волекс» для стимулирования труда работников руководство использует экономические, социальные и административные методы мотивации. Значимым экономическим методом мотивации на предприятии является зарплата, начисляемая по повременно-премиальной и сдельной системам оплаты труда, вся выработка рабочего оплачивается по одной постоянной сдельной расценке, для руководителей, специалистов и служащих используется система должностных окладов, выплачивается премия.