数学建模社区-数学中国
标题:
这个插值代码到底哪里错了嘛
[打印本页]
作者:
gaoyuan_1993
时间:
2013-8-26 10:51
标题:
这个插值代码到底哪里错了嘛
A=load('频率.txt');
\/ W8 {8 v9 C+ m
x=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. m
axis 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 `+ @& o
Error 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! r
Error in ==> aa at 7
: p2 H6 N/ T$ O) B. ^4 ]: N: m$ _7 J
Z=griddata(x,y,z,xlin,ylin,'cubic');
9 E0 m( k4 S/ W u
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5