Equation différentielle & Matlab

[invite50429c49]

Bonjour,
j'aimerai réaliser un graphe du pic maximal d'un circuit RLC en fonction de la résistance utilisée. Je sais déja avoir le graphe du circuit, mais je ne sais comment faire varier la résistance dans cette équation, je me plante complètement avec les martices...
Voici mon petit programme matlab:
function dy = rlc2(t,y) %Equa diff RLC avec antenne
dy = zeros(2,1); %Un vecteur colonne
L=1.5/10000
R=100;
C=2.5/10000000000
Uc=1 %Courant induit par l'antenne
w=1/(L*C)^(1/2) %la fréquence résonnante
dy(1) = y(2);
dy(2) = Uc*sin(w*t)/(L*C)-y(1)/(L*C)-(R/L*y(2) )
end

Avec comme commande:
[T,Y] = ode45('rlc2',[0 0.00006],[-0.01 0.01]);plot(T,Y(:,2),'-')
end
-----

[invite551c2897]

Bonjour.
Tu peux utiliser global :

clear
global R
n=0;
for R=10:10:50
  n=n+1;
[T,Y] = ode45('rlc2',[0 0.00006],[-0.01 0.01]);
figure(n)
plot(T,Y(:,2),'-')
 grid
 title(['R = ' num2str(R)])
end

[invite50429c49]

Ah merci beaucoup!! J'aurais jamais trouvé sinon =) Ils sont fous ces profs de nous demander de modéliser avec matlab alors qu'on a jamais eut un seul cours d'aucune notion de matlab ^^'

[invite551c2897]

Bonjour.

du pic maximal d'un circuit RLC en fonction de la résistance utilisée

Pour avoir la courbe pic=f(R) :

clear
global R
n=0;
for R=10:10:100
n=n+1;
[T,Y] = ode45('rlc2',[0 0.00006],[-0.01 0.01]);
 maxR(n)=max(Y(:,2));
 resist(n)=R;
end
plot(resist,maxR)
grid

[A voir en vidéo sur Futura]