Creazione di una matrice 3D (Matlab)

[Ianero]

Ho 3 vettori di lunghezze diverse:







e vorrei costruire tre distinte matrici, tutte della stessa dimensione mxnxl (3D), fatte così (per ogni posizione, riporto 3 elementi, uno per ogni matrice che voglio ottenere):

primo strato


secondo strato


...

ultimo strato


Vorrei ovviamente evitare di crearla con cicli for o cose simili, che in Matlab sono quasi sempre di troppo.

Esiste qualche scorciatoia che me lo fa fare in un paio di righe?

[RLC]

Ciao
Sei sicuro che le componenti delle varie matrici siano giuste? perché trovare delle correlazioni/serie dei pedici (trovarle servirà per automatizzare il codice) sembra un po' un rompicapo in stile "test QI"

Tuttavia, comincerei guardando qui:
https://it.mathworks.com/help/matlab/math/multidimensional-arrays.html

[Ianero]

Non è che esista un giusto.
A me servono così, chiedevo se ci fosse un modo standard di farlo con le potenzialità di Matlab o se in effetti bisogna ricorrere al for.

[MarkyMark]

Il comando repmat è utile in questo caso

Codice: Seleziona tutto

theta = (1:2)';
v = (3:5)';
k = (6:9)';
l = length(theta);
m = length(v);
n = length(k);

A1(1:m,1:n,1:l) = theta(1);
A1(:,1,:) = repmat(theta',m,1);

A2(1:m,1:n) = repmat(v,1,n);
A2(:,:,1:l) = repmat(A2,1,1,l);

A3(1:m,1:n) = repmat(k',m,1);
A3(:,:,1:l) = repmat(A3,1,1,l);

Spero di aver interpretato correttamente la tua richiesta.

[Ianero]

A1(:,1,:) = repmat(theta',m,1);

Questa mi mancava, grande.
Non sapevo che Matlab capisse da solo di dover distribuire le colonne.
Ottimo, grazie mille.

[Ianero]

No però aspetta c'è un problema con A1, le altre due vanno bene.
La matrice è corretta solo nella prima colonna di ogni strato. Ora bisognerebbe dirgli di replicare la matrice A1(:,1,:) lungo la seconda dimensione per n volte.

Dovrebbe andare così:

Codice: Seleziona tutto

A1(1:m,1,1:l) = repmat(theta',m,1);
A1(:,1:n,:)=repmat(A1,1,n,1);

Concordi?

EDIT: avevo indicato male in [1] la richiesta sulla matrice dei vari , colpa mia, ho corretto adesso.

[MarkyMark]

Sì, hai ragione. Ho sbagliato a leggere il primo post

Sorgente FidoCadJ | Ingrandisci | Cos'è Fidocad

Come hai detto, si può costruire prima uno degli strati mxl e poi lo si replica n volte

[Ianero]

Grazie MarkyMark

[Ianero]

Ho appena scoperto che si può fare tutto in una sola riga:

Codice: Seleziona tutto

theta = (1:2)';
v = (3:5)';
k = (6:9)';

[A1,A2,A3] = ndgrid(v,k,theta);


Quant'è potente questo linguaggio.

[MarkyMark]

Fantastico! Grazie