Problème d'allocation de valeur dans un for
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:
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) |