%% ========
clear;clc;
%% ========
% 假设系数
a=1;
b=0.8;
c=1;
% 假设初始条件
C=[1:7]';
% 其他假设
t=10;% 计算次数
%% ========
% 建立F矩阵
F=sparse(7,7);% 建立稀疏矩阵(节省储存空间)
F(1,1:2)=[b,c];% 头
for i=2:6
F(i,[i-1:i+1])=[a,b,c];% 中
end
F(7,6:7)=[a,b];% 尾
%% ========
% 开始计算
for k=1:t
D=F*C% 计算D矩阵
C=D;% 用D矩阵取代C矩阵
D=nan(7,1);% 将D矩阵清除(这行可有可无)
end
S=full(F)% 打印F
以上是我之前写数值的写法
其中我将D与C矩阵的关系写为D=FC
（为了将中间的矩阵命名）
而使用稀疏矩阵的目的在此可能看不出来
但当在点数多时稀疏矩阵的优势就很明显了
以上供您参考～