4

私はSymPyにかなり慣れていないので、基本的な質問があります。あるいは、SymPyがどのように使用されるのかを単に誤解しているかもしれません。

アトムではなく、他の式の組み合わせで表される式を作成する方法はありますか？

例：

>>> from sympy.physics.units import *
>>> expr1 = m/s
>>> expr2 = mile/hour
>>> expr1
m/s
>>> expr2
1397*m/(3125*s)
>>> expr1.in_terms_of([mile,hour]) #in my dreams?
3125*mile/(1397*hour)
>>>

補足：完全なSymPyドキュメントの「公式」PDF（または他の印刷可能な）バージョンを見つけることができますか？（私は職場でのインターネットの使用制限が厳しく、週末に自宅で仕事をするのにうんざりしています。）

アップデート：

これは私がプレリュードの提案に従うことになったものですが、それは汚れているように感じるので、このようにあまり使用されない可能性があります。コメントとWTFの歓迎。

def in_terms_of(self, terms):
    expr2 = eval(str(self), physics.units.__dict__)
    converter = {}
    for term in terms:
        term_expr = physics.units.__dict__[term]
        coeff = term_expr.as_coeff_terms()[0]
        unit = physics.units.Unit('my_'+term, term)
        converter[term_expr/coeff] = unit/coeff
    return str(expr2.subs(converter))

使用法：

>>> x = in_terms_of('J',['gram','mile','hour'])
>>> x
'9765625000*mile**2*gram/(1951609*hour**2)'
>>> in_terms_of(x,['J'])
'J'
4

3 に答える 3

2

sympy.physics.unitsのソースを確認すると、すべての単位がメートル、キログラム、秒、アンペア、ケルビン、モル、カンデラで定義されていることがわかります。これらは基本単位です。

次に、1マイルは5280フィートとして定義され、1フィートは0.3048メートルとして定義されます。

したがって、非基本単位を使用するすべての式では、非基本単位が基本単位に置き換えられます。

特定の単位を使用するために式が必要な場合に、式に置き換えることができる独自の単位を定義できます。

import sympy.physics.units as units
from sympy import Rational

my_mile = units.Unit('my_mile', 'mile')
my_hour = units.Unit('my_hour', 'hour')

次に、新しい単位の代わりに基本単位を使用する辞書を定義します。

converter = {units.m: my_mile/Rational('1609.344'),
             units.s: my_hour/Rational('3600')}

基本単位を使用してすべての計算を実行します。次に、マイルと時間を使用して値が必要な場合は、新しい単位を式に代入できます。

v = 10*units.miles/units.hour
print v # = 2794*m/(625*s)

print v.subs(converter) # = 10*mile/hour

ドキュメントを取得するには、arsanswerを使用してください。Sympyのリポジトリはここにあります：https ：//github.com/sympy/sympy

docsフォルダーには、htmlドキュメントの作成方法を説明するREADMEファイルがあります。

于 2010-01-11T11:16:48.493 に答える
2

数量もあります：

In [1]: from quantities import *
In [2]: v = 1397*m / (3125*s)
In [3]: v
Out[3]: array(0.44704) * m/s
In [8]: v.units = mile/hour
In [9]: v
Out[9]: array(1.0) * mi/h
于 2010-01-11T19:55:09.203 に答える
1

sympy単位では、次のように除算してexpr1からexpr2に変換します。

In [120]: import sympy.physics.units as units

In [121]: expr1 = units.m / units.s

In [122]: expr2 = units.miles / units.hour

In [123]: a = (1397/3125) * expr1 / expr2

In [124]: a
Out[124]: 1

唯一の問題は、無次元量になってしまうことです。このタイプのもののQuantitiesパッケージを検討することもできます。

In [125]: import quantities as pq

In [126]: a = pq.Quantity(1397/3125, 'm/s')

In [127]: a
Out[127]: array(0.44703999999999999) * m/s

In [128]: a.units = pq.mile / pq.hour

In [129]: a
Out[129]: array(1.0) * mi/h

Unumもあります。

sympyドキュメントのPDFバージョンはわかりませんが、リポジトリからディストリビューションをチェックアウトし、Sphinxを使用して自分用のPDFを生成できる可能性があります。

于 2010-01-11T09:43:51.097 に答える