Mercurial > octave-nkf

comparison scripts/general/interp2.m @ 6448:2110cc251779

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2007-03-24 02:47:36 by jwe]
author jwe
date Sat, 24 Mar 2007 02:47:36 +0000
parents 7fad1fad19e1
children b2391d403ed2
comparison
equal deleted inserted replaced
6447:3f79532415b5 6448:2110cc251779
237 %! x=[0,1,4]; y=[10,11,12]; 237 %! x=[0,1,4]; y=[10,11,12];
238 %! xi=linspace(min(x),max(x),17); 238 %! xi=linspace(min(x),max(x),17);
239 %! yi=linspace(min(y),max(y),26)'; 239 %! yi=linspace(min(y),max(y),26)';
240 %! mesh(xi,yi,interp2(x,y,A,xi,yi,'linear')); 240 %! mesh(xi,yi,interp2(x,y,A,xi,yi,'linear'));
241 %! [x,y] = meshgrid(x,y); 241 %! [x,y] = meshgrid(x,y);
242 %! __gnuplot_raw__ ("set nohidden3d;\n")
243 %! hold on; plot3(x(:),y(:),A(:),"b*"); hold off; 242 %! hold on; plot3(x(:),y(:),A(:),"b*"); hold off;
244 243
245 %!demo 244 %!demo
246 %! A=[13,-1,12;5,4,3;1,6,2]; 245 %! A=[13,-1,12;5,4,3;1,6,2];
247 %! x=[0,1,4]; y=[10,11,12]; 246 %! x=[0,1,4]; y=[10,11,12];
248 %! xi=linspace(min(x),max(x),17); 247 %! xi=linspace(min(x),max(x),17);
249 %! yi=linspace(min(y),max(y),26)'; 248 %! yi=linspace(min(y),max(y),26)';
250 %! mesh(xi,yi,interp2(x,y,A,xi,yi,'nearest')); 249 %! mesh(xi,yi,interp2(x,y,A,xi,yi,'nearest'));
251 %! [x,y] = meshgrid(x,y); 250 %! [x,y] = meshgrid(x,y);
252 %! __gnuplot_raw__ ("set nohidden3d;\n")
253 %! hold on; plot3(x(:),y(:),A(:),"b*"); hold off; 251 %! hold on; plot3(x(:),y(:),A(:),"b*"); hold off;
254 252
255 %!#demo 253 %!#demo
256 %! A=[13,-1,12;5,4,3;1,6,2]; 254 %! A=[13,-1,12;5,4,3;1,6,2];
257 %! x=[0,1,2]; y=[10,11,12]; 255 %! x=[0,1,2]; y=[10,11,12];
258 %! xi=linspace(min(x),max(x),17); 256 %! xi=linspace(min(x),max(x),17);
259 %! yi=linspace(min(y),max(y),26); 257 %! yi=linspace(min(y),max(y),26);
260 %! mesh(xi,yi,interp2(x,y,A,xi,yi,'cubic')); 258 %! mesh(xi,yi,interp2(x,y,A,xi,yi,'cubic'));
261 %! [x,y] = meshgrid(x,y); 259 %! [x,y] = meshgrid(x,y);
262 %! __gnuplot_raw__ ("set nohidden3d;\n")
263 %! hold on; plot3(x(:),y(:),A(:),"b*"); hold off; 260 %! hold on; plot3(x(:),y(:),A(:),"b*"); hold off;
264 261
265 %!test % simple test 262 %!test % simple test
266 %! x = [1,2,3]; 263 %! x = [1,2,3];
267 %! y = [4,5,6,7]; 264 %! y = [4,5,6,7];