Discussion :

[glenngrauwels]

Bonjour, alors voila j'ai un gros problème. Je dois mettre l'équation différentielle (pièce jointe) sous Matlab via Euler, je n'y arrive pas... Quelqu'un pourrait-il m'aider?
Merci d'avance
ps: petite erreur: il faut changer chaque signe de la 2ème équation du système, mettre les signes opposés

[glenngrauwels]

Bonjour, j'ai un problème avec mon code, je dois trouver une fréquence de résonance. J'ai basé ce code sur un autre qui fonctionnait, mais dans celui-ci il y a un problème mtimes... Quelqu'un pourrait-il m'éclairer svp
Merci d'avance

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
% L'Èquation ‡ rÈsoudre est:
 
%v''(t)+1/(R*C)*v'(t)+1/(L*C)*v(t)=A*W*cos(wt))
 
h=0.0001;		
t=0:h:2;        
A=5;
L=10^-4
C=10^-9
W=2*pi*500000:1000:2*pi*1000000
R=100:100:1000000
 
x(1,:)=[0 0];	
 
for i=2:length(t),	
 
   k1=[x(i-1,2) A*W*(cos(W*t(i)))-R^(-1)*C^(-1)*x(i-1,2)-L^(-1)*C^(-1)*x(i-1,1)];	%Calculate derivatives;
   x(i,:)=x(i-1,:)+k1*h;		
end
 
plot(t,x(:,1),t,(cos(W*t(i)));

[christophe_halgand]

L'erreur se trouve dans ta première ligne de code après le "for" :

En effet, quand tu calcules R*C, il faut que tu prennes en compte que ce sont tous deux des vecteurs...

Code :

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

1
2
3
size(L^(-1)*C^(-1)*x(i-1,1))  
size(R^(-1)*C^(-1)*x(i-1,2))
size(A*W*(cos(W*t(i))))

je te conseille dans un premier temps de réaliser ta boucle sans faire varier trop de paramètres en même temps. Valides d'abord ton code avec des calculs de scalaires puis tu pourras passer au vecteur...

Christophe