Это не трехдиагональный, а вариант пятиугольной матрицы.Вы можете использовать diag, или sparse, или spdiags, и все это может создать его.Вы можете прочитать мой blktridiag код, который можно найти на File Exchange , чтобы узнать, как эффективно построить такую матрицу.что ваша матрица имеет особую форму, матрицу Теплица.
>> toeplitz([1:5,zeros(1,4)]',[1, zeros(1,4)])
ans =
1 0 0 0 0
2 1 0 0 0
3 2 1 0 0
4 3 2 1 0
5 4 3 2 1
0 5 4 3 2
0 0 5 4 3
0 0 0 5 4
0 0 0 0 5