如何使用matlab生成高斯随机过程?

时间:2013-12-03 17:32:45

标签: matlab gaussian stochastic-process

如何使用具有零均值和单位方差的Matlab生成高斯随机过程?

高斯随机变量可以通过

实现

W =(1 / SQRT(2 * PI))* EXP( - (叔^ 2)/ 2);

但是高斯随机过程呢?

2 个答案:

答案 0 :(得分:2)

如果高斯过程为白色(不同时刻的样本之间没有相关性),请使用

w = randn(1,n);

其中n是所需的样本数。

如果需要在样本之间引入相关(也就是说,不同时刻的值是相关的),通常的方法是生成白高斯过程,然后应用低通滤波器(使用convfilter)。过程的自相关由过滤器形状决定。

例如,

w = randn(1,500);
y = conv(w,ones(1,100)/10,'same'); %// apply a simple low-pass filter
plot(w)
hold on
plot(y,'r')

由于滤波器引入的(自动)相关性,您可以看到滤波后的信号(红色)具有更平滑的时间变化。

enter image description here

答案 1 :(得分:1)

通过将平均值为0的白噪声和标准差hRMSE通过高斯滤波器g=exp(-(x.^2)/(cl^2/2)),可以生成具有指定相关长度(cl)和RMSE-高度(hRMSE)的随机高斯过程。

此外,您可以在以下链接下找到Matlab代码:http://www.mysimlabs.com/matlab/surfgen/rsgeng1D.m

以下内容已转录:

function [f,x] = rsgeng1D(N,rL,h,cl)
%
% [f,x] = rsgeng1D(N,rL,h,cl) 
%
% generates a 1-dimensional random rough surface f(x) with N surface points. 
% The surface has a Gaussian height distribution function and a Gaussian 
% autocovariance function, where rL is the length of the surface, h is the 
% RMS height and cl is the correlation length.
%
% Input:    N   - number of surface points
%           rL  - length of surface
%           h   - rms height
%           cl  - correlation length
%
% Output:   f   - surface heights
%           x   - surface points
%
% Last updated: 2010-07-26 (David Bergström).  
%

format long;

x = linspace(-rL/2,rL/2,N);

Z = h.*randn(1,N); % uncorrelated Gaussian random rough surface distribution
                     % with mean 0 and standard deviation h

% Gaussian filter
F = exp(-x.^2/(cl^2/2));

% correlation of surface using convolution (faltung), inverse
% Fourier transform and normalizing prefactors
f = sqrt(2/sqrt(pi))*sqrt(rL/N/cl)*ifft(fft(Z).*fft(F));