pkg load signal
[y,fs] = audioread('modulator22.wav');
audiowrite('outmodulator22.wav',y,fs/2);
n = length(y); 
t=0:1/fs:(n-1)/fs; % time range
figure;
plot(t,y);
xlabel('time (s)');
ylabel('amplitude');


[power,duration]=frequencySpectrum(y, fs, false);
figure;
plot(duration,power);

i=1;
ldft=[];
while (i<6)
  [power,duration]=frequencySpectrum(y, fs, false);
  disp(i);
  disp (duration);
  ldft(i)=duration;
  i=i+1;
end
disp (ldft);
lfft=[];
i=1;
while (i<6)
  [power,duration]=frequencySpectrum(y, fs, true);
  disp(i);
  disp (duration);
  lfft(i)=duration;
  i=i+1;
end
disp (lfft);
i=1;
ratio=[];
while (i<6)
  ratio(i)=(lfft(i))/(ldft(i));
  i=i+1;
end
disp('ratio');
disp(ratio);
avgratio=mean(ratio);
disp('avgratio');
disp(avgratio);

stdratio= std(ratio);
disp('stdratio');
disp(stdratio);
>>>>>>> estevan