Лабы / ОКМвТС лабы / лаб3_Построение спектральной функции сигнала

Пример построения спектральной функции через быстрое Фурье-преобразование

получается сигнал, поврежденный шумовым сигналом со средним значением, равным нулю). Исследовать временной интервал от 0 до 0.6. 1) Построить зашумленный сигнал. 2) Построить спектральную функцию этого сигнала.

t = 0:0.001:0.6; % временной интервал

x = sin(2*pi*50*t)+sin(2*pi*120*t); % несущий сигнал y = x + 2*randn(size(t)); % зашумленный сигнал figure(1)

plot(1000*t(1:100),y(1:100)) % 1000 – это коэффициент, вводится, чтобы увеличить количество точек

title ('Signal Corrupted with Zero-Mean Random Noise') xlabel('time (milliseconds)')

% Начинается выполнение второй части задания. Для построения спектра условно выбрано 512 точек. Строим спектр сигнала y .

Y = fft(y,512); figure(2)

Pyy = Y.* conj(Y) / 512; % учет комплексной части решения

%Строим первые 257 точек

%Формируется частота (256 – это половина от 512 точек) f = 1000*(0:256)/512;

plot(f,Pyy(1:257)) % Так и используется 256 точек, число 257 вводится потому, что отсчет

%начинается не с «0», а с «1».

title('Frequency content of y') xlabel('frequency (Hz)')

Результат:

Задание: постройте спектральную функцию А) синусоидального сигнала; В) зашумленного синусоидального сигнала;

С) одного из сигналов, рассмотренного в лабораторной работе 1.