Inscription / Connexion Nouveau Sujet
Posté par
joug
programme maple 23-06-09 à 16:48

slt vous avez trouver votre erreur  svp trouver mon erreur voici le programme restart;
> U:=unapply(Diff( u,u$2)+r^(-1)*Diff(u,r)=c^(-2)*Diff(u,t$2),r,t);
> #ceci est un comentaire
> print(`les condition initiales`);
> F:=f(r)=u(r,0);
> V:=Diff(u(r,0),t)=c*g(r);
> #ceci est un comentaire
> print(`opérateur différentielles linéaires`);
> R1:=L[r]=Diff(u(r,t),t);
> T:=L[t,t]=Diff(u(r,t),t$2);
> R:=L[r,r]=Diff(u(r,t),r$2);
> #ceci est un comentaire
> print(`opérateur d'intégration double`);
> H:=L[t,t]^(-1)=int(int(u(r,t),t=0..t),t=0..t);
> Q := u(r,0)+c*g(r)*t;
> u1(r,t):=unapply(Q+c^2*((L[t,t]^(-1)*L[r,r]*u(r,t))+(L[t,t]^(-1)*(r^(-1))*L[r]*u(r,t))),r,t);u=Sum('u[n](r,t)','n'=0..infinity);
> E:=u[0]=Q;
> Z:=u[1]=unapply(c^2*((L[t,t]^(-1)*L[r,r]*E)+(L[t,t]^(-1)*(r^(-1))*L[r]*E)));
> Y:=u[2]=c^2*((L[t,t]^(-1)*L[r,r]*Z)+(L[t,t]^(-1)*(r^(-1))*L[r]*Z));
> sol:= u[k+1]=unapply(a^2*((L[t,t]^(-1)*L[r,r]*u(k))+L[t,t]^(-1)*(r^(-1)*L[r]*u(k))));
quand je donne des valeur à F et g(r) et l'es réaffiche
et je veut tracer en prenent différente valeur de c=4,5,6,7 en dim 3d avec F=r; et g(r)=1;

*** message déplacé ***

Posté par
joug
programme de l'équation de la corde vibrante avec la méthode déc 23-06-09 à 16:57


slt j'ai aussi une erreur  svp trouver mon erreur voici le programme
restart;
> U:=unapply(Diff( u,u$2)+r^(-1)*Diff(u,r)=c^(-2)*Diff(u,t$2),r,t);
> #ceci est un comentaire
> print(`les condition initiales`);
> F:=f(r)=u(r,0);
> V:=Diff(u(r,0),t)=c*g(r);
> #ceci est un comentaire
> print(`opérateur différentielles linéaires`);
> R1:=L[r]=Diff(u(r,t),t);
> T:=L[t,t]=Diff(u(r,t),t$2);
> R:=L[r,r]=Diff(u(r,t),r$2);
> #ceci est un comentaire
> print(`opérateur d'intégration double`);
> H:=L[t,t]^(-1)=int(int(u(r,t),t=0..t),t=0..t);
> Q := u(r,0)+c*g(r)*t;
> u1(r,t):=unapply(Q+c^2*((L[t,t]^(-1)*L[r,r]*u(r,t))+(L[t,t]^(-1)*(r^(-1))*L[r]*u(r,t))),r,t);u=Sum('u[n](r,t)','n'=0..infinity);
> E:=u[0]=Q;
> Z:=u[1]=unapply(c^2*((L[t,t]^(-1)*L[r,r]*E)+(L[t,t]^(-1)*(r^(-1))*L[r]*E)));
> Y:=u[2]=c^2*((L[t,t]^(-1)*L[r,r]*Z)+(L[t,t]^(-1)*(r^(-1))*L[r]*Z));
> sol:= u[k+1]=unapply(a^2*((L[t,t]^(-1)*L[r,r]*u(k))+L[t,t]^(-1)*(r^(-1)*L[r]*u(k))));
quand je donne des valeur à F et g(r) et l'es réaffiche
et je veut tracer en prenent différente valeur de c=4,5,6,7 en dim 3d avec F=r; et g(r)=1;

*** message déplacé ***

Niveau algorithmique
Partager :

programme de l'équation de la corde vibrante avec maple

Posté par
joug
23-06-09 à 17:06



slt j'ai une erreur  svp trouver mon erreur voici le programme  par la méthode décompositionnellle d'Adomian  merci
restart;
> U:=unapply(Diff( u,u$2)+r^(-1)*Diff(u,r)=c^(-2)*Diff(u,t$2),r,t);
> #ceci est un comentaire
> print(`les condition initiales`);
> F:=f(r)=u(r,0);
> V:=Diff(u(r,0),t)=c*g(r);
> #ceci est un comentaire
> print(`opérateur différentielles linéaires`);
> R1:=L[r]=Diff(u(r,t),t);
> T:=L[t,t]=Diff(u(r,t),t$2);
> R:=L[r,r]=Diff(u(r,t),r$2);
> #ceci est un comentaire
> print(`opérateur d'intégration double`);
> H:=L[t,t]^(-1)=int(int(u(r,t),t=0..t),t=0..t);
> Q := u(r,0)+c*g(r)*t;
> u1(r,t):=unapply(Q+c^2*((L[t,t]^(-1)*L[r,r]*u(r,t))+(L[t,t]^(-1)*(r^(-1))*L[r]*u(r,t))),r,t);u=Sum('u[n](r,t)','n'=0..infinity);
> E:=u[0]=Q;
> Z:=u[1]=unapply(c^2*((L[t,t]^(-1)*L[r,r]*E)+(L[t,t]^(-1)*(r^(-1))*L[r]*E)));
> Y:=u[2]=c^2*((L[t,t]^(-1)*L[r,r]*Z)+(L[t,t]^(-1)*(r^(-1))*L[r]*Z));
> sol:= u[k+1]=unapply(a^2*((L[t,t]^(-1)*L[r,r]*u(k))+L[t,t]^(-1)*(r^(-1)*L[r]*u(k))));
quand je donne des valeur à F et g(r) et l'es réaffiche
et je veut tracer en prenent différente valeur de c=4,5,6,7 en dim 3d avec F=r; et g(r)=1;

Posté par
jandri Correcteur
re : programme de l'équation de la corde vibrante avec maple 24-06-09 à 08:57

Bonjour,

Il y a une erreur dès la première ligne.
Tu dérives u par rapport à u.

Posté par
joug
re : programme de l'équation de la corde vibrante avec maple 24-06-09 à 12:33

Bonjour,
merci de vore réponse
mais le programme ne veut pas s'executer car je veut donner des valeur a F et g(r) pour calculer la somme

Posté par
jandri Correcteur
re : programme de l'équation de la corde vibrante avec maple 24-06-09 à 16:16

Le programme est incompréhensible.
Quand on écrit F:=f(r)=u(r,0) cela affecte à la variable F l'égalité f(r)=u(r,0) mais cela ne définit pas une fonction f.
De plus F n'intervient même pas dans la suite du programme!

Posté par
joug
re : programme de l'équation de la corde vibrante avec maple 24-06-09 à 18:38

je voudrai faire le programme suivant :
((d²u)/(dr²))+((du)/(rdr))=((d²u)/(c²dt²))            r≥0,t≥0
    les conditions initiales suivantes :
u(r,0)=f(r)  
((d )/(dt))u(r,0)=cg(r)  
L_{tt}u(r,t)=c²[L_{rr}u(r,t)+(1/r)L_{r}u(r,t)] tq L_{tt}=d²u/dt² et  ainsi pour L_{rr} et L_{t}
et je multiplie PAR l'integrale d'inverse double L_{tt}⁻¹=∫∫(.)dtdt de 0 à t des deux coté
ça domme u(r,t)=Q+c²[L_{tt}⁻¹L_{rr}u(r,t)+L_{tt}⁻¹(1/r)L_{r}u(r,t)]
tq Q=f(r)+c*t*g(r)
pour avoir
uO = Q
u₁ = c²[L_{tt}⁻¹(L_{rr}u₀)+L_{tt}⁻¹((1/r)L_{r}u₀)]
u₂ = c²[L_{tt}⁻¹(L_{rr}u₁)+L_{tt}⁻¹((1/r)L_{r}u₁)]
u_{k+1} = c²[L_{tt}⁻¹(L_{rr}u_{k})+L_{tt}⁻¹((1/r)L_{r}u_{k})]       ,k≥0  
aprés je donne pour f(r) la valeur r et g(r)=1

Posté par
joug
problème de procedure ds maple 25-06-09 à 13:25

Bonjour à tous aider moi pour ce programme
moder:=proc(F,G)
> local U,V,U1,R1,T,H,E,R,Q,k,u,sol;
> restart;
> U:=unapply(Diff( u,t$2)+r^(-1)*Diff((u,r),r),r,t)=c^(-2)*Diff(u,t$2);
> #ceci est un comentaire
> print(`les condition initiales`);
> F:=f(r)=u(r,0);
> G:=g(r);
> V:=Diff(u(r,0),t)=c*g(r);
> #ceci est un comentaire
> print(`opérateur différentielles linéaires`);
> R1:=L[r]=Diff(u(r,t),t);
> T:=L[t,t]=Diff(u(r,t),t$2);
> R:=L[r,r]=Diff(u(r,t),r$2);
> #ceci est un comentaire
> print(`opérateur d'intégration double`);
> H:=L[t,t]^(-1)=int(int(u(r,t),t=0..t),t=0..t);
> Q := F+c*G*t;
> U1:=unapply(Q+c^2*((L[t,t]^(-1)*L[r,r]*u)+(L[t,t]^(-1)*(r^(-1))*L[r]*u)),r,t);
> u=Sum('u[n](r,t)','n'=0..infinity);
> E:=u[0]=Q;
> for k from 1to K+1 do
> sol:= u[k+1]=unapply(a^2*((L[t,t]^(-1)*L[r,r]*u(k))+L[t,t]^(-1)*(r^(-1)*L[r]*u(k))));
> end do;
> end proc ;
comment faire executer cette procedure je veut donner des valeur à F=r et G=1 merci

Posté par
gui_tou
re : programme de l'équation de la corde vibrante avec maple 25-06-09 à 14:04

salut,

et si tu donnais un schéma et des équations ?

Posté par
joug
re : programme de l'équation de la corde vibrante avec maple 25-06-09 à 15:51

voici l'équation de la corde vibrante:
^2u/t^2+(1/r)*u^2/r=c^2*u/t^2 (3.1)
les condition initiales sont:
u(r,0)=f(r)
u(r,0)/t=c*g(r)
Où u(r,t) représente le déplacement de trouver une particule au point  r dans le temps t instantané, c est  la vitesse de propagation de l'onde u. Nous considérons l'éq  (3.1) telque:
L[tt]=c^2*(L(r,r)+(1/r)*L(r))
ou L  =^2u/t^2
   L(r,r)  =^2u/r^2
   L(r)   =^2u/r
en multiplie l'inversed'integration double
L[t,t]^(-1)=((.)dtdt de 0 àt pour les deux
on aura u(r,t)=[t]+c^2*(L[r,r]L[t,t]^(-1)u(r,t)+(1/r)*L[r]*L[t,]^(-1)u(r,t))
[t]=u(r,0)+u[r]*t
                          =f(r)+c*t*g(r)
on aboutit a u(r,t)=u(r,t)n de 0à l'infinie
u[0]=[t]
u[1]=[t](c^2*(L[r,r]L[t,t]^(-1)u[0]+(1/r)*L[r]*L[t,]^(-1)u[0]))
u[2]=c^2*(L[r,r]L[t,t]^(-1)u[1]+(1/r)*L[r]*L[t,]^(-1)u[1])
u[k+1]=c^2*(L[r,r]L[t,t]^(-1)u[k]+(1/r)*L[r]*L[t,]^(-1)u[k])

quand je met L[t,t]c Len desous tt
est L[t,t]^(-1) c'est l'inverse
merci

Répondre à ce sujet

Seuls les membres peuvent poster sur le forum !

Vous devez être connecté pour poster :

Connexion / Inscription Poster un nouveau sujet
Une question ?
Besoin d'aide ?
(Gratuit)
Un modérateur est susceptible de supprimer toute contribution qui ne serait pas en relation avec le thème de discussion abordé, la ligne éditoriale du site, ou qui serait contraire à la loi.


Vous devez être membre accéder à ce service...

Pas encore inscrit ?

1 compte par personne, multi-compte interdit !

Ou identifiez-vous :


Rester sur la page

Désolé, votre version d'Internet Explorer est plus que périmée ! Merci de le mettre à jour ou de télécharger Firefox ou Google Chrome pour utiliser le site. Votre ordinateur vous remerciera !