插值与拟合_MATLAB智能算法模型+课件讲义代码.zip

得用拟合或插值。

常用的拟合有多项式拟合 POLYFIT

插值有 INTERP1，SPLINE，LAGR1 等。。。

在 Matlab 中，用于曲线和曲面平滑的方法与函数很多，曲线平滑可用 smooth 和 smoothts

连接。

如：[X,Y,Z] = peaks(30);surf(X,Y,Z)

mesh，和 surf 一样，只不过生成的是网格。

surface 用法也一样。

smooth3 作用于体数据，使光滑

isosurface

X,Y,Z 如 meshgrid 的定义。

V 中元素为

1 则表示存在，即要显示。但要连成片的 1 才会显示。

V 中元素如 a>1 时，表示要显示的这个点离上方的网格距离是单位距离的 1/a

圆滑程度由 isovalue 决定，0.9999 是最硬，越接近 0 越圆滑。可同时配合 isocaps.

附录：

一. 样条函数的建立

cs* 三次样条

pp* 分段多项式样条，系数为 t^n 的系数

csapi 插值生成三次样条函数

csape 生成给定约束条件下的三次样条函数

csaps 平滑生成三次样条函数

cscvn 生成一条内插参数的三次样条曲线

getcurve 动态生成三次样条曲线

2. 分段多项式样条函数

ppmak 生成分段多项式样条函数

ppual 计算在给定点处的分段多项式样条函数值

3. B 样条函数

fnval 计算在给定点处的样条函数值

fmbrk 返回样条函数的某一部分（如断点或系数等）

fncmb 对样条函数进行算术运算

fn2fm 把一种形式的样条函数转化成另一种形式的样条函数

fnder 求样条函数的微分(即求导数)

fndir 求样条函数的方向导数

6. 样条曲线端点和节点处理函数

augknt 在已知节点数组中添加一个或多个节点

aveknt 求出节点数组元素的平均值

chbpnt 求出用于生成样条曲线的合适节点数组

2 Responses to “Matlab 曲面插值和拟合”

1.

#

waiting

on 14 Nov 2007 at 2:22 pm

非常之高手，想请教一下，如果我只有三，四个已知三维坐标的点，能不能做到拟

合呢？

t=csvread(’c:/www/dat/1chang.csv’);

x=[t(1,1,1),t(1,4,1),t(1,7,1),t(1,10,1)];

y=[t(1,2,1),t(1,5,1),t(1,8,1),t(1,11,1)];

z=[t(1,3,1),t(1,6,1),t(1,9,1),t(1,12,1)];

xlabel(’X’), ylabel(’Y’), zlabel(’Z’),

colormap, colorbar;

axis equal

如果你想画曲面，要把数据整理成固定的格式，下面一个简单的例子，希望对你有

x=[1:1:2]

y=[1:1:2]

[xx,yy]=meshgrid(x,y)

z=[4.1,4.2;4.3,4.4]

surf(x,y,z)

v(1,:,:) = [1 0 s2 1 s2 0 -1].'*[1 0 -1 0 1];

sph = csape({x,y},v,{'clamped','periodic'});

values = fnval(sph,{0:.1:4,-2:.1:2});

sphere=surf(squeeze(values(1,:,:)),squeeze(values(2,:,:)),squeeze(values(3,:,:)));

axis equal, axis off

[xx,yy]=ndgrid(-1:.5:1,-1:.5:1);

请教:如何求拟合球面sph(图形sphere)上对应点[xx,yy]的z坐标?