Discussion :

[TiMeWaLk]

Bonjour à tous,

Je sais qu'une matrice de fonction H peut être défini de la façon suivante sous matlab:

H=@(z) [R(z) G(z)], R et G étant des matrices 2*2 par exemple défini telle que:

R=@(z) [z+1 z+2;z+3 z+4];
G=@(z) [z+5 z+6;z+7 z+8];

Mon problème survient lorsque H, R ou G sont beaucoup plus grosses et que je veux définir chaque coefficient. Hors il m'est impossible d'écrire avec matlab quelque chose comme:

Code :

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

1
2
3
4
5
6
7
8
 
for x=1:100
for y=1:100
 
R(x,y)=@(z) z+x+y;
 
end
end

Connaissez vous une méthode pour obtenir tout de même une matrice de fonctions?

Je précise que je ne peux pas créer une function (sous forme d'un fichier .m) en effet, je désire avoir ma matrice en quelque sorte précalculée car elle sera utilisée dans une équation différentielle du type: v'(z)=H(z)=v(z) et qu'étant très grosse je n'ai pas envie qu'elle soit recalculée pour chaque z utilisé par ode45

Merci d'avance.

[duf42]

Bonjour,

Je suis pas sûr que ca répond exactement à ton besoin, mais tu pourrais créer quelquechose du genre:

Code :

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

1
2
3
4
5
6
7
8
 
fcn = @(a,b,c) (@(x,y,z) (a*x+b*y+c*z));
 
for i=1:100
   for j=1:100
       M{i,j} = fcn(2*i, 4*j, i+j);
   end
end

M contiendra les function handles de (a*x+b*y+c*z) où a,b et c sont déjà définis (dans la boucle).
De plus, il est conseillé dans l'aide de MATLAB d'utiliser un cell array pour les function handles (ce qui ne change pas grand chose au final).

Bonne journée,
Duf

[TiMeWaLk]

Merci,

Je vais voir ce que je peux faire avec ça. Cependant je ne peux pas utiliser les cells pour résoudre une équation différentielle si je ne me trompe pas.

[duf42]

Cependant je ne peux pas utiliser les cells pour résoudre une équation différentielle si je ne me trompe pas.

Je comprends pas ce que tu veux dire par là

En fait, je sais pas quelle est ta version de MATLAB mais dans la version 7.1, il y a un avertissement dans l'aide informant que les matrices de function handle ne seront plus supportées dans les versions suivantes. Je ne sais pas ce qu'il en est exactement mais tu risques de ne pas avoir trop le choix...

[mihaispr]

Une autre solution:

Code :

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

1
2
3
4
 
[I,J]=ndgrid(1:100,1:100);
 
M=@(x,y,z) 2*I*x + 4*J*y + (I+J)*z

Code :

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

1
2
 
>>help ndgrid

[TiMeWaLk]

Je voulais juste dire qu'on ne peut pas faire d'opérations comme ceci:

erty =

[1] [2] [3]
[4] [5] [6]

>> ytre={6 5 4;3 2 1}

ytre =

[6] [5] [4]
[3] [2] [1]

>> erty.*ytre

En effet je comptait me servir de cette matrice H.

[mihaispr]

OK.

Donc ca t'aide un peu ou non? (la solution avec ndgrid?)

Poste la solution ou si tu as des nouveaux problemes nous sommes ici pour essayer de t'aider!

Donc s'il y a d'autres problemes continue la discussion!

Amicalment,

Michel

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
>> erty ={1 2 3 ; 4 5 6}
 
erty = 
 
    [1]    [2]    [3]
    [4]    [5]    [6]
 
>> ytre={6 5 4 ; 3 2 1}
 
ytre = 
 
    [6]    [5]    [4]
    [3]    [2]    [1]
 
>> cellfun(@mtimes, ytre,erty,'uniformoutput',false)
 
ans = 
 
    [ 6]    [10]    [12]
    [12]    [10]    [ 6]
 
>> cellfun(@mtimes, ytre,erty,'uniformoutput',true)
 
ans =
 
     6    10    12
    12    10     6

[TiMeWaLk]

En fait j'ai abandonnée l'idée d'utiliser une matrice de fonctions. J'ai améliorée la façon dont je calculais H via l'indexage linéaire. Maintenant je cherche à créer quelquechose du type de ode45 (en moins compliqué) pour que ce solver gère les équations du type V'=H*V avec V une matrice, ce qui me permettrait de calculer H, 5 fois uniquement par palier du solver mais ceci pour toutes les colonnes. On va dire que le sujet dérive un peu donc...