- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
# necessary imports
import numpy as np
import matplotlib.pyplot as plt
def binet(n):
phi = (1 + 5 ** 0.5) / 2
return ((phi**n) - (-1/phi)**n) / (5**0.5)
对于
phi
值,我用了
this .
binet(n)
为
n=[0.1,0.2,0.3,0.4,0.5,...,4.9,5.0]
:
[binet(x/10) for x in range(1,51)]
让我们绘制它:
# our results
plt.plot([n.real for n in binetn],[n.imag for n in binetn])
# classic fibonacci numbers
plt.scatter([1,1,3,5],[0,0,0,0],c='r')
binetn=[binet(x) for x in np.arange(0.1,5.1,0.1)]
然而,事实并非如此。
binetn
变成:
[nan,nan,nan,nan,nan,nan,nan,nan,nan,1.0,...,nan,nan,5.000000000000001]
即它是
nan
除非
binet(n)
是真实的。
/usr/local/lib/python3.6/dist-packages/ipykernel_launcher.py:2:RuntimeWarning: invalid value encountered in double_scalars
range()
生成的数字列表& 得到复杂的结果,而我不能用 np.arange()
做同样的事情?
最佳答案
副本集中在 np.arange
.但这不是这里的问题。
比较将浮点数传递给 binet
与 np.float64
:
In [52]: binet(0.1)
Out[52]: (0.06391735396852471-0.13170388861716523j)
In [53]: binet(np.float64(0.1))
/usr/local/bin/ipython3:3: RuntimeWarning: invalid value encountered in double_scalars
# -*- coding: utf-8 -*-
Out[53]: nan
或这些迭代:
In [54]: [binet(float(x)) for x in np.arange(0.1, 5.1, 0.1)];
In [55]: [binet(x) for x in np.arange(0.1, 5.1, 0.1).tolist()];
In [56]: [binet(x) for x in np.arange(0.1, 5.1, 0.1)];
/usr/local/bin/ipython3:3: RuntimeWarning: invalid value encountered in double_scalars
# -*- coding: utf-8 -*-
进一步挖掘,它似乎与
complex
有关。 .如果
numpy
参数复杂,计算没问题:
In [58]: binet(np.complex(.1))
Out[58]: (0.06391735396852471-0.13170388861716523j)
In [59]: [binet(x) for x in np.arange(0.1, 5.1, 0.1).astype(complex)];
事实上,我们不需要迭代。一个复杂的数组工作得很好:
In [60]: x = np.arange(0.1, 5.1, 0.1).astype(complex)
In [61]: binet(x)
Out[61]:
array([0.06391735-1.31703889e-01j, 0.16378803-2.38746028e-01j,
0.28913834-3.13167258e-01j, 0.42813892-3.50853867e-01j,
0.56886448-3.51577584e-01j, 0.70044744-3.18660871e-01j,
0.81402504-2.58333837e-01j, 0.9034083 -1.78872498e-01j,
...
4.32034432-3.76903988e-02j, 4.54051382-2.60971457e-02j,
4.76690299-1.30754885e-02j, 5. +0.00000000e+00j])
关于python - 为什么 numpy.arange() 和 range() 生成的相同数字的(复杂)函数不同?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/63569585/
我在最近更新了运行 Ubuntu 的计算机并且 Python 的默认版本更改为 2.7 时注意到了这个问题。 import json import numpy as np json.dumps(lis
我需要将一些 Matlab 代码转换成 P。我被困在了 numpy.arange我使用给定角度(以弧度为单位)在圆弧上连续设置点。 我做到了这一点(例如 x 轴上的点): def sensor_dat
试试这个: import numpy as np np.arange(0,3*0.1,0.1) 输出将是: 数组([ 0., 0.1, 0.2, 0.3]) 这令人难以置信,因为对于 np.arang
我原以为 numpy 的 arange(start,end) 会生成 [start,end] 范围内的值。以下示例演示了情况并非总是如此(最终值大于 end): import numpy as n
这个问题已经有答案了: Is floating point math broken? (33 个回答) 已关闭 5 年前。 我正在尝试为 N 个离散分数的分布创建一个 N+1 个 bin 的数组。 我
我想创建每月间隔而不使用 np.arange 进行迭代。作为一个简单的示例,我想要一个包含 np.datetime64 对象的数组,表示 1990-2000 范围内的元旦。我想是这样的: np.ara
我正在教自己一些关于 numpy 的知识,我已经整理了一些旧的本科教材以用作示例。因此,我编写了一个 without numpy 函数来计算悬臂梁在任意点的单点载荷引起的挠度。非常简单,除了偏转方程会
有人可以向我解释这里发生了什么吗? 为什么 0.3 和 0.7 值有更多的小数点。我只想要 1 个小数点值。 threshold_range = np.arange(0.1,1,0.1) thresh
我有一个函数要应用于排列: import math from numpy import arange x = arange(7.0,39.0,0.0001) fx = math.exp(-2.0 /
np.arange 接受三个参数:开始、停止和步骤。 我想要的步长是-0.3048。我还有单独的数组用于启动和停止。 首先,我有一个完整的元素数组: array([5.000, 5.000, 5.00
我正在尝试用基于像素的图像的所有 x 坐标的列表填充一个 numpy 数组。所以基本上,n x m 图像将有 n 1,n 2的,依此类推,直到 n m。有没有一种简单的方法来填充一个 numpy 数组
以下是我使用的代码: import numpy as np import pandas as pd from pandas import DataFrame, Series animals = Dat
这个问题在这里已经有了答案: Concatenate range arrays given start, stop numbers in a vectorized way - NumPy (1 个回
有没有办法带... >>> x = np.array([0, 8, 10, 15, 50]).reshape((-1, 1)); ncols = 5 ……然后把它变成…… array([[ 0, 1
我已经使用 numpy 的 arange 函数来制作以下范围: a = n.arange(0,5,1/2) 这个变量本身可以正常工作,但是当我尝试将它放在我的脚本中的任何位置时,我会收到一条错误消息:
我正在使用 arange 函数来定义我的 for 循环迭代并得到意想不到的结果。 i = arange(7.8,8.4,0.05) print i 产生以下结果: [ 7.8 7.85 7.9
我在另一个使用 Numpy 函数 arange 的文件中有一个函数,即 import numpy as np def plot_2D_boundary(plot_range, points, deci
考虑以下 numpy 数组: import numpy as np arr = np.array([np.random.permutation(4) for _ in range(4)]) array
假设我想生成一个介于 0 和 1 之间、间距为 0.1 的数组。在 R 中,我们可以这样做 > seq(0, 1, 0.1) [1] 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
我想做我上传的 png 插图所具有的功能,主要有两列,其中有数字。带星号的代码为我提供了类和“SquareMiles”。百分比栏更难获得。我将系列号数组除以该列的总系列数,但没有成功。这个部门如何才能
我是一名优秀的程序员,十分优秀!