Witam, mam problem, mam za zadanie wygenerowac szum bialy, wyswietlic jego widmo, nastepnie przefiltrowac filtrem IIR. Problem w tym, ze widmo takiego szumu powinno byc plaskie, tymczasem mi dla czestotliwosci 0 pojawia się prążek wysoki, a przy reszcie dośc równe, płaskie. Prosiłbym o skorygowanie błedu:
clear all;
close all;
%DANE
'Częstotliwosc próbkowania [Hz]'
fs=input('fs= ');
'Czas obserwacji [s]'
to=input('to= ')
N=ceil(to*fs);
'Czas dyskretny'
t=N/fs
%Generacja szumu białego
szum=rand(1,N);
'PARAMETRY WYGENEROWANEGO SZUMU'
'Wartosc srednia'
m=mean(szum)
'Standardowa dewiacja szumu'
sigma=std(szum)
'ANALIZA W DZIEDZINIE CZĘSTOTLIWO?CI'
PSD=abs(fft(szum));
figure
f=(1:N/2)*fs/N;
plot(f,PSD(1,(1:N/2)) )
hold on
