diff --git a/modifiedwithdecimate.wav b/modifiedwithdecimate.wav new file mode 100644 index 0000000..19c3f7c Binary files /dev/null and b/modifiedwithdecimate.wav differ diff --git a/modifiedwithdownsample.wav b/modifiedwithdownsample.wav new file mode 100644 index 0000000..a5baab8 Binary files /dev/null and b/modifiedwithdownsample.wav differ diff --git a/modulator22.wav b/modulator22.wav new file mode 100644 index 0000000..a6ad482 Binary files /dev/null and b/modulator22.wav differ diff --git a/speech_analysis.m b/speech_analysis.m index b499136..7a73938 100644 --- a/speech_analysis.m +++ b/speech_analysis.m @@ -1 +1,42 @@ -pkg load signal \ No newline at end of file +pkg load signal; + +[y, fs] = audioread("modulator22.wav"); +range = 0:1/fs:(1/fs)*(size(y)-1); +#plot (range, y); +#xlabel ("Time (s)"); +#ylabel ("Amplitude (arbitrary unit)"); +b=fir1(30,1000/(fs/2)); +yfiltered=filter(b,1,y); + +a=butter(8,1000/(fs/2)); +ybuttered=filter(a,1,y); + +[h,w]=freqz(a,1) +freqz_plot(w,h); + + + + + +newfsdwn = 4000; +n=fs/newfsdwn; +ydwn = downsample(ybuttered, n); +audiowrite("modifiedwithdownsample.wav",ydwn,newfsdwn); + +ydec = decimate(ybuttered,6); +fsndec = fs*size(ydec)/size(y); +audiowrite("modifiedwithdecimate.wav",ydec,fsndec); + + + + + +y_min = find(range >= 0,1); +y_max = find(range <= 0.3,1,"last"); +yf= y(y_min:y_max); + +[power, duration] = frequencySpectrum(ybuttered, fs, true); + + + +spectrogram(ybuttered, fs, 5, 30); \ No newline at end of file