numerical-methods

    1

    1答えて

    現在、C++の有限差分法を使って微分方程式を解く小さなプログラムを書いています。問題は私の三重対角行列を解くためにlapackを使用していることですが、ルーチンが呼び出された後にSegmentation Faultが発生しています。ルーチンはinfo = 0で終了しますので、問題の内容がわかりません。ここにプログラムがあります。 #include <iostream> #include <cma

    0

    1答えて

    ラプラス変換をMathematica、Matlab、Pythonなどの言語で1次元(例:TalbotまたはEulerアルゴリズム)に数値的に逆転させるアルゴリズムには、多くの実装があります。 私は多次元の場合や、自分自身をコード化することができるものを(任意の言語で)実装したいと考えています。もしあなたがこれを知っているなら、私に知らせてください。

    0

    2答えて

    最近C++を使い始めたばかりで、ユーザー定義のodeのシステムの統合を可能にするクラスを作成しました。その性能を比較するために2つの異なるインテグレータを使用します。ここで、コードの一般的なレイアウトは、次のとおり class integrators { private: double ti; // initial time double *xi; // ini

    -1

    2答えて

    数値のxとyデータをnumpyの配列として取り、それをxの領域に(数値的に)積分するプログラムを書いています。私はこれをして、今は精度について疑問に思っています。 私は、math.fsum()がsum()でない特定のエラーを修正していることに気がつきました。numpy.sumもこれらの修正を行っていますか? sum = numpy.sum(array)などと書いても、結果(合計)はできるだけ正確で

    1

    1答えて

    私はプロット可能なベクトルを持っており、その積分を計算したいと思います。私は下の総面積を意味するのではなく、統合の領域をどのように進化させているのでしょうか。基本的に、その「不確定」積分。これは可能ですか?多分補間によって? Chebyshevの微分行列を扱っているので、積分に相当するものがあるかどうか知っていますか? はおそらくcumtrapzたいあなたに

    1

    1答えて

    MATLABを使用して次の式の誤差の確率をプロットしようとしています。数値積分にtrapzというコマンドを使用したいのですが、問題は、 y軸の値は間違っていますが、全体の曲線は0と1.2の間でなければなりませんが、0.492と0.5の間です!誰かが自分のコードで何が間違っているかを教えてくれますか、それとも私にヒントを与えてもらえますか?私は本当に助けが必要です。ここで私は(Maketexを使用し

    1

    2答えて

    これは、ANSI Cコードです:このコードの float x = 3.14159264351134890172; double y = 3.14159264351134890172; long double z = 3.14159264351134890172; printf("%f\n",x); printf("%f\n",y); prin

    0

    1答えて

    私はMATLABのFMINCONの "interior-point"アルゴリズムを実行しています。このアルゴリズムはBFGSを使用しており、数値ヘッセ行列の情報を考慮しています。 私は、一意の解法に関して非凸であることが知られているパラメータ空間に対して最適化を行っています。 arctan関数で一意でないパラメータの組み合わせにペナルティを課します。 私は最適化手順が高い(約10回)、いくつかの反

    1

    1答えて

    いくつかの積分方法を検討しています。私は長方形の統合方法(すなわち、オイラーの方法)に遭遇しました。私が読んでいる本によると、実際のアルゴリズムは コード自明です。この例では、dxdt=cos(t)の初期値はx(0)=0です。分析解はx=sin(t)です。私たちは、t=1で誤差を計算することができますし、実際、著者は次の表にそれを提供しています。 問題は、私の次のコードでは、誤差は9.1%であると