Calcul et tracé d'une surface 3D

Salut à tous
J'ai un petit problème d'implémentation.
Matlab refuse de prendre en compte ma fonction avec le message

??? Subscript indices must either be real positive integers or logicals

Pour le code suivant:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

clear all l=0.1; m=0.01; n=48; a=15; p=4; vpm=20600; x=0.1:0.01:5; y=0.1:0.01:4.9; marge=zeros(length(x),1); r2=((-vpm).*((((1+(y+l+m)).^(-n))-1)/(y+l+m))); r1=((-vpm).*((((1+(x+l+m)).^(-n))-1)/(x+l+m))); d=((a)-(n/p)); marge(x,y)=1-(((r2-r1)/r1)/(d)); [xx,yy]=meshgrid(x,y); figure; surf(xx,yy,marge(x,y));

---

Je pense que l'erreur est basique mais je suis débutant.
Le but étant de tracer un plan en trois dimensions avec comme variables x,y (y<x car y décroissance de x) et marge(x,y)

Bonjour.

Ton problème vient de cette ligne :

Code:

---

marge(x,y)=1-(((r2-r1)/r1)/(d));

---

X et y sont des float, selon la définition que tu leurs a donné. Lorsque tu crées une matrice (ce qui est le cas de "marge"), les numéros de ligne et de colonnes doivent être entiers positifs.

Une solution serait la suivante :

Code:

---

marge((x-0.09)*100,(y-0.09)*100)=...

---

Reste que lorsque tu fais ça (comme dans ce que tu essayais de faire avant), tu crées une matrice de taille (approximativement) 500*500, alors qu'avant tu as initialisé cette variable comme un vecteur de taille 500*1.
Assures toi bien de ce que tu veux obtenir, un vecteur colonne ou une matrice ;)

je propose:

Code:

---

1 2 3 4 5 6

[xx,yy]=meshgrid(x,y); Z=zeros(size(xx)); Z(:)=1-(((r2-r1)/r1)/(d)); figure; surf(xx,yy,Z)

---

(d'autant plus que tu crées marge comme une matrice n*1 et ensuite tu veux le manipuler comme une matrice n*p.)

Bonjour,

:faq: Subscript indices must either be real positive integers or logicals qui concerne ta ligne marge(x,y) = ... ainsi que celle du plot.
J'aurais quelques remarques:
Tu sembles pre-allouer ta variable marge comme un tableau 1D, mais tu sembles vouloir marge(x,y), donc que veux-tu faire exactement?
De plus r1 et r2 doivent être des scalaires ou plutôt des vecteurs? Dans ce dernier cas, ce serait plutôt:

Code:

---

1 2	r2=-vpm*(((1+(y+l+m)).^(-n)-1)./(y+l+m)); r1=-vpm*(((1+(x+l+m)).^(-n)-1)./(x+l+m));

---

Une dernière chose: certaines paranthèses que tu as mises sont inutiles.
exemple: d=a-n/p;

Citation:

---

Envoyé par Winjerome

Bonjour,

:faq: Subscript indices must either be real positive integers or logicals qui concerne ta ligne marge(x,y) = ... ainsi que celle du plot.
J'aurais quelques remarques:
Tu sembles pre-allouer ta variable marge comme un tableau 1D, mais tu sembles vouloir marge(x,y), donc que veux-tu faire exactement?
De plus r1 et r2 doivent être des scalaires ou plutôt des vecteurs? Dans ce dernier cas, ce serait plutôt:

Code:

---

1 2	r2=-vpm*(((1+(y+l+m)).^(-n)-1)./(y+l+m)); r1=-vpm*(((1+(x+l+m)).^(-n)-1)./(x+l+m));

---

Une dernière chose: certaines paranthèses que tu as mises sont inutiles.
exemple: d=a-n/p;

---

Merci beaucoup!
Je souhaite en fait que r1 soit un vecteur calculé à partir de x pour des valeurs de 0.1 à 5; de même poour r2 pour des y de 0.1 à 4.9.
Pour ensuite faire un taux de variation au prorata de d. Après ces taux diminueront un seuil de 1 marge=1-...%

thxsmuch

Tu as alors au moins deux solutions:
Soit le faire à l'aide de boucles for:

Code:

---

1 2 3 4 5 6 7 8

N1 = length(r1); N2 = length(r2); marge = zeros(length(r1),length(r2)); for i = 1:N1 for j = 1:N2 marge(i,j) = 1 - (r2(j) - r1(i))/r1(i); end end

---

Soit de façon plus optimisée en utilisant la fonction repmat:

Code:

---

1 2 3

R1 = repmat(r1(:),1,N2); R2 = repmat(r2,N1,1); marge = 1-(R2-R1)./R1;

---

Salut à tous

Encore une fois sur le même code :-(

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23

l=0.1; m=0.01; n=48; a=15; p=4; vpm=20600; x=0.001:0.001:0.1; y=0.001:0.0001:0.1; N1=length(x); N2=length(y); [xx,yy]=meshgrid(x,y); marge=zeros(length(x),length(y)); for i=1:N1 for j=1:N2 x(i)=-vpm*(((1+((i./1000)+l+m)).^(-n)-1)./((i./1000)+l+m)); y(j)=-vpm*(((1+(((j+9)./10000)+l+m)).^(-n)-1)./((j+9)./10000)+l+m); d=a-(n/p); marge(i,j)=1-(((y(j)-x(i))./x(i))./d); end end x=Points(:,0.001:0.001:0.1); y=Points(:,0.001:0.0001:0.1); dx=1; dy=1; x_edge=[floor(min(x)):dx:ceil(max(x))]; y_edge=[floor(min(y)):dy:ceil(max(y))]; [X,Y]=meshgrid(x_edge,y_edge); Z=griddata(x,y,marge(i,j),X,Y);figure;

---

Le but étant de calculer une marge à partir de deux variable, x et y, qui sont eux-mêmes deux taux.
Je pense que le code est bon après avoir eu l'aide hier de monsieurs sympathiques.
Il me reste un problème

Code:

---

1 2 3 4 5

??? Error using ==> surf at 75 Z must be a matrix, not a scalar or vector.   Error in ==> mger at 17 figure; surf(xx,yy,Z);

---

J'ai tenté de le résoudre avec des Points et griddata sans succès, je veux juste faire un plan trois dimensions, la marge(i,j) en fonction de x et y. (pas x(i) et pas y(j))

Code:

---

1 2 3 4

??? Undefined variable Points.   Error in ==> mger at 16 x=Points(:,0.001:0.001:0.1);

---

Molto grazie!

Rémi

Bonjour,

qu'est censé faire cette ligne :

Code:

---

x=Points(:,0.001:0.001:0.1);

---

?

Tu n'as pas créer de variables Points, et de toute manière la sélection par des indices soit se faire avec des entiers, pas des nombres décimaux.

Quand à la première erreur, apparemment tu ne nous as pas donné le code associé...

Alors déjà première remarque:

Code:

---

1 2	x=Points(:,0.001:0.001:0.1); y=Points(:,0.001:0.0001:0.1);

---

La variable Points que tu veux utiliser n'est pas définie, de plus même si elle l'était, tu ferais la même erreur que dans ton premier message: les indices de tableaux doivent être des entiers positifs non nuls.
Ensuite toujours dans ton premier message, tu calculais des variables intermédiaires r1 et r2 pour obtenir le contenu de marge, or ici tu remplaces le contenu des variables x et y.
Pour finir, tu as une erreur sur l'utilisation de la fonction surf, mais je ne vois pas cette ligne dans le code que tu nous donnes (je présume que c'est la ligne suivante), de plus je te conseille de revoir les résultat de ces lignes qui selon moi n'est pas ce que tu cherches:

Code:

---

1 2 3 4

x_edge=[floor(min(x)):dx:ceil(max(x))]; y_edge=[floor(min(y)):dy:ceil(max(y))]; [X,Y]=meshgrid(x_edge,y_edge); Z=griddata(x,y,marge(i,j),X,Y);figure;

---

Bonjour à tous,

Voilà le code qui me pose problème (encore!)=

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

l=0.0; m=0.01; n=48; a=15; p=4; vpm=25393; x=0.001:0.001:0.1; y=0.001:0.0001:0.1; N1=length(x); N2=length(y); z=zeros(length(x),length(y)); for i=1:N1 for j=1:N2 z(i,j)=1-((((vpm*((1-((1+(((j+9)./10000)+l+m)./4).^(-n)))./(((j+9)./10000)+l+m)./4))-(vpm*((1-((1+((i./1000)+l+m)./4).^(-n)))./((i./1000)+l+m)./4)))/(vpm*((1-((1+((i./1000)+l+m)./4).^(-n)))./((i./1000)+l+m)./4)))/a-(n/p)); end end plot3(i/1000,(j+9)/10000,z(i,j),'-r.');

---

Pièce jointe 78740

Je ne parviens pas à dessiner une surface, pourtant z(i,j) varie en fonction de deux variables et il n'y pas de soucis sur l'integer... Ce point rouge me pique les yeux :cry:

Merci à vous!

Rémi

Citation:

---

Envoyé par mantus

Code:

---

plot3(i/1000,(j+9)/10000,z(i,j),'-r.');

---

---

Ici tu ne traces qu'un seul point :

• ta variable i est égale à sa dernière valeur d'itération: N1
• De même pour j = N2
• Toujours pareil z(i,j) représente une seule valeur: z(N1,N2)

Alors que toi tu veux tracer tous les points. Tu dois donc fournir toutes les coordonnées x et y et leurs valeurs correspondantes.
Ce que tu avais très bien essayé de faire au début en utilisant les fonctions meshgrid et surf (ou plot3),

Salut

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

figure; l=0.0; m=0.01; n=48; a=15; p=4; vpm=25393; x=0.001:0.001:0.1; y=0.001:0.0001:0.1; N1=length(x); N2=length(y); z=zeros(length(x),length(y)); for i=1:N1 for j=1:N2 z(i,j)=1-((((vpm*((1-((1+(((j+9)./10000)+l+m)./4).^(-n)))./(((j+9)./10000)+l+m)./4))-(vpm*((1-((1+((i./1000)+l+m)./4).^(-n)))./((i./1000)+l+m)./4)))/(vpm*((1-((1+((i./1000)+l+m)./4).^(-n)))./((i./1000)+l+m)./4)))/a-(n/p)); end end [i,j]=meshgrid(1:N1,1:N2); surf(i/1000,(j+9)/10000,z(i,j));

---

Que j'utilise "surf" ou "plot3" rien n'y fait, la surface est un point rouge. J'ai pourtant tenté de déclarer meshgrid avec i et j ou même x et y pour que surf le prenne en compte; même indiqué dans un autre endroit.

Il m'indique aussi

Code:

---

??? Maximum variable size allowed by the program is exceeded

---

Ou invariablement:

??? Error using ==> surf at 75 
Data dimension must agree

Merci

Bonjour,

Par exemple, tu peux faire :

Code:

---

1 2	[i,j]=meshgrid(1:N1,1:N2); surf(i,j,z)

---

ou :

Code:

---

surf(x,y,z,'EdgeColor','none')

---

Pour l'erreur, c'est un problème d'espace mémoire disponible. Tu peux lire ceci :
1107 - Avoiding 'Out of Memory' Errors

Thxsmuch
it runs well!

merci bcp

Rémi