数学建模社区-数学中国

标题: 这个插值代码到底哪里错了嘛 [打印本页]

作者: gaoyuan_1993    时间: 2013-8-26 10:51
标题: 这个插值代码到底哪里错了嘛
A=load('频率.txt');
  \/ W8 {8 v9 C+ mx=A(:,1);y = A(:,2);, h. h" \" d9 z0 X9 W8 M
z = A(:,3);4 F+ J. V; O! B2 m2 h( F& D/ K
xlin = linspace(min(x),max(x),30);( O- M, ~# ?7 x8 L9 `- C' U+ {( z8 l7 S6 u
ylin = linspace(min(y),max(y),30);* N4 e  m7 D8 z0 ]
[X,Y] = meshgrid(xlin,ylin);
( k4 P: n. C2 B1 F( g% R& {Z=griddata(x,y,z,xlin,ylin,'cubic'); 0 p0 Y1 J( p" q; c
mesh(X,Y,Z) %interpolated
) @2 R. A' T- D) t# {3 P# k. maxis tight; hold on( q5 S% P1 P# k% l
plot3(x,y,z,'.','MarkerSize',15) %nonuniform
0 d! x% L  `+ M  Y* {2 Q, w
, E7 q" ]$ o! |6 f- g0 D) D输入后会出现如下提示:! y( w7 c& c( U! Z
Error in ==> delaunayn at 81+ h9 S, A0 u7 C- r, R% y
t = qhullmx(x', 'd ', opt);

+ O6 Z2 \! P8 `+ @& oError in ==> griddata>cubic at 168
! f  a1 T7 Q8 u+ Y7 a* ?$ R& L    tri = delaunayn([x(:) y(:)]);
  {. X; r/ P5 S1 b
Error in ==> griddata at 92
* j1 t( v6 [* t7 q/ D3 E" x    zi = cubic(x,y,z,xi,yi,opt);

& h/ z% N/ O* K" b$ t! rError in ==> aa at 7
: p2 H6 N/ T$ O) B. ^4 ]: N: m$ _7 JZ=griddata(x,y,z,xlin,ylin,'cubic');
9 E0 m( k4 S/ W  u





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