Discussion :

[Jefeh]

Bonjour,
voilà je travaille pour le moment sur un projet dans lequel je suis amené à travailler avec des matrices.

A un moment donné, j'utilise une boucle dans laquelle je fais un od45, jusqu'à aucun problème.
Mais dans cette même boucle, après le ode45, j'incrémente un des éléments de mon vecteur de la manière suivante :

Code :

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

x0v(6) = x0v(6) + 1/t^2

Le problème, c'est que lorsque je lance mon programme, j'ai l'erreur suivante

Code :

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

1
2
3
4
5
??? Error using ==> mpower
Matrix must be square.
 
Error in ==> Question5 at 34
        x0v(6) = x0v(6) + (1/t^2);

Et j'aimerais donc bien savoir ce qu'est la fonction "mpower" et savoir aussi quelle opération l'utilisée.

[tug83]

IL suffit de faire

Code :

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

docsearch mpower

pour trouver de la doc sur cette fonction depuis le prompt MATLAB. Une recherche sur le site de TMW te renvoie sur les Arithmetic Operators

[Jefeh]

Ah ben merci je ne connassais que la commande doc .

Je sais donc ce qui en va pas mais je ne sais pas pourquoi ça ne va pas...

J'ai plus qu'à chercher merci

Ce qui est bizarre c'est que quand je vais les mêmes opérations mais dans la "console" de MatLab, je n'ai pas de problème...

[sjrd]

Il semble que ta variable t ne soit pas une valeur simple, mais bien une matrice (non carrée, d'ailleurs). Du coup, l'opération t^2 ne peut être calculée.

Vérifie que t ne devrait pas être remplacé par une valeur simple à cet endroit.

[tug83]

à mon avis ajoute le . devant la puissance:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 >>t=rand(2,1)
t =
   0.814723686393179
   0.905791937075619
>> 1/t^2
??? Error using ==> mpower
Matrix must be square.

>> 1/t.^2
ans =
                   0   1.218829898200963

[LordPeterPan2]

à mon avis ajoute le . devant la puissance:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 >>t=rand(2,1)
t =
   0.814723686393179
   0.905791937075619
>> 1/t^2
??? Error using ==> mpower
Matrix must be square.

>> 1/t.^2
ans =
                   0   1.218829898200963

Le 0 me choque, je ferais plutôt :

Code :

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

1./t.^2

Enfin çà dépend de ce qu'il veut faire une inversion terme à terme ou une version tel que

Code :

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

(1/t.^2)*(t.^2) = 1

[Jefeh]

Voilà merci pour vos réponses, j'ai trouvé la solution ...
En faite le problème que j'avais été vraiment stupide, j'ai vraiment pas fait attention à ce que je faisais puisque j'utilisais 2 fois la variable t dans mon programme pour désigner 2 choses différentes.
Il m'a donc suffit de changer le nom de celle-ci.
Un grand merci pour vos réponses qui m'ont obligés à faire attention à ce que j'écrivais lol