利用python 解对称正定矩阵方程组的平方根法
大二在读,代码水平有限,轻喷
import math
a=[]
L=[]
U=[]
ling=0
jieshu=int(input("矩阵的阶数:"))
print("请输入方程的系数矩阵")
for i in range(jieshu):onerow=[]Lonerow=[]Uonerow=[]for j in range(jieshu):x=float(input())onerow.append(x)Lonerow.append(ling)Uonerow.append(ling)a.append(onerow)L.append(Lonerow)U.append(Uonerow)
print("请输入方程等式右边的常数项")
b=[]
for i in range(jieshu):onerow1=[]for i in range(1):xx=float(input())onerow1.append(xx)b.append(onerow1)
for i in range(jieshu):if i == 0:L[0][0]=math.sqrt(a[0][0])else:L[i][0]=a[i][0]/L[0][0]
for j in range(1,jieshu):he=0for k in range(j):he+=(L[j][k])**2L[j][j]=math.sqrt(a[j][j]-he)for i in range(j,jieshu):he1=0for o in range(i):he1+=L[i][o]*L[j][o]L[i][j]=(a[i][j]-he1)/L[j][j]
print('L=',L)
print()
Y=[]
X=[]
for i in range(jieshu):onerow2=[]onerow3=[]for i in range(1):onerow2.append(0)onerow3.append(0)X.append(onerow2)Y.append(onerow3)
for i in range(jieshu):if i == 0:Y[0][0]=b[0][0]/L[0][0]else:he=0for k in range(i):he+=L[i][k]*Y[k][0]Y[i][0]=(b[i][0]-he)/L[i][i]
print('Y=',Y)
print()
for i in range(jieshu):wuhu=jieshu-iif wuhu==jieshu :X[wuhu-1][0]=Y[jieshu-1][0]/L[jieshu-1][jieshu-1]else:he=0for j in range(wuhu-1,jieshu):he+=L[j][wuhu-1]*X[j][0]X[wuhu-1][0]=(Y[wuhu-1][0]-he)/L[wuhu-1][wuhu-1]
print("X=",X)
输入案例
利用python 解对称正定矩阵方程组的平方根法相关推荐
- 解对称正定矩阵线性方程组的平方根方法
对称正定矩阵线性方程组 0. 引言 1. 对称正定矩阵线性方程组 1.1 对称正定矩阵及其三角分解法 1.2 平方根法 1.3 改进平方根法 1.4 代码实现 2. 应用举例:最小二乘法多项式拟合函数 ...
- 如何利用python解方程_Python 解方程的三种方法
# 首发于我的博客 The North. 新年第一篇,搞起. 这回写一个好久之前想做,一直搁着没做的东西-- Python 解方程(其实是放假回家,趁着家里电脑重装 LOL 的时间过来写一篇). 咱这 ...
- python解常微分方程组
美国气象学家洛伦兹(E.N.Lorenz)是混沌理论的奠基者之一.20世纪50年代末到60年代初,他的主要工作目标是从理论上进行长期天气预报研究.他在使用计算机模拟天气时意外发现,对于天气系统,哪怕初 ...
- 如何利用python解方程_如何用python解方程
用Python解数学方程,需要用到Python的一个库--SymPy库.SymPy是符号数学的Python库,它的目标是成为一个全功能的计算机代数系统,同时保持代码简洁.易于理解和扩展. 如果你的电脑 ...
- python解压gz_python解压gz
广告关闭 云服务器1核2G首年99年,还有多款热门云产品满足您的上云需求 执行程序,会在当前目录下生成aa.tar.gz文件二.解压defuntar(fname,dirs):解压tar.gz文件:pa ...
- python ks值计算_利用Python计算KS的实例详解
在金融领域中,我们的y值和预测得到的违约概率刚好是两个分布未知的两个分布.好的信用风控模型一般从准确性.稳定性和可解释性来评估模型.sOf免费资源网 一般来说.好人样本的分布同坏人样本的分布应该是有很 ...
- python能处理nc文件吗_利用python如何处理nc数据详解
前言 这两天帮一个朋友处理了些 nc 数据,本以为很简单的事情,没想到里面涉及到了很多的细节和坑,无论是"知难行易"还是"知易行难"都不能充分的说明问题,还是& ...
- python zxing 识别条码_详解利用python识别图片中的条码(pyzbar)及条码图片矫正和增强...
前言 这周和大家分享如何用python识别图像里的条码.用到的库可以是zbar.希望西瓜6辛苦码的代码不要被盗了.(zxing的话,我一直没有装好,等装好之后再写一篇) 具体步骤 前期准备 用open ...
- python解多元多次方程组_Python求解多重或非线性方程,python,多元,多次,方程组,线性方程组...
背景: 如何使用python求解多元多次方程组或者非线性方程组. 原创内容,转载注明出处!请勿用于商业用途! (上篇用python拟合2019nCov感染人数的文章被不少博主转载了,发的比较早,不少博 ...
最新文章
- c语言程序与设计苏小红,c语言程序设计苏小红
- np.c_和np.r_用法
- 乐视手机权限开启方法
- ListBox combobox的常用功能
- Python教程:json中encode与decode区别
- OpenGL创建hello Window窗口
- 解决Failed to load resource: net::ERR_CONNECTION_TIMED_OUT问题
- Unity 2017 Game Optimization 读书笔记(2)Scripting Strategies Part 2
- 《那些年啊,那些事——一个程序员的奋斗史》——78
- android view滑动到顶部并固定在顶部
- myeclipse2017安装与破解
- 模型融合之stacking和blending
- Machine Learning and Data Science 教授大师
- 2021-06-07java继承
- 安装redis 5.0.5版本 真香警告
- 使用VMware创建一个虚拟机,并安装乌班图系统
- pde中微元分析法的主要思想_有限元法(FEM)
- python读取大智慧数据_大智慧数据格式
- 集群搭建——SGE与PBS
- 深入理解Java虚拟机之走进Java之Java的发展历程