clc;
clear;

%% read audio data (amplitude vector of the signal)
data = mydataread('Aretha Franklin - Think.wav');

%% compute average
Average = mean(data)

%% compute variance
Variance = var(data)

%% compute autocorrelation
x = data';  % original signal
k = 1       % amount of signal delay, change k to see different results
x_k = [zeros(1, k) x(1:(length(x)-k))];  % delayed signal
Autocorrelation = x * x_k' / length(x)