sample_rate = 44100;


split1 = 100 / (sample_rate/2);
split2 = 200 / (sample_rate/2);
split3 = 300 / (sample_rate/2);
split4 = 400 / (sample_rate/2);

signal = zeros(1, 5000);

signal_power_sectrum = abs(fft(signal)).^2;

band1 = sum(signal_power_sectrum(1:split1));
band2 = sum(signal_power_sectrum(split1:split2));
band3 = sum(signal_power_sectrum(split2:split3));
band4 = sum(signal_power_sectrum(split3:split4));