主要内容

预测

预测高斯过程回归模型的响应

语法

ypred =预测(gprMdl Xnew)
[ypred,ysd]=预测(gprMdl,Xnew)
[ypred, ysd yint] =预测(gprMdl Xnew)
[ypred,ysd,yint]=predict(gprMdl,Xnew,Name,Value)

描述

伊普雷德=预测(gprMdl,Xnew)返回预测的响应伊普雷德对于完全或紧凑高斯过程回归(GPR)模型,gprMdl的预测值Xnew

[伊普雷德,ysd) =预测(gprMdl,Xnew)还返回中预测值处新响应的估计标准偏差Xnew来自训练有素的探地雷达模型。

[伊普雷德,ysd,英特) =预测(gprMdl,Xnew)也返回95%的预测区间,英特的每一行对应的真实响应Xnew

[伊普雷德,ysd,英特) =预测(gprMdl,Xnew,名称、值)还返回带有一个或多个指定的附加选项的预测间隔名称、值成对参数。例如,您可以指定预测间隔的置信水平。

输入参数

全部展开

高斯过程回归模型,指定为RegressionGP(全部)或CompactRegressionGP(紧凑的)对象。

预测器的新值fitrgp用于训练探地雷达模型,指定为桌子或者一个M——- - - - - -D矩阵。M是观察次数和D是训练数据中预测变量的数量。

如果你受过训练gprMdl在一个桌子那么Xnew一定是个桌子包含用于训练的所有预测变量gprMdl

如果你受过训练gprMdl那么,在一个矩阵上Xnew必须是一个数字矩阵D柱。

数据类型:||桌子

名称-值对参数

指定可选的逗号分隔的字符对名称、值参数。的名字是参数名和价值为对应值。的名字必须出现在引号内。您可以按任意顺序指定多个名称和值对参数,如下所示:Name1, Value1,…,的家

预测间隔的显著性水平,指定为逗号分隔对,包括“α”一个标量的取值范围是0到1。

例子:“阿尔法”,0.01指定99%的预测间隔。

数据类型:|

输出参数

全部展开

预测响应值,作为N-by-1向量。

新响应值的估计标准偏差,以N-乘1向量ysd(),= 1, 2, ...,N,包含与预测值对应的新响应的估计标准偏差th一排Xnew从训练有素的探地雷达模型。

的每一行对应的真实响应值的预测区间Xnew,返回为N-by-2矩阵。矩阵的第一列英特包含预测区间的下限,第二列包含预测区间的上限。

例子

全部展开

生成示例数据。

n = 10000;rng (1)%的再现性x = linspace (0.5, 2.5, n) ';Y = sin(10* x) / (2. x)+(x-1)^ 4 + 1.5 *兰德(n, 1);

使用Matern 3/2核函数拟合探地雷达模型,每个预测器具有单独的长度尺度,有效集大小为100。使用回归子集近似方法进行参数估计,使用完全独立的条件方法进行预测。

gprMdl=fitgp(x,y,“内核函数”,“ardmatern32”,...“ActiveSetSize”,100,“FitMethod”,“sr”,“预测方法”,膜集成电路的);

计算预测。

[ypred, ~,开办]=预测(gprMdl x);

绘制数据以及预测和预测间隔。

图(x,y,“r.”);持有图(x,ypred);图(x,yci(:,1),“k——”);情节(x,开办(:,2),“k——”); xlabel(“x”);ylabel (“是的”);

图中包含一个坐标轴。轴包含4个线型对象。

加载示例数据并存储在桌子

负载鱼腥草台=表(量(:1),量(:,2),(:,3),(4):,,,...“VariableNames”, {“我是1”,“meas2”,“meas3”,“meas4”,“物种”});

使用第一个测量值作为响应,其他变量作为预测值,拟合探地雷达模型。

mdl=fitrgp(待定,“我是1”);

计算预测和99%置信区间。

[ypred,~,yci]=预测(mdl,tbl,“α”,0.01);

绘制真实响应和预测以及预测间隔。

图();情节(mdl。Y,“r.”);持有;情节(ypred);情节(开办(:1),凯西:”);情节(开办(:,2),凯西:”);传奇(“真正的反应”,“探地雷达预测”,...“预测下限”,“预测上限”,...“位置”,“最好的”);

图中包含一个坐标轴。轴包含4个线型对象。这些对象代表真实响应、探地雷达预测、预测下限、预测上限。

加载示例数据。

负载(“gprdata.mat”);

数据包含训练和测试数据。训练数据中有500个观测值,测试数据中有100个观测值。该数据有8个预测变量。这是模拟数据。

拟合一个GPR模型使用平方指数核函数与每个预测器的长度尺度。标准化训练数据中的预测器。使用精确的拟合和预测方法。

gprMdl = fitrgp (Xtrain ytrain,“基础”,“不变”,...“FitMethod”,“准确”,“预测方法”,“准确”,...“内核函数”,“ardsquaredexponential”,“标准化”,1);

预测测试数据的响应。

[ytestpred,~,ytestci]=预测(gprMdl,Xtest);

将测试响应和预测绘制出来。

图;情节(欧美,“r”);持有;情节(ytestpred“b”); 图(ytestci(:,1),凯西:”);绘图(ytestci(:,2),凯西:”);传奇(“实际响应”,“探地雷达预测”,...‘下降95%’,‘95%以上’,“位置”,“最好的”);持有

图中包含一个坐标轴。轴包含4个线型对象。这些目标代表实际响应,GPR预测,95%低,95%高。

提示

  • 在训练GPR模型时,可以选择预测方法预测方法中的名称-值对参数fitrgp. 默认的预测方法是“准确”对于N≤10000,N是训练数据中的观察数,和“bcd”(块坐标下降),否则。

  • 标准差的计算,ysd,预测区间,英特,在以下情况下不受金宝app支持:预测方法“bcd”

  • 如果gprMdl是一个CompactRegressionGP对象,则无法计算标准偏差,ysd,或预测间隔,英特对于预测方法等于“sr”膜集成电路的.来计算ysd英特对于预测方法等于“sr”膜集成电路的,使用完全回归(RegressionGP)对象。

选择

你可以用再预测计算训练后的GPR模型在训练数据观测时的预测响应。

扩展功能

介绍了R2015b