Inner matrix dimensions must agree

Bonjour, j'ai un petit programme en matlab

Code:

---

1 2 3 4 5 6 7 8 9

n=2; h=20; d=10; L=300; f=8; x=-150; H=(f*((2*x)./L).^2+h) hc=[H*ones(n,1)] xc=[d*ones(n,1)]

---

si je veux faire varier le point x comme suit

Code:

---

1 2 3 4 5 6 7 8 9

n=2; h=20; d=10; L=400; f=8; x=-150:50:150; H=(f*((2*x)./L).^2+h) hc=[H*ones(n,1)] xc=[d*ones(n,1)]

---

j'obtiens

??? Error using ==> mtimes 
Inner matrix dimensions must agree. 
Error in ==> Untitled at 8 
hc=[H*ones(n,1)]

Merci d'avance

Bonsoir,

Je vois deux problèmes dans ton code:
- Pour commencer, tes vecteurs x et H comptent length(x) et length(H) 7 valeurs, donc tu ne peux multiplier H par ones(n,1) pour n égal à 2: il te faut le même nombre de valeurs.
- Ensuite, il faut faire attention que x et H sont des vecteurs 7x1 (à différencier de 1x7), donc la multiplication se fera plutôt avec ones(1,N) avec N = length(x) (égal à 7 ici).

PS: pense à regarder la :faq: du site, elle contient un bon nombre d'erreurs courantes expliquées.

merci jerome pour ta belle reponse et pour l'explication exactement que dois je faire pour avoir la reponse suivante:

Code:

---

1 2 3 4 5 6 7 8 9

hc =   24.5000 22 20.5000 20 20.5000 22 24.5000 24.5000 22 20.5000 20 20.5000 22 24.5000   xc =   10 10 10 10 10 10 10 10 10 10 10 10 10 10

---

merci pour ton aide

Comme je te l'ai dis dans mon dernier point, tu peux déjà introduire une variable N représentant la longueur de x.
Ensuite, je présume que tu veux que le nombre de ligne de hc dépende de n et ne soit pas simplement de 2, donc tu peux pour cela utiliser la fonction repmat (doc repmat pour obtenir sa documentation sur Matlab).
Ensuite pour xc, il te faudra juste rajouter la longueur N dans ton ones.

merci ,c'est-à-dire n*length(x) 2*7
est-ce que si possible comme suivant

Code:

---

1 2	hc=[H*repmat(n,length(x))] xc=[d*repmat(n,length(x))]

---

Pas tout à fait, il faut rentrer d et H dans la fonction

Code:

---

1 2 3

N = length(x); hc = repmat(H,n,1); H est déjà de longueur 7 xc = repmat(d,n,N); % ou d*ones(n,N); mieux pour un scalaire

---

Il suffisait d'essayer pour le voir ;)

merci ca marche bien .
si je dois executer l'operation suivante

Code:

---

1 2 3 4 5

for i=1:14; for j=1:14; x1(i)=xc(i)+0.05; y1(j)=hc(j)+0.025; R1(i,j)=sqrt((x1(i)+1).^2+y1(j).^2);

---

je dois executer seulement
i=1 et 2 pour j=1 et 2
puis i=3 et 4 pour j=3 et 4 et ainsi de suite jusqu' i=13 et 14 pour j=13 et 14 .
c'est un peu compliqué.