Discussion :

[mostafaziad]

bonjour à tous

j’ai résolu un système des équations différentielles . mais quand je fait un zoome dans la figure qui représente la solution je trouve que la courbe n’est pas tout à fait courbé .je ne comprend pas ou sa viens se problème /voila le code:

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
22
23
24
25
26
27
28
29
30
31
32
function dy = eq_bc_ouv(t,y,Ex)
 
global tstar
 
dy = zeros(size(y,1),1);
 
imput=0*10^-7;
index=max(find(Ex(:,1)*3600<=t));
Ta = Ex(real(index),2);
dy(1)=(4.403*10^-5)*(y(2)-y(1))+(1.260*10^-4)*(Ta-y(1));
dy(2)=(218.414*10^-5)*(y(1)-y(2))+93.289*10^-5*(y(3)-y(2));
dy(3)=(408.173*10^-7)*(y(2)-y(3))+(33.583*10^-7)*(10-y(3))+ imput;
 
if t >3*10^5
    imput=3000*10^-7;
    index=max(find(Ex(:,1)*3600<=t));
    Ta = Ex(real(index),2);
    dy(1)=(4.403*10^-5)*(y(2)-y(1))+(1.260*10^-4)*(Ta-y(1));
    dy(2)=(218.414*10^-5)*(y(1)-y(2))+93.289*10^-5*(y(3)-y(2));
    dy(3)=(408.173*10^-7)*(y(2)-y(3))+(33.583*10^-7)*(10-y(3))+ imput;
end
if t > 8*10^5
 
    tstar = 8*10^5;
 
    imput2=3300*10^-7;
    index=max(find(Ex(:,1)*3600<=t));
    Ta = Ex(real(index),2);
    dy(1)=(4.403*10^-5)*(y(2)-y(1))+(1.260*10^-4)*(Ta-y(1));
    dy(2)=(218.414*10^-5)*(y(1)-y(2))+93.289*10^-5*(y(3)-y(2));
    dy(3)=(408.173*10^-7)*(y(2)-y(3))+(33.583*10^-7)*(10-y(3))+ imput2;
end

la fonction main :

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
% global tstar
%resolution
Ex=xlsread('test');
x0=[10;20;15];
tspan=[0:3600:3600*400];
 
[t,y]=ode15s(@eq_bc_ouv,tspan,x0,[],Ex)
tnouv=t(223:end);
ynouv=y(223:end,2);
 
% % affichage
plot((tnouv-222*3600)/3600,ynouv-13.7779)
axis([-1 400 0 0.4]);
 
 
tnouv2=tnouv-(222*3600);
ynouv2=ynouv-13.7779;

[ancrou]

le continu n'existe pas en programmation
Tu as forcément un échantillonnage, des points.
Il doit y avoir un moins de réduire le pas, avoir plus de points et donc une représentation plus proche de la réalité
Et avoir une belle courbure