径向基函数插值（4）三维数据的插值

直接上代码：过程跟一维二维的一样

clear all; figure;
%**************************************************************************
% 3D Interpolation (takes some time)
%**************************************************************************
%Matlab
[x,y,z,v] = flow(10);
[xi,yi,zi] = meshgrid(.1:.25:10, -3:.25:3, -3:.25:3);
vi = interp3(x,y,z,v,xi,yi,zi);
subplot(2,2,1); slice(xi,yi,zi,vi,[6 9.5],2,[-2 .2]), shading flat; ; title('Interpolation using Matlab function interp3');
subplot(2,2,2); slice(x,y,z,v,[6 9.5],2,[-2 .2]), shading flat; title('data used for interpolation');%RBF
op=rbfcreate([x(:)'; y(:)'; z(:)'], v(:)','RBFFunction', 'multiquadric', 'Stats', 'on');
rbfcheck(op);
rbfvi = rbfinterp([xi(:)'; yi(:)'; zi(:)'], op);
rbfvi = reshape(rbfvi, size(xi));subplot(2,2,3); slice(xi,yi,zi,rbfvi,[6 9.5],2,[-2 .2]), shading flat; ; title('RBF Interpolation');[xr,yr,zr,vr] = flow(100);
subplot(2,2,4); slice(xr,yr,zr,vr,[6 9.5],2,[-2 .2]), shading flat; title('Original 3D data');clear all; figure;
%**************************************************************************
% 3D Interpolation using isosurface vizualization
%**************************************************************************
%Matlab
[x,y,z,v] = flow(10);
[xi,yi,zi] = meshgrid(.1:.55:10, -3:.55:3, -3:.55:3);vi = interp3(x,y,z,v,xi,yi,zi);
subplot(2,2,1); p = patch(isosurface(xi,yi,zi,vi,-3));  title('Interpolation using Matlab function interp3');
isonormals(xi,yi,zi,vi,p)
set(p,'FaceColor','red','EdgeColor','none');
daspect([1 1 1])
view(3); axis tight
camlight
lighting gouraud
% show data used for interpolation
subplot(2,2,2); p = patch(isosurface(x,y,z,v,-3));  title('data used for interpolation');
isonormals(x,y,z,v,p)
set(p,'FaceColor','red','EdgeColor','none');
daspect([1 1 1])
view(3); axis tight
camlight
lighting gouraud%RBF
op=rbfcreate([x(:)'; y(:)'; z(:)'], v(:)','RBFFunction', 'multiquadric', 'Stats', 'on');
rbfcheck(op);
rbfvi = rbfinterp([xi(:)'; yi(:)'; zi(:)'], op);
rbfvi = reshape(rbfvi, size(xi));subplot(2,2,3); p = patch(isosurface(xi,yi,zi,rbfvi,-3)); title('RBF Interpolation');
isonormals(xi,yi,zi,rbfvi,p)
set(p,'FaceColor','red','EdgeColor','none');
daspect([1 1 1])
view(3); axis tight
camlight
lighting gouraud[x,y,z,v] = flow(100);
subplot(2,2,4); p = patch(isosurface(x,y,z,v,-3)); title('Original 3D data');
isonormals(x,y,z,v,p)
set(p,'FaceColor','red','EdgeColor','none');
daspect([1 1 1])
view(3); axis tight
camlight
lighting gouraud

代码下载地址：

点击打开链接

径向基函数插值（4）三维数据的插值相关推荐

matlab径向基函数插值,径向基函数（Radial Basis Function）插值
将RBF用于插值标签(空格分隔):径向基函数插值算法 RBF 曲面重构当高维数据稀疏,需要预测一些数据,需要使用曲面重构的方法. 曲面重构一般可以分为: 插值重构曲面插值里我们一般使用径向基 ...
python给出数据点进行插值_Python对数据进行插值和下采样的方法
使用Python进行插值非常方便,可以直接使用scipy中的interpolate import numpy as np x1 = np.linspace(1, 4096, 1024) x_new = ...
径向基函数（RBF）插值
RBF函数插值径向基函数(Radial Basis Function, RBF)插值的基本形式为 F(r)=∑i=1Nwiφ(∥r−ri∥)F(\boldsymbol{r})=\sum_{i=1}^ ...
matlab rbf函数_基于径向基函数（RBF）的无网格伪谱法与程序实现（2）——微分矩阵...
参考资料 Gregory E. Fasshauer. Meshfree Approximation Methods with MATLAB. P.387 P401 数值实现 Matlab 2019a ...
径向基函数插值（3）二维数据的插值
二维数据的插值过程跟一维数据的过程一样, 只不过在显示二维的数据插值的过程中,需要meshgrid函数产生数据,进行显示插入数据的输出值. 产生二维数据: clear all; figure; %** ...
径向基函数插值（2）一维数据的插值
假设我们有N组数据(xi,yi),,,,,,,,这时我们根据径向基函数我们的目的主要是找到径向基函数中的位置参数的值, 我们要用这些已知数据的值用最小二乘法得到这些参数.. 现在我们用一般的方法mat ...
径向基函数插值（1）
当我们遇到高维的数据,二维或者三维的数据时,由于数据的不完整,需要预测一些数值,一般我们的方法是用曲面重构, 曲面重构可分为:1.插值2.逼近曲面插值我们一般使用径向基插值的方法: RBF (Ra ...
基于径向基函数(RBF)的函数插值
基于径向基函数的函数插值 1. 函数插值 2. RBF函数插值代码实现 1. 函数插值函数插值问题: 用形式简单的插值函数 f^(x)\hat f(x)f^(x) 近似原函数 (1)\qquad ...
基于MATLAB的三维数据插值拟合与三次样条拟合算法（附完整代码）
目录一. 三维插值例题1 二. 高维度插值拟合格式一格式二格式三格式四格式五例题2 三. 单变量三次样条插值例题3 例题4 四. 多变量三次样条插值例题6 一. 三维插值首先三维 ...

径向基函数插值（4）三维数据的插值

径向基函数插值（4）三维数据的插值相关推荐

最新文章

热门文章