乗算

2017-01-01 20 views
2

機能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) 
    dotprod2 = np.dot(vec1,vec2) 
print(dotprod1) 
print(dotprod2) 

答えて

7

vec1vec2はnumpyのアレイの代わりにPythonのなければなりませんリスト:

import numpy as np 
from numpy import dot, exp, log, sqrt 
a, b = 1.2, 0.4 
x, y = 0., 0. 

vec1 = np.array([[1], [0]]) 
vec2 = np.array([[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 = dot(vec1,vec1) 
    dotprod2 = dot(vec1,vec2) 
print(dotprod1) 
print(dotprod2) 

出力:

[[ 2.32511745e+112 1.52483358e+056] 
[ 0.00000000e+000 0.00000000e+000]] 
[[ 1.07150861e-199 0.00000000e+000] 
[ 0.00000000e+000 0.00000000e+000]]