Discussion :

[Udani]

Bonjour,

J' ai programme une fonction "john" qui doit effectuer des calculs en fonction de plusieurs paramétres fixes et variables "t" et "s".

En analysant mon programme, on vois bient que le vecteur "t" est de longuer 200, "s" aussi. le "i" de "for", varie bien de 1 jusqu'à 200, donc pourquoi le message erreur en bas ?? j'ai passé des heures à essayer de trouver une solution, en vain. Merci de votre aide!!!
le CODE:

Code :

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

1
2
3
4
5
6
7
8
9
 
n=200;
t=linspace(0,5/12,n)';
s=linspace(0,100,n)';
J=0;
 
for i=1:n
  J(i)=john(s(i),50,0.1,0.4,t(i),5/12)
end

In an assignment A(I) = B, the number of elements in B and I must be the same.

Merci

[christophe_halgand]

Salut,

Désolé, je ne peux pas tester ton code mais il faut que tu initialises "J" comme ceci :

Code :

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

J=zeros(n,1);

De manière à ce que l'allocation en mémoire fasse gagner du temps à matlab.

Ensuite pour ton message d'erreur, généralement tu le résous en utilisant la fonction "size" sur tout les paramètres avant de faire le calcul afin d'être sûr que tu n'aies pas de problème de dimension comme ici.

Voilà, j'espère que ça t'aidera

Christophe

[Jerome Briot]

In an assignment A(I) = B, the number of elements in B and I must be the same.

Voir la : In an assignment A(I) = B, the number of elements in B and I must be the same.

[Udani]

Merci , maintenant ca marche

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
function johnson2=john_plot(N)
 
n=N;
 
t=linspace(0,5/12,n)';
s=linspace(0,100,n)';
J=zeros(n,1);
 
for i=1:n
  J(i)=john(s(i),50,0.1,0.4,t(i),5/12)
end
 
Z1=J;
save johnson;
Aus=[t s Z1 ];
save john.dat Aus -ascii;
plot(s,Z1,'b');
axis([0 101 0 50]);  
legend('Methode de John');

[Udani]

Bonjour,

problème résolu, nouveau problème!!
donc j'ai reussi à faire marche la function "john" comme fonction de temps 't" et d'espace "s". La j'iamerais faire une representation à 3 Dimension.

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
function plot3d=plot3d_joh(N)
 
n=N;
sigma = 0.4; 
K = 50; 
r = 0.1; 
T = 5/12; 
Svals = linspace(0,100,n);
tvals = linspace(0,5/12,n);
 
[Smat,tmat]=meshgrid(Svals,tvals);
 
P=johnson(Smat,K,r,sigma,tmat,T); %la fonction john(So, K, r, sig, t, T)
 
mesh(Smat,tmat,P);
 
xlabel('S'),ylabel('t'),zlabel('P(S,t)','rot',0);
box off;

Message erreur:

Code :

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

1
2
3
4
5
6
7
8
9
EDU» plot3d_john(20)
??? Error using ==> surface
Matrix dimensions must agree.
 
Error in ==> C:\MATLAB_SE_5.3\toolbox\matlab\graph3d\mesh.m
On line 93  ==>       hh = surface(x,y,z,'FaceColor',fc,'EdgeColor','flat', 'FaceLighting', 'none', 'EdgeLighting', 'flat');
 
Error in ==> C:\MATLAB_SE_5.3\work\plot3D_john.m
On line 17  ==> mesh(Smat,tmat,P);

je ne sais plus quoi faire, car j'ai testé un autre fonction et ca marche!
Merci de votre aide.

[Jean Dumoncel]

Salut,

quelles sont les dimensions de Smat,tmat et P avant l'appel de mesh? (puisque apparemment c'est un problème de dimension)

Code :

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

1
2
3
4
5
P=johnson(Smat,K,r,sigma,tmat,T); %la fonction john(So, K, r, sig, t, T)
size(Smat)
size(tmat)
size(P)
mesh(Smat,tmat,P);

[Udani]

bonjour,

es dimensions sont:

size(Smat) = 20 20
size(tmat) = 20 20
size(P) = 1 1

[Jean Dumoncel]

Donc, là tu devrais te rendre compte qu'il y a un problème : P devrait être une matrice de dimension 20 par 20...


Dans ta fonction Johnson, as-tu pris en compte qu'elle devait te retourner une variable de la taille de Smat et tmat?

[Udani]

Donc, là tu devrais te rendre compte qu'il y a un problème : P devrait être une matrice de dimension 20 par 20...


Dans ta fonction Johnson, as-tu pris en compte qu'elle devait te retourner une variable de la taille de Smat et tmat?

c'est vrais ! sauf que d'habitude, la fonction johnson, peux prendre des valeurs à partir d'un vecteur. Comment peut-elle prendre des valeurs d'une matrice (deux matrices Smat et tmat , les deux de dimensions 20x20)

j'ai déja essayé ceci pour plot de dimension 2, mais ca marche pas pour la dimension 3:

Code :

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

1
2
3
4
P=zeros(n,1);
for i=1:n
   P=johnson(Smat(i),K,r,sigma,tmat(i),T);
end

??

[Jean Dumoncel]

A partir de là, c'est un problème de conception mais ce n'est pas directement lié à matlab.

Mais ne sachant pas ce que tu fais dans ta fonction johnson, c'est à toi de réfléchir à ce que tu calcules et à ce que tu veux afficher.

[Udani]

oui, je crois aussi que je devrais revoir la conception de mes programmes.
Merci pour votre aide et vos conseils!
Je vous tiens au courant.