You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: run_validation_FS_fmod.m
    path of file: ./repos/SQAT/validation/FluctuationStrength_Osses2016/1_AM_tones_fmod
    
    the code of the file until where you have to start completion: % Script run_validation_FS_fmod
%
% This routine plot the validation of fluctuation strength (Osses et al. model)
%  for AM tones as a function of the modulation frequency
%
% - Inputs, signals: AM tones, fc=1 kHz, m=100%, SPL=70 dB, fmod=[1 2 4 8 16 32]
%
% - Reference values taken from: 
%   Osses, Alejandro, Rodrigo García, and Armin Kohlrausch (2016). "Modelling 
%   the sensation of fluctuation strength." Proceedings of Meetings on 
%   Acoustics. Vol. 28. 050005. doi: 10.1121/2.0000410
%
% Fluctuation strength computed using:
%   OUT = FluctuationStrength_Osses2016(insig,fs,method,time_skip,show)
%   type <help FluctuationStrength_Osses2016> for more info
%
% In order to run this code, the user needs to download the dataset of 
%  sound files from zenodo (https://doi.org/10.5281/zenodo.7933206).
%  The obtained folder called `validation_SQAT_v1_0` has to be included in 
%  the `sound_files` folder of the toolbox. 
%
% Author: Gil Felix Greco, Braunschweig, 02/03/2020 (updated in 13.03.2023)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc; clear all; close all;

save_figs = 0; %% save figs flag

dir_out = [fileparts(mfilename('fullpath')) filesep];

%% reference data from Osses et al 2019

ref=[1    2    4    8    16    32    ; % fmod (Hz)
     0.39 0.84 1.25 1.30  0.36  0.06]; % FS values (vacil)

%% compute FS from signals

levelOut = 70;  % target signal's SPL
j=1; % init counter

res=cell(1,size(ref,2));  % declaring for memory allocation

SQAT_version=1; % v1.0
dir_sounds = get_dir_validation_sounds('FluctuationStrength_Osses2016',SQAT_version);

dBFS_in  = 100; % dB full scale convention from the input sounds
dBFS_out =  94; % dB full scale convention in SQAT: amplitude of 1 = 1 Pa, or 94 dB SPL
dB_correction = dBFS_in - dBFS_out;

tic
for i=1:size(ref,2)
    
    fname = sprintf('%sAM-tone-fc-1000_fmod-%.0f_mdept-100-SPL-70-dB.wav',dir_sounds,ref(1,j));
    [insig,fs]=audioread(fname);
   
    insig  = insig * 10^(dB_correction/20); % correct rms SPL to desired levelOut
    SPL(i) = 20.*log10(rms(insig)/2e-5); % verify final SPL of the signal
    
    res{i} = FluctuationStrength_Osses2016(insig,fs,...  % input signal and sampling freq.
                                                       0,...  % method, stationary analysis =0 - window size=length(insig), time_varying analysis - window size=2s
                                                       0,...  % time_skip, in seconds for statistical calculations

    
function of the modulation frequency
%
% - Inputs, signals: AM tones, fc=1 kHz, m=100%, SPL=70 dB, fmod=[1 2 4 8 16 32]
%
% - Reference values taken from: 
%   Osses, Alejandro, Rodrigo García, and Armin Kohlrausch (2016). "Modelling 
%   the sensation of fluctuation strength." Proceedings of Meetings on 
%   Acoustics. Vol. 28. 050005. doi: 10.1121/2.0000410
%
% Fluctuation strength computed using:
%   OUT = FluctuationStrength_Osses2016(insig,fs,method,time_skip,show)
%   type <help FluctuationStrength_Osses2016> for more info
%
% In order to run this code, the user needs to download the dataset of 
%  sound files from zenodo (https://doi.org/10.5281/zenodo.7933206).
%  The obtained folder called `validation_SQAT_v1_0` has to be included in 
%  the `sound_files` folder of the toolbox. 
%
% Author: Gil Felix Greco, Braunschweig, 02/03/2020 (updated in 13.03.2023)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc; clear all; close all;

save_figs = 0; %% save figs flag

dir_out = [fileparts(mfilename('fullpath')) filesep];

%% reference data from Osses et al 2019

ref=[1    2    4    8    16    32    ; % fmod (Hz)
     0.39 0.84 1.25 1.30  0.36  0.06]; % FS values (vacil)

%% compute FS from signals

levelOut = 70;  % target signal's SPL
j=1; % init counter

res=cell(1,size(ref,2));  % declaring for memory allocation

SQAT_version=1; % v1.0
dir_sounds = get_dir_validation_sounds('FluctuationStrength_Osses2016',SQAT_version);

dBFS_in  = 100; % dB full scale convention from the input sounds
dBFS_out =  94; % dB full scale convention in SQAT: amplitude of 1 = 1 Pa, or 94 dB SPL
dB_correction = dBFS_in - dBFS_out;

tic
for i=1:size(ref,2)
    
    fname = sprintf('%sAM-tone-fc-1000_fmod-%.0f_mdept-100-SPL-70-dB.wav',dir_sounds,ref(1,j));
    [insig,fs]=audioread(fname);
   
    insig  = insig * 10^(dB_correction/20); % correct rms SPL to desired levelOut
    SPL(i) = 20.*log10(rms(insig)/2e-5); % verify final SPL of the signal
    
    res{i} = FluctuationStrength_Osses2016(insig,fs,...  % input signal and sampling freq.
                                                       0,...  % method, stationary analysis =0 - window size=length(insig), time_varying analysis - window size=2s
                                                       0,...  % time_skip, in seconds for statistical calculations
                                                       0);    % show results, 'false' (disable, default value) or 'true' (enable)
    
    results(i)=res{1,i}.FSmean; % store mean FS value in vector results[nfmod,nFc] 
    
    j=j+1;
    
end