Po przeglądnięciu kilku plików pdf, do edukowałem się trochę i przyłożyłem się do tego zadania bardziej. Oto efekt moich wypocin:
clear; %wyczyszczenie wszystkich dotychczasowych zmiennych
[x, fs] = wavread("c:/audio.wav"); %wczytanie rzeczywistego dźwieku z pliku (x - dane, fs - czestotliwość próbkowania
n = length(x); %ilość próbek
%t = (0:length(x)-1)/fs; %wyznaczenie wektora czasowego
%x = x + rand(size(t)); %tutaj chciałem dodać szum biały ale coś nie chce działać...
figure(1);
plot(x);
title("Wczytany dźwięk rzeczywisty");
xlabel("Częstotliwość(Hz)");
ylabel("Amplituda(Db)");
legend("Widmo sygnału");
dt = 1/fs; %przyrost wartośći
fx = fft(x); %oblicznie szybkiej transforamty Fourieara sygnału
nx = length(fx); %ilość próbek
base = inv(dt)*(0:(n/2-1))/n; %wyznaczenie osi częstotliwości
powerx = abs(fx(1:nx/2)); %wyznaczenie widma
powerxn = 2*powerx./nx; %normalizacja odpowiedzi
figure(2);
plot(base, powerxn);
title("Charakterystyka amplitudowa wczytanego sygnału");
xlabel("Częstotliwość(Hz)");
ylabel("Amplituda(Db)");
legend("Charakterystyka amplitudowa");
Coś nie mogę dodać tam tego szumu białego, pewnie robię coś źle... Jakbyś mógł jeszcze zerknąć na poprawność tego rozwiązania, to byłbym wdzięczny.