Hypothèse basée sur Fourier - page 10

Nouveau commentaire
 

Voici mon code C, daté du 03.1995 (quel désordre poussiéreux dans la mezzanine !) S'il n'a pas été effacé, il doit fonctionner.


int Invert (double * t)
/*---------------------------------------------------------*
 *   Вычисление обратной матрицы путем ее перевода   	   *
 *   к треугольному виду  с последующим умножением 	   *
 *      на преобразованную единичную матрицу.        	   *
 *   Обратная матрица распологается на месте исходной.	   *
 *   Требуется <buf> область размерности t[1..NM*NM] 	   *
 *---------------------------------------------------------*
 *  t[0] =  NM  - размерность матрицы.               	   *
 *  t[1..NM*NM] - исходная (обратная) матрица по столбцам. *
 *  e[1..NM*NM] - элементы единичной матрицы.       	   *
 *  ВОЗВРАТ :  0 - OK                           	   *
 *---------------------------------------------------------*/
{   int    js, jk, jj, jp, n= t[0];
    double a1, a2, * e=(double *) buf;
//  -----------  Заполнение Е матрицы  -----------
    for ( jk = 1; jk <= n; jk++) {
       jp = n*( jk-1);
       for ( js = 1; js <= n; js++) {
	  e[ js+ jp] = t[ js+ jp];  t[ js+ jp] = 0.0; }
       t[ jk+ jp] = 1.0; }
// ----------------  Проход по столбцам ----------
    for ( jk = 1; jk < n; jk++) {
//  ---------------  Обнуление столбца  ----------
       a1 = e[( jp = n*( jk-1))+ jk];  //  Диагональный элемент
       for ( js = jk+1; js <= n; js++) {
	    a2 = e[ jp+ js]/ a1;   e[ jp+ js] = 0.0;
	    for ( jj = jk+1; jj <= n; jj++)
	       e[ n*( jj-1)+ js] -= a2* e[ n*( jj-1)+ jk];
	    for ( jj = 1; jj <= n; jj++)
	       t[ n*( jj-1)+ js] -= a2* t[ n*( jj-1)+ jk]; }
    } //  END for (JK)
//  -------  Определение обратной матрицы  -------
    for ( js = n; js > 0; js--)
       for ( jk = 1; jk <= n; jk++) {
	  for ( a1 = 0.0, jj = js+1; jj <= n; jj++)
	     a1 += e[ n*( jj-1)+ js]* t[ n*( jk-1)+ jj];
	  t[ n*( jk-1)+ js] = ( t[ n*( jk-1)+ js]- a1)/ e[ n*( js-1)+ js]; }
    return 0;
}
/**********************************************************************/
[Supprimé]  
Mathemat >> :

Bien sûr, c'est une méthode lente. Je me demande comment vous avez obtenu quelque chose pour une matrice de 100x100.

>> Je me joins à eux. Je suis aussi un peu surpris.

montre (de manière réfléchie) la grande efficacité du compilateur MQL4.

 

à Mathemat, AlexEro

Chers collègues, ne faites pas la fine bouche, le premier algorithme que nous avons vu a été mis en œuvre. Nous allons l'accélérer, l'approfondir et l'améliorer.

[Supprimé]  
grasn >> :

PS : Une correction. Si c'est 1:1 et que la ligne d'entrée est citée, c'est assez bizarre. Si l'image est stable, c'est déjà très bizarre.

Oui, 1:1. J'ai trouvé une coquille. Voici le résultat.


Hélas, on est encore loin d'une prédiction fiable...

[Supprimé]  
grasn >> :

à Mathemat, AlexEro

Chers collègues, ne faites pas la fine bouche, le premier algorithme que nous avons vu a été mis en œuvre. Nous allons l'accélérer, l'approfondir et l'améliorer.



Vraiment, collègue ? ! On ne s'en prend pas à vous ! Où avez-vous vu ça ? Nous ne faisons que souligner le râteau que nous avons nous-mêmes touché il y a longtemps. Et nous vous montrons comment les contourner. Tu préfères ne rien dire ?

 

à equantis

Je vous ai prévenu sur la page précédente en grosses lettres. Vous pouvez passer des années et des années sans trouver de solution acceptable, ou vous pouvez tout résoudre en 10 minutes. Mais je recommande quand même de le laisser pour le dessert :o).


à AlexEro

Je vous rappelle juste que jusqu'à ce moment-là, il n'y avait pas de bibliothèque du tout, sauf le travail d'Ilnur:o)

[Supprimé]  
grasn >>:А я предупреждал на прошлой странице крупными буквами. Можно годы потрать, но не найти приемлемого решения, а можно и за 10 мин во всем разобраться. Все же рекомендую оставить на десерт :о)

C'est vrai, ça vaut le coup de garder...

Afin de clore la question, je voudrais montrer une image avec des résultats typiques :

Bleu - prix

Rouge - prévision avec transformation cosinus avec départ à 0

violet - la même courbe, mais calculée à partir du point de départ de la prévision (100)

Vert - prévision simple basée sur la courbe des prix (j'ai utilisé la fonction predict intégrée)

[Supprimé]  
grasn >> :

à AlexEro

Je rappelle juste que jusqu'à présent, il n'y avait pas de bibliothèque du tout, sauf le travail d'Ilnur :o)

vous ne l'aviez pas....

 
AlexEro >> :

vous ne l'aviez pas....

On peut se demander, et où étiez-vous ? :о))))))

1...34567891011
Nouveau commentaire