лаб4 / lab4_3

clc; clear; close all;

Fd = 44100; % Частота дискретизации
Ts = 60; % Продолжительность сигнала в секундах
N = Fd * Ts; % Количество отсчетов

% Генерация белого шума
x = randn(N, 1);

% Частоты для фильтра нижних частот
lowFreq = 20; % Нижняя граница
highFreq = 300; % Верхняя граница

% Применяем ФНЧ к шуму
y = lowpass(x, highFreq, Fd);

% Добавляем эффект "биения" для имитации повторяющихся звуковых сигнатур
for k = 1:N
y(k) = y(k) * (1 + 0.25 * sin(2 * pi * 2 * k / Fd)); % 2 Гц для имитации биения
end

% Нормализация
y = y / max(abs(y));

% Преобразуем в INT16
Am = 8192; % Максимальная амплитуда
y = int16(y * Am);

% Запись в аудиофайл
audiowrite('march.wav', y, Fd);