Résolution analytique sur matlab

Version imprimable

Bonjour,
J'essaye de résoudre une EDP sur matlab. J'ai trouvé (par calcul) que U(x,t)=(exp(-t)+8/(pi)^3)*sin(pi*x)+8/(pi)^3*sin(x)*Somme(...)

J'ai calculé ma somme mais mon programme ne marche pas:
Code:

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 clear all close all clc N=1000; % Initialisation somme S: S=1/27; % Calcul somme: for i=2:N S=S+1/(2*i+1)^3; end S for t=0:0.1:1 for x=0:0.1:1 U=(exp(-t)+8/(pi)^3)*sin(pi*x)+8/(pi)^3*sin(x)*S; end end % Affichage solution en 3D: imagesc(x,U)
Est-ce que quelqu'un saurait d'où vient le problème ? Merci !

23/02/2010, 09h06
Jerome Briot

Citation:

Envoyé par dodo1245

J'ai calculé ma somme mais mon programme ne marche pas:

Peux-tu être plus précis ?
MATLAB renvoie-t-il un message d'erreur ?
Si oui, lequel ?
23/02/2010, 13h31
dodo1245

Non matlab ne renvoie pas de message d'erreur, c'est juste que j'ai un fond vert comme graphique quand j'écris imagesc(U) (c'est une équation de chaleur normalement, donc je devrai observer la source et la température autour, donc un graphique qui prend différentes teintes).

Bonjour.
Peut_être :

Code:

1
2
3
4
5
6
7
8
9
10
i=0;j=0;
for t=0:0.1:1
    i=i+1;j=0;
    for x=0:0.1:1
        j=j+1;
        U(i,j)=(exp(-t)+8/(pi)^3)*sin(pi*x)+8/(pi)^3*sin(x)*S;
    end 
end
% Affichage solution en 3D:
imagesc(U)