ElectroYou

Ho 3 vettori di lunghezze diverse:

$\underline{\theta}=\begin{bmatrix} \theta_1\\ \theta_2\\... \\ \theta_l \end{bmatrix}$

$\underline{v}=\begin{bmatrix} v_1\\ v_2\\... \\ v_m \end{bmatrix}$

$\underline{k}=\begin{bmatrix} k_1\\ k_2\\... \\ k_n \end{bmatrix}$

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
$\begin{bmatrix} \theta_1/v_1/k_1 & \theta_1/v_1/k_2 & ... & \theta_1/v_1/k_n \\ \theta_1/v_2/k_1 & \theta_1/v_2/k_2 & ... & \theta_1/v_2/k_n \\ &...&\\ \theta_1/v_m/k_1 & \theta_1/v_m/k_2 & ... & \theta_1/v_m/k_n \end{bmatrix}$

secondo strato
$\begin{bmatrix} \theta_2/v_1/k_1 & \theta_2/v_1/k_2 & ... & \theta_2/v_1/k_n \\ \theta_2/v_2/k_1 & \theta_2/v_2/k_2 & ... & \theta_2/v_2/k_n \\ &...&\\ \theta_2/v_m/k_1 & \theta_2/v_m/k_2 & ... & \theta_2/v_m/k_n \end{bmatrix}$

...

ultimo strato
$\begin{bmatrix} \theta_l/v_1/k_1 & \theta_l/v_1/k_2 & ... & \theta_l/v_1/k_n \\ \theta_l/v_2/k_1 & \theta_l/v_2/k_2 & ... & \theta_l/v_2/k_n \\ &...&\\ \theta_l/v_m/k_1 & \theta_l/v_m/k_2 & ... & \theta_l/v_m/k_n \end{bmatrix}$

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?

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

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.

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.

MarkyMark ha scritto: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. :-)

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ì: