薄板样条差值

最新推荐文章于 2020-09-09 10:16:52 发布

林中漫

最新推荐文章于 2020-09-09 10:16:52 发布

阅读量3.5k

点赞数 1

本文链接：https://blog.csdn.net/caoniyadeniniang/article/details/78107057

版权

本文详细介绍了薄板样条插值的概念，包括问题的引入、差值方法的阐述，并提供了MATLAB实现代码及证明过程。

摘要由CSDN通过智能技术生成

薄板样条插值

一。问题引入

已知三维空间存在一些点集（x,y,z）,现想通过差值的方法插入一些点形成一个光滑的面，问题在于怎么确定免得方程，得到面的方程的方法很多，介绍下薄板样条差值。

二。差值方法

设已知的点集为{X1,X2,....Xj,.......Xn},其中X=（x,y,z）,下面不加证明的给出任意待插入点（x,y）的差值方程

上述公式中共有N+3个系数，系数有下面的公式确定

联立求解上述方程，即可求出插值面的系数，而后带入任意的待插值点（x,y）,即可求出其z值。样条差值的求解步骤如下：

3. MATLAB 代码

clear;
clc;
%% 
%已知点x/3+y/4+y/5=1
point_group=[3 0 0;0 4 0;0 0 5;3 -4 5;-3 4 5;3 4 -5;6 -4 0;6 0 -5;0 8 -5;-3 8 0;0 -4 10;-3 0 10];  % 已知点
%差值范围
x=-5:0.1:5;
y=-5:0.1:5;
[xx,yy]=meshgrid(x,y);
[m,n]=size(xx);
%待插值点
point_est =zeros(m*n,3);
point_est(:,1)=xx(:);
point_est(:,2)=yy(:);

l=length(point_group);
distance=zeros(l); %构建已知点距离矩阵
for i=1:l
    for j=i+1:l
        tmp=norm(point_group(i,1:2)-point_group(j,1:2));
        distance(i,j)=tmp^2*log(tmp);
    end
end
distance=distance+distance';

distance2=zeros(m*n,l); % 带插入点到已知点的距离矩阵
for i=1:m*n
    for j=1:l
        tmp=norm(point_est(i,1:2)-point_group(j,1:2))^2;
        distance2(i,j)=tmp^2*log(tmp);
    end
end

%%
R=zeros(l+3);   % 距离矩阵R
   R1=distance;
   R2=[ones(l,1),point_group(:,1:2)];
   R3=R2';
   R4=zeros(3);
R=[R1,R2;R3,R4];

Z=zeros(l+3,1);
Z(1:l)=point_group(:,3);  % z值矩阵
coff=R\Z;    %系数矩阵A
%% 
R2=zeros(m*n+3,l+3);
     R21=distance2;
     R22=[ones(m*n,1),point_est(:,1:2)];
R2=[R21,R22];
point_est(:,3)=R2*coff;  %插值点的Z值


scatter3(point_est(:,1),point_est(:,2),point_est(:,3))
hold on
scatter3(3,0,0,'red')
scatter3(0,4,0,'red')
scatter3(0,0,5,'red')