51130

Вейвлет-аналіз сигналів

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

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

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

Украинкский

2014-02-06

914.01 KB

7 чел.

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

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

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

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

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

«Вейвлет-аналіз сигналів»

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

групи ДП-92

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

Київ – 2011

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

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

      1. Обрати материнські вейвлет-функції за варіантом з таблиці 1. Отримати та уяснити інформацію щодо обраної вейвлет-родини (функція waveinfo).Розрахувати значення вейвлет-функцій та масштабуючих функцій (якщо вони існують) побудувати їх графіки (функція wavefun). Сформувати вектор масштабних коефіцієнтів a=[.01:.02:.11 .2:.2:1 2:2:30].

      2. Сформувати вектор відліків часу тривалістю 5 с для частоти дискретизації 128 Гц. Сформувати сигнали ділянки синусоїди частотою 2, 2.5, 40, 100, 600, 600.5 Гц. З використанням материнських функцій згідно варіанту побудувати графік модулів вейвлет-коефіцієнтів (функція cwt), скейлограмму та тривимірний графік вейвлет-коефіцієнтів. Зробити висновки щодо скейлограми періодичних синусоїдальних функцій.

      3. Сформувати вектор відліків часу тривалістю 5 с для частоти дискретизації 128 Гц. Сформувати сигнали ділянки синусоїди частотою 40 Гц амплітудою 2 В, зашумленої випадковим сигналом з нульовим середнім

значенням, амплітудою 1 В, 2 В, 10 В. Побудувати спектрограми сигналів,зробити висновки щодо можливості визначити наявність синусоїдального сигналу в шумі.

      4. Сформувати вектор відліків часу тривалістю 10 с для частоти дискретизації 128 Гц. Сформувати сигнали ділянки синусоїди частотою 5 Гц амплітудою 3 В (S1) та 30 Гц амплітудою 1 В (S2). Сформувати на їх основі сигнал(тривалістю 5 с), що дорівнює сумі цих двох сигналів, та інший сигнал, який спочатку містить сигнал 2*S1, а потім сигнал 2*S2 (матиме тривалість 10 с). З використанням материнських функцій згідно варіанту побудувати скейлограму.Повторити для частот 15 та 16 Гц. Зробити висновки щодо можливості визначити

момент зміни частоти сигналу за скейлограмою.

      5. Сформувати вектор відліків часу тривалістю 10 с для частоти дискретизації 128 Гц. Сформувати сигнал ділянки синусоїди частотою 20 Гц. Створити розрив (вставити п’ять нульових відліків) в сигналі в момент часу 5 с. З використанням материнських функцій згідно варіанту побудувати скейлограмму. Зробити висновки щодо можливості часової локалізації моменту розриву в сигналі за скейлограмою.

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

     7. Сформувати вектор відліків часу тривалістю 30 с для частоти дискретизації 128 Гц. Сформувати сигнал одиночного прямокутного імпульсу (функція rectpuls) для тривалості імпульсу 0.1, 1, 5 сек. (для величин зсуву

відносно початку відліку часу 0 та 5 с). З використанням материнських функцій згідно варіанту побудувати скейлограму. Зробити висновки.

    8. Побудувати скейлограму для сигналу, який складається з суми одиночного імпульсу та синусоїдального сигналу, а також для суми одиночного імпульсу, синусоїдального сигналу та випадкового сигналу. Зробити висновки.

    9. Сформувати випадкові сигнали з ненульовим та нульовим середніми значеннями тривалістю 10 с. З використанням материнських функцій згідно варіанту побудувати скейлограму.

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

та плетизмограми.

   11. Побудувати скейлограму та тривимірний графік вейвлет-коефіцієнтів звукових сигналів, які отримані з різною частотою дискретизації.__

close all;

clear all;

clc;

 

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

waveinfo('mexh');

waveinfo('sym');

[psi,x]=wavefun('mexh');

figure

plot(psi);

title('Материнська функція, Mexican hat')

[phi1,psi1,x]=wavefun('sym5');

figure

subplot(2,1,1);

plot(psi1);

title('Материнська функція, Sym')

subplot(2,1,2);

plot(phi1);

title('Функція масштабування, Sym')

 

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

a=[.01:.02:.11 .2:.2:1 2:2:30];

T=5;

Fs=128;

t=[0:T/Fs:T-T/Fs];

fs = [2 2.5 40 100 600 600.5];

for i = 1:length(fs)

   s=sin(2*pi*fs(i)*t);

   figure

   c1=CWT(s, a, 'mexh','scal');

   figure

   c2=CWT(s, a, 'sym5','scal');

   figure

   subplot(1,2,1);

   c11=CWT(s, a, 'mexh','abslvl');colorbar

   title('Розклад модулів вейвлет-коефіцієнтів Mexican hat');

   subplot(1,2,2);

   c22=CWT(s, a, 'sym5','abslvl');colorbar

   title('Розклад модулів вейвлет-коефіцієнтів sym');

   figure

   subplot(1,2,1);

   cwt (s,a,'mexh','3Dplot');

   subplot(1,2,2);

   cwt (s,a,'sym5','3Dplot');

end;

 

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

T=5;

a=[.01:.02:.11 .2:.2:1 2:2:30];

Fs=128;

t=[0:T/Fs:T-T/Fs];

s=2*sin(2*pi*40*t);

A=[1 2 10];

for i=1:length(A)

   N=A(i)*rand(size(s))-0.5;

   sN=s+N;

   figure

   c1=CWT(sN, a, 'mexh','scal');

   figure

   c2=CWT(sN, a, 'sym5','scal');

end;

 

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

T=10;

a=[.01:.02:.11 .2:.2:1 2:2:30];

Fs=128;

t=[0:T/Fs:T-T/Fs];

F=[5 15 30 16];

for i=1:2

   s1=sin(2*pi*F(i)*t);

   s2=sin(2*pi*F(i+2)*t);

   s3=s1+s2;

   s4=[s1 s2];

   figure

   c1=CWT(s3, a, 'mexh','scal');

   figure

   c2=CWT(s3, a, 'sym5','scal');

   figure

   c11=CWT(s4, a, 'mexh','scal');

   figure

   c22=CWT(s4, a, 'sym5','scal');

end;

 

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

T=10;

Fs=128;

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

f=[0:1/T:Fs-1/T];

s=sin(2*pi*20*t);

a=[.01:.02:.11 .2:.2:1 2:2:30];

figure

cwt(s,a,'sym5','plot'); colorbar

figure

cwt(s,a,'mexh','plot'); colorbar

for i=2:length(t)

   t(i)=t(i-1)+1/Fs;

   if(fix(t(i))==5)

       for k=i:i+5

       s(k)=0;

       end;

        break;

   end;

end;

figure

cwt(s,a,'sym5','plot'); colorbar

figure

cwt(s,a,'mexh','plot'); colorbar

 

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

Fs=128;

T=30;

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

F=[0.1 1 10 ];

a=[.01:.02:.11 .2:.2:1 2:2:30];

ts=0;

for i=1:length(F)

   s=rectpuls(t-ts,F(i));

   figure;

   cwt(double(s), a, 'sym5', 'plot'); colorbar

   figure;

   cwt(double(s), a, 'mexh', 'plot'); colorbar

end;

 

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

Fs=128;

a=[.01:.02:.11 .2:.2:1 2:2:30];

T=30;

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

s1=sin(2*pi*t*20);

s2=rectpuls(t,5);

s3=randn(1,length(t));

s4=s1+s2;

s5=s1+s2+s3;

figure

cwt(s4,a,'sym5','plot'); colorbar

figure

cwt(s5,a,'sym5','plot'); colorbar

figure

cwt(s4,a,'mexh','plot'); colorbar

figure

cwt(s5,a,'mexh','plot'); colorbar

 

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

Fs=128;

a=[.01:.02:.11 .2:.2:1 2:2:30];

T=10;

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

s1=randn(1,length(t));

s2=rand(size(t))-0.5;

figure

cwt(s1,a,'sym5','plot'); colorbar

figure

cwt(s2,a,'sym5','plot'); colorbar

figure

cwt(s1,a,'mexh','plot'); colorbar

figure

cwt(s2,a,'mexh','plot'); colorbar

 

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

fs=128;

T=10;

a=[.01:.02:.11 .2:.2:1 2:2:30];

%--------------------------------------------------

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

R1=fread(EKG,'int16');

figure

cwt(R1,a,'sym5','plot'); colorbar

figure

cwt(R1,a,'sym5','3Dplot');

%--------------------------------------------------

load eeg_healthy_10

figure

cwt(sig,a,'sym5','plot'); colorbar

figure

cwt(sig,a,'sym5','3Dplot');

%--------------------------------------------------

load eeg_sick_10

figure

cwt(sig,a,'sym5','plot'); colorbar

figure

cwt(sig,a,'sym5','3Dplot');

%--------------------------------------------------

EKG=fopen('rec_20110922_162251_7_Walenko.bin','r');

R=fread(EKG,[9,inf],'int16');

figure

cwt(R,a,'sym5','plot'); colorbar

figure

cwt(R,a,'sym5','3Dplot');

 

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

a=[.01:.02:.11 .2:.2:1 2:2:30];

[zapz1,fs1,bits1]=wavread('44.1kHz.wav');

[zapz2,fs2,bits2]=wavread('8kHz.wav');

figure

cwt(zapz1,a,'sym5','plot'); colorbar

figure

cwt(zapz1,a,'sym5','3Dplot');

figure

cwt(zapz2,a,'sym5','plot'); colorbar

figure

cwt(zapz2,a,'sym5','3Dplot');

 

Графіки

1)Пункт

2)Пункт

3)Пункт

4)Пункт

5)Пункт

7)Пункт

8)Пункт

9)Пункт

10)Пункт

11)Пункт


 

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

83626. МОЛНИЕЗАЩИТА ПОДСТАНЦИЙ 32.29 KB
  Молниезащита Iкатегории Защита от прямых ударов молнии зданий и сооружений относимых по устройству молниезащиты к I категории должна выполняться отдельно стоящими стержневыми или тросовыми молниеотводами. Защита от прямых ударов молнии зданий и сооружений II категории с неметаллической кровлей должна быть выполнена отдельно стоящими или установленными на защищаемом объекте стержневыми или тросовыми молниеотводами обеспечивающими зону защиты в соответствии с требованиями табл. Установка молниеприемников или наложение молниеприемной сетки не...
83627. Условия и требования норм проектирования по выбору трансформаторов тока (встроенные или отдельно стоящие, 10% погрешность, чувствительность продольной дифференциальной защиты) 39.83 KB
  Трансформаторы тока предназначены для понижения первичного тока до стандартной величины и для отделения цепей измерения и защиты от первичных цепей высокого напряжения. Основные номинальные параметры трансформаторов тока: Номинальное напряжение линейное Uном кВ Номинальный первичный ток I1ном А Номинальный вторичный ток I2ном А 1 или 5 Номинальная вторичная нагрузка с коэффициентом мощности cosφ2=0.8 ВА Номинальный класс точности для измерений Номинальный класс точности для защиты Коэффициент трансформации...
83628. Требования нормами технологического проектирования и САНПИНом к городским подстанциям и электрическим сетям 32.97 KB
  Нормами технологического проектирования к городским ПС и электрическим сетям рассматриваются и регламентируются следующие разделы: 1 Общие положения общие указания обьем и состав проектной документации 2 Расчетные электрические нагрузки расчетные электрические нагрузки жилых зданий электрические нагрузки общественных зданий и промышленных предприятий электрические нагрузки распределительных линий до 1 кВ электрические нагрузки сетей 106 кВ и ЦП укрупненные показатели расхода электроэнергии коммунальнобытовых потребителей 3...
83629. Очерёдность выполнения чертежей «План и разрезы подстанции и плана фундаментов» и что отражено на этих чертежах 28.85 KB
  В рабочие чертежи включают: 1 общие данные по рабочим чертежам; 2 принципиальную схему главных цепей; 3 принципиальные полные схемы релейной защиты управления измерения сигнализации и т.; 4 планы расположения электрооборудования ошиновки и прокладки сетей заземления; 5 планы прокладки электрических сетей; 6 схемы таблицы подключения; 7 кабельный журнал; 8 рабочую документацию задания МЭЗ; 9 эскизные чертежи общего вида НКУ. На схеме указывают: 1 номинальное напряжение сборных шин; 2 типы номинальные токи и сопротивление...
83630. Перечислить основные виды спецификаций и что отражает экспликация на чертеже 30.3 KB
  Описание спецификации дается в ГОСТ 2. Над основной надписью помещаются графы спецификации. В основной надписи спецификации указывают наименование сборочной единицы масштаб ее изображения. Спецификации первого типа чаще используются в конструкторскотехнологических подразделениях а второго на сборочных участках и при работе по заказам.
83631. Перечислить мероприятия предотвращающие электромагнитную наводку на кабели в ОРУ и устройства в ЗРУ 29.45 KB
  Должны выполняться мероприятия исключающие электростатические и электромагнитные наводки на металлических элементах расположенных в помещениях аккумуляторных батарей а также заносы туда высоких потенциалов. Для защиты от электростатической индукции на указанных элементах достаточно надежно присоединить к общему заземляющему устройству электростанций и подстанций гладкие трубы в помещениях аккумуляторных батарей предназначенные для отопления и выдержать расстояния от токоведущих шин до частей здания и других заземленных элементов не менее...
83632. Назначение кабельного журнала и что отражено на листах кабельного журнала 30.27 KB
  В кабельном журнале описывается маркировка каждого кабеля откуда и куда он идет его марка длина кабельной линии и его назначение. Назначение: Всю необходимую информацию о кабелях можно представить непосредственно на схемах: можно указать марку длину способ прокладки кабеля. Однако при построении достаточно большой системы во избежание перегруженности чертежей удобнее вынести эти данные в отдельную таблицу оставив на схемах лишь присвоенные кабелям обозначения. Во время монтажа в кабельный журнал заносятся следующие сведения: номер...
83633. Выполнение чертежей оперативной блокировки на ПС. Что должно быть отражено на чертеже. Какие виды блокировки коммутационных оборудований используются на ПС 30.11 KB
  Основные требования к оперативной блокировке: Блокзамки блокировки должны запирать приводы разъединителей только в крайних положениях включено и отключено; они не должны запирать привод разъединителя в промежуточном положении; Оперативная блокировка не должна давать ложное разрешение на операции с разъединителями при исчезновении напряжения оперативного тока или неисправностях самой оперативной блокировки. Механическая блокировка это блокировка непосредственного действия которая может быть выполнена на близко...
83634. Нелинейные магнитные цепи при постоянных потоках 161 KB
  Для концентрации магнитного поля и придания ему желаемой конфигурации отдельные части электротехнических устройств выполняются из ферромагнитных материалов. Векторные величины характеризующие магнитное поле Наименование Обозначение Единицы измерения Определение Вектор магнитной индукции Тл тесла Векторная величина характеризующая силовое действие магнитного поля на ток по закону Ампера Вектор намагниченности А м Магнитный момент единицы объема вещества Вектор напряженности магнитного поля А м где Гн м магнитная постоянная Основные...