LR1 / Analyse

% Считывание данных аудиофайла
[signal, fs] = audioread('model_audio_signal.wav');

% Временная ось для визуализации волновой формы
t = (0:length(signal)-1)/fs;

% Визуализация волновой формы
figure;
subplot(2,1,1);
plot(t, signal);
xlabel('Время (с)');
ylabel('Амплитуда');
title('Волновая форма');

% Вычисление и визуализация амплитудного спектра
% Переход к одному каналу, если сигнал стерео
if size(signal,2) == 2
signal = signal(:,1);
end

% Вычисление FFT
Y = fft(signal);
L = length(signal);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:(L/2))/L;

% Визуализация амплитудного спектра
subplot(2,1,2);
plot(f, P1);
xlabel('Частота (Гц)');
ylabel('|P1(f)|');
title('Амплитудный спектр');

% Дополнительная настройка для улучшения визуализации
xlim([0 2000]); % Ограничение оси частот для лучшей визуализации