[Home] [Lexikon] [Aufgaben] [Tests] [Kurse] [Begleitmaterial] [Hinweise] [Mitwirkende] [Publikationen]
	Mathematik-Online-Aufgabensammlung: Lösung zu
	Interaktive Aufgabe 1093: Programm zur linearen Interpolation

function B = refine(A)

% function B = refine(A)
% refinement of matrix via linear interpolation 


[m,n] = size(A);
% expand grid
B(1:2:2*m-1,1:2:2*n-1) = A;

% interpolating columns
ind = [2:2:2*m-1];
B(ind,:) = (B(ind-1,:) + B(ind+1,:))/2;

% interpolating rows
ind = [2:2:2*n-1];
B(:,ind) = (B(:,ind-1) + B(:,ind+1))/2;

>> sum(sum(refine(magic(5))))

ans =

        1053