数学建模社区-数学中国

标题: 这个插值代码到底哪里错了嘛 [打印本页]
作者: gaoyuan_1993    时间: 2013-8-26 10:51
标题: 这个插值代码到底哪里错了嘛
A=load('频率.txt');
  _' I7 S/ a( K% \: Y0 I! @" px=A(:,1);y = A(:,2);$ j$ e* O5 r8 T. J
z = A(:,3);  |# [; e  M0 h6 `. p* k! I: ?
xlin = linspace(min(x),max(x),30);
  s7 W0 v4 ?7 Hylin = linspace(min(y),max(y),30);( Q' q6 ]+ y7 Q: ]; O/ f& l
[X,Y] = meshgrid(xlin,ylin);
) h3 d' `6 P# f# F. xZ=griddata(x,y,z,xlin,ylin,'cubic');
3 ~  }+ o( G# Z2 }4 Fmesh(X,Y,Z) %interpolated& m; L# E) o6 Y
axis tight; hold on% e- {( ]5 s" h
plot3(x,y,z,'.','MarkerSize',15) %nonuniform. r) C8 ^6 h* f# G; s* N& C
) I/ A0 }! `; ~, w1 W5 M# n
输入后会出现如下提示：  j! v; x0 S5 Q- ^
Error in ==> delaunayn at 81
4 {- x. ], ~0 ?0 N1 At = qhullmx(x', 'd ', opt);6 D& b4 |3 B# y4 k
Error in ==> griddata>cubic at 168
# d( u# W. n! `+ m  }$ P( F$ G+ z    tri = delaunayn([x(:) y(:)]);
/ m! b: \9 d5 m8 i% Q# fError in ==> griddata at 92
. @* K$ A  A  G! ^" s5 p) s- x    zi = cubic(x,y,z,xi,yi,opt);
4 @: K2 \2 L3 u4 V/ LError in ==> aa at 7
0 A5 W0 Y; ]6 @' [Z=griddata(x,y,z,xlin,ylin,'cubic');9 w  H9 {7 q: Q9 ]% f




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5