問題タブ [chaos]

ダイナミクス用の Fortran プログラム (基本的には verlet アルゴ) があります。速度をより速く計算するために、アルゴリズムを MPI で並列化しました。私が緊張しているのは、4 つのプロセッサがある場合、各プロセッサが Verlet を実行し、並列化のポイントに到達すると情報を共有することです。ただし、わずかな数値の違い (たとえば、各ノードでコンパイルされた LAPACK) により、各 Verlet 軌跡は長期的には完全に異なる方向に展開する可能性があります。軌跡。したがって、発散を防ぐためにすべての時間ステップで情報を同期することにしましたが、これは明らかに障壁になります。

この問題 (ノードの発散) は通常どのように解決されますか? 参照はありますか？

私は混沌としたアトラクタを扱っており、いくつかの連続 - > 離散等価性をテストしています。この方法で、ロスラー システムの連続シミュレーションを作成しました。

さて、RSolve で離散等価システムを評価しようとすると、Mma は何もせず、エラーも出さず、解けないだけです。

NDSolve が DSolve の場合と同様に、RSolve の数値関数があるかどうかを知りたいです。For[] サイクルで計算できることはわかっていますが、そのような関数が存在するかどうかを知りたいだけです。

おそらく混沌とした時系列の予測可能性を推定したいと思います。このため、リアプノフ指数が良い候補になると思いました。私はインターネットでいくつかの調査を行いましたが、私が見つけた説明は (少なくとも私にとっては) Java で実装するのに十分ではありません。

時系列のリアプノフ指数推定の詳細なアルゴリズム/疑似コードをどこで入手できるか知っていますか? あるいは、時系列の予測可能性の他の測定値も歓迎します。

ありがとうございました！

このアドレスに投稿された DNA シーケンスのカオス ゲームを作成するための Mathematica コードを試しました: http://facstaff.unca.edu/mcmcclur/blog/GeneCGR.html

これは次のようなものです:

私が持っている fasta シーケンスは AACCTTTGATCAAA のような文字のシーケンスであり、生成されるグラフは次のようになります。

コードは小さなシーケンスでは問題なく動作しますが、たとえば約 40Mb の染色体など、巨大なシーケンスを配置したい場合、プログラムは多くの時間を要し、黒い四角が表示されるだけなので分析できません。前述のコードを改善して、表示される正方形が大きくなるようにすることはできますか?ところで、正方形は正方形単位のみである必要があります。事前にご協力いただきありがとうございます

ローレンツシステムの最大のLyaponov指数を計算するためにこの手順ode45()を使用して、ダイナミクスのいくつかの作業にmatlabで使用しています。

このプロセスでは、から始まる微分方程式のシステムを解き、x0これをに非常に近い軌道と比較しx0ます。

ode45() 各タイムステップで、タイムステップを進める前に2番目の軌道を再調整する必要があるので、一度だけ呼び出すことができるようにしたいと思います-これは可能ですか？

代替ソリューションの開始の試みがここに示されていますが、機能しません。私が見る限り、以下で得られる行列rとR取得は類似しているはずです。

...しかし、そうではありません：

の最後の行rがゼロであることは問題ではありません。

何か案は？乾杯！\ T

キーワードを取得し、CSS で強調表示/色付けするプログラムがあります。

最終的に強調表示される単語を変更したいので、jQueryを使用しています。すべてユーザーが行います。

現時点での問題は、何も強調表示されていないことではなく(yay)、代わりにすべてが強調表示されていない(ああ、いや)ということではありません。いくつかの単語を通過して強調表示するように設定していますが(バルーン、ドロシー、サイクロン)、特定の単語のすべてが色付けされるわけではありません。

このjsFiddleのことが、私ができるよりもうまく説明できることを願っています。2 行目を見ると、最初の単語とキーワードである Dorothy が思ったほど赤くないことがわかります。

ありがとう。支離滅裂で申し訳ありませんが、フィードバックをお待ちしております

乱数ジェネレーターの正直さを保証するために、ユーザーは、必要に応じて、その数値が実際にエントロピーの公開ソースから生成されていることを確認できるという考え方です。これにより、システムは、サーバーが乱数を選択できなかったことをユーザーに確認できます。

私の質問は次のとおりです。文字列に連結してmd5にフィードできる、エントロピー（好ましくは不変で混沌とした）で絶対的に参照可能ないくつかの優れた公開ソースは何ですか？いくつかのアイデアは、特定の株価、温度（正直な情報源から）、ビットコインブロックチェーンに含まれるハッシュです...

プログラムは次のことを行うことになっています。

1. 0に等しい値のグリッド「ActualGridValues」（X行×Y列）を初期化します。
2. 別の値のグリッド「RandomGridValues」（X行×Y列）を作成し、グリッド内の各セルを-ZからZの範囲のN個の乱数の平均に等しいものとして定義します。
3. RandomGridValuesの数値をActualGridValuesの対応するセルに追加します。
4. ActualGridValuesのセルごとに、隣接するセルの値も追加されます。
5. 手順2〜4を無限に繰り返します。

その結果、ランダムな構造がグリッド内で合体します。高い正の値のグループ化は近くのセルに影響を与え、それらを高く歪ませます。低い負の値のグループ化は同じことを行いますが、低く歪めます。

すべてのセルのすべての値の合計積はゼロである必要があります。高スキュー数と低スキュー数のローカライズされたグループがありますが、十分に大きいサンプルサイズでは、すべてのセルの合計はゼロである必要があります。

問題は、たとえば1000回の反復の後、値が一貫して負に歪むことです。構造があり、高点と低点がローカライズされていますが、合計値は常に負に偏っています。（これは、時間の経過とともにグリッド全体が負の数のみで満たされることを意味します）。

シミュレーションが実行されるたびに、値は負に偏ります。理由について何か考えはありますか？

編集：私は次の機能に問題を切り分けました。RandomGridValueのすべての数値の平均値は、ほとんどの場合負になります。