Youngist的博客

本人是湖南大学的一名力学学生，2020年上数值分析课时写下的作业，希望对大家有所帮助，一起进步。1、下面是自己整理的LU解线性方程组的基本理论。2、接下来是python实现代码import numpy as npdef LU(A): ''' 生成值全位0的U矩阵，和单位矩阵L ''' L = np.eye(len(A)) U = np.zeros(np.shape(A)) for r in range(1, len(A)): # 求U的第一行和

def Jacobi(mx,mr,n=100,c=0.0001): if len(mx) == len(mr): # mx为系数矩阵，mr x = [] for i in range(len(mr)): x.append([0]) # x矩阵初始化 count = 0 # 迭代次数的计算 while count < n: nx = [] #保存单次迭代后的值的集合

2020数值分析作业，已成功实现功能，可直接复制代码就能实现1、先给出大概的理论推导过程这里取高斯点坐标为Ɛ=-Ƞ=1/，权重为W1=W2=12、边界条件和解题思路如下3、源代码双手奉上from math import *import numpy as np# from my_inv import invmat # 自己在网上看的矩阵求逆方法# from LU import my_LUsolve # LU解线性方程# 1个单元四个节点def shapefunction

2020年数值分析作业，已成功实现，可直接复制代码运行！！1、大概的理论讲解。这里主要讲总刚度矩阵的组装原理，如何得到单元刚度矩阵请看