¿Cómo puedo indexar una matriz MATLAB devuelta por una función sin asignarla primero a una variable local?

Por ejemplo, si quiero leer el valor medio de magic(5), puedo hacerlo así:

M = magic(5);
value = M(3,3);

para llegar value == 13. Me gustaría poder hacer algo como uno de estos:

value = magic(5)(3,3);
value = (magic(5))(3,3);

prescindir de la variable intermedia. Sin embargo, MATLAB se queja sobre Unbalanced or unexpected parenthesis or bracketel primer paréntesis antes del 3.

¿Es posible leer valores de una matriz / matriz sin asignarlo primero a una variable?

En realidad, es posible hacer lo que quiera, pero debe usar la forma funcional del operador de indexación. Cuando realiza una operación de indexación utilizando (), en realidad está haciendo una llamada a la subsreffunción. Entonces, aunque no puedas hacer esto:

value = magic(5)(3, 3);

Usted puede hacer esto:

value = subsref(magic(5), struct('type', '()', 'subs', {{3, 3}}));

Feo, pero posible. ;)

En general, solo tiene que cambiar el paso de indexación a una llamada de función para que no tenga dos conjuntos de paréntesis inmediatamente seguidos. Otra forma de hacerlo sería definir su propia función anónima para hacer la indexación con subíndice. Por ejemplo:

subindex = @(A, r, c) A(r, c);     % An anonymous function for 2-D indexing
value = subindex(magic(5), 3, 3);  % Use the function to index the matrix

Sin embargo, cuando todo está dicho y hecho, la solución de variable local temporal es mucho más legible, y definitivamente lo que sugeriría.

Hubo una buena publicación de blog en Loren sobre el Arte de Matlab hace un par de días con un par de gemas que podrían ayudar. En particular, el uso de funciones auxiliares como:

paren = @(x, varargin) x(varargin{:});
curly = @(x, varargin) x{varargin{:}};

donde paren()se puede usar como

paren(magic(5), 3, 3);

volvería

ans = 16

También supongo que esto será más rápido que la respuesta de gnovice, pero no lo he verificado (¡¡¡Usa el generador de perfiles !!!). Dicho esto, también debe incluir estas definiciones de funciones en alguna parte. Personalmente, los he convertido en funciones independientes en mi camino, porque son súper útiles.

Estas funciones y otras ahora están disponibles en el complemento Construcciones de programación funcional que está disponible a través del Explorador de complementos de MATLAB o en el Intercambio de archivos .

¿Cómo te sientes al usar características no documentadas?

>> builtin('_paren', magic(5), 3, 3)               %# M(3,3)
ans =
    13

o para matrices de células:

>> builtin('_brace', num2cell(magic(5)), 3, 3)     %# C{3,3}
ans =
    13

Justo como magia :)

ACTUALIZAR:

Malas noticias, ¡el truco anterior ya no funciona en R2015b ! Está bien, era una funcionalidad no documentada y no podemos confiar en ella como una característica compatible :)

Para aquellos que se preguntan dónde encontrar este tipo de cosas, busquen en la carpeta fullfile(matlabroot,'bin','registry'). Hay un montón de archivos XML que enumeran todo tipo de golosinas. Tenga en cuenta que llamar a algunas de estas funciones directamente puede bloquear fácilmente su sesión de MATLAB.

Al menos en MATLAB 2013a puedes usar getfieldcomo:

a=rand(5);
getfield(a,{1,2}) % etc

para obtener el elemento en (1,2)

desafortunadamente la sintaxis como magic(5)(3,3)no es compatible con matlab. necesita usar variables intermedias temporales. puede liberar la memoria después de su uso, p. ej.

tmp = magic(3);
myVar = tmp(3,3);
clear tmp

Tenga en cuenta que si compara los tiempos de ejecución con la forma estándar (asigne el resultado y luego acceda a las entradas), son exactamente iguales.

subs=@(M,i,j) M(i,j);
>> for nit=1:10;tic;subs(magic(100),1:10,1:10);tlap(nit)=toc;end;mean(tlap)

ans =

0.0103

>> for nit=1:10,tic;M=magic(100); M(1:10,1:10);tlap(nit)=toc;end;mean(tlap)

ans =

0.0101

En mi opinión, la conclusión es: MATLAB no tiene punteros, hay que vivir con ellos.

Podría ser más simple si realiza una nueva función:

function [ element ] = getElem( matrix, index1, index2 )
    element = matrix(index1, index2);
end

y luego úsalo:

value = getElem(magic(5), 3, 3);

Su notación inicial es la forma más concisa de hacer esto:

M = magic(5);  %create
value = M(3,3);  % extract useful data
clear M;  %free memory

Si está haciendo esto en un bucle, puede reasignar M cada vez e ignorar la declaración clara también.

Para complementar la respuesta de Amro, puede usar en fevallugar de builtin. Realmente no hay diferencia, a menos que intentes sobrecargar la función del operador:

BUILTIN (...) es lo mismo que FEVAL (...) excepto que llamará a la versión incorporada original de la función incluso si existe una sobrecargada (para que esto funcione, nunca debe sobrecargar BUILTIN).

>> feval('_paren', magic(5), 3, 3)               % M(3,3)
ans =
    13

>> feval('_brace', num2cell(magic(5)), 3, 3)     % C{3,3}
ans =
    13

Lo interesante es que fevalparece ser un poco más rápido que builtin(en ~ 3.5%), al menos en Matlab 2013b, lo cual es extraño dado que fevalnecesita verificar si la función está sobrecargada, a diferencia de builtin:

>> tic; for i=1:1e6, feval('_paren', magic(5), 3, 3); end; toc;
Elapsed time is 49.904117 seconds.
>> tic; for i=1:1e6, builtin('_paren', magic(5), 3, 3); end; toc;
Elapsed time is 51.485339 seconds.