[signal_in, fs]= audioread ('Siren_DSP.wav'); % Input from wave file 
vs= 30; % Function parameter setting
Final_output = Doppler_effect_total_new(signal_in, fs, vs); % Function call 
Final_output = Final_output./max(abs(Final_output)); % Normalise the output file
audiowrite ('final_output_doppler_30.wav', Final_output, fs); % Output as wave file