Absorption et Emission de photons

Ok, mais le problème n'est toujours pas résolu...   

J'ai relu et comparé caractère par caractère le code sur le forum et sur la calculatrice, je ne trouve aucune différence cette fois. Et le test sur calculatrice fonctionne tel que mes deux exemples.

code de test SEULEMENT l'absorption:

EDIT 3 :
If 0.001 > Ans
Une erreur et un bug d'affichage du forum corrigé.

EDIT 2 : j'ajoute ça pour aider quelqu'un qui passerait plus tard et aurait les mêmes difficultés que Pito2901.

affs:

EDIT 1 :

résultat:

Pourquoi "seulement l'absorption, de toute façon le résultat reste le même, et même avec ce bout de code le problème persiste...

(j'ai recopié le programme en entier 3 fois...)

Ce bout de code aurais dû permettre de chercher une erreur dans la zone concernée au lieu de scruter d'autres parties du programme qui ne sont pas en cause.

J'ajoute avec le programme test de l'absorption quelques tests montrant toutes les valeurs, en espérant cibler l'origine du problème.

Ok, merci.

avec 266, j'obtiens:

Résultats:

Avec 609:

Résultats:

Voilà, qu'est ce qu'on peut faire avec ça ???

EDIT: En fait je viens de m'apercevoir que j'obtiens les mêmes résultats, comment faire pour afficher juste le résultat voulu maintenant ?

La première valeur: 4.667400572 correspond à celle recherchée.
Ensuite il y a:
4.67 //La valeur en eV dans ce cas
5.566227856E-4 //Plus cette valeur est proche de 0 et plus le résultat est probable.
{0 1} //Le cas en question

Donc tu as 2 solutions: soit tu compares les résultats, soit tu cherche la plus petite valeur avant l'affichage du cas (ici il s'agit de 5.566227856E-4 donc du cas {0 1}).

Pour 609 il s'agit de 6.692488141E-4, soit le cas {1 4}


Pour afficher juste le cas voulu il faut revenir aux algo précédents, avec juste l'affichage de L2-1 à la fin.

Dac, et donc que faut-il changer dans le dernier bout de code de Linka pour obtenir le L2-1 ?

Essaye donc:

Code:

1241,5285520433/L->E
For(B,1,Z-1
For(C,B+1,Z
If 0.001>abs(1-E/abs(L1(C)-L1(B
Pause {B-1,C-1
End
End

Ca affichera tous les résultats plausibles normalement.
Tu peux aussi changer le 0.001 par la valeur qui te convient pour faire une recherche plus ou moins sélective.

Ca devrait marcher non?

ça marche du feu de dieu !

Par contre je viens de le réadapter sur le code initial, et ben pour 609, il m'affiche rien ! ! !    HHHAAAAA ! ! ! !

J'ai ce code là, avec 266 et 2177 ça marche mais pour 608.59... NON ! ! !

Code:

EffEcr
Input "NB ETATS EXC:",Z
If Z<6
Then
If Z>=1
Input "ETAT 0:",A
Input "ETAT 1:",B
If Z>=2
Input "ETAT 2:",C
If Z>=3
Input "ETAT 3:",D
If Z>=4
Input "ETAT 4:",E
If Z>=5
Input "ETAT 5:",F
If Z>=6
Input "ETAT 6:",G
Else
Disp "TROP !"
Stop
End
{A,B,C,D,E,F,G->L1
Lbl M
Menu("LAMBDA","CALCUL",A,"DEMANDER",B
Lbl A
EffEcr
Disp "POUR CHOISIR","L'ETAT 1, TAPEZ 1"
Input "1ER ETAT:",G
Input "2EME ETAT:",H
abs(L1(G+1)-L1(H+1->E
1241,5285520433/E->L
Disp "E(eV)=",E
Disp "LAMBDA(NM)=",L

Pause
Goto M
Lbl B
Input"LAMBDA(NM)=",L
1241,5285520433/L->E
Disp E
For(B,1,Z-1
For(C,B+1,Z
If .001>abs(1-E/abs(L1(C)-L1(B))
Pause {B-1,C-1
End
End
Pause
Goto M

Je viens de faire des essais, j'ai fais les calculs avec:
état 1 vers état 4 → 608.5924275
2 vers 4 → 682.1585451
3 vers 4 → 993.2228416
0 vers 4 → 185.0266098

Et quand je rentre les valeurs dans le input Lambda du lbl B, et bien ça ne marche pas, donc j'en conclus que dès que c'est un calcul vers l'état 4 (soit dans ce cas, c'est le dernier état excité) et bien ça ne marche pas.

EDIT: chose bizarre, il y a aussi de 0 vers 3, qui donne 227.3861817.

Tu n'as pas pris le bon code...
Remplace cette ligne:
If .001>abs(E-L1(C)+L1(B))
par:
If .001>abs(1-E/abs(L1(C)-L1(B

Ca devrait marcher...


(PS: Évites le triple post par contre.)

Désolé faute de recopiage (j'édite), j'ai bien mis ça dans ma Ti...  

PS: Désolé...

Et ça marche du coup ou pas?

non

Rajoute cette ligne qu'on puisse voir ce qui ne va pas (parce que ça devrait marcher je ne comprends pas...).

Code:

For(B,1,Z-1
For(C,B+1,Z
Pause abs(1-E/abs(L1(C)-L1(B //Ligne de test
If .001>abs(1-E/abs(L1(C)-L1(B
Pause {B-1,C-1
End
End

Ok donc voici les résultats:
état 1 vers état 4 → 608.5924275
.5631691649
.5828220859
.6263736264
8.272727272
1.582278481
2.578947368

2 vers 4 → 682.1585451
.6102783726
.627811861
.6666666667
7.272727272
1.303797468
2.192982456

3 vers 4 → 993.2228416
.7323340471
.7443762781
.7710622711
4.681818182
.5822784811
1.192982456

0 vers 4 → 185.0266098
.4368308354
.3721881394
.2289377292
29.50000001
7.493670888
10.77192983

0 vers 3 → 227.3861817 (quand je le mets en entier ça marche mais quand je rentre 227, non)
.1711539134
.1184639623
.0017012409
23.86040353
5.923150349
8.595243466

Voilà...

Mais il manque des valeurs dans ta liste d'affichages !

Test : Ajoute un affichage systématique de la paire. Pour ne pas perturber la lecture de tout le reste, je propose de le placer en haut à gauche de l'écran avec un Output.
Puis tu verras si toutes les paires sont analysées. (en tout il doit y en avoir 10)
Pense à affecter Z, car ce sont les paires d'indices finissant par 4 qui te manquent

code test:

=========

J'avais préparé un pavé sur les tolérances, mais comme il semble y avoir un autre problème, je le poserai plus tard.

=========

Confusion entre les types de comparaison...
Je résume les différentes opérations de comparaison pour ne plus les confondre :

Spoiler:

ah les bugs:

état 1 vers état 4 → 608.5924275
.5631691649
.5828220859
.6263736264
.695976155
8.272727272
1.582278481
4.569E-11
{1 4}

Puis 2.578947368 avec en haut à gauche : {2,3}
Puis .1208791208 avec en haut à gauche : {2,4}
Puis .6319999999 avec en haut à gauche : {3,4}

Tout ça juste pour 608.5924275

Est-ce que je fais les autres ?

Comment se fait-il que tu ais les output seulement à la fin ???
Si le problème venait de la boucle et non du test, ce ne serait pas la peine de refaire ce test pour d'autres valeurs de L.

Non tu as toutes les valeurs, maintenant essayes d'enlever les lignes de test:

Code:

5->Z
For(B,1,Z-1
For(C,B+1,Z
If .001>abs(1-E/abs(L1(C)-L1(B
Pause {B-1,C-1
End
End

Edit: J'ai saisis le problème:

Code:

If Z>=1
Input "ETAT 0:",A
Input "ETAT 1:",B
If Z>=2
Input "ETAT 2:",C
If Z>=3
Input "ETAT 3:",D
If Z>=4
Input "ETAT 4:",E
If Z>=5
Input "ETAT 5:",F
If Z>=6
Input "ETAT 6:",G
Else

Si on rentre 6 états, en réalité on en demande 7...

Donc corrige:

Code:

//Enlève le 5->Z
For(B,1,Z //Ici enlève le -1
For(C,B+1,Z+1 //Et ici rajoute le +1
If .001>abs(1-E/abs(L1(C)-L1(B
Pause {B-1,C-1
End
End

T'es un boss ! Tu as enfin résolu le problème. Félicitations à vous deux ! ! !

Sinon j'ai toujours le problème de 0 à 3 quand je mets seulement 227, si j'écris en entier le problème cesse. cela doit donc venir du 0.001 que dois je mettre pour qu'il soit plus large ?

EDIT: J'ai mis 0.01 et ça marche ! Youpi !

MERCI MERCI MERCI BEAUCOUP ! ! !

petite rétrospective : tu as eu à la fois :
-variable N au lieu de Z
-tolérance inadaptée
-mes erreurs dans les formules de tolérance
-la valeur Z inadaptée

Maintenant je sais pourquoi mes test du programme complet fonctionnaient et pas les tiens : je ne saisissais pas avec les Input A,B,C ... c'est de ma faute.

Toutefois j'avais signalé le problème de la saisie au tout début, et l'absence de réponse à ce propos, l'a faite tomber dans l'oubli.

============

Quant à la tolérance à choisir, je signale que dans certaines situations la tolérance sera trop grande devant les données, dans d'autres elle sera trop petite.
Il n'existe pas de valeur miracle de tolérance.

Si la longueur d'onde est un résultat de calcul, et que tu veux l'inclure dans le programme, tu as intérêt à stocker le résultat dans une variable directement, et saisir son nom dans l'input "LAMDA(NM):" lors de l'exécution de programme.
Et si tu as calculé la longueur d'onde avec l'émission dans le programme, tu peux simplement saisir "L" dans le champ Lamda de l'absorption.

Ce sont les chiffres significatifs qui comptent, pas le nombre de décimales, donc tu n'as pas forcément besoin de saisir tout ce que tu vois.

incertitudes et abandon de l'automate:

=================

Je poste des corrections complètes (ou je l'espère), et en adaptant Z lors de la saisie plutôt que dans la boucle de recherche.
Il manque quelques cas mais je fais confiance à ceux qui pourraient le lire, ils recouperont les cas pour adapter chaque morceau.

début : variables VS liste
début et fin : Z est SOIT le nombre d'états, SOIT l'indice maximal
fin : automatique VS comparaison manuelle

variables, nombre_etats, auto:

variables, indice_max, auto:

liste, nombre_etats, auto:

liste, indice_max, auto:

liste, nombre_etat, manuel:

Je crois que tu t'es trompé sur les 2 premiers codes:
Le premier il faut prendre
For(B,1,Z
For(C,B+1,Z+1
Et le second
For(B,1,Z-1
For(C,B+1,Z
(tu as fait l'inverse)

Euh... j'ai un problème mais pas celui là. Ou bien on ne compte pas la même chose. Je regarde encore ça.
Je voulais adapter le début et donner un sens à Z.
EDIT: maintenant ce devrais être mieux. Je teste.

En revanche j'avais mal étiqueté le premier spoiler, il s'agit d'une saisie de variables, pas d'une liste. Modifié.

Sinon, à un moment j'ai évoqué si on pouvait rendre l'affichage de {0 1} par:

ETAT 0 à ETAT 1 ?

et s'il n'y a rien, on afficherai: PAS DE TRANSLATION

Maintenant que tu as édité c'est bon.

Pour afficher état 0 à etat 1:
Au lieu de

Code:

Pause {B-1,C-1

Tu peux mettre:

Code:

//Avant les 2 for:
0->I

//Après le If (à la place du pause actuel):
Then
I+1->I
Output(I,1,"ETAT   A ETAT //3 espaces entre ETAT et A
Output(I,6,B-1
Output(I,15,C-1
If I=8
Then
Pause
EffEcr
0->I
End
End

quand je mets 266 il affiche:
ETAT 0 A ETAT61

En fait il faut mettre un espace après le 2ème ETAT

Il peut y avoir des caractères sur la 16ième colonne, et dans ce cas il serait préférable d'écrire 3 espaces après le second ETAT, ou bien d'effacer l'écran avant la boucle de recherche.

Je pencherai plutôt pour la deuxième option...