Discussion :

[Invité]

Bonjour, j'ai un petit programme en matlab

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Invité]

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 du site, elle contient un bon nombre d'erreurs courantes expliquées.

[Invité]

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Invité]

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.

[Invité]

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Invité]

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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