数学建模社区-数学中国
标题:
这个插值代码到底哪里错了嘛
[打印本页]
作者:
gaoyuan_1993
时间:
2013-8-26 10:51
标题:
这个插值代码到底哪里错了嘛
A=load('频率.txt');
- W" z2 G# p: K) S3 X6 f
x=A(:,1);y = A(:,2);
4 @" Z" Y/ W/ J) Q
z = A(:,3);
7 W) G& L* E! @, @
xlin = linspace(min(x),max(x),30);
. Y0 J# g. |4 R$ e
ylin = linspace(min(y),max(y),30);
) P6 N# I2 i* w6 z
[X,Y] = meshgrid(xlin,ylin);
, N6 }1 Y) l0 a Q' ]( a. ^- h; E
Z=griddata(x,y,z,xlin,ylin,'cubic');
0 J, ^1 ^+ T4 j/ v* [! e
mesh(X,Y,Z) %interpolated
% |8 K0 b3 ]' {4 O0 U7 P
axis tight; hold on
1 q3 t$ {+ }0 o5 K9 B0 F) L! [, U
plot3(x,y,z,'.','MarkerSize',15) %nonuniform
`7 [" j/ B n$ X7 M! v1 P
! V" i1 e }! r$ D
输入后会出现如下提示：
1 t/ o$ W" e6 V& U
Error in ==> delaunayn at 81
% H5 k0 Q% N' ^* ]4 D. q
t = qhullmx(x', 'd ', opt);
4 H8 m- k8 I3 d+ I
Error in ==> griddata>cubic at 168
% h t5 Z0 Z7 A
tri = delaunayn([x(:) y(:)]);
) P( ^) r7 h" ^8 F/ o5 Q+ e \
Error in ==> griddata at 92
8 w: r, G- g5 J7 k% \
zi = cubic(x,y,z,xi,yi,opt);
1 C$ B; @$ \; ~ N k7 x7 h8 X
Error in ==> aa at 7
9 e- n) `9 \( d, D9 }
Z=griddata(x,y,z,xlin,ylin,'cubic');
* @ G% L* r6 P/ v9 m# F! C
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5