Частина тексту файла (без зображень, графіків і формул):
МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА»
Кафедра ЗІ
/
Звіт
до лабораторної роботи №1
з курсу «Методи та засоби технічного захисту інформації»
№п/п
Амплітуда, В
Період коливання, с
Кількість спектральних коефіцієнтів
Роздільча здатність по частоті ΔF, Гц
15
1
1,5
10
1/15
Частина “a”
Навести аналітичний вираз та обчислити спектральні коефіцієнти періодичного сигналу, одержаного шляхом двонапівперіодного випрямлення гармонічного коливання, із параметрами в таблиці 1.
Добрати параметри ДПФ для спектрального аналізу періодичного сигналу, щоб забезпечити вимоги в таблиці 1. Показати графіки часової функції сигналу і його спектра.
Лістинг програми в середовищі Matlab:
clc; % очищення командного вікна
kmax=10; % к-ть спектральних коефіцієнтів
Tk=1.5; % період коливання
dF=1/15; % роздільча здатність по частоті
T0=Tk/2; % період півколивання
Fmax=kmax*1/T0; % максимальна частота
T=1/dF;
Ts=1/(2*Fmax); % період дискретизації
t=0:Ts:(T-Ts); % формування вектора часу
x=abs(sin(2*pi*t/T0)); % аналітичний вираз сигналу
y=fft(x); % застосування ШПФ
f=-Fmax:dF:Fmax-dF; % формування вектора частоти
yy=fftshift(y); % перестановка
xx=abs(yy); % взяття абсолютних значень
figure(1);
t1=0:Ts:T/4; % формування вектора часу для графіка
x1=abs(sin(2*pi*t1/T0)); % аналітичний вираз сигналу для графіка
plot(t1,x1); % виведення графіка сигналу
figure (2);
plot (f,xx); % виведення графіка спектру
Результати роботи:
/
/
Навести аналітичний вираз спектральної густини експоненціального імпульсу s(t)=A×exp(-|a×t|), параметри якого наведено в таблиці 2.
Добрати параметри ДПФ для спектрального аналізу імпульсного сигналу, щоб забезпечити вимоги в таблиці 2. Показати графіки часової функції сигналу і його спектра.
№п/п
Амплітуда, В
Стала згасання a, с-1
Частотний інтервал, Гц
Роздільча здатність по частоті ΔF, Гц
15
1
2,5
5
0,25
Лістинг програми в середовищі Matlab:
clc; % очищення командного вікна
a=2.5; % стала згасання
Am=1; % амплітуда сигналу
dF=0.25; % роздільча здатність по частоті
Fmax=2; % частотний інтервал
Ts=1/(2*Fmax); % частота дискретизації
T=1/dF;
t=0:Ts:(T-Ts); % формування вектора часу
x=Am*exp(-abs(a*t)); % аналітичний вираз сигналу
y=fft (x); % застосування ШПФ
f=(-1/(2*Ts)):dF:1/(2*Ts)-dF; % формування вектора частоти
yy=fftshift (y); % перестановка
xx=abs(yy); % взяття абсолютних значень
figure (1);
plot (t,x); % виведення графіка сигналу
figure (2);
plot (f,xx);% виведення графіка спектру
Результати роботи:
/
/
Навести аналітичний вираз, що описує спектр дискретних сигналів.
Добрати параметри ДПФ для спектрального аналізу дискретизованого трикутного вікна, щоб забезпечити вимоги в таблиці 3. Показати графіки часової функції сигналу і його спектра.
№п/п
Амплітуда, В
Тривалість імпульсу, с
Кількість спектральних пелюсток
Роздільча здатність по частоті ΔF, Гц
15
1
1,5
2
1/15
Лістинг програми в середовищі Matlab:
clc;
tau=1.5; % тривалість імпульса
Am=1; % амплітуда сигналу
dF=1/15; % роздільча здатність по частоті
k=2; % к-ть спектральних пелюсток
Fmax=k/tau; % максимальна частота
N=Fmax/dF; % к-ть відліків
Ts=1/(2*Fmax); % період дискретизації
T=1/dF;
% блок побудови трикутного імпульса
i=1;
for t=0:Ts:(tau/2)-Ts
x(i)=(2*Am*t)/tau;
i=i+1;
end
for t=tau/2:Ts:tau
x(i)=(-2*Am*t)/tau +2*Am;
i=i+1;
end
for t=tau+Ts:Ts:(T-Ts)
x(i)=t*0;
i=i+1;
end
t=0:Ts:(T-Ts); % формування вектора часу
figure (1);
plot (t,x); % виведення графіка сигналу
y=fft(x); % застосування ШПФ
f=-Fmax:dF:Fmax-dF; % формування вектора частоти
yy=fftshift(y); % перестановка
xx=abs(yy); % взяття абсолютних значень
figure (2);
plot (f,xx); % виведення графіка спектру
Результати роботи:
/
/
Лістинг програми в середовищі Matlab:
clc;
load Lab_1_15.mat
Fs=8000;
N=1024;
Ts=1/Fs;
df=1/(N*Ts);
y=Signal;
F0=(N-1)*df;
f=(-F0/2):df:F0/2;
w=blackman (N);
x=fft(y.*w);
xx=fftshift (x);
xx=abs(xx);
t=0:Ts:(N-1)*Ts;
plot (t,Signal);
figure (2);
stem (f,xx);
Результати роботи:
/
/
Висновок: