3

私はいくつかの異なる関数をグラフ化するコードを書いていますが、その多くは指数関数です。これは私が持っているものです:

>>> from fractions import Fraction
>>> f13=Fraction('1/3')
>>> f23=Fraction('2/3')
>>> f43=Fraction('4/3')
>>> f53=Fraction('5/3')
>>> f12=Fraction('1/2')
>>> f32=Fraction('3/2')
>>> f56=Fraction('5/6')
>>> fm1=Fraction('-1')
>>> fm23=Fraction('-2/3')
>>> fm32=Fraction('-3/2')
>>> from matplotlib import pyplot as plt
>>> import numpy as np
>>> x=np.arange(0.01,2,0.01)
>>> xa=x/(1+0.1071*x)
>>> rate=(4.817e+6)*(x**fm23)*np.exp(-14.964/(x**f13))*(1+0.0325*(x**f13)-(1.04e-3)*   (x**f23)-(2.37e-4)*x-(8.11e-5)*(x**f43)-(4.69e-5)*(x**f53))+(5.938e+6)*(xa**f56)*(x**fm32)*np.exp(-12.859/(xa**f13))

そして、それを入力すると、次のエラーが表示されます。

Traceback (most recent call last):
 File "<pyshell#22>", line 1, in <module>
  rate=rate=(4.817e+6)*(x**fm23)*np.exp(-14.964/(x**f13))*(1+0.0325*(x**f13)-(1.04e-3)*   (x**f23)-(2.37e-4)*x-(8.11e-5)*(x**f43)-(4.69e-5)*(x**f53))+(5.938e+6)*(xa**f56)*(x**fm32)*np.exp(-12.859/(xa**f13))
AttributeError: exp

ただし、この関数に対して同じコードを作成しました。

>>> rate=(7.29e+2)+2.40*((10**3)*(x**fm32)*np.exp(-0.223/x))

そして、それはうまく機能します。誰が何が間違っているのか知っていますか?

4

1 に答える 1

3

NumPy 型は とうまく混ざりませんfractions.Fraction。式x**f13の ,xは dtype の通常の NumPy 配列ですがfloat64、オブジェクトです。その結果、式は dtype を持つ配列になり、そのすべての要素は NumPy float ではなく (通常の) Python float になります。これは基本的に、NumPy が結果の型を特定できなかったことを意味するため、パックされた効率的な同種型の配列ではなく、Python オブジェクトのコレクションを格納することになります。f13Fractionx**f13object

あなたが与える2番目の式で-0.223/xは、 dtype のNumPy配列であるfloat64ため、それに適用しても問題はありませんnp.exp。

np.expdtype の配列を適用すると、各要素objectのメソッドが個別に検索され、表示されているものが得られます。次に例を示します。expAttributeError

Python 2.6.9 (unknown, Nov 18 2013, 14:53:18) 
[GCC 4.2.1 Compatible Apple LLVM 5.0 (clang-500.2.79)] on darwin
Type "help", "copyright", "credits" or "license" for more information.
>>> import numpy as np
>>> x = np.array([0.1, 0.2], dtype=object)
>>> np.exp(x)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
AttributeError: exp

Python 2.7.8 では、少なくともexp属性を見つけることができなかったオブジェクトのタイプについて言及している、わずかに優れたエラー メッセージが表示されます。

AttributeError: 'float' object has no attribute 'exp'

このエラー メッセージは、Python 2.7.5 と Python 2.7.8 の間のある時点で改善されたようですが、いつになるかはわかりません。

推奨される解決策: NumPy 配列と結合する前に、Fractionインスタンスを次のように変換します。float

>>> import numpy as np
>>> from fractions import Fraction
>>> np.exp(np.arange(10) ** Fraction(1, 3))  # Fails as above
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
AttributeError: exp
>>> np.exp(np.arange(10) ** float(Fraction(1, 3)))  # Works.
array([ 1.        ,  2.71828183,  3.52514317,  4.23020126,  4.89102089,
        5.52882849,  6.15411272,  6.77291238,  7.3890561 ,  8.0051399 ])

Fractionまたは、タイプを完全に避けます。ここでは、実際には何も購入しません。いずれにしても、フロートと組み合わせるとすぐに正確さが失われます。

>>> f13 = 1.0 / 3.0  # Look Ma, no fractions!
>>> np.exp(np.arange(10) ** f13)
array([ 1.        ,  2.71828183,  3.52514317,  4.23020126,  4.89102089,
        5.52882849,  6.15411272,  6.77291238,  7.3890561 ,  8.0051399 ])
于 2014-07-24T18:27:44.470 に答える