scipy

2熱

1答えて

私は、配列Aを持っている、と私はスライスインデックス（S、T）のリストを持っている、のは L.このリストを呼び出してみましょう私は[S1：T1]の85パーセンタイルを見つけたい、A [S2 ：t2] ... numpyでこれらの演算をベクトル化する方法はありますか？ ans = [] for (s,t) in L: ans.append(numpy.percentile(A[s:t],

6熱

1答えて

シフィー補間の勾配を直接求める方法は？

私はスカラー値の大規模な3D numpy配列を持っています（必要ならばOKと呼んでください）。私は、不均一な連続でこの上にスムーススカラー場を補間したいと思っています。すべてのの非整数xyz座標が知られているわけではありません。今、このためscipyのダウンロードのサポートは、単に優れている：私は（ filtered_volume = scipy.ndimage.interpolation.s

4熱

3答えて

冗長要素を持つPython配列を使用したPython配列のインデックス

配列のインデックス付けに問題があります。 b要素が指し示す位置にある値を代入するために、配列aとbのインデックスとして使う別の配列bがあるとします。私は [ 0. 2. 0. 0. 0.] 配列bには何の冗長性（すべての値を持っていない場合は何も問題はありませんが期待しながら、 a=numpy.zeros(5) print a [ 0. 0. 0. 0. 0.] 今、私は二回 b=n

13熱

3答えて

PythonでのMATLABプログラマー向けのサイエンスコンピューティング

MATLABプログラマー向けのPython/SciPyの科学計算の入門書や優れたチュートリアルを誰かが知っているのだろうかと思っていました。私は、David Hiebeler's textのRの行に沿って、両方の言語間の並列参照を提供するMATLABプログラマーのための何かを考えていました。私の現在の環境はPython(x,y)でSpyderですが、MATLABとSciPy +サードパーティエ

2熱

3答えて

Pythonの最速の階乗階乗（Pochhammer関数）

大きな数字のrising factorialを計算する必要があります。これまでsympyパッケージsympy packageからの階乗階乗関数が見つかりました。もっと早く。は私がまさに必要となることは、こののは本当に速いバージョンです： from itertools import combinations from numpy import prod def my_rising_facto

22熱

2答えて

scipy.spatial.KDTreeとscipy.spatial.cKDTreeの相違点

これら2つのアルゴリズムの違いは何ですか？

0熱

2答えて

numpyの配列を回転するには？

私はいくつかの気が/ scipyの問題があります。私は2進法[楕円体から0]で楕円体を表す3D配列を持っています。私はある程度の形を回転させたいと思っています。それが可能だと思いますか？回転を使って楕円体方程式を直接書き込む効率的な方法はありますか？

2熱

1答えて

Scipy Leastsqオプションの出力変数（メッシュ）

SciPyのLeastSqに関するドキュメントはhereです。 ierは整数フラグです。 1、2、3、または4に等しい場合、解が見つかりました。それ以外の場合は、解決策が見つかりませんでした。どちらの場合でも、オプションの出力変数 'mesg'はより多くの情報を提供します。しかし、オプション変数mesgを取得するにはどうすればよいですか？ x,mesg,ier=leastsq(residual

4熱

1答えて

SciPy LeastSq Dfun Usage

私はJacobianをSciPyのOptimizeライブラリのleastsq関数で動作させようとしています。私は、次のコードを持っている：私はこれを実行したとき、私は私のリターンとしてplsq=[10,10,10]を取得し、今 #!/usr/bin/python import scipy import numpy from scipy.optimize import leastsq #

1熱

1答えて

係数の大部分がゼロのときにIIRフィルタを行う方法

私はPythonで何らかのオーディオエフェクトを実行する必要があります。例えば、最も単純なエコー効果式：これはIIRフィルタであり、scipy.signalにlfilter（）により算出することができる y[n] = x[n] + k*y[n-1000] ： import numpy as np import time import scipy.signal as signal pul