フーリエベース仮説 - ページ 10

新しいコメント
 

これは私のCコードで、1995年03月のものです（なんと中二階で埃をかぶっています！）消されていないのなら、きっと動いているのでしょう。


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 >> :

もちろん、ゆっくりした方法です。100×100のマトリックスで、どうやって何かを手に入れたのか不思議です。

>> 私も参加します。私も少し驚いています。

(思わせぶりに)MQL4コンパイラの効率の良さをアピールしています。

 

to数学、AlexEro

同僚よ、甘えるな！最初に見たアルゴリズムが実装されて いた。スピードアップし、深化させ、向上させる。

 
grasn >> :

追記：訂正です。1:1で入力列が引用されると、かなり変なことになります。絵が安定していると、もう本当に気持ち悪いんですよ。

はい、1：1です。 誤字がありました。その結果がこちらです。


残念ながら、信頼できる予測にはまだ程遠いのですが...。

 
grasn >> :

to 数学、AlexEro

同僚よ、甘えるな！最初に見たアルゴリズムが実装されていた。スピードアップし、深化させ、向上させていきます。



本当に、同僚か！？いじめてるんじゃないよ！どこで見たんだ？私たちは、自分たちが昔打った熊手を指摘しているだけです。そして、それらを回避する方法を紹介します。何も言わない方が良いですか？

 

toequantis

前のページで大きな文字で警告 しました。何年も何年もかけても納得のいく解決策が見つからないかもしれないし、10分ですべてを理解できるかもしれない。でも、やっぱりデザートに残しておくのがおすすめです:o)


toAlexEro

その瞬間まで、イルヌールの 作品以外、図書館は全くなかったことを思い出してほしいんです :o)

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

そうです、残す価値があるのです...。

この問題を解決するために、典型的な結果を写真でお見せしたいと思います。

ブルー - 価格

赤 - 開始点を0とするコサイン変換による予測

紫 - 同じ曲線だが、予測の開始点(100)から計算される。

緑色：価格カーブに基づく単純な予測（内蔵の予測機能を使用しました）

 
grasn >> :

to AlexEro

私は、この瞬間まで、イルヌールの作品以外のライブラリーは全くなかったことを思い出しただけです :o)

持っていなかったのか......。

 
AlexEro >> :

持っていなかったのか......。

という問いかけが残されており、どこに行っていたのか？:о))))))

1...34567891011
新しいコメント