回 帖 发 新 帖 刷新版面

主题:请问如何为三对角矩阵赋值

我在编程时要用到三对角矩阵(主对角线元素及左右相临的两条对角线元素不为零,其余的矩阵元素全部为零)这是不为零的元素 U(i,i)=-(1+b*cos((i+1/2)*a))**3-(1+b*cos((i-1/2)*a))**3
                   U(i+1,i)=(1+b*cos((i-1/2)*a))**3,
                   U(i,i+1)=(1+b*cos((i+1/2)*a))**3。
        请问这个矩阵如何表示,提前谢谢了。

回复列表 (共2个回复)

沙发

今天问了老师,他说0可以不用赋值。

板凳

跟你的存储方式有关的吧。如果化成一维数组来存储的话,记住对角线元素位置,然后再处理细部了。要是存储整个矩阵的话,应该是按行存储,共三列吧。如
[X 0 0]
[X X 0]
[X X X]
.
.
.
.
.
.
[X X X]
[0 X X]
[0 0 X]
这样存储,零不管了。
不知道我这个思路,你用不用的上了。

我来回复

您尚未登录,请登录后再回复。点此登录或注册