問題タブ [numerical-analysis]
For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.
python - Python の数値精度
階乗の計算など、Python を使用した簡単な数値実験を行ってきました。たとえば、32 の階乗を計算します。
私のルーティーン:
scipy.misc から:
私のルーチンは階乗の対数を計算し、1 から 32 (この場合) までの対数の合計を計算し、次に exp 関数を使用するだけであることを指摘したいと思います ( Fortran 90)。
正解は驚きです
pari/gpによると。
誰かがPythonで正しい数値の答えを探すことができる参考文献を教えてくれたらとてもうれしいです. ドキュメントは問題ありませんが、「短い」番号が必要な場合のみです。
c++ - If ステートメントが TRUE を評価すべきときに TRUE を評価しない
「x」値の配列 (PDE ソルバーのグリッド) があり、これらの x 値に基づいて別の配列を満たす関数に渡すと、1 つの x 値を含む特定の式が正しく評価されません。x 値の範囲は -1:1 で 0.0125 ずつ増加し、x = -0.5 と x = 0.5 では、これらのケースを他の値とは異なる方法で処理する必要があります。ただし、以下のブロックは点 x = 0.5 で TRUE と評価されません (x = -0.5 では問題ありません)。以下は、問題のあるブロックのスニペットを削除したもので、詳細は次のとおりです。
x = 0.5 の出力は次のとおりです。
したがって、if ステートメントの式は、x = 0.5 の場合に TRUE を返すように見えますが、範囲の「リム」内にあるため、正確には 0.5 ではありません。何かご意見は??
c - C の newton-raph func ポインター
単一の関数ポインターを使用して Newton-Raphson メソッドを実装しようとしています。関数には、方程式とその導関数の両方が含まれている必要があります。テスト関数内でこれら 2 つの別個の関数を通過するのが困難です。
ありがとうございました
javascript - Javascript の bigInt の
数値解析を始めたばかりです。
整数を扱う必要があるJavascriptでいくつかの幾何学的操作を実行しているとしましょう。
扱うことができる最大の整数サイズは?
JS が 32 ビットに対応できることは知っていますが、次のようないくつかのライブラリを使用すると、http ://www-cs-students.stanford.edu/~tjw/jsbn/ 64 ビットに対応できます。
ここでのロジックは、64 ビットを 2 組の 32 ビットに分割することだと思います。
同じロジックに従って、256 ビット数で計算を実行できるようにすることはできますか? ライブラリは存在しますか?
c++ - 記号演算から C まで、カスタムの非線形最小化を実装
コンピューター ビジョンの問題を解決しようとすると、非線形エネルギー関数を最小化し、C++ で実装する必要があります。特定の機能に役立つライブラリは見つかりませんでしたが、そのための数学はあります。では、記号演算から C++ コードに移行する最善の方法は何でしょうか?
例: 関数 g(x):=x^2 および f(x):=x+2 が与えられた場合、f(g(x)) を C コードに変換することに関心があるとします。明らかな C コードは y=x^2+2; です。ただし、ヤコビアンなどを含む複雑な数学の場合、操作のページとページに変換するのはそれほど簡単ではありません。
私はすでにMatlabとCコードへの変換モジュールを試しましたが、コードは最適化されていません(例:結果を再利用する代わりに同じ操作が何度も繰り返されます)。
matlab - x^b+a = a となる最大の x を求めます
この質問で見た答えを適用しようとすると、 でa+x=a
うまくいきましたa+eps(a)/2
。があるとしますx^b+a=a
。ここb
で、 は小さな整数で、たとえば 3 とa=2000
です。その後a+(eps(a))^3
、またはa+(eps(a)/2)^3
常に number を返しますa
。誰かがの測定を手伝ってくれますかx
? いずれにせよ、eps とは異なっていても問題ありません。
ps1938+(eps(1938)/0.00000000469)^3
は、 を返す最後の数値ですans = 1.9380e+003
。
1938+(eps(1938)/0.0000000047)^3
戻りますa=1938
。それは何か関係がありますか?
scipy - 発散の補間
均一に配置された 2D グリッド上に所定の速度フィールドがあります。グリッド ポイントまたは中間点の位置と一致しないいくつかの任意のポイントで速度場の発散を計算する必要があります。
グリッド点の発散を計算し、未知の点でバイキュービック補間を使用して補間することが可能です。チェッカーボード効果が発生する可能性があるため、これが良いアイデアかどうかはわかりません。
それを行うためのより良い方法があるかどうか疑問に思います。
私はpython/scipyを使用しています。
matlab - MATLAB: 最小二乗法を使用して 2 次、3 次、および 4 次多項式で明日の気温を概算する
以下は、数値解析タスクの演習 3 で、このテーマに関する大学のコースの一環として行う必要があります。
あなたの都市の過去 23 時間の気温の値に基づいて、明日の気温の概算を見つけます (テッサロニキ {私の大学の都市} の気象履歴は、ここで見つけることができます: http://freemeteo.com )
最小二乗法を使用して、温度関数を 2 次、3 次、および 4 次の多項式で近似します。それに続いて、興味のあるポイントで関数の値を見つけます。概算を定性的に比較し、概算を行っている日時をメモします。
最初の 2 つの作業を休むことなくこなしたことによる疲労のためか、数値解析の経験がないためか、完全に困惑しています。どこから始めればよいかさえわかりません。
努力の兆しさえ見せずに解決策を求めるのはうんざりすることだと私は知っていますが、何でも感謝します. リード、チュートリアル、取り組まなければならないことの概要、次々と。
大変お世話になりました。
注: MATLAB の組み込みの近似関数を使用できません。