lyapnov 指數計算

2021-09-26 12:04:12 字數 2964 閱讀 4599

from sympy import

*import numpy as np

import copy

n =2000

#控制迭代次數

defhenon

(x,y,n)

:for i in

range

(n):

x1 =1-

1.4* x **

2+ y

y1 =

0.3* x

x = x1

y = y1

return x,y

defle_calculate()

: sum_lambda1 =

0 sum_lambda2 =

0 a =

0.123456789

b =0.123456789

# 使用符號方式求解

x, y = symbols(

"x,y"

) f_mat = matrix([1

-1.4

* x **

2+ y,

0.3* x]

)# 求解雅各比矩陣

jacobi_mat = f_mat.jacobian(

[x, y]

)# 帶變數的雅各比矩陣形式是固定的

a, b = henon(a, b,

1001

)# 先迭代5000次,消除初始影響.以第5001次的值作為初始值

u1 = matrix([1

,0])

# 初始列向量

u2 = matrix([0

,1])

for i in

range

(n):

j = jacobi_mat.subs(

)# 將變數替換為當前迭代值,得到當前的雅各比矩陣(數字)

column_vector1 = u1#初始列向量為上一次的u1和u2

column_vector2 = u2

v1 = j * column_vector1 # 初始列向量乘上雅各比矩陣之後得到的向量

v2 = j * column_vector2

u1 = v1 /

(v1.norm(2)

) v2 = v2 -

(v2.dot(u1)

)* u1

u2 = v2 /

(v2.norm(2)

) lambda1 = ln(v1.norm(2)

) lambda2 = ln(v2.norm(2)

) sum_lambda1 = sum_lambda1 + lambda1

sum_lambda2 = sum_lambda2 + lambda2

a, b = henon(a,b,1)

#進行下一次迭代

le1=sum_lambda1/n

le2=sum_lambda2/n

print

(le1)

print

(le2)

if __name__ ==

'__main__'

: le_calculate(

)

0.414873360027928

-1.61884616435387

#-*-coding:utf-8-*-

from sympy import

*import numpy as np

# np.set_printoptions(suppress=true)

import copy

n =3000

#控制迭代次數

deflogistic

(x,n)

:for i in

range

(n):

x =4* x *(1

- x)

return x

defle_calculate()

: sum_lambda1 =

0 a =

0.123456789

# 使用符號方式求解

x = symbols(

"x")

f_mat = matrix([4

* x *(1

- x)])

# 求解雅各比矩陣

jacobi_mat = f_mat.jacobian(

[x])

# 帶變數的雅各比矩陣形式是固定的

a = logistic(a,

5001

)# 先迭代5000次,消除初始影響.以第5001次的值作為初始值

u1 = matrix([1

])# 初始列向量

for i in

range

(n):

j = jacobi_mat.subs(x,a)

# 將變數替換為當前迭代值,得到當前的雅各比矩陣(數字)

column_vector1 = u1#初始列向量為上一次的u1

vector1 = j * column_vector1 # 初始列向量乘上雅各比矩陣之後得到的向量

v1 = vector1 # 將vector1複製給v1

u1 = v1 /

(v1.norm(2)

)# 向量u1等於向量v1除以向量v1的模(2範數)

lambda1 = ln(v1.norm(2)

) sum_lambda1 = sum_lambda1 + lambda1

a= logistic(a,1)

#進行下一次迭代

le1=sum_lambda1/n

print

(le1)

if __name__ ==

'__main__'

: le_calculate(

)

0.692823875867068

計算BMI指數

coding utf 8 height float input 請輸入您的身高 m weight float input 請輸入您的體重 kg bmi weight height height if bmi 29.9 print 您的bmi指數為 str bmi print 肥胖 else if b...

計算BMI指數(身體質量指數)

題目描述 問題 計算bmi指數 身體質量指數 bmi指數 即身體質量指數,簡稱體質指數又稱體重,英文為body mass index,簡稱bmi 是用體重公斤數除以身高公尺數平方得出的數字,是目前國際上常用的衡量人體胖瘦程度以及是否健康的乙個標準。主要用於統計用途,當我們需要比較及分析乙個人的體重對...

PAT L1 012 計算指數

題目描述 真的沒騙你,這道才是簡單題 對任意給定的不超過10的正整數n,要求你輸出2n。不難吧?輸入格式 輸入在一行中給出乙個不超過10的正整數n。輸出格式 在一行中按照格式 2 n 計算結果 輸出2n的值。輸入樣例 5輸出樣例 2 5 32哈 pow 底數,次數 如下 include includ...