%% ========
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
（為了將中間的矩陣命名）
而使用稀疏矩陣的目的在此可能看不出來
但當在點數多時稀疏矩陣的優勢就很明顯了
以上供您參考～