diff --git a/semestralka2/analyza.m b/semestralka2/analyza.m new file mode 100644 index 0000000..170a652 --- /dev/null +++ b/semestralka2/analyza.m @@ -0,0 +1,72 @@ +clc; +close all; +clear all; + +[y, Fs] = audioread('sem2.wav'); +t = (0:length(y)-1) / Fs; + +figure; %1 +plot(t, y); +xlabel('Time'); +ylabel('Amplitude'); +title('Waveform sem2.wav'); +grid on; + +[y1, Fs1] = audioread('flac.wav'); +[y2, Fs2] = audioread('flac2.wav'); + +t1 = (0:length(y1)-1) / Fs1; +t2 = (0:length(y2)-1) / Fs2; + +figure; %2 +plot(t1, y1); +xlabel('Time [s]'); +ylabel('Amplitude'); +title('Waveform flac.wav'); +grid on; + +figure; %3 +plot(t2, y2); +xlabel('Time [s]'); +ylabel('Amplitude'); +title('Waveform flac2.wav'); +grid on; + +N1 = length(y1); +N2 = length(y2); + +X1 = fft(y1); +X2 = fft(y2); + +freq_shift1 = (-N1/2 : N1/2 - 1) * (Fs1 / N1); +freq_shift2 = (-N2/2 : N2/2 - 1) * (Fs2 / N2); + +figure; %4 +subplot(2,1,1); +plot(freq_shift1, real(fftshift(X1)), 'b'); +title('Reálna časť posunutého spektra flac.wav'); +xlabel('Frekvencia [Hz]'); +ylabel('Reálna hodnota'); +grid on; +subplot(2,1,2); +plot(freq_shift1, imag(fftshift(X1)), 'r'); +title('Imaginárna časť posunutého spektra flac.wav'); +xlabel('Frekvencia [Hz]'); +ylabel('Imaginárna hodnota'); +grid on; +set(gcf, 'Color', 'none'); set(gca, 'Color', 'none'); + +figure; %5 +subplot(2,1,1); +plot(freq_shift2, real(fftshift(X2)), 'b'); +title('Reálna časť posunutého spektra flac2.wav'); +xlabel('Frekvencia [Hz]'); +ylabel('Reálna hodnota'); +grid on; +subplot(2,1,2); +plot(freq_shift2, imag(fftshift(X2)), 'r'); +title('Imaginárna časť posunutého spektra flac2.wav'); +xlabel('Frekvencia [Hz]'); +ylabel('Imaginárna hodnota'); +grid on; +set(gcf, 'Color', 'none'); set(gca, 'Color', 'none');