Fourier-basierte Hypothese - Seite 10

Neuer Kommentar
 

Dies ist mein C-Code, datiert auf 03.1995 (was für ein staubiges Durcheinander im Zwischengeschoss!) Wenn er nicht gelöscht wurde, muss er noch funktionieren.


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;
}
/**********************************************************************/
 
Mathemat >> :

Das ist natürlich eine langsame Methode. Ich frage mich, wie Sie etwas für eine 100x100-Matrix bekommen haben.

>> Ich schließe mich an. Ich bin auch ein wenig überrascht.

zeigt (nachdenklich) die hohe Effizienz des MQL4-Compilers.

 

an Mathemat, AlexEro

Kolleginnen und Kollegen, seien Sie nicht pingelig, der erste Algorithmus, den wir gesehen haben, wurde umgesetzt. Wir werden sie beschleunigen, vertiefen und verbessern.

 
grasn >> :

PS: Eine Korrektur. Wenn 1:1 und die Eingabezeile in Anführungszeichen steht, ist das ziemlich seltsam. Wenn das Bild stabil ist, dann ist es schon sehr seltsam.

Ja, 1 zu 1. Ich habe einen Tippfehler gefunden. Hier ist das Ergebnis.


Von einer zuverlässigen Vorhersage sind wir allerdings noch weit entfernt...

 
grasn >> :

an Mathemat, AlexEro

Kolleginnen und Kollegen, seien Sie nicht pingelig, der erste Algorithmus, den wir gesehen haben, wurde umgesetzt. Wir werden sie beschleunigen, vertiefen und verbessern.



Wirklich, Kollege?! Wir hacken nicht auf Ihnen herum! Wo haben Sie das gesehen? Wir weisen nur auf die Harke hin, die wir selbst schon vor langer Zeit geschlagen haben. Und wir zeigen Ihnen, wie Sie sie umgehen können. Möchten Sie lieber gar nichts sagen?

 

zu equantis

Ich habe Sie auf der vorherigen Seite in großen Buchstaben gewarnt. Sie können Jahre damit verbringen, keine akzeptable Lösung zu finden, oder Sie können alles in 10 Minuten herausfinden. Aber ich empfehle trotzdem, es beim Dessert zu belassen :o).


an AlexEro

Ich erinnere Sie nur daran, dass es bis zu diesem Zeitpunkt überhaupt keine Bibliothek gab, außer Ilnurs Arbeit :o)

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

Das stimmt, es lohnt sich, es zu behalten...

Um das Thema abzuschließen, möchte ich ein Bild mit typischen Ergebnissen zeigen:

Blau - Preis

Rot - Vorhersage mit Kosinustransformation mit Start bei 0

violett - dieselbe Kurve, aber berechnet vom Ausgangspunkt der Prognose (100)

Grün - einfache Vorhersage auf der Grundlage der Preiskurve (ich habe die integrierte Vorhersagefunktion verwendet)

 
grasn >> :

an AlexEro

Ich erinnere nur daran, dass es bis zu diesem Zeitpunkt überhaupt keine Bibliothek gab, außer Ilnurs Arbeit :o)

du hattest es nicht....

 
AlexEro >> :

du hattest es nicht....

Es bleibt die Frage: Wo warst du? :о))))))

1...34567891011
Neuer Kommentar
Grund der Beschwerde: