Port-R-Expression in MATLAB

Question

Wenn ich eine NXP Matrix G und und PX1 Vektor Lambda ich dies in R tun können

G = matrix(c(1/2,1/4,3/4,5/2),nrow=2,ncol=2,byrow=TRUE) 
Lambda = c(1,2) 
M = t(G)%*%(c((1/(1+G%*%Lambda)^(2)))*G) 
M 
      [,1]  [,2] 
[1,] -0.07484568 -0.07240226 
[2,] -0.07240226 -0.15279921 

Und ich bekomme eine PxP Matrix, das ist, was ich will. Aber in Matlab, ist dies das Beste, was ich Figur könnte aus so etwas wie

b=1./((1+G*Lambda).^2) 
transpose(G(i,:))*G(i,:)*b(i) 

tut und Schleifen über diese Idee. Aber ich möchte eine Simulationsstudie machen und würde einen vektorisierten Weg bevorzugen. Wenn mir jemand helfen könnte, das herauszufinden, wäre ich sehr glücklich.

Answer 1

Ich glaube, dass Matlab Vektorelemente nicht wiederverwenden kann, wenn Dimensionen nicht übereinstimmen, daher ist die Lösung etwas komplizierter.

n = 2; 
p = 4; 

G = reshape(1:(n*p), [n, p]); 
Lambda = (1:p).'; 

A = (1./(1+G*Lambda).^(2)); 
% since Matlab cannot reuse the vector elements if dimensions mismatch as R 
B = reshape(repmat(A,p,1).*G(:), size(G)); 
G.'*B 

Sie können auch einen Einstrich aus dem obigen Code erstellen. Ich finde immer noch den Weg, wie Matlab Operationen intuitiver behandelt, im Einklang mit der Mathematik.

EDIT: basierend auf @ Dan Vorschlägen:

bessere Lösung B zu erhalten:

B = bsxfun(@times, A, G);