¿Cómo utilizar MATLAB para realizar ajustes de mínimos cuadrados de funciones multivariadas?

Código de referencia

% Generar datos de prueba

a0=1;a1=2;a2=3;

x1=rand( 10,10);

x2=rand(10,10);

Y=aa1*x1+a2*x2 )

% Agregar aleatorio Ruido

Y=Y+0.1*randn(size(y));

% Forma de función de ajuste

f=@(k,x)k ( 1)+k(2)*x(:,1)+k(3)*x(:,2);

x = [x1(:) x2(:)];

k=lsqcurvefit(f,[1 1 1 1],x,Y(:))

n = tamaño(x,1);

plot(1:n, Y (:),'bo',1:n,f(k,x),'r:.')

Resultados en ejecución

k =

1.0309 1.9782 2.9800

corresponden a a0~a2 respectivamente. Se puede ver que son más consistentes con los coeficientes originales.