Discussion :

[Alecks]

Bonjours à tous,

J'ai un plotting à faire qui me cause des problèmes. L'axe des x est une valeur incrémentielle 0:0.01:6.25 et l'axe des y est un vecteur de la même taille que je remplit avec une boucle for.

Le problème est que seule la première valeur veut bien s'insérer dans mon vecteur. Les loops suivantes sont incapables d'assigner la valeur. Voici le code et merci pour votre aide.

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
28
29
30
31
 
clear all
close all
clc
 
%% cste
rad = pi/180; %% 1 rad = 180/pi degré
m2 = 600;
g = 9.81;
B = 8000;
gamma = 2.91 * rad;
 
%% init. table
Fa1_table = ones(1,626);
 
%% FA - phase 1
for th1 = 0:0.01:6.25
    theta1 = th1 * rad;
    i = 1;
    beta1 = ( 2.6 - 2.5*cos(theta1) ) / ( 1.5*sin(theta1) + 1.5 );
    alpha1 = 2 * rad; %% rad/s²
    w1 = 0.8*theta1 * rad; %% rad/s
    %% store
    Fa1_table(i) = ( 4*m2*alpha1 + 6*m2*g*sin(pi-theta1-gamma) + B*w1 ) / ( 1.5*cos(beta1)*sin(theta1) + 1.5*sin(beta1)*cos(theta1) );
    i = i + 1;  
end
 
%% graph
figure(1)
Th1 = 0 : 0.01 : 6.25;
plot(Th1,Fa1_table)

[Merel]

Salut,

tu as juste fait une erreur à l'intérieur de ta boucle for.

Tu as

Code :

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

i=1

et ensuite dans la même boucle

Code :

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

i=i+1

Du coup, à chaque itération, tu écrases i.

suffit donc de sortir l'initialisation du i de la boucle ;-)

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
28
29
30
clear all
close all
clc

%% cste
rad = pi/180; %% 1 rad = 180/pi degré
m2 = 600;
g = 9.81;
B = 8000;
gamma = 2.91 * rad;

%% init. table
Fa1_table = ones(1,626);
  i = 1;
%% FA - phase 1
for th1 = 0:0.01:6.25
    theta1 = th1 * rad;
  
    beta1 = ( 2.6 - 2.5*cos(theta1) ) / ( 1.5*sin(theta1) + 1.5 );
    alpha1 = 2 * rad; %% rad/s²
    w1 = 0.8*theta1 * rad; %% rad/s
    %% store
    Fa1_table(i) = ( 4*m2*alpha1 + 6*m2*g*sin(pi-theta1-gamma) + B*w1 ) / ( 1.5*cos(beta1)*sin(theta1) + 1.5*sin(beta1)*cos(theta1) );
    i = i + 1;  
end

%% graph
figure(1)
Th1 = 0 : 0.01 : 6.25;
plot(Th1,Fa1_table)

[Alecks]

Mais quel idiot!

Je cherchais une erreur trop compliqué.

Merci

[Merel]

Quand on a la tête dans le truc, on ne voit plus trop ce qu'on a écrit...

Et ce sont en général les plus petites boulettes qui sont les plus pénibles .