Построить график автокорреляционной функции гармонического колебания; изобразить спектр полигармонического сигнала с тремя составляющими

Задания по дисциплине

«Теоретические основы информационно-измерительной техники»

Для студентов специальности 1-360402 «Промышленная электроника» заочного отделения

Задание №1.

Синтезировать:

a. гармонический сигналс частотой F= 2 Гц;

b. полигармонический сигнал с двумя составляющими F1= 2 Гц, F₂=8 Гц;

c. полигармонический сигнал с тремя составляющими F1= 2 Гц, F₂=8 Гц, F₃=16 Гц;

Амплитуды и начальные фазы составляющих выбираются произвольно. Частота дискретизации F_s=512 Гц, Интервал наблюдения T = 2 c.

Построить график автокорреляционной функции гармонического колебания; изобразить спектр полигармонического сигнала с тремя составляющими.

%-------------------------inputdata---------------------------%

T = 2; % время реализации

Fs = 512; % частота дискретизации

n = T *Fs; % длина реализации

t = (0:n-1)/Fs; % шаг дискретизации

Rn = (0:n-1)'; % вектор целых чисел

%%%---------------------signal’ parametrs--------------------%%%

f1 = 2; A1 = 40; phi1 = 0;

f2 = 8; A2 = 30; phi2 = 0;

f3 = 16; A3 = 15; phi3 = 0;

%%%------------------signals’ generation---------------------%%%

s1 = A1*sin(2*pi*t*f1+phi1);

s2 = A1*sin(2*pi*t*f1+phi1) + A2*sin(2*pi*t*f2+phi2);

s3 = A1*sin(2*pi*t*f1+phi1) + A2*sin(2*pi*t*f2+phi2) + A3*sin(2*pi*t*f3+phi3);

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

subplot(3,1,1)

plot(t,s1);title('Signal s1');ylabel('Amplitude');xlabel('Time')

subplot(3,1,2)

plot(t,s2);title('Signal s2');ylabel('Amplitude');xlabel('Time')

subplot(3,1,3)

plot(t,s3);title('Signal s3');ylabel('Amplitude');xlabel('Time')

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

Рис.1. Графикисигналов сигнала

%%%-----------------------Correlation------------------------%%%

[c1,intervals11]= xcorr(s1,'coeff'); % АКФ гармонического колебания

plot(intervals11/n,c1)

title('Autocorrelation');xlabel('Time');ylabel('ACF');grid on;

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

Рис.2. График автокорреляционной функции гармонического колебания.

%%%------------------------Spectrum--------------------------%%%

FT=fft(bartlett(n).*s3');

FT=FT(1:n/2);

df = Rn(1:(n/2))*(1/T); % разрешениепо частоте

plot(df,abs((FT)/max(FT)),'b')

gridon

ylabel('Magnitude')

xlabel('Frequency (Hertz)')

title('Power spectrum')

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

Рис. 3. Спектр полигармонического сигнала с тремя составляющими.

Поиск по сайту: