Hello, can you please help me with this MATLAB script, especially with the code and the commands. I have no idea what to do. Please, do NOT use any code involving FFT. This must literally use the .wav files given in the instructions and follow the instructions given in the comments; otherwise, it's of no use. (It's actually very simple--if you know what you're doing--which I don't. lol) Thank you so much!

Convolution reverb in audio processing. For this assignment, take the short recording of an opera singer in a sound-proof room (singing.wav file) and see how it would sound if recorded at the Milan Opera House, St. Nicolas Church, and the Grand Canyon. Use the following impulse responses from each place. Download the sound file singing.wav for the opera singer, as well as the impulse responses, ScalaMilanOperaHall.wav, StNicolasChurch.wav and GrandCanyon.wav. Submit the Matlab plots for each scenario/location (three in total opera house, church, and Grand Canyon). On each plot there should be three subplots: the song as the first subplot, the impulse response as the second output, and the convoluted song as the third plot. The impulse response of a system/venue/environment is recorded while a pseudo audio impulse is played, such as a loud clap or any loud short bang. The impulse responses of the Milan Opera Hall and St Nicolas Church are downloaded from http://www.voxengo.com/files/impulses/ Their file names are: ScalaMilanOperaHall.wav and StNicolasChurc.wav. The impulse responses of the Grand Canyon is downloaded from http://eleceng.dit.ie/dorran/matlab/impulse_responses/ Its file name is: GrandCanyon.wav

close all; clear all;

% make sure all .wav files are in the same directory of this script

% all audio files are sampled at 44100Hz

% import, plot and play singing recording

[song, fs] = audioread('singing.wav'); % import the song

t = [1:length(song)]/fs; % change from index number to time

subplot(3, 1, 1) % first plot of the three plots in the plot window

plot(t, song) % plot the song

xlabel('t (second)')

ylabel('Relative signal strength')

title('Song')

soundsc(song, fs) % play the song

pause % need to press spacebar to continue the progam to listen to next tone

% USE the code segment above as an example

% now import, plot and play the impulse response

subplot(3, 1, 2)

% CONVOLUTE the audio with the system impulse responses

% using the command conv(sig1, sig2)

% the convolution could take a couple of seconds to complete

% plot and play the output,

% it should sounds as if it was actually recorded in those spaces subplot(3, 1, 3)