![](/user_photo/71224_AGCDL.jpg)
Добавил:
KaFaka
t.me
Инфо для ГУАП студентов от меня: https://kafaka.notion.site/99e6d9b70ca74f7baef3daea17839e5a
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:LR4 / LR4_create
.m % Параметры симуляции
Fs = 44100; % Частота дискретизации, Гц
duration = 10; % Длительность сигнала, секунды
t = linspace(0, duration, Fs*duration); % Вектор времени
x = randn(1, length(t)); % Генерация белого шума
% Формирующий фильтр: комбинация полосового и низкочастотного фильтров
% Полосовой фильтр с центром на 2000 Гц и шириной 1000 Гц
y_band = bandpass(x, [1500, 2500], Fs);
% Низкочастотный фильтр с частотой среза 300 Гц для модуляции
y_low = lowpass(x, 300, Fs);
% Подготовка параметров для динамической модуляции
a = 1; % Скорость изменения сигмоиды
b = 0.2; % Смещение сигмоиды
mod_signal = 0.5 + 0.5 * sin(2 * pi * 0.1 * t) .* sigmoid(t - 5, a, b);
% Динамическое изменение интенсивности
y_modulated = y_band .* mod_signal;
% Сохранение и воспроизведение результата
audiowrite('dynamic_intensity_noise.wav', y_modulated, Fs);
sound(y_modulated, Fs);
pause(duration + 1);
% Локальная функция сигмоиды для плавного перехода интенсивности
function y = sigmoid(x, a, b)
y = 1 ./ (1 + exp(-a * (x - b)));
end
Fs = 44100; % Частота дискретизации, Гц
duration = 10; % Длительность сигнала, секунды
t = linspace(0, duration, Fs*duration); % Вектор времени
x = randn(1, length(t)); % Генерация белого шума
% Формирующий фильтр: комбинация полосового и низкочастотного фильтров
% Полосовой фильтр с центром на 2000 Гц и шириной 1000 Гц
y_band = bandpass(x, [1500, 2500], Fs);
% Низкочастотный фильтр с частотой среза 300 Гц для модуляции
y_low = lowpass(x, 300, Fs);
% Подготовка параметров для динамической модуляции
a = 1; % Скорость изменения сигмоиды
b = 0.2; % Смещение сигмоиды
mod_signal = 0.5 + 0.5 * sin(2 * pi * 0.1 * t) .* sigmoid(t - 5, a, b);
% Динамическое изменение интенсивности
y_modulated = y_band .* mod_signal;
% Сохранение и воспроизведение результата
audiowrite('dynamic_intensity_noise.wav', y_modulated, Fs);
sound(y_modulated, Fs);
pause(duration + 1);
% Локальная функция сигмоиды для плавного перехода интенсивности
function y = sigmoid(x, a, b)
y = 1 ./ (1 + exp(-a * (x - b)));
end
Соседние файлы в папке LR4