高斯过程(Gaussian Process。GP)是概率论和数理统计中随机过程(stochastic process)的一种,是一系列服从正态分布的随机变量(random variable)在一指数集(index set)内的组合。
高斯过程中任意随机变量的线性组合都服从正态分布,每个有限维分布都是联合正态分布,且其本身在连续指数集上的概率密度函数即是所有随机变量的高斯测度,因此被视为联合正态分布的无限维广义延伸。高斯过程由其数学期望和协方差函数完全决定,并继承了正态分布的诸多性质。
高斯过程的例子
高斯过程的例子包括维纳过程、奥恩斯坦-乌伦贝克过程等。
对高斯过程进行建模和预测是机器学习、信号处理等领域的重要内容,其中常见的模型包括高斯过程回归(Gaussian Process Regression,GPR)和高斯过程分类(Gaussian Process Classification,GPC)。
高斯过程的命名来自德国数学家卡尔·弗里德里希·高斯(Carl Friedrich Gauss)以纪念其提出正态分布概念。
傻不傻的clc,clear
N=250
alp=1
sig=1
delt=1
a=1
mm=zeros(1,N)
x1=rand(1,N)
x2=rand(1,N)
x3=a.*sqrt(-2*log(x1)).*cos(2*pi*x2)+mm %产生随高斯分布的随机数
y(1)=sig*x3(1)
for n=2:N%产生高斯过程一个样本
y(n)=exp(-alp)*y(n-1)+sig*sqrt(1-exp(-2*alp*delt))*x3(n)
end
i=1:N
plot(i,y)
hold on
plot(i,mm,'-')
title('高斯过程样本函数')
M=0
for i=1:N
M=M+y(i)
end
M=M/N
D=0
for i=1:N
D=D+(y(i)-M)^2
end
D=D/N
for m=1:N %计算自相关函数正半轴
for n=1:N-m+1
rr(n)=y(n)*y(n+m-1)
end
r2(m)=sum(rr)/N
end
for m=1:N-1 %由对称性计算负半轴
r1(N-m)=r2(m+1)
end
j=-N+1:N-1
plot(j,[r1,r2])
title('自相关函数')
ss=fft(y,512)
s=(ss.*conj(ss))/512 %计算功率谱密度估值
f=1000*(0:255)/512
plot(f,s(1:256))
title('功率谱密度')
clc,clear
t=-250:250
r=exp(-abs(t))
plot(t,r)
title('自相关函数')
clc,clear
syms t
f=exp(-abs(t))
g=fourier(f)
clc,clear
w=0:0.01:20
h=2./(w.^2+1)
plot(w,h)
title('实际功率谱密度')
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)