102 lines
2.9 KiB
Matlab
102 lines
2.9 KiB
Matlab
pkg load signal
|
|
|
|
#[signal, Fs] = audioread("modulator22.wav");
|
|
#t=0:1/Fs:length(signal)/Fs - 1/Fs;
|
|
#figure; % Create a new figure
|
|
#plot(t,signal);
|
|
#xlabel('Time(s)');
|
|
#ylabel('Signal Amplitude (normalized unit)');
|
|
|
|
|
|
#audiowrite("modifiedmodulator.wav",signal,Fs/2);
|
|
[signal, Fs] = audioread("Sound/modulator22.wav");
|
|
t=0:1/Fs:length(signal)/Fs - 1/Fs;
|
|
#player=audioplayer(signal, Fs, 8)
|
|
#play(player);
|
|
##figure; % Create a new figure
|
|
##plot(t,signal);
|
|
##xlabel('Time(s)');
|
|
##ylabel('Signal Amplitude (normalized unit)');
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
##% Parameters for measurements
|
|
##num_measurements = 100; % Number of measurements
|
|
##durations_dft = zeros(1, num_measurements);
|
|
##durations_fft = zeros(1, num_measurements);
|
|
##
|
|
##for i = 1:num_measurements
|
|
## % Measure time taken for DFT
|
|
## tic;
|
|
## [power_dft, duration_dft] = frequencySpectrum(signal, Fs, false, false, [100 2700]);
|
|
## durations_dft(i) = duration_dft;
|
|
##
|
|
## % Measure time taken for FFT with padding
|
|
## tic;
|
|
## [power_fft, duration_fft] = frequencySpectrum(signal, Fs, true, false, [100 2700]);
|
|
## durations_fft(i) = duration_fft;
|
|
##end
|
|
##
|
|
##% Calculate average and standard deviation
|
|
##avg_duration_dft = mean(durations_dft);
|
|
##std_dev_dft = std(durations_dft);
|
|
##
|
|
##avg_duration_fft = mean(durations_fft);
|
|
##std_dev_fft = std(durations_fft);
|
|
##
|
|
##fprintf('Average duration for DFT: %f seconds\n', avg_duration_dft);
|
|
##fprintf('Standard deviation for DFT: %f seconds\n', std_dev_dft);
|
|
##fprintf('\n');
|
|
##fprintf('Average duration for FFT (with padding): %f seconds\n', avg_duration_fft);
|
|
##fprintf('Standard deviation for FFT (with padding): %f seconds\n', std_dev_fft);
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%Searching for formants without low pass filter
|
|
|
|
##start_time = 0.7; % start time in seconds
|
|
##end_time = 1.3; % end time in seconds
|
|
##
|
|
##start_index = round(start_time * Fs) + 1; % start index
|
|
##end_index = round(end_time * Fs) + 1; % end index
|
|
##
|
|
##cropped_signal = signal(start_index:end_index); % cropped signal
|
|
##
|
|
##[power_dft, duration_dft] = frequencySpectrum(cropped_signal, Fs, false, true, [0 160]);
|
|
|
|
#spectrogram(signal, Fs, 5, 30);
|
|
#spectrogram(signal, Fs, 5, 5);
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%formants with low pass filter
|
|
|
|
N=8;
|
|
fc=2700;
|
|
[b,a]= butter(N,fc/(Fs/2));
|
|
%freqz(b,a);
|
|
|
|
signal_filtered=filter(b,a,signal);
|
|
|
|
%One : 0.75-0.95 / 0.7-1.2
|
|
%Two : 1.3-1.7
|
|
%Three : 1.85-2.2
|
|
|
|
start_time = 2.0; % start time in seconds
|
|
end_time = 2.2; % end time in seconds
|
|
|
|
start_index = round(start_time * Fs) + 1; % start index
|
|
end_index = round(end_time * Fs) + 1; % end index
|
|
|
|
cropped_signal = signal_filtered(start_index:end_index); % cropped signal
|
|
|
|
[power_dft, duration_dft] = frequencySpectrum(cropped_signal, Fs, false, true, [100 3000]);
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|