Matlab, Cross Produkt Multi Dim Matrix

Question

Ich habe u1 v1 w1 (Geschwindigkeitskomponenten) volumetrische Daten. Die Größe von u1 ist NxNxN (N ist 100) und dasselbe gilt für v1 und w1.

Ich habe noch einen Fall von u2, v2, w2, jeder der Größe NxNxN. Beide Datensätze stammen aus TriScatteredInterp und Meshgrid. Grundsätzlich habe ich Stromlinien für diese beiden Fälle gezeichnet und möchte nun sehen, wie groß die Abweichung zwischen diesen Stromlinien ist (von Fall 1 zu 2). Ich muss den Winkel zwischen ihnen an den durch meshgrid gegebenen Punkten finden. Aber Cross-Produkt wird nicht funktionieren (A und B müssen mindestens eine Dimension der Länge 3 haben). Vielen Dank

Answer 1

der Winkel gefunden

, lassen sich mit

angle = atan2(norm(cross(a,b)),dot(a,b)) 

[mit Vier-Quadranten-Arkustangens]; aber für meinen Fall funktioniert der Kreuzproduktteil nicht als die Größe jeder Geschwindigkeitskomponente in 100x100x100.

Mithilfe der Concatenate-Funktion von Matlab können die drei Komponenten Velocity1 und Velocity2 kombiniert werden, um eine Matrix mit der Größe NxNxNx3 zu erhalten. Nun kann für jedes Tripel/Kombination von V1 (u1, v1, w1) und V2 (u2, v2, w2) ein Kreuzprodukt (nennen wir es C) in der vierten Dimension genommen werden, was zu einer Matrix mit der Größe NxNxNx3 führt. Um in der oben genannten Formel die Größe/Norm des Ergebnisses des Kreuzprodukts zu erhalten, müssen wir das Skalarprodukt von C mit C selbst in der 4. Dimension und der Quadratwurzel nehmen, um die Größe zu erhalten und seine Größe wieder zu NxNxN zu machen. so schließlich kann es like--

A = cat(4,u1,v1,w1); % Combine the three components in the 4th dimension 
    B = cat(4,u2,v2,w2); % Ditto 
    C = cross(A,B,4); % Take the cross products there. 
    ang = atan2(sqrt(dot(C,C,4)),dot(A,B,4)); 

Diese Lösung wird von Herrn Roger Stafford (Matlab-Website) und alle Kredite geht an ihn gegeben erfolgen.

Antwort finden Sie hier .. gefunden http://fr.mathworks.com/matlabcentral/answers/281635-cross-product-multi-dim-data