数学建模社区-数学中国
标题:
这个插值代码到底哪里错了嘛
[打印本页]
作者:
gaoyuan_1993
时间:
2013-8-26 10:51
标题:
这个插值代码到底哪里错了嘛
A=load('频率.txt');
* F% j0 `/ d6 f4 m/ n2 W7 s, n* l
x=A(:,1);y = A(:,2);
' L. q6 e7 A6 b# ]
z = A(:,3);
2 {9 L5 M( s! o, R
xlin = linspace(min(x),max(x),30);
7 [1 J5 x% d( u7 q. O8 h2 e
ylin = linspace(min(y),max(y),30);
* |: H& O3 i1 K D& Z8 l: w) P
[X,Y] = meshgrid(xlin,ylin);
! H: K+ a* [3 W2 c) v- u3 r
Z=griddata(x,y,z,xlin,ylin,'cubic');
* P" J+ z& @+ |2 U3 V2 {0 D/ p
mesh(X,Y,Z) %interpolated
; J% J" D) W6 r' b) [
axis tight; hold on
' E3 l3 ]! x7 @" G- Q+ M
plot3(x,y,z,'.','MarkerSize',15) %nonuniform
8 ^ t0 K& x% T
4 J9 T9 }5 L% U* [5 k$ h
输入后会出现如下提示:
# p4 ?) p8 o0 k/ [, l" b
Error in ==> delaunayn at 81
9 Z5 l" i5 E* M1 M
t = qhullmx(x', 'd ', opt);
' Y! ]& J; [5 x) e! b1 S1 ~
Error in ==> griddata>cubic at 168
4 q: W. y7 z: H) A0 O
tri = delaunayn([x(:) y(:)]);
, \1 n$ x8 U: v, Z6 C$ Y8 q
Error in ==> griddata at 92
" F4 Z+ V- q& {' O7 v$ P8 }
zi = cubic(x,y,z,xi,yi,opt);
1 p+ w+ [0 Q/ i
Error in ==> aa at 7
* h8 h3 h) o+ @% P0 t* z4 j0 }: t
Z=griddata(x,y,z,xlin,ylin,'cubic');
$ M, e1 `0 n. |' F
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5