python三对角矩阵_用NumPy实现三对角矩阵算法(TDMA)

我正在使用NumPy在Python中实现TDMA.三对角矩阵存储在三个数组中：

a = array([...])

b = array([...])

c = array([...])

我想有效地计算α系数.算法如下：

# n = size of the given matrix - 1

alpha = zeros(n)

alpha[0] = b[0] / c[0]

for i in range(n-1):

alpha[i+1] = b[i] / (c[i] - a[i] * alpha[i])

但是,由于Python的for循环,这不是很有效.我想要的是这样的方法：

# n = size of the given matrix - 1

alpha = zeros(n)

alpha[0] = b[0] / c[0]

alpha[1:] = b[1:] / (c[1:] - a[1:] * alpha[:-1])

在后一种情况下,结果是不正确的,因为NumPy将最后一个表达式的右侧部分存储在temprorary数组中,然后将其元素的引用分配给alpha [1：].因此,[1：] * alpha [： – 1]只是一个零数组.

有没有办法告诉NumPy使用在其内部循环中的先前步骤计算的alpha值？

谢谢.

解决方法:

如果你要解决它的三对角系统,那么numpy.linalg中就有solve_banded().不确定这是不是你想要的.

标签：python,numpy,numerical-methods

python三对角矩阵_用NumPy实现三对角矩阵算法(TDMA)相关推荐

python 函数作用于矩阵_图解NumPy：常用函数的内在机制
原标题:图解NumPy:常用函数的内在机制来源: 机器之心支持大量多维数组和矩阵运算的 NumPy 软件库是许多机器学习开发者和研究者的必备工具,本文将通过直观易懂的图示解析常用的 NumPy 功 ...
python计算图解_图解NumPy，看这一篇就够了！
晓查编译整理量子位报道 | 公众号 QbitAI NumPy是Python的最重要的扩展程序库之一,也是入门机器学习编程的必备工具.然而对初学者来说,NumPy的大量运算方法非常难记. 最近,国 ...
python numpy矩阵索引_python – Numpy中的矩阵索引
索引后查看形状: In [295]: A=np.matrix([1,2,3]) In [296]: A.shape Out[296]: (1, 3) In [297]: A[0] Out[297]: ...
java 矩阵_(JAVA)计算机图形学与矩阵
(JAVA)计算机图形学与矩阵平移旋转变比一)矩阵 1)矩阵的乘法运算矩阵是计算机图形学的重要工具. 计算机图形学要常常用到矩阵乘法(矩阵与矩阵相乘). 矩阵(Aij)m*n 与矩阵(Bi ...
python 空矩阵_在NumPy中创建空数组/矩阵的方法
如何在NumPy中创建空数组/矩阵? 在添加行的情况下,你最好的选择是创建一个与数据集最终一样大的数组,然后向它添加数据 row-by-row: >>> import numpy & ...
python生成正态分布矩阵_使用Numpy生成正态分布数据
如何使用Python生成正态分布的人员身高体重信息? 下面是生成的数据的分布情况,身高数据在各个区间大致呈正态分布. 对应的代码如下,先用pd.cut进行各个区间分桶,然后使用sort_index按照 ...
python创建列向量_关于Numpy中的行向量和列向量详解
关于Numpy中的行向量和列向量详解行向量方式1 import numpy as np b=np.array([1,2,3]).reshape((1,-1)) print(b,b.shape) 结 ...
python数组赋值_基本Numpy数组值赋值
在我开始使用python中的数字代码之前,作为一个小练习,我正在尝试创建一个LDLT算法.只是为了"把脚弄湿".在不过,我似乎对纽比阵法缺乏基本的了解.请参见以下示例:def l ...
python维度扩展_在TensorFlow中实现矩阵维度扩展
一般TensorFlow中扩展维度可以使用tf.expand_dims().近来发现另一种可以直接运用取数据操作符[]就能扩展维度的方法. 用法很简单,在要扩展的维度上加上tf.newaxis就行了. ...
python 数列筛选_对numpy中的数组条件筛选功能详解
在程序设计中,时常会遇到数据的唯一化.相同.相异信息的提取等工作,在格式化的向量存储矩阵中南,numpy能够提供比较不错的快速处理功能. 1,唯一化的实现: In [63]: data = np.ar ...

python三对角矩阵_用NumPy实现三对角矩阵算法(TDMA)

python三对角矩阵_用NumPy实现三对角矩阵算法(TDMA)相关推荐

最新文章

热门文章