|
6721
|
1 function interpimages (nm, typ) |
|
|
2 bury_output (); |
|
|
3 if (strcmp (nm, "interpft")) |
|
|
4 t = 0 : 0.3 : pi; dt = t(2)-t(1); |
|
|
5 n = length (t); k = 100; |
|
|
6 ti = t(1) + [0 : k-1]*dt*n/k; |
|
|
7 y = sin (4*t + 0.3) .* cos (3*t - 0.1); |
|
|
8 yp = sin (4*ti + 0.3) .* cos (3*ti - 0.1); |
|
|
9 plot (ti, yp, 'g', ti, interp1(t, y, ti, 'spline'), 'b', ... |
|
|
10 ti, interpft (y, k), 'c', t, y, 'r+'); |
|
|
11 legend ('sin(4t+0.3)cos(3t-0.1','spline','interpft','data'); |
|
|
12 print (strcat (nm, ".", typ), strcat ("-d", typ)) |
|
|
13 elseif (strcmp (nm, "interpn")) |
|
|
14 x = y = z = -1:1; |
|
|
15 f = @(x,y,z) x.^2 - y - z.^2; |
|
|
16 [xx, yy, zz] = meshgrid (x, y, z); |
|
|
17 v = f (xx,yy,zz); |
|
|
18 xi = yi = zi = -1:0.1:1; |
|
|
19 [xxi, yyi, zzi] = ndgrid (xi, yi, zi); |
|
|
20 vi = interpn(x, y, z, v, xxi, yyi, zzi, 'spline'); |
|
6723
|
21 mesh (zi, yi, squeeze (vi(1,:,:))); |
|
6721
|
22 print (strcat (nm, ".", typ), strcat ("-d", typ)) |
|
6743
|
23 elseif (strcmp (nm, "interpderiv1")) |
|
|
24 t = -2:2; |
|
|
25 dt = 1; |
|
|
26 ti =-2:0.025:2; |
|
|
27 dti = 0.025; |
|
|
28 y = sign(t); |
|
|
29 ys = interp1(t,y,ti,'spline'); |
|
|
30 yp = interp1(t,y,ti,'pchip'); |
|
|
31 plot (ti, ys,'r-', ti, yp,'g-'); |
|
|
32 legend('spline','pchip', 4); |
|
|
33 print (strcat (nm, ".", typ), strcat ("-d", typ)) |
|
|
34 elseif (strcmp (nm, "interpderiv2")) |
|
|
35 t = -2:2; |
|
|
36 dt = 1; |
|
|
37 ti =-2:0.025:2; |
|
|
38 dti = 0.025; |
|
|
39 y = sign(t); |
|
6721
|
40 ddys = diff(diff(interp1(t,y,ti,'spline'))./dti)./dti; |
|
|
41 ddyp = diff(diff(interp1(t,y,ti,'pchip'))./dti)./dti; |
|
6743
|
42 plot (ti(2:end-1),ddys,'r*', ti(2:end-1),ddyp,'g+'); |
|
|
43 legend('spline','pchip'); |
|
6721
|
44 print (strcat (nm, ".", typ), strcat ("-d", typ)) |
|
|
45 endif |
|
|
46 bury_output (); |
|
|
47 endfunction |
|
|
48 |
|
|
49 ## Use this function before plotting commands and after every call to |
|
|
50 ## print since print() resets output to stdout (unfortunately, gnpulot |
|
|
51 ## can't pop output as it can the terminal type). |
|
|
52 function bury_output () |
|
|
53 f = figure (1); |
|
|
54 set (f, "visible", "off"); |
|
|
55 endfunction |
