linear-predictive-speech-synth/func/get_impulse_train.m

28 lines
799 B
Mathematica
Raw Normal View History

2020-11-06 19:08:42 +00:00
%% get_impulse_train.m
%%
%% Generate periodic impulse train for use in speech synth
%%
%% Signal of pitch fundamental_freq sampled at sampling_freq
%% for time length_ms
function signal = get_impulse_train(fundamental_freq, sampling_freq, length_ms)
if fundamental_freq > sampling_freq
disp('Fundamental frequency greater than sampling_freq')
signal = [];
return
end
required_samples = ms_to_samples(length_ms, sampling_freq);
pitch_period = 1 / fundamental_freq;
sample_period = 1 / sampling_freq;
cell_length = round(pitch_period / sample_period);
% cell to be repeated into periodic signal
pitch_cell = [1 zeros(1, cell_length - 1)];
required_cells = ceil(required_samples / cell_length);
signal = repmat(pitch_cell, 1, required_cells);
signal = signal(1:required_samples);
end