問題タブ [sympy]

nsolve最初の推測を与えるいくつかの関数の解決策を見つけるのが難しいという問題がありました。次に、numpy/scipy ソルバーを試してみたいと思いました。

これはsympyを使用したプログラムで、このソリューションを提供すると非常にうまく機能します:[0.0, -9.05567e-72, 9.42477, 3.14159]

これは同じプログラムですが、numpy と互換性のある関数を使用しています。

しかし、私がプログラムを実行したとき：

x0 に値 0 を指定しようとしましたが、2 番目のプログラム (numpy を使用) は 0 に近い数値を指定して動作しましたが、pi/4 から開始するとエラー メッセージが表示されます。numpy で何かを見逃しましたか?

配列変数を記号的に扱う方法があるかどうか疑問に思っています（sympyのようなものですが、数値変数ではなく配列用です）。次のような配列式を持つことができるように

または、ネストされた配列関係など、さらに高度なもの

A、B を記号的に扱うことで、substitution を使用して A、B の内容を変更できます。たとえば、次のようになります。

次に、f1 と f2 などを追加できます。

上記の構造のようなものを使用して、引数のy切片を返す関数を作成するにはどうすればよいですか？

編集：方程式を取得した参照には、いくつかのエラーが含まれていました。ここで修正しました。解決策は、今では実際に理にかなっているかもしれません!

2 層の流体が地形上を流れる場合、流体内の流速と波の速度の相対的なサイズに応じて、さまざまな解が存在します。

これらは「超臨界」、「亜臨界」、「臨界」と呼ばれます (ここでは最初の 2 つを「超臨界」と呼びます)。

次の方程式は、(h, U0) パラメーター空間での臨界動作と超臨界動作の間の境界線を定義します。

を消去してd_1c(つまり、それが何であるかは気にしません)、これらの方程式の解を で見つけたいと思い(h, U_0)ます。

単純化要因:

• 与えられた答えだけが必要です d_0
• 正確な解は必要ありません。解曲線の概要だけが必要なので、これは解析的または数値的に解決できます。
• 領域 (h, U0) = (0,0) から (0.5, 1) だけをプロットしたいだけです。

Enthoughtディストリビューション（numpy、scipy、sympy）で利用可能なモジュールを使用してこれを解決したいのですが、どこから始めればよいのか本当にわかりません。私を本当に混乱させるのは、変数 d1c の削除です。

Python の方程式は次のとおりです。

私は、多くのソリューション ブランチ (必ずしも物理的ではありませんが、心配する必要はありません) を持ち、おおよそ次のようなソリューションを期待しています。

これを実装するにはどうすればよいですか？

sympy の解決機能は、解決策を思い付くまでに時間がかかります。必要な精度を下げることで高速化することは可能ですか (小数点以下 15 桁は必要ありません!)。

Boost-Pythonプログラムを作成して、ユーザーからシンボリックpython関数を取得し、その派生物をプログラムで評価したいと思います。

たとえば、ユーザーはF = sin（x）* cos（x）のような関数を定義するPythonファイル（Function.py）を提供します。

次に、Sympyのシンボリック微分能力を使用してF'（x）（F（x）の導関数）にアクセスしたいと思います。数値微分は使いたくない。

Boost-Pythonを使用してC++でそのような関数F'（x）にアクセスできるようにする方法はありますか？

重複の可能性:
SymPy - 任意の数のシンボル

反復によって形成される連立方程式を作成しようとしています。各方程式は次のように生成されます。

summation(binomial(n,i) * (1-t) ** (n-i) * t ** i * P[i]) - V[j]

すべてV[j]が入力として与えられます。アイデアは、作成された連立方程式を使用しP[i]てすべてを解決することP[i]ですが、そうするためには、sympy シンボル オブジェクトでなければならず、そのようなリストを生成する方法がわかりません。

名前のないリストを作成したり、シンボルのリストを生成したりすることは可能ですか?

質問が明確化を求めている場合はお知らせください。言葉にするのに苦労しました。

前もって感謝します ：）

今、私はエラーが発生しています：アプリケーションのサーバーエラー。入力された名前をインポートできません 説明: 現在の Web 要求の実行中に未処理の例外が発生しました。エラーの詳細とコード内のどこでエラーが発生したかについては、スタック トレースを確認してください。例外の詳細: IronPython.Runtime.Exceptions.ImportException: 入力された名前をインポートできません ソース エラー: 44 行目: expr)); 45行目: 46行目: script.Execute(scope); 48行目: return scope.GetVariable("result");

可変数の引数を持つpython関数があります：

各引数に関して F の導関数を表す N 個の関数を自動的に生成したいと考えています。

たとえば、F(x1) = sin(x1) と F(x1, x2) = sin(x1) * cos(x2) の両方を与えて、すべての導関数を自動的に取得できます。

Edit2：関数Fが2変数（引数の固定数）の場合、使用できます

次のコードの答えがなぜ 10 ではなく 10 なのか、私は困惑してlambdifyいます。

が より小さいのに、答えがk11 でないのはなぜですか?x = 0loc = 0.5

ただし、そうすると正しい答えが返されます

kfun後で積分の引数の一部として使用するため、x と y の関数として持つ必要があります。また、最終的には y に依存します。

Mac 10.6.xでpython 2.6.8を使用しています。