A=load('频率.txt'); ) G, [ |2 r9 ] G; G. T# tx=A(:,1);y = A(:,2);3 L: ?$ N. m% O- F! J \
z = A(:,3); [; U# }0 }1 f9 P4 Y( S
xlin = linspace(min(x),max(x),30); + e2 P5 g* V* R8 N1 s. P& K( Cylin = linspace(min(y),max(y),30);) D5 w& c2 z! t) [. a) [
[X,Y] = meshgrid(xlin,ylin); 3 F0 r, E7 @+ f6 N4 KZ=griddata(x,y,z,xlin,ylin,'cubic'); ) \6 D3 |1 \4 l, ]9 k. U5 ]+ d
mesh(X,Y,Z) %interpolated 8 A; S% p& E3 T q. M Faxis tight; hold on/ T$ X H# ~+ ?4 c
plot3(x,y,z,'.','MarkerSize',15) %nonuniform2 `( r7 y5 e8 I! r. v) D# ~
0 |0 c2 L2 D8 j; p6 `2 B输入后会出现如下提示:- F$ x& I3 O- g6 g- g/ v
Error in ==> delaunayn at 81 * @6 q/ `7 ~7 V! mt = qhullmx(x', 'd ', opt); 9 g9 I8 d9 L9 Z+ w. zError in ==> griddata>cubic at 168- t9 k$ h9 f" w7 F. d. D
tri = delaunayn([x(:) y(:)]);7 s' q7 g$ C1 G$ H# p" I
Error in ==> griddata at 92. I9 o, j# K& {, g# B
zi = cubic(x,y,z,xi,yi,opt); ( b6 @) k; q" J2 h/ v3 mError in ==> aa at 7/ k8 ]6 Y" r* t: X/ {
Z=griddata(x,y,z,xlin,ylin,'cubic');# \! `& C/ S' d# C. Q