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

[duf42]

Je présume que l'erreur que remonte MATLAB correspond à des dimensions qui ne correspondent pas?

Je pense que dans ton calcul, tu devrais utiliser l'opérateur '.*' au lieu de '*' notamment pour les opérations sur R et W.

Sinon pour débugger ton équation, essaye dans un premier temps de la "découper" en parties "élémentaires", ca devrait te permettre de cerner le problème plus facilement