如何用python计算行列式_如何将Python数组转换为SymPy矩阵来计算行列式?

根据文档示例,只需执行以下操作:from sympy import *

import math

n=3

n_atoms=8

MainMatrix=Matrix([[0 for x in range(n)] for y in range(n)])

KappaMatrix=Matrix([0 for x in range(n-1)])

MassMatrix=Matrix([0 for x in range(n)])

Kappa=3

ka_by_pi=4

for i in range(1,n):

KappaMatrix[i-1]=3

for i in range(1,n+1):

MassMatrix[i-1]=3

for i in range(2,n):

MainMatrix[i-1]= -KappaMatrix[i-2]/MassMatrix[i-1]

MainMatrix[i-1]= (KappaMatrix[i-2] + KappaMatrix[i-1])/MassMatrix[i-1]

MainMatrix[i-1]= -KappaMatrix[i-1]/MassMatrix[i-1]

MainMatrix[0]=(KappaMatrix[0]+Kappa)/MassMatrix[0]

MainMatrix[0]=-(KappaMatrix[0])/MassMatrix[0]

i=symbols('i')

MainMatrix[0]= -(Kappa/MassMatrix[0])*(math.cos(math.pi*ka_by_pi)+(i*math.sin(math.pi*ka_by_pi)))

MainMatrix[n-1]= -(Kappa/MassMatrix[n-1])*(math.cos(math.pi*ka_by_pi)-(i*math.sin(math.pi*ka_by_pi)))

MainMatrix[n-1]=-(KappaMatrix[n-2])/MassMatrix[n-1]

MainMatrix[n-1]=(KappaMatrix[n-2]+Kappa)/MassMatrix[n-1]

t=symbols('t')

for j in range(0,n_atoms):

MainMatrix[j]=MainMatrix[j]-(t**2)

print(MainMatrix)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值