問題タブ [determinants]

このラプラス展開を Matlab で行列式に対して実行すると、「未定義の関数または変数 'A'」が表示されます。誰かがそれを修正して実行可能にするのを手伝ってくれれば幸いです。

ありがとう

行ベクトルがnあり、ベクトルを選択して行列kを形成できます。k*k次に、各行列の行列式の絶対値を取ります。したがって、(n choose k)絶対値決定要因の可能性があります。

最高の行列式を生成する行ベクトルを選択したいと考えています。これをすばやく行うために検索を最適化することは可能ですか?

私は現在 と で作業してn ~= 10,000いk = 31ますが、さまざまなセットアップで作業する予定です。

私はPythonとC ++でコードを書いています

次数 N < 10 の行列の Cramer の規則 (2 つの行列式の商) を使用して線形方程式系を解くと、LAPACK ソリューションと比較してかなり大きな残差エラーが発生することがわかりました。

次に例を示します。

これをlinsolvingすると、LAPACK (Intel MKL から) が生成されます。

x = [-0.000314947 -0.000589154 -0.00587876 0.0184799 0.01738 -0.0170484]

Cramer の規則 (独自の実装) により、次の結果が得られます。

x = [-0.000314933 -0.000798058 -0.00587888 0.0184808 0.017381 -0.0170508]

違いに注意x[1]。

私の行列式計算が正しいことを保証できます。誰かが同様の観察をしたか、これについて何か言うことができますか?

行列 nxn の行列式の計算の複雑さを実験的に決定する助けが必要です

私のコード:

しかし、すべての n に対して同じ時間を取得するため、複雑さを計算します: O(N) は、O(N^3) であることが意図されているため、正しくありません。どんな助けでも大歓迎です。

上三角行列の行列式を見つける関数を機能させるのに問題があります。私のコードは明らかに正しくない値を返すようです。通常はゼロです。これは、関数を何らかの方法で誤って定義したことが原因であると確信しています。それは私の側の基本的なエラーだと思いますが、しばらくそれを見つめた後、私はそれを理解することができませんでした. 関数と印刷コードは次のとおりです。

および印刷機能：

詳細を提供するために、これまでに作成した完全なコードを含めます。現時点での基本的なアプリケーションは、n x n 行列 (2