sympy

2熱

2答えて

私は import sympy as sp x = sp.Symbol('x') sp.simplify(sp.log(sp.exp(x))) を入力すると、私は log(e^x) 代わりのxを得ます。私はこの機能に「保証はありません」と知っています。質問。指数の対数をアイデンティティ関数に変換するいくつかの特定の簡素化（級数展開などによる）がありますか？ import sympy

1熱

1答えて

Sympyでx/2を1/2 xとして印刷

たとえば、x/2を1/2 xとして印刷するようにSympyに指示するにはどうすればよいですか？私。分子を分母に1を乗じて分数を出力しますか？

1熱

1答えて

SymPyの再帰的区分関数と係数のリスト/配列を合計する

ラテックスエディタがないため、ここではSympyを使用してプロットしたい区分的な関数の画像を示します。私は係数の2つの配列とxの値を渡し、それを評価して関数をプロットしたいと思います。（編集：アルファがあるよりも、正確に一つ以上のpがある、画像更新）これは、これまで（アルファおよびpはリスト/配列で、tは数である）私の試みです： def getf(alpha,p,t): #Create t

3熱

1答えて

Sympyが複素双曲線割符（Sech）関数のフーリエ変換を見つけられない

複雑な引数を持つ双曲線割線関数（ "Sech"）のフーリエ変換を見つけるためにSympyを使用しようとしています。 import sympy as sy C, t, T0, f, w = sy.symbols('C, T, T_0, f, omega', real=True) Ut = sy.sech(t/T0) * sy.exp(-sy.sqrt(-1) * C/2 * t * t/(T0

1熱

1答えて

同じ入力に対して異なる式を返すSympify

私はSymPy式に解析する必要がある文字列のリストを持っています。問題は、そのうちの1人がSympifyError例外を発生させた場合、次の文字列は例外の前に同じ文字列が返された同じSymPy式を返しません。例えば、以下のコード： import traceback from sympy.core import sympify from sympy.core.evaluate import

1熱

1答えて

sympyを使ってITEをピース単位に変換する

sympyを使ってITEをpiecewiseに変換するにはどうすればよいですか？ ITEクラスは_eval_rewrite_as_Piecewise(self, *args)ですが、どのように使用すればよいかわかりません。は、私が試した： ITE(x,y,y)._eval_rewrite_as_Piecewise(x,y) が、それは実際にエラーが発生しました。だから私はそれをどのように使用

0熱

1答えて

PythonのLHS値の方程式を解くには

私はsympyライブラリを使って次の方程式を解くことを試みています。 I = -0.53V + 65.88、私は次のコード行使用されている： from sympy import * V = Symbol('vi') I = -0.5*vi + 65.88 solve(V - 5, I) をしかし、それは私に空のリストを返します。私は、このコードがRHS上で未知のものを見つけるために設計

1熱

1答えて

SymPyで特定のシンボルを取得する

ギリシャ文字nuの後にプライマリの上付き文字が続くシンボルを達成しようとしています。これはLaTeXのを使用して容易に実現することができますいくつかを言及する nuprime = symbols('{\nu}\'') nuprime = symbols('{nu}{\'}') nuprime = symbols('nu\'') nuprime = symbols('$nu\'$') ：

2熱

1答えて

大規模な行列を使ったsympyコンパイル関数

私はsympyを使って微分方程式のシステムを扱っています。私は式を象徴的に書いて、autowrapを使ってcythonでそれらをコンパイルし、その結果の関数をscipy ODEソルバに渡します。これを行うことの大きな利点の1つは、sympy jacobian関数を使用して象徴的にジャコビアを解くことができ、それをコンパイルしてODEソルバーにも渡せることです。これは約30変数のシステムでうまく機

2熱

1答えて

sympyでの条件付き期待

sympyでランダム変数の条件付き期待値を計算するにはどうすればよいですか？私はthisを読んでみました： from sympy.stats import * v = Uniform("v",0,1) E(v) これは、正しく1/2を返しますが、： E(v, v>1/2) をNaNを返し。私も試してみました： E(v, where(v > 1/2)) これは正しくありません（3/