小编aja*_*eep的帖子

在 Julia 中访问数组中的一对行元素

假设我有以下数组

A = [1 2; 1 3; 1 5; 2 1; 2 3; 2 4; 3 1; 3 3; 3 4; 3 1; 4 2; 4 3; 4 5; 5 2; 5 3; 5 5]
Run Code Online (Sandbox Code Playgroud)

在 Julia 中查找给定数组 A 中是否存在一对行元素的一个 liner 是什么。例如,一个 liner 必须为 [1 4] 返回 false,因为 [1 4] 行在 A 中不存在,并为 [ 返回 true 3 4]。我尝试了以下方式

[3 4] .== A
Run Code Online (Sandbox Code Playgroud)

返回

0  1

1  0

1  0

1  1

1  0
Run Code Online (Sandbox Code Playgroud)

而我希望我的输出仅为“真实”

1 1
Run Code Online (Sandbox Code Playgroud)

julia

3
推荐指数
1
解决办法
153
查看次数

python中矩阵的乘法导致类型错误

在尝试从以下函数的jacobian矩阵计算点积时x(n+1) = 1 -a*x(n)**2 + y(n); y(n+1) = b*x(n),我得到以下错误:vec1 = J*vec1 TypeError: can't multiply sequence by non-int of type 'list'

import numpy as np
a, b = 1.2, 0.4
x, y = 0., 0.

vec1 = [[1], [0]]
vec2 = [[0], [1]]

for i in range(500):
    x1 = 1. - a*x**2 + y
    y1 = b*x
    x = x1
    y = y1
    J = [[-2.*a*x, 1], [b, 0]]
    vec1 = J*vec1
    vec2 = J*vec2
    dotprod1 = np.dot(vec1,vec1) …
Run Code Online (Sandbox Code Playgroud)

python numpy python-2.7

2
推荐指数
1
解决办法
74
查看次数

标签 统计

julia ×1

numpy ×1

python ×1

python-2.7 ×1