matlab - 在matlab中制作三对角矩阵

Question

我想用matlab制作一个三对角矩阵，使用

full(gallery('tridiag', 10,  1, -4, 6, -4, 1)) 

我认为我在函数中有太多参数。还有其他方法吗？

我正在尝试制作以下矩阵:

 6    -4     1     0     0
-4     6    -4     1     0
 1    -4     6    -4     1
 0     1    -4     6    -4
 0     0     1    -4     6

Answer 1

因为你的矩阵是 pentadiagonal , 我认为最好的解决方案是使用 spdiags :

>> n = 5;
>> full(spdiags(ones(n,1)*[1,-4,6,-4,1],[-2,-1,0,1,2],n,n));
ans =
     6    -4     1     0     0
    -4     6    -4     1     0
     1    -4     6    -4     1
     0     1    -4     6    -4
     0     0     1    -4     6

full是可选的，不推荐用于较大的 n。