克里金插值(Kriging)在MATLAB中的实现【优化】
该部分是基于克里金插值(Kriging)在MATLAB中的实现(克里金工具箱),由于在运行过程中有部分问题,基于此做的一些理解+优化。
工具箱的下载见上面的链接,其提供了工具箱。
clc
clearload('data_kriging.mat') %载入数据;我在这里将工具箱中的data1替换了名字%模型参数设置,无特殊情况不需修改,见说明书
theta = [10 10]; lob = [1e-1 1e-1]; upb = [20 20];[dmodel, perf] = dacefit(S, Y, @regpoly0, @corrgauss, theta, lob, upb);
% corrgauss表示高斯相关函数
% lob:相关函数参数的下界向量
% upb:相关函数参数的上界向量%S存储了点位坐标值,Y为观测值
X = gridsamp([0 0;100 100], 40); %创建一个40*40的格网,标注范围为0-100,即格网间距为2.5% X=[83.731 32.36]; %单点预测的实现
%格网点的预测值返回在矩阵YX中,预测点的均方根误差返回在矩阵MSE中
[YX,MSE] = predictor(X, dmodel);
X1 = reshape(X(:,1),40,40); X2 = reshape(X(:,2),40,40);
YX = reshape(YX, size(X1)); %size(X1)=40*40figure(1), mesh(X1, X2, YX) %绘制预测表面
hold on,
plot3(S(:,1),S(:,2),Y,'.k', 'MarkerSize',10) %绘制原始散点数据
hold offfigure(2),mesh(X1, X2, reshape(MSE,size(X1))); %绘制每个点的插值误差大小
在dacefit函数中,参数的含义如下:
- S:输入变量的样本数据矩阵,每一行代表一个样本点,每一列代表一个输入变量。
- Y:响应变量的样本数据矩阵,每一行代表一个样本点,每一列代表一个响应变量。
- @regpoly0:回归多项式函数的句柄,用于拟合输入变量和响应变量之间的回归关系。regpoly0表示零阶多项式,即常数回归模型。
- @corrgauss:相关函数的句柄,用于描述输入变量之间的相关性。corrgauss表示高斯相关函数。
- theta:相关函数的参数向量,用于调整相关函数的形状和范围。具体含义根据相关函数的定义而定。
- lob:相关函数参数的下界向量,用于限制参数的范围。
- upb:相关函数参数的上界向量,用于限制参数的范围。
函数返回值包括:
- dmodel:拟合好的DACE模型,可以用于进行预测和插值。
- perf:拟合模型的性能指标,如均方根误差(RMSE)、均方误差(MSE)等。
通过调整相关函数的参数和范围,可以对DACE模型进行优化,以最好地拟合输入变量和响应变量之间的关系。
备注:如果缺失数据不是规整的,这种方法还是比较难应用。
相关的原理可借鉴:
克里金插值MATLAB程序
克里金插值(Kriging)在MATLAB中的实现【优化】相关推荐
- 克里金插值详细步骤_openlayers4 入门开发系列之前端动态渲染克里金插值 kriging 篇(附源码下载)...
前言 openlayers4 官网的 api 文档介绍地址 openlayers4 api,里面详细的介绍 openlayers4 各个类的介绍,还有就是在线例子:openlayers4 官网在线例子 ...
- 克里金插值参数设置Matlab,克里金插值 调用matlab工具箱
克里金插值 克里金插值是依据协方差函数对随机过程或随机场进行空间建模和插值的回归算法. 克里金插值法的公式为: 式中为待插入的各点的重金属污染值,为已知点的重金属污染值,为每个点的权重值. 用BLUP ...
- 克里金插值 Kriging
https://blog.csdn.net/sinat_23619409/article/details/79974990
- matlab 克里金插值,克里金插值(arcgis克里金插值步骤)
1. 克里格方法概述 克里格方法(Kriging)又称空间局部插值法,是以变异函数理论和结构分析为基础, 在有限区域内对区域化变量进行无偏最优估计的一种方法,是地. 克里金差值最后的出来的克里金误差有 ...
- 克里金插值(Kriging)在MATLAB中的实现(克里金工具箱)
一,直接献上克里金插值MATLAB工具箱 链接:https://pan.baidu.com/s/1p0lt2G8KQ-els7Dyw5jPFg 提取码:wcss 下载后将该程序添加到MATLAB安装文 ...
- ArcGIS中使用协同克里金插值(co-kriging interplotation )对气象数据插值
ArcGIS中如何使用协同克里金插值(co-kriging interplotation )对气象数据插值 ANUSPLIN气象站点数据插值局限性 百度搜索ArcGIS 克里金插值 搭建梯子搜索Arc ...
- 克里金插值---MATLAB程序
最近在研究克里金插值,拜读了@lanainluv的笔记,备受启发, 在这里做一些补充,并分享自己的代码,希望对各位有所帮助,有误的地方请批评指正!(有帮助的话点个赞吧~) @lanainluv的原文链 ...
- GNSS速度场简易MATLAB克里金插值
引言 由于GPS观测点分布离散且不均匀,在进行应变计算和分析前,需要对速度场进行插值,获得均匀分布的速度场.一般采用Kriging 法估计在均匀网格节点上的速度值,需要下载克里金MATLAB工具箱 ...
- Arcgis pro中Arcpy对Excel文件进行克里金插值
具体流程:excel→table(投影)→shp→克里金插值(掩膜) import arcpy from arcpy import env from arcpy.sa import * import ...
最新文章
- 用C语言解“打印九九乘法表”题
- 团队冲刺第二阶段-2
- python无限次输入_在Raspberry Pi上的Python用户输入无限循环内部输入时很多输入错过了输入...
- linux shell 数值计算let和expr
- CrashFinder,找到崩溃代码行
- 值类型和引用类型的区别
- Flash 与数学:圆的切线(2)
- Visio2019安装
- 计算机网络五层结构要点以及功能,试述五层协议的网络体系结构的要点,包括各层的主要功能...
- java基于安卓Android微信小程序的音乐论坛uniAPP小程序
- 祭奠信任Nokia的那些日子
- 计算机主板电杆,嵌入式主板的常见故障解决办法
- 王者荣耀服务器什么时候增加人数,王者荣耀服务器连续两天崩!春节每人游戏时间暴涨75%,玩家要背锅?...
- CSS预处理之less,sass
- 列表的增删改查,嵌套等,及元组的查
- 单点登录(一)| LDAP 协议
- python做三维图片挑战眼力_挑战眼力!三维立体图片你能看出几个
- PTA题目集python判断题
- 小米笔记本Air12.5与Deepin系统体验
- 如何更快地学会任何事情?
热门文章
- PassFab iPhone Backup Unlocker(itunes备份密码找回神器)官方中文版V5.2.10.2 | iphone备份密码解除工具下载 | itunes备份密码忘记怎么办
- 写给大一同学的C语言--顺序结构
- Spring Cloud Alibaba-服务雪崩效应
- 北师大计算机专业硕士提供住宿吗,京研教育:北师北大部分专业硕士需要自己解决住宿,但学校提供奖学金和一些补助...
- Function echartR with package recharts
- mysql主键generated_MySQL之虚拟列(generated-columns)
- Create_generated_clock
- idea 禁止生成generated
- js正则 6-20位只能包含 字母 或者字母+数字 必须字母开头
- 1.hacking-lab基础关[详细writeup]