Bonjour,
Je réalise sous matlab un programme permettant de résoudre une EDO de Runge Kutta d'ordre 4. Il s'agit d'une EDO de thermique.
Mais Matlab détecte des erreurs que je ne comprends pas
Les erreurs sont:
-la fonction ' f ' n'est pas défini
-il y a une erreur à la ligne a=feval(f,t,T)
-de même pour la ligne [t,T]=f_rangeK(' f ',t0,T0,tf,h)
Ci-dessous le fichier de commande et le fichier de la fonction
Pouvez vous me dire que signifie ces erreurs
merci beaucoup
%*************************************************************************%
% PROJET MATLAB
%
% Problème 1 : Etude thermique d'un fusible en plomb
%*************************************************************************%
%*************************************************************************%
% FICHIER DE COMMANDE c_probleme1.m
%*************************************************************************%
% paramètres du fusible
I=1;
Cp = 130; % capacité calorifique massique du plomb (J/kg/K)
u = 11300; % masse volumique du plomb (kg/m^3)
k = 2.5; % paramètre du terme de convection
o = 5.67*10^(-8); % paramètre du terme de rayonnement
rho = 20.7*10^(-8); % resistivité electrique ( ohm.m)
r = 10^(-3); % rayon du fusible (m)
Pv = rho*(I/(pi*r.^2)).^2; % puissance volumique
% condition initiale
T0 = 298; % Température initiale du fusible (K)
Tf = 600; % température de fusion du plomb (K)
t0= 0; % borne inférieure de l'intervalle de temps
tf = 80; % borne supérieure de l'intervalle de temps
f = 'Pv./(Cp.*u)-2*(k.*(T-T0).^(5/4)+o.*(T.^4-T0.^4))./(r.*Cp.*u)';
h=1;
[t,T] = f_rangeK('f',t0,T0,tf,h);
% tracé de T(t)
figure(1);
hold on
plot(t,T)
title('evolution de la température en fonction du temps');
xlabel('Temps');
ylabel('Température (K)');
%*************************************************************************%
% PROJET MATLAB
%
% Problème 1 : Etude thermique d'un fusible en plomb
%*************************************************************************%
%*************************************************************************%
% FICHIER DE FONCTION f_RangeK.m
%*************************************************************************%
% RESOLUTION D'EDO DU 1er ORDRE PAR LA METHODE DE RANGE KUTA D'ORDRE '
function [t,T] = f_rangeK(f,t0,T0,tf,h)
% [t,T] : tableau de valeurs de la fonction T(t)
% t0, tf : bornes de l'étude de la fonction
% T0 : conditions initiales : T(t0) = T0
% h : pas de discrétisation
% a,b,c,d : élément de calcul pour la méthode
tx = t0:h:tf;
Ty = T0;
for i=1:floor((tf-t0)/h)
t = tx(i);
T = Ty(i);
a = feval(f,t,T);
t = tx(i)+h/2;
T = Ty(i)+a*h/2;
b = feval(f,t,T);
t = tx(i)+h/2;
T = Ty(i)+b*h/2;
c = feval(f,t,T);
t = tx(i)+h;
T = Ty(i)+c*h;
d = feval(f,t,T);
Ty(i+1) = Ty(i)+h*(a+2*b+2*c+d)/6;
end
t = tx;
T = Ty;
... moi aussi ca marche pa