Jak mogę zindeksować tablicę MATLAB zwróconą przez funkcję bez uprzedniego przypisania jej do zmiennej lokalnej?

Na przykład, jeśli chcę odczytać środkową wartość magic(5), mogę to zrobić w następujący sposób:

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

dostać value == 13. Chciałbym móc zrobić coś takiego:

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

zrezygnować ze zmiennej pośredniej. Jednak MATLAB narzeka Unbalanced or unexpected parenthesis or bracketna pierwszy nawias przed 3.

Czy możliwe jest odczytanie wartości z tablicy / macierzy bez uprzedniego przypisania jej do zmiennej?

Możliwe jest robienie tego, co chcesz, ale musisz użyć funkcjonalnej formy operatora indeksowania. Podczas wykonywania operacji indeksowania za pomocą ()faktycznie wywołujesz subsreffunkcję. Więc nawet jeśli nie możesz tego zrobić:

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

Możesz to zrobić:

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

Brzydkie, ale możliwe. ;)

Ogólnie rzecz biorąc, wystarczy zmienić krok indeksowania na wywołanie funkcji, aby nie było dwóch zestawów nawiasów następujących po sobie. Innym sposobem na to byłoby zdefiniowanie własnej anonimowej funkcji do indeksowania w indeksie dolnym. Na przykład:

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

Jednak kiedy wszystko zostanie powiedziane i zrobione, rozwiązanie tymczasowej zmiennej lokalnej jest znacznie bardziej czytelne i zdecydowanie to, co sugerowałbym.

Kilka dni temu był tylko dobry post na blogu o Lorenie na temat Art of Matlab z kilkoma klejnotami, które mogą pomóc. W szczególności za pomocą funkcji pomocniczych, takich jak:

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

gdzie paren()można użyć jak

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

wróciłby

ans = 16

Przypuszczam również, że będzie to szybsze niż odpowiedź gnovice, ale nie sprawdziłem (użyj profilera !!!). To powiedziawszy, musisz gdzieś dołączyć te definicje funkcji. Osobiście uczyniłem je niezależnymi funkcjami na mojej ścieżce, ponieważ są one bardzo przydatne.

Te i inne funkcje są teraz dostępne w dodatku Functional Programming Constructs, który jest dostępny za pośrednictwem Eksploratora dodatków MATLAB lub wymiany plików .

Co sądzisz o korzystaniu z nieudokumentowanych funkcji:

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

lub w przypadku tablic komórkowych:

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

Tak jak magia :)

AKTUALIZACJA:

Złe wieści, powyższy hack już nie działa w R2015b ! W porządku, to była nieudokumentowana funkcjonalność i nie możemy polegać na niej jako obsługiwanej funkcji :)

Jeśli zastanawiasz się, gdzie znaleźć tego rodzaju rzeczy, zajrzyj do folderu fullfile(matlabroot,'bin','registry'). Jest tam kilka plików XML z listą różnego rodzaju gadżetów. Ostrzegamy, że bezpośrednie wywołanie niektórych z tych funkcji może łatwo zawiesić sesję MATLAB.

Przynajmniej w MATLAB 2013a możesz używać getfieldtakich jak:

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

uzyskać element w (1,2)

niestety składnia like magic(5)(3,3)nie jest obsługiwana przez matlab. musisz użyć tymczasowych zmiennych pośrednich. możesz zwolnić pamięć po użyciu, np

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

Zauważ, że jeśli porównasz czasy działania ze standardowym sposobem (przypisz wynik, a następnie uzyskaj dostęp do wpisów), będą one dokładnie takie same.

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

Moim zdaniem, podstawową kwestią jest: MATLAB nie ma wskaźników, musisz z tym żyć.

Może być prostsze, jeśli utworzysz nową funkcję:

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

a następnie użyj go:

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

Twoja początkowa notacja jest najbardziej zwięzłym sposobem:

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

Jeśli robisz to w pętli, możesz po prostu ponownie przypisać M za każdym razem i zignorować również jasne polecenie.

Aby uzupełnić odpowiedź Amro, możesz użyć fevalzamiast builtin. Tak naprawdę nie ma różnicy, chyba że spróbujesz przeciążyć funkcję operatora:

BUILTIN (...) jest taki sam jak FEVAL (...) poza tym, że wywoła oryginalną wbudowaną wersję funkcji, nawet jeśli istnieje przeciążona (aby działało, nigdy nie należy przeciążać BUILTIN).

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

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

Interesujące jest to, że fevalwydaje się to nieco szybsze niż builtin(o ~ 3,5%), przynajmniej w Matlabie 2013b, co jest dziwne, biorąc pod uwagę, że fevalnależy sprawdzić, czy funkcja jest przeciążona, w przeciwieństwie do 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.