問題タブ [equation-solving]

n*log2(n) = A、A既知の値。nMatlabでどのように解決しますか? n必ずしも整数ではないことに注意してください。

C++のLapackパッケージを使用して線形方程式系を解きたい。ここからのルーチン、つまりdgesvを使用して、このように実装する予定です。

これは私のコードです：

私は次の問題を抱えています：

dgesvが要素{3、3}を使用してベクトルxを計算するのはどうしてですか？解は{0.6、0}である必要があります（matlabで確認）。

ご挨拶

編集：dgesvは正方行列に対して機能する可能性があります。以下の私の解決策は、dgelsを使用して過剰決定系を解決する方法を示しています。

Cで次の方程式を解く必要があります。

x の値はわかっているので、この方程式を満たす 'a' の値を見つける必要があります。この方程式を解くための効率的な方法を提案できる人はいますか? ニュートン・ラフソンか、それとももっと効率的なものでしょうか?

プロジェクトの一環として、C++ で閉じた形式で 4 次多項式を解く必要があります。

A*x ⁴ + B*x ³ + C*x ² + D*x + E = 0

この目的へのリンクがいくつか見つかりました。それらの 1 つがここにあります。しかし、すべての根を計算しますが、実際の根だけが必要です。アルゴリズムは、主にフェラーリの方法を使用して次数を減らします。

これは、実数と虚数の両方の解を与える solveCubic() を使用します。

コードを説明するリンクをいくつか示します。solveCubicとsolveQuartic

実根の 4 次多項式を解くためにコードを変更できる人はいますか?

できるだけ効率よく導入したい。ところで、誰かがこの目的のために LAPACK のような便利なライブラリを導入してくれたら幸いです (4 次多項式の根を直接計算することはできないようです)。

各変数が配置されている列ごとに、次の行から変数を取得して繰り返したい一連の方程式があります。私は初心者なので、コーディングはおそらくスクラッチではありません。

これは、各変数値に相対的な範囲として設定したいくつかの変数に対して行いました。そして、次の行では、方程式全体をやり直して、次の行の変数を再度設定する必要がありました。次の計算のために次の行のセルにスキップする変数を（おそらく相対セル参照を使用して）設定する方法はありますか？計算用に36行あり、約9つの変数があることも考慮してください。

私は、A、B、Cの3つの値を与える物理的な測定器（ロードセルを備えた力プラットフォーム）を持っています。ただし、これらの値（直交している必要があります）は、力とトルクの適用値と戻り値の間にクロストークを引き起こす測定装置。

次に、次のように、キャリブレーションマトリックスを使用して、測定値を実際の値のより適切な推定値に変換することをお勧めします。

問題は、一連の測定を実行する必要があることです。これにより、システム全体に最適なC行列を取得するために、異なる値measured(Fz, Mx, My)と最小二乗法が得られます。actual(Fz, Mx, My)

Ax = B1つの測定で問題を解決することscipy.linalg.lststqも、scipy.linalg.solve（正確な解決策を提供することで）解決することもできますが、それぞれが潜在的に異なる3x3行列を与える独自の方程式を持つ、一連の異なる測定を検討するにはどうすればよいですか？

読んでくれてありがとう、どんな助けでも大歓迎です。

重複の可能性:
A と B の両方が大きな行列の場合、MATLAB で AX=B の X を解く効率的な方法

X=A\b を実行して、matlab で方程式 Ax=b を解こうとしています。ただし、問題は、A が巨大な行列 (2048000x2048000) であり、実行中に matlab がクラッシュすることです。matlab がクラッシュしないように、この方程式を解く別の効率的な方法はありますか?

私は非常に下手なプログラマーなので、非常にばかげた間違いを犯している可能性が高いことに注意してください。私は C を書いていて、 gsl_multiroot_fsolver_hybridを間違って使用しています。

私のコードには、次の 2 つの関数が含まれています。

関数int COMPUTE_QUASILOGNORMAL_PARAMETERSは、 int quasilognormal_fで定義されている 3 つの連立方程式を解く必要があります。2 番目の (iter ==1) 反復が開始されるとすぐに、パラメーター N、A、および S がint quasilognormal_f内でNANになるため、これはうまくいきません。これは、プログラムの出力から次のように結論付けました。

ソルバーは私にとってブラックボックスであり、GSL マニュアルに従ってみました。

固定次元 (N=9) (行列は対称、半正定値) の密な線形システムの高速なソリューションに推奨するアルゴリズムはどれですか?

• ガウス消去法
• LU分解
• コレスキー分解
• 等？

型は 32 ビットと 64 ビットの浮動小数点です。

このようなシステムは何百万回も解決されるため、アルゴリズムは次元 (n=9) に関してかなり高速である必要があります。

提案されたアルゴリズムの堅牢なC++ 実装の PS の例は高く評価されます。

1) 「何百万回も解決した」とはどういう意味ですか? 100 万の異なる右辺項を持つ同じ係数行列、または 100 万の異なる行列ですか?

何百万もの異なるマトリックス。

2) 正の _semi_defined は、行列が (マシンの精度で) 特異になる可能性があることを意味します。この場合、どのように対処したいですか？エラーを発生させるだけですか、それとも賢明な回答を返そうとしますか?

上げ間違いOK。

実装固有の理由により、Java 1.2 を使用する必要があります。String数値 (そのステップを抽象化するために事前に変数を置き換えます) と演算子 (PDMAS) のみでオブジェクトを解析しようとしています。私はこれをうまく行う多くのライブラリを見つけましたが、残念ながら Java 1.2 と互換性のあるライブラリはありません (いじっていても、それらはすべてジェネリックなどに依存しています)。明らかに私はこれを自分で作ることができますが、ホイールを作り直したくないのは確かです. これを行うことができる、まだ見つけていないライブラリはありますか? ありがとう。

(要件: 二項演算子と括弧)

編集: 要求に応じて、入力と出力のいくつかの例:

• 「(10 / 5) + 4.5 - (8)」は -1.5 になります。
• "(1/3) * 4" は 1.3333333 になります...
• 「5^3 + 4 * 2」は 133 になります。
• 「-10 + 5」は-5になります

うまくいけば、それは理にかなっています。