Equation de récurrence

Version imprimable

Voir 40 message(s) de cette discussion en une page

26/11/2007, 17h56
Combomaster

Equation de récurrence
Bonjour je prenome Nicolas, Voila j'aimerais realiser une equations de reccurence avec matlab, et je ne vois pas du tout comment faire. C'est mes tout debuts de programmation, donc pas trop doué :)

voila ce que j'ai pour le moment :
Code:

1 2 3 4 5 6 7 8 9 10 clear all;close all;echo off k=[1,2,3,4,5,6,7,8]; y=[0,0,0,2,2,2,2,2]; figure(1);plot(k,y); Yref=0 Sh=zeros(8,1) for k=1:8 Sh(k)=Sh(k-1)+(y(k)-Yref(k)) end
si quelqu'un pouvait m'aider merci d'avance
je ne sais pas du tout si c'est juste je sais que mes 3 premieres lignes sont bonnes ainsi que ma formule :s
26/11/2007, 20h00
Jerome Briot

L'indexage de Sh dans la boucle n'est pas bon

Pour k=1 => Sh(k-1) => Sh(1-1) => Sh(0)

Or avec MATLAB, l'indexage commence toujours à 1. Tu dois donc très certainement commencer les itérations à partir de k=2, non ?

Ensuite, pour Yref(k), tu indexes un vecteur, alors que tu as définis un scalaire. Yref(1) existe bien puisque tu as initialisé Yref=0, mais Yref(2) n'existe pas
28/11/2007, 09h02
Combomaster

merci beaucoup j'ai reussi :)
28/11/2007, 10h10
Jerome Briot

Tu pourrais nous montrer la solution que tu as codé ?

Code:

1
2
3
4
5
6
7
8
9
10
k=[1,2,3,4,5,6,7,8];
y=[0,0,0,2,2,2,2,2];
figure(1);plot(k,y);
 
yref=zeros(8,1);
Sh=zeros(8,1);
for i=2:8;
   Sh(i)=Sh(i-1)+(y(i)-yref(i));
end
Sh

28/11/2007, 11h56
Jerome Briot

Es-tu sûr du résultat ?

yref vaut toujours 0 dans la boucle :?
28/11/2007, 11h58
Combomaster
oui je suis sur du resultat j'ai verifié en calculant la récurrence moi meme

par contre j'aurais une autre question

voila j'ai programmé cela afin d'avoir 4 graphe sur un seul
Code:

1 2 3 4 5 g = plot(P2,Vb); g1 = plot(P2,xbarbar1); g3 = plot(P2,LSCXbar); g4 = plot(P2,LICXbar); g5 = graph_union(g,g1,g3,g4);
et je voudrais mettre un autre graphe, mais en dessous de celui la donc un subgraph, hors je n'ai pas trouver le moyen de le faire :s

Citation:

Envoyé par Combomaster

oui je suis sur du resultat j'ai verifié en calculant la récurrence moi meme

Permet moi d'être sceptique :
Code:

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 28 29 30 31 32 k=[1,2,3,4,5,6,7,8]; y=[0,0,0,2,2,2,2,2]; % Sh yref=zeros(8,1); Sh=zeros(8,1); for i=2:8; Sh(i)=Sh(i-1)+(y(i)-yref(i)); end % Shh Shh=zeros(8,1); for i=2:8; Shh(i)=Shh(i-1)+y(i); end % Shhh Shhh=cumsum(y); %Tracés figure subplot(3,1,1) plot(Sh,'r-*') title('Sh') subplot(3,1,2) plot(Shh,'r-*') title('Shh') subplot(3,1,3) plot(Shhh,'r-*') title('Shhh')
:koi:

28/11/2007, 12h09
Combomaster

moi en faisant mon calcule j'obtiens la meme chose que ce que me donne matlab

Citation:

Sh =

0.
0.
0.
2.
4.
6.
8.
10.
28/11/2007, 12h12
Jerome Briot

yref ne sert à rien dans ton équation puisque cette variable vaut toujours 0 et que tu ne fais que l'ajouter à Sh
28/11/2007, 12h16
Combomaster

oui je sais, mais mon prof voullais que sa soit dedans
28/11/2007, 12h18
Jerome Briot

OK, donc pour obtenir plusieurs graphiques, utilise la fonction SUBPLOT comme dans le code que j'ai donné précédemment
28/11/2007, 12h30
Combomaster

merci beaucoup

voila ca que j'ai programmé pour avoir les 2 graphes hors il me met que le premier :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
 
figure
subplot(2,1,1)
g = plot(P2,Vb);
g1 = plot(P2,xbarbar1);
g3 = plot(P2,LSCXbar);
g4 = plot(P2,LICXbar);
g5 = graph_union(g,g1,g3,g4);
plot(g5,'r-*')
title('moyenne')
 
subplot(2,1,2)
g6=plot(P2,Rb);
g7=plot(P2,LSCRbar);
g8=plot(P2,LICRbar);
g9=graph_union(g6,g7,g8);
plot(g9,'r-*')
title('etendue')

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
figure
subplot(2,1,1)
hold on
g = plot(P2,Vb);
g1 = plot(P2,xbarbar1);
g3 = plot(P2,LSCXbar);
g4 = plot(P2,LICXbar);
g5 = graph_union(g,g1,g3,g4);
plot(g5,'r-*')
title('moyenne')

subplot(2,1,2)
hold on
g6=plot(P2,Rb);
g7=plot(P2,LSCRbar);
g8=plot(P2,LICRbar);
g9=graph_union(g6,g7,g8);
plot(g9,'r-*')
title('etendue')

28/11/2007, 12h51
Combomaster

est ce que sous scilab il y a une fonction que ressemble a hold on sur matlab?
28/11/2007, 12h59
Jerome Briot

Pour les équivalences des fonctions MATLAB/Scilab : http://www.scilab.org/product/dic-mat-sci/M2SCI_doc.htm

Cherche "hold" dans la liste et n'oublie pas de cliquer sur la petite flèche à coté du nom de la fonction ;)

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
figure
subplot(2,1,1)
set(gca(1),"auto_clear","off")                     //a remplacer par 'hold on' sous matlab
g = plot(P2,Vb);
g1 = plot(P2,xbarbar1);
g3 = plot(P2,LSCXbar);
g4 = plot(P2,LICXbar);
g5 = graph_union(g,g1,g3,g4);
plot(g5,'r-*')
title('moyenne')
 
subplot(2,1,2)
set(gca(1),"auto_clear","off")                     //a remplacer par 'hold on' sous matlab
g6=plot(P2,Rb);
g7=plot(P2,LSCRbar);
g8=plot(P2,LICRbar);
g9=graph_union(g6,g7,g8);
plot(g9,'r-*')
title('etendue')

toujours le meme problème il fait bien les 2 emplacements mais pas de 2ème graphe et il me dit que ma variable g n'est pas definit hors g est le nom de mon graphe

28/11/2007, 13h22
Jerome Briot

Il faudrait être un peu plus précis dès le début dans tes questions.

Utilises-tu MATLAB ou Scilab ? C'est quoi la fonction graph_union ?
Quel message d'erreur est exactement retourné ?
28/11/2007, 13h27
Combomaster

j'utilise scilab
le graphe_union permet : l'union de plusieurs graphes en un seul

Citation:

g5 = graph_union(g,g1,g3,g4);
!--error 4
undefined variable : g
at line 72 of exec file called by :
exec("C:/DOCUME~1/COMBOM~1/LOCALS~1/Temp/SCI_TMP_3828_/Untitled1.sce");
in execstr instruction called by :

Voir 40 message(s) de cette discussion en une page