matlab griddata？请教Matlab的griddata的用法
大家好，今天来为大家解答matlab griddata这个问题的一些问题点，包括请教Matlab的griddata的用法也一样很多人还不知道，因此呢，今天就来为大家分析分析，现在让我们一起来看看吧！如果解决了您的问题，还望您关注下本站哦，谢谢~
Matlab常用命令
Matlab常用命令汇总
记住Matlab中一些常用的命令，对于我们学习EDA来说用处很大!下面我为大家整理了关于Matlab的常用命令，希望对你有所帮助。
一、常用对象操作：除了一般windows窗口的常用功能键外。
1、!dir可以查看当前工作目录的文件。!dir&可以在dos状态下查看。
2、who可以查看当前工作空间变量名， whos可以查看变量名细节。
3、功能键：
功能键快捷键说明
方向上键 Ctrl+P返回前一行输入
方向下键 Ctrl+N返回下一行输入
方向左键 Ctrl+B光标向后移一个字符
方向右键 Ctrl+F光标向前移一个字符
Ctrl+方向右键 Ctrl+R光标向右移一个字符
Ctrl+方向左键 Ctrl+L光标向左移一个字符
home Ctrl+A光标移到行首
End Ctrl+E光标移到行尾
Esc Ctrl+U清除一行
Del Ctrl+D清除光标所在的字符
Backspace Ctrl+H删除光标前一个字符
Ctrl+K删除到行尾
Ctrl+C中断正在执行的命令
4、clc可以命令窗口显示的内容，但并不清除工作空间。
二、函数及运算
1、运算符：
+：加，-：减，*：乘，/：除，：左除 ^：幂，‘：复数的共轭转置，()：制定运算顺序。
2、常用函数表：
sin()正弦(变量为弧度)
Cot()余切(变量为弧度)
sind()正弦(变量为度数)
Cotd()余切(变量为度数)
asin()反正弦(返回弧度)
acot()反余切(返回弧度)
Asind()反正弦(返回度数)
acotd()反余切(返回度数)
cos()余弦(变量为弧度)
exp()指数
cosd()余弦(变量为度数)
log()对数
acos()余正弦(返回弧度)
log10()以10为底对数
acosd()余正弦(返回度数)
sqrt()开方
tan()正切(变量为弧度)
realsqrt()返回非负根
tand()正切(变量为度数)
abs()取绝对值
atan()反正切(返回弧度)
angle()返回复数的相位角
atand()反正切(返回度数)
mod(x,y)返回x/y的余数
sum()向量元素求和
3、其余函数可以用help elfun和help specfun命令获得。
4、常用常数的值：
pi 3.1415926…….
realmin最小浮点数，2^-1022
i虚数单位
realmax最大浮点数，(2-eps)2^1022
j虚数单位
Inf无限值
eps浮点相对经度=2^-52
NaN空值
三、数组和矩阵：
1、构造数组的方法：增量发和linspace(first,last,num)first和last为起始和终止数，num为需要的数组元素个数。
2、构造矩阵的方法：可以直接用[ ]来输入数组，也可以用以下提供的函数来生成矩阵。
ones()创建一个所有元素都为1的矩阵，其中可以制定维数，1，2….个变量
zeros()创建一个所有元素都为0的矩阵
eye()创建对角元素为1，其他元素为0的矩阵
diag()根据向量创建对角矩阵，即以向量的元素为对角元素
magic()创建魔方矩阵
rand()创建随机矩阵，服从均匀分布
randn()创建随机矩阵，服从正态分布
randperm()创建随机行向量
horcat C=[A,B]，水平聚合矩阵，还可以用cat(1,A,B)
vercat C=[A;B]，垂直聚合矩阵,还可以用cat(2,A,B)
repmat(M,v,h)将矩阵M在垂直方向上聚合v次，在水平方向上聚合h次
blkdiag(A，B)以A，和B为块创建块对角矩阵
length返回矩阵最长维的的长度
ndims返回维数
numel返回矩阵元素个数
size返回每一维的长度，[rows,cols]=size(A)
reshape重塑矩阵，reshape(A,2,6),将A变为2×6的矩阵，按列排列。
rot90旋转矩阵90度，逆时针方向
fliplr沿垂轴翻转矩阵
flipud沿水平轴翻转矩阵
transpose沿主对角线翻转矩阵
ctranspose转置矩阵，也可用A’或A.’，这仅当矩阵为复数矩阵时才有区别
inv矩阵的逆
det矩阵的行列式值
trace矩阵对角元素的和
norm矩阵或矢量的范数，norm(a，1)，norm(a，Inf)…….
normest估计矩阵的最大范数矢量
chol矩阵的cholesky分解
cholinc不完全cholesky分解
lu LU分解
luinc不完全LU分解
qr正交分解
kron(A，B) A为m×n，B为p×q，则生成mp×nq的矩阵，A的每一个元素都会乘上B，并占据p×q大小的空间
rank求出矩阵的刺
pinv求伪逆矩阵
A^p对A进行操作
A.^P对A中的每一个元素进行操作
四、数值计算
1、线性方程组求解
(1)AX=B的解可以用X=AB求。XA=B的解可以用X=A/B求。如果A是m×n的矩阵，当m=n时可以找到唯一解，mn，超定系统，至少找到一组解。如果A是奇异的，且AX=B有解，可以用X=pinv(A)×B返回最小二乘解
(2)AX=b, A=L×U，[L,U]=lu(A), X=U(Lb),即用LU分解求解。
(3)QR(正交)分解是将一矩阵表示为一正交矩阵和一上三角矩阵之积，A=Q×R[Q,R]=chol(A), X=Q(Ub)
(4)cholesky分解类似。
2、特征值
D=eig(A)返回A的所有特征值组成的矩阵。[V,D]=eig(A),还返回特征向量矩阵。
3、A=U×S×UT，[U,S]=schur(A).其中S的对角线元素为A的特征值。
4、多项式Matlab里面的多项式是以向量来表示的，其具体操作函数如下：
conv多项式的乘法
deconv多项式的除法，【a，b】=deconv(s)，返回商和余数
poly求多项式的系数(由已知根求多项式的系数)
polyeig求多项式的特征值
Polyfit(x，y，n)多项式的.曲线拟合，x，y为被拟合的向量，n为拟合多项式阶数。
polyder求多项式的一阶导数，polyder(a，b)返回ab的导数
[a,b]=polyder(a，b)返回a/b的导数。
polyint多项式的积分
polyval求多项式的值
polyvalm以矩阵为变量求多项式的值
residue部分分式展开式
roots求多项式的根(返回所有根组成的向量)
注：用ploy(A)求出矩阵的特征多项式，然后再求其根，即为矩阵的特征值。
5、插值常用的插值函数如下：
griddata数据网格化合曲面拟合
Griddata3三维数据网格化合超曲面拟合
interp1一维插值(yi=interp1(x,y,xi,’method’)Method=nearest/linear/spline/pchip/cubic
Interp2二维插值zi=interp1(x,y,z,xi,yi’method’),bilinear
Interp3三维插值
interpft用快速傅立叶变换进行一维插值，help fft。
mkpp使用分段多项式
spline三次样条插值
pchip分段hermit插值
6、函数最值的求解
fminbnd(‘f’，x1，x2，optiset(,))求f在x1和x2之间的最小值。Optiset选项可以有‘Display’+‘iter’/’off’/’final’,分别表示显示计算过程/不显示/只显示最后结果。fminsearch求多元函数的最小值。fzero(‘f’，x1)求一元函数的零点。X1为起始点。同样可以用上面的选项。
五、图像绘制：
1、基本绘图函数
plot绘制二维线性图形和两个坐标轴
plot3绘制三维线性图形和两个坐标轴
fplot在制定区间绘制某函数的图像。fplot(‘f’，区域，线型，颜色)
loglog绘制对数图形及两个坐标轴(两个坐标都为对数坐标)semilogx绘制半对数坐标图形
semilogy绘制半对数坐标图形
2、线型：颜色线型
y黄色.圆点线 v向下箭头
g绿色-.组合>向右箭头
b蓝色+点为加号形<向左箭头
m红紫色 o空心圆形 p五角星形
c蓝紫色*星号 h六角星形
w白色.实心小点 hold on添加图形
r红色 x叉号形状 grid on添加网格
k黑色 s方形-实线
d菱形--虚线 ^向上箭头
3、可以用subplot(3，3，1)表示将绘图区域分为三行三列，目前使用第一区域。此时如要画不同的图形在一个窗口里，需要hold on。
;MATLAB中的插值函数griddata()运行出现NAN怎么办
同遇到这个隐蔽的问题，解决了，所以顺便来挖坟答一下题。
========================================
没有griddata的具体实现算法方面的信息，但是插值原理应该是基于散乱数据点生成局部区域的插值查询。而且这个插值似乎是要求“内插”的，即查询点必须处于输入样本XY的”包围“状态中，否则就会报NaN查询结果。
matlab 2012a中的帮助是这样说的：
The method defines the type of surface fit to the data. The'cubic' and'v4' methods produce smooth surfaces while'linear' and'nearest' have discontinuities in the first and zero'th derivatives, respectively. All the methods except'v4' are based on a Delaunay triangulation of the data. If method is [], then the default'linear' method is used.
Occasionally, griddata might return points on or very near the convex hull of the data as NaNs. This is because roundoff in the computations sometimes makes it difficult to determine if a point near the boundary is in the convex hull.
========================================
所以这个问题无法避免。（虽然实测使用nearest方法没有产生NaN，但因没有解读其算法，不确定是否绝对不出现NaN结果）
如果只是为了绘出有效数据，把结果中的NaN数据删掉就行了。
如果想得到所有的查询值，把NaN结果全部取出来，相应的X,Y重新用nearest方法查询一次。或者自己写一个允许用近邻点外推插值的算法对其特殊处理。但还是会与griddata内部方法产生较大偏差，影响结果的”平滑性“。
请教Matlab的griddata的用法
MATLAB散乱点插值函数
griddata函数
语法：
ZI= griddata(x,y,z,XI,YI)
[XI,YI,ZI]= griddata(x,y,z,XI,YI)
[...]= griddata(...,method)
[...]= griddata(...,method,options)
说明：
ZI= griddata(x,y,z,XI,YI)调整形如z= f(x,y)的曲面，使之与非等间距矢量(x,y,z)中的数据吻合。griddata函数在指定的(XI,YI)点处插补此曲面，生成ZI.此曲面一定通过这些数据点。 XI和 YI通常构成均匀网格(与meshgrid函数生成的相同). XI可以是行矢量,这种情况下该矢量确定一个具有固定列数的矩阵。与之类似，YI可以是列矢量，确定一个具有固定行数的矩阵。
[XI,YI,ZI]= griddata(x,y,z,XI,YI)函数返回与上述矩阵相同的插补后的矩阵ZI。并返回由行矢量XI和列矢量yi形成的矩阵XI和YI.后者与meshgrid函数返回的矩阵相同。
[...]= griddata(...,method)使用规定的插补方法:
'linear'基于三角形的线性插补法(缺省)
'cubic'基于三角形的三次插补法
'nearest'最近邻居插补法
'v4' MATLAB 4 griddata方法。
这些方法定义了匹配数据点的曲面类型。'cubic'和'v4'方法生成平滑曲面，而'linear'和'nearest'分别具有一阶导数和零阶导数不连续。除'v4'外所有方法基于数据的三角化。如果方法为[],则使用缺省的'linear'方法。
[...]= griddata(...,method,options)指定一串将通过delaunayn函数在Qhull中使用的单元阵列选项.如果选项为 [],则使用缺省的三角化选项。如果选项为{''},不使用任何选项，包括缺省选项。有时，griddata函数可能将位于数据凸壳上或靠近凸壳的点返回为NaNs。这是因为有时计算中的圆整处理使得很难确定一个靠近边界的点是否处于凸壳内。
算法：
griddata(...,'v4')命令使用 [3]中文档化的方法。其他griddata方法基于使用Qhull [2]的数据三角化。
实例：
对一个函数在±2.0范围内随机采样100点:
rand('seed',0)
x= rand(100,1)*4-2; y= rand(100,1)*4-2;
z= x.*exp(-x.^2-y.^2);
x, y,和 z均为包含非均匀采样数据的矢量。定义一个规范的网格，将数据与网格匹配:
ti=-2:.25:2;
[XI,YI]= meshgrid(ti,ti);
ZI= griddata(x,y,z,XI,YI);
Plot the gridded data along with the nonuniform data points used to generate it: mesh(XI,YI,ZI), hold
plot3(x,y,z,'o'), hold
matlab如何画指定等值线
clabel可以指定标注等值线的数值。clabel(C,'manual')，是要标注的内容，'manual'就是人工标注。一般结合contour使用，contour(X,Y,Z,Vector)，XYZ为对应的横纵坐标，Z为高度，Vector如果输入向量[1 3 7]，则显示Z等于1 3 7的三圈等值线，如果是一个数6，会将最大值和最小值之间均匀的划分六分，画6条等值线。具体代码如下：
官方解释：clabel(C,'manual') placescontour labels at locations you select with a mouse. Click the mouseor press the space bar to label the contour closest to the centerof the crosshair. Press the Return key while thecursor is within the figure window to terminate labeling.
举例说明：
x=0:0.07:1;y=0:0.07:1;[X1,Y1]=meshgrid(x,y);Z1=0.001./((X1-0.5).^2+(Y1-0.5).^2);
[C,h]=contour(X1,Y1,Z1,[1 3 7]);axis([0.4 0.6 0.4 0.6]);clabel(C,'manual');
结果运行上面代码可以看到手工标注（最近百度知道出问题了，不能上传图片）
感谢您花时间阅读本文！我们希望通过对matlab griddata和请教Matlab的griddata的用法的问题进行探讨，为您提供了一些有用的见解和解决方案。如果您需要更多帮助或者有其他疑问，请不要犹豫与我们联系。