A=load('频率.txt'); * L! ?5 H( a' G! B
x=A(:,1);y = A(:,2);' ~. e4 X$ b4 |3 e+ f% B
z = A(:,3);& i! \# m2 |8 `2 Y
xlin = linspace(min(x),max(x),30);5 p. z: t- n0 P7 l5 M
ylin = linspace(min(y),max(y),30); * F6 {& o0 x7 L8 R& B: x[X,Y] = meshgrid(xlin,ylin);: J+ W4 _. l6 z3 d% |! z9 E
Z=griddata(x,y,z,xlin,ylin,'cubic'); % r- s% ^7 s/ E. L% t4 c; S: ]mesh(X,Y,Z) %interpolated , F: k5 U+ U0 Aaxis tight; hold on 9 d$ p5 X- v8 H2 f. Kplot3(x,y,z,'.','MarkerSize',15) %nonuniform* n8 r5 I5 M9 A3 O+ e% c' J9 ?
5 x6 `2 t6 y. x输入后会出现如下提示: ) G( B* ~4 A% AError in ==> delaunayn at 81+ h7 h1 r! P! U, ?
t = qhullmx(x', 'd ', opt);5 V" ~9 `# H/ s' V
Error in ==> griddata>cubic at 168 9 f. y" p$ P9 ~8 S U) i tri = delaunayn([x(:) y(:)]); $ ~, ^1 S) l: |( FError in ==> griddata at 924 W# l# G% `; C
zi = cubic(x,y,z,xi,yi,opt); 4 y& D# B2 j) `/ ]9 pError in ==> aa at 72 k+ H% Z' t& ^& t; B
Z=griddata(x,y,z,xlin,ylin,'cubic'); * \0 l- Q6 v' v& f' x: |9 B8 z