Discussion :

[the_reward]

bonjour

je veux écrire un programme en utilisant les fichiers pour tracer la fonction gaussienne (c la solution initiale de mon problème ) quand je plote la fonction je trouve pas la bonne allure par contre c quand j'écris la formule de la fonction f(x) = exp((-(x+l)**2)/(2.D0))/sqrt((2.D0)*pi) je trouve la bonne allure je vois pas c quoi le problème.

Si quelqu'un pourrait m'aider je serai vraiment reconnaissante

merci d'avance!

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
 
program solution
  implicit none
 
  real, dimension(:),allocatable :: x
  integer :: i,n
  real::l,xdeb,xfin,dx,pi
 
  open(unit=10, file='parametre')
  read(10,*) n,xdeb,xfin
  close(10)
 
  allocate(x(n))
 
  l=xfin-xdeb
  dx = l*1.D0/(n+1)
  pi = 4.D0*datan(1.D0)
  do i=1,n
     x(i) = (1.D0/sqrt(2.D0*pi))*exp((-((i*dx)-(l/2.D0))*((i*dx)-(l/2.D0)))/2.D0)
  end do
 
  open(unit=10, file='solution_ini.dat')
  do i = 1,n
     write(10,*) x(i)
  end do
  close(10)
 
end program solution

et pour les parametres j'ai pris 100 0.D0 1.D0

[__dardanos__]

Salut.
C'est parce que les abscisses sont différents. Je retrouve la même courbe avec :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

n=100 ; xdeb=-10.d0 ; xfin = 10.d0

Vérifie si l'argument de l'exponentiel est correct (facteur pi manquant ?)

[the_reward]

salut __dardanos__

j'ai mis pi = 4.D0*datan(1.D0) (ligne 16 du prog) mais pour quoi ça marche pas avec xdeb = 0 et xfin = 1?!!!

[__dardanos__]

On peut réécrire la boucle :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
    do i = 1, n
        absc = i * dx - l / 2.D0
        x(i) = 1.D0 / sqrt(2.D0*pi) * exp(-.5d0 * absc**2)
    end do

La forme générale d'une gaussienne étant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

 1.D0 / sqrt(2.D0*pi) / sigma * exp(-.5d0 * (absc-mu)**2 / sigma**2)

Tu as pris une gaussienne avec un écart type sigma de 1. Elle est centrée, mu=0 (elle est maximale en zéro). Rappelle toi que l'écart type est lié à la largeur à mi-hauteur.
Avec xdeb=0 et xfin=1, x varie de -0,49 à +0,49. Comme sigma=1, on a toujours |x|<sigma/2. D'où la courbe qui varie peu.
Avec xdeb=-10 et xfin=+10, x varie de -9,80 à +9,80. x prend des valeurs telles que |x|>>sigma. On retrouve la forme en cloche présentée dans les manuels.

[the_reward]

merci __dardanos__