Forum de mathématiques - Bibm@th.net

Bonjour Yoshi,

Il y a pas de différence entre les deux écritures mais j'ai écrit exp(-(x-l/2)) pour poser une question sur la périodicité de la gaussienne

math@

pour le 1er point, oui je souhaite que mon abcs couvre l'intervalle [xbed,xfin] et qu'elle soit centrée au milieu de cet intervalle.

et pour la gaussienne dans mon cas elle n'est pas paire (j'ai exp(-(x-l/2))!!!

Pour la condition initiale je peux pas mettre que u=0? car j'ai un système de 2 équations et 2 inconnus W([tex]\zeta[/tex],u)!

Bonjour,

Voilà le code que j'ai implémenté. je travaille sur les équation de green naghdi, le but c'est résoudre un système de type tAAX=tACY ou tA est la transposée de A et  C est une matrice aussi
l'erreur c dans le calcul de la solution finale et je comprends pas d'où vient le problème

http://pastebin.com/5kkyec1t
http://pastebin.com/YhbynDPi
http://pastebin.com/WyAmsexR
http://pastebin.com/7MdG9w6p
http://pastebin.com/LnZeLRLG
http://pastebin.com/9iidjXCi
http://pastebin.com/1twEPtSt
le système que je veux résoudre
[tex]
$\left\{\begin{array}{ll}
\partial_{t}\zeta + h\partial_{x}u = 0\\\\
\left( I-\dfrac{\alpha}{3}h^{2}\partial_{xx}\right) \left( \partial_{t}u+\left( 1-\dfrac{1}{\alpha}\right) g\partial_{x}\zeta\right) +\dfrac{g}{\alpha}\partial_{x}\zeta = 0\\
\end{array}\right.$
[/tex]
le schéma numérique
[tex]

$\left\{\begin{array}{ll}
\zeta_{i}^{n+1} = \zeta_{i}^{n} - h\dfrac{\delta t}{4\delta x}(u_{i+1}^{n+1}-u_{i-1}^{n+1} + u_{i+1}^{n}-u_{i-1}^{n})\\\\

\dfrac{\alpha}{12}\left(1-\dfrac{1}{\alpha}\right)g \zeta^{n+1}_{i-2}-\left(\left(1-\dfrac{1}{\alpha}\right)g\dfrac{\delta x^2}{4h^2}+\dfrac{\alpha}{6}\left(1-\dfrac{1}{\alpha}\right)g+\dfrac{g\delta x^2}{4\alpha h^2}\right)\zeta^{n+1}_{i-1}+\\\\ \left(\left(1-\dfrac{1}{\alpha}\right)g\dfrac{\delta x^2}{4h^2}+\dfrac{\alpha}{6}\left(1-\dfrac{1}{\alpha}\right)g+\dfrac{g\delta x^2}{4\alpha h^2}\right)\zeta^{n+1}_{i+1}-\dfrac{\alpha}
{12}\left(1-\dfrac{1}{\alpha}\right)g \zeta^{n+1}_{i+2}\\\\-\dfrac{\alpha\delta x}{3\delta t}U^{n+1}_{i-1}+\left(\dfrac{\delta x^3}{\delta t h^2}+\dfrac{2\alpha\delta x}{3\delta t}\right)U^{n+1}_{i}-\dfrac{\alpha\delta x}{3\delta t}U^{n+1}_{i+1} =

-\dfrac{\alpha}{12}\left(1-\dfrac{1}{\alpha}\right)g \zeta^{n}_{i-2}+\\\\
\left(\left(1-\dfrac{1}{\alpha}\right)g\dfrac{\delta x^2}{4h^2}+\dfrac{\alpha}{6}\left(1-\dfrac{1}{\alpha}\right)g+\dfrac{g\delta x^2}{4\alpha h^2}\right)\zeta^{n}_{i-1}-
\left(\left(1-\dfrac{1}{\alpha}\right)g\dfrac{\delta x^2}{4h^2}+\dfrac{\alpha}{6}\left(1-\dfrac{1}{\alpha}\right)g+\dfrac{g\delta x^2}{4\alpha h^2}\right)\zeta^{n}_{i+1}\\\\-\dfrac{\alpha}{12}\left(1-\dfrac{1}{\alpha}\right)g \zeta^{n}_{i+2}-\dfrac{\alpha\delta x}{3\delta t}U^{n}_{i-1}+\left(\dfrac{\delta x^3}{\delta t h^2}+\dfrac{2\alpha\delta x}{3\delta t}\right)U^{n}_{i}-\dfrac{\alpha\delta x}{3\delta t}U^{n}_{i+1}

\end{array}\right.$   

[/tex]

j'ai multiplié la 2ème équation par dx^3/h_0^2
les coeffs des matrices

[tex]
$a_1 = \dfrac{\alpha\delta t}{12}$\\\\
$a_2 = \left(1-\dfrac{1}{\alpha}\right)g $\\\\
$a = a_1*a_2$\\\\
$b_1 = \dfrac{g\delta t\delta x^2}{4h_0^2}$\\\\
$b_2 = \dfrac{\alpha\delta t}{6}$\\\\
$b = b_1+(b_2*a_2)$\\\\
$c = \dfrac{\alpha\delta x}{3}$\\\\
$d_1 = \dfrac{\delta x^3}{h_0^2}$\\\\
$d_2 = \dfrac{2\alpha\delta x}{3}$\\\\
$e = \dfrac{\delta t h_0}{4\delta x}$\\\\
[/tex]

la matrice A

[tex]
$ \left(\begin{array}{c|c}

\begin{matrix}
1 &0 &0 &0 &0 &\cdots &0\\
0 &1 &0 &0  &0 &\cdots &0 \\
0 &0 &1 &0 &0 &\cdots &0\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &0 &\cdots &0 &1 &0 \\
0 &0 &0 &0 &\cdots &0 &1

\end{matrix}

&
  \begin{matrix}
0 &e &0 &0 &\cdots &0 &-e\\
-e &0 &e &0  &0 &\cdots &0\\
0 &-e &0 &e &0 &0 &\cdots\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &\cdots &0 &-e &0 &e\\
e &0 &0 &\cdots &0 &-e &0

\end{matrix}\\

\hline

\begin{matrix}
0 &b &-a &0 &\cdots &a &-b\\
-b &0 &b &-a  &0 &\cdots &a\\
a &-b &0 &b &-a &0 &\cdots\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
-a &0 &\cdots &a &-b &0 &b\\
b &-a &0 &\cdots &a &-b &0

\end{matrix}
&

\begin{matrix}
d &-c &0 &0 &0 &\cdots &-c\\
-c &d &-c &0  &0 &0 &\cdots\\
0 &-c &d &-c &0 &\cdots &0\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &0 &\cdots &-c &d &-c \\
-c &0 &0 &0 &\cdots &-c &d

\end{matrix}

\end{array}\right)$
[/tex]

la transposée de A
[tex]
$ \left(\begin{array}{c|c}

\begin{matrix}
1 &0 &0 &0 &0 &\cdots &0\\
0 &1 &0 &0  &0 &\cdots &0 \\
0 &0 &1 &0 &0 &\cdots &0\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &0 &\cdots &0 &1 &0 \\
0 &0 &0 &0 &\cdots &0 &1

\end{matrix}

&
  \begin{matrix}
0 &e &0 &0 &\cdots &0 &-e\\
-e &0 &e &0  &0 &\cdots &0\\
0 &-e &0 &e &0 &0 &\cdots\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &\cdots &0 &-e &0 &e\\
e &0 &0 &\cdots &0 &-e &0

\end{matrix}\\

\hline

\begin{matrix}
0 &b &-a &0 &\cdots &a &-b\\
-b &0 &b &-a  &0 &\cdots &a\\
a &-b &0 &b &-a &0 &\cdots\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
-a &0 &\cdots &a &-b &0 &b\\
b &-a &0 &\cdots &a &-b &0

\end{matrix}
&

\begin{matrix}
d &-c &0 &0 &0 &\cdots &-c\\
-c &d &-c &0  &0 &0 &\cdots\\
0 &-c &d &-c &0 &\cdots &0\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &0 &\cdots &-c &d &-c \\
-c &0 &0 &0 &\cdots &-c &d

\end{matrix}

\end{array}\right)$
[/tex]

la matrice C

[tex]
$ \left(\begin{array}{c|c}
      \begin{matrix}
1 &0 &0 &0 &0 &\cdots &0\\
0 &1&0 &0  &0 &\cdots &0 \\
0 &0 &1 &0 &0 &\cdots &0\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &0 &\cdots &0 &1 &0 \\
0 &0 &0 &0 &\cdots &0 &1

\end{matrix}
 
&
\begin{matrix}
0 &-e &0 &0 &\cdots &0 &e\\
e &0 &-e &0  &0 &\cdots &0\\
0 &e &0 &-e &0 &0 &\cdots\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &\cdots &0 &e &0 &-e\\
-e &0 &0 &\cdots &0 &e &0

\end{matrix}\\

\hline
\begin{matrix}
0 &-b &a &0 &\cdots &-a &b\\
b&0 &-b &a  &0 &\cdots &-a\\
-a &b &0 &-b &a &0 &\cdots\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
a &0 &\cdots &-a &b &0 &-b\\
-b &a &0 &\cdots &-a &b &0
\end{matrix}

&
  \begin{matrix}
d &-c &0 &0 &0 &\cdots &-c\\
-c &d &-c &0  &0 &0 &\cdots\\
0 &-c &d &-c &0 &\cdots &0\\
\vdots &\ddots &\ddots &\vdots\\
\vdots &\ddots &\ddots &\vdots\\
0 &0 &0 &\cdots &-c &d &-c \\
-c &0 &0 &0 &\cdots &-c &d

\end{matrix}

\end{array}\right)$
[/tex]
Si quelqu'un peut m'aider je serai vraiment reconnaissante.
Merci d'avance

Cordialement
math@

Bonjour tout le monde,

Je suis vraiment désolée de ne pas pouvoir répondre plus tôt à vos message, ma matrice est un peu "spéciale" je peux pas prendre une taille petite le problème qui n’embête c'est la condition initiale qui est une gaussienne normalement je dois avoir une gaussienne sur un intervalle [-5,5] centré en 0 mais je trouve pas ça et même pour les valeurs des abscisse je trouve x= 50, 100 ...
Volà la partie du code de la solution initiale j'espère avoir des réponse merci d'avance

Bonjour tout le monde,
merci beaucoup pour votre aide et excusez moi car j'ai pas répondu j'étais en pleine révision pour les examens.
je suis débutante en programmation et en python mais je dois l'apprendre car la programmation est un outil indéspensable pour ma spécialité. En fait pour la fonction polyno je dois donner des arguments arbitrairessoit (1) des coefficients d'un polynome (2) une chaine de caractère( un polynome) et le résultat sera tjr la liste des coefficients j'ai réussi à faire la 1ère partie le problème c que comment je vais faire ces deux chose en une seule fonction ce que je veux dire c que
je défini une fonction dont les arguments sont arbitraie :def poly(*arg):
pour le 1er cas j'ai réussi à le faire mais pour le second je ne sais pas comment dire que 'arg' est un polynome et faire le teste pour avoir seulement les coefficients de ce polynome.
Je m'excuse une autre fois et merci beaucoup pour votre aide à la prochainne