问题标签 [scipy]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
5779 浏览

python - 如何使用 SciPy 运行统计累积分布函数和概率密度函数?

我是 Python 新手,也是 SciPy 库的新手。在深入 SciPy 世界之前,我想从列表中的专家那里得到一些问题。

我想知道是否有人可以提供有关如何运行两个统计函数的粗略指南:累积分布函数 (CDF) 和概率分布函数 (PDF)。

我的用例如下:我有一个有 1000 个浮点值的 sampleSpaceList []。当在我的程序中生成一个新的浮点值时,我想在 sampleList 上为其运行 CDF 和 PDF,并获得 CDF 的值小于或等于的概率和 PDF 的概率分布。

更多信息

基本上,在我的程序中有可能成功或失败的事件。如果他们成功了,那么我计算那个事件的事件比率并添加到我的 sampleSpaceList 直到它达到 1000 的阈值。一旦达到阈值,那么对于任何下一个事件比率;我想知道该事件比率在我的系统中是否成功的概率。

我基本上想得到的是特定事件比率的成功概率。

我不太确定 CDF 或 PDF 是否与我的问题相关,这就是为什么我想学习如何同时使用这两种方法,但在任何给定时刻,我只会使用 CDF 或 PDF 来获得事件的概率-成功率。

0 投票
3 回答
3226 浏览

math - 使用 octave/matlab ODE 求解器时拟合 ODE 的参数

我在 Octave 中使用 OdePkg 来解决僵硬的 ODE 系统,例如通过 ode5r:

求解器将因变量 Y 存储在关于时间 t(向量)的矩阵中:

我想在 param 中拟合参数,以便生成的变量 Y 最适合我的参考值,例如:

哪个 Octave/Matlab(欢迎使用其他语言)例程可以执行多参数(最小二乘/样条)拟合?如何在拟合中组合不同初始值Y0 的参数集?如果你能给我一些提示和可能性,我会很高兴。

最好的问候,西蒙

0 投票
2 回答
594 浏览

python - 使用 Enthought

我需要计算一个问题的互补误差函数 (erfc^(1)) 的倒数。

我正在研究它的 Python 工具,许多线程说 Enthought 拥有大部分所需的数学工具,所以我下载并安装在我的本地用户帐户中。但是我不太清楚如何使用它?

有任何想法吗?

0 投票
1 回答
1277 浏览

python - 多维列表(数组)重新分配问题

美好的一天编码员和编码员,

我正在编写一段代码,它通过一堆统计数据并返回我从中询问的内容。为了完成其任务,该方法从一个多维数组中读取并写入另一个数组。给我带来问题的代码是:

两个切片的大小都是,27:1:6:1但会引发异常:

我大吃一惊。

谢谢你。

0 投票
2 回答
692 浏览

python - 严格限制参数值在范围内

我正在尝试使用 scipy 中的 l_bfgs 约束优化例程来优化函数。但是优化例程将值传递给函数,这些值不在 Bounds 范围内。

我的完整代码看起来像,

基于初始值(初始)和边界(间隔) opt = optimize.fmin_l_bfgs() 将值传递给我的软件进行模拟,但传递的值应该在“边界”内。情况并非如此..请参阅下面在各种迭代中传递的值

在本次迭代中,-4.4139714 被传递给我的第二个参数,但它应该与 (10.0, 50000.0) 不同,我不知道 -4.4139714 从哪里来?

我应该在哪里更改代码?以便它传递应该在范围内的值

0 投票
2 回答
7167 浏览

python - 如何使用 csv dictreader、matplotlib 和 numpy 创建概率密度函数图?

我正在尝试使用 csv dictreader、matplotlib 和 numpy 使用 csv 文件的一列中的数据创建一个简单的概率密度函数(pdf)图...

有没有一种简单的方法可以将 CSV DictReader 与 numpy 数组结合使用?下面是不起作用的代码。错误消息是 TypeError: len() of unsized object,我猜这与我的数据不是 numpy 数组格式有关吗?我的数据也有负数和正数。提前致谢!

0 投票
1 回答
639 浏览

scipy - 在数据集的端点处插值数据点

从 scipy 我使用 interpolate.splrep 和 interpolate.splev 函数来插入我的数据集。不出所料,如果我尝试在数据集的边缘附近获得一个插值,这不会很好。

我想出了一个解决方法(通过两个附加条目扩展数据集,这些条目与数据集的最后一个“常规”条目具有相同的值;这似乎可行),但我想知道是否应该使用另一个插值函数,或者如果有一个通用的方法来处理这个。

提前致谢!

0 投票
2 回答
3179 浏览

python - 在 scipy 中为 fmin_cobyla 指定约束

我使用 Python 2.5。

我将界限传递给 cobyla 优化:

基于初始值Initial并根据/在范围内b1b10值传递给opt(). 但是这些值是有偏差的,尤其是b9. 这是我的问题的一个非常重要的边界条件!

“在每次迭代中x[2]传递给我的函数的值opt()必须始终大于x[3]”——如何实现这一点?

我的 bounds ( b1to b9) 定义有什么问题吗?

还是有更好的方法来定义我的界限?

请帮我。

0 投票
7 回答
11913 浏览

python - 分配问题,一个 NumPy 函数?

由于分配问题可以以单个矩阵的形式提出,我想知道 NumPy 是否具有解决这样一个矩阵的功能。到目前为止,我还没有找到。也许你们中的一个人知道 NumPy/SciPy 是否具有分配问题解决功能?

编辑:与此同时,我在http://software.clapper.org/munkres/找到了一个 Python(不是 NumPy/SciPy)实现。我仍然认为 NumPy/SciPy 实现可能会快得多,对吧?

0 投票
3 回答
5744 浏览

python - scipy '最小化一组方程的平方和'

我在 scipy 'leastsq' 优化例程中遇到问题,如果我执行以下程序,它会说

有时index out of range for an array...

我不知道我的 func optimize.leastsq() 调用有什么问题,请帮助我