МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА» Кафедра ЗІ / Звіт до лабораторної роботи №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); Результати роботи: / / Висновок: