Octave: Wie Zufallsfunktion mit exponentiell abfallendem Spektrum erzeugen

Hallo,

ich will eine Zufallsfunktion mit exponentiell abfallendem Spektrum 
erzeugen. Mein Ansatz: "Weißes Rauschen erzeugen und anschließend 
filtern" klappt aber anscheinend nicht.

-------------------------------------------------
fs = 1000;    % Abtastfrequenz
fg = 1;       % Grenzfrequenz des Filters
ord=4;        % Ordnung des Tiefpaßfilters
N=10000;      % Anzahl der Messwerte
f = linspace(0,fs,N);
t = 0:1/fs:(N-1)/N*fs;
v = zeros(1,N);
[b,a]=butter(ord,fg/(fs/2),'low');
v = v + 30*filter(b,a,wgn(1,N,0,'dBW','real'));
-------------------------------------------------

Das Amplitudenspektrum müßte eigentlich mit 4*20dB/Dekade abfallen.

Das Ergebnis von
loglog(f,abs(fft(v)))
hat mit diesem Abfall allerdings nicht viel zu tun.

Was mache ich falsch, oder kann Octave das bloß nicht?



Gruß,
  Michael