推 GreenBeret: 3Q 感謝這篇的大大+上一篇推文的大大 等一下來跑跑看 06/22 21:30
→ GreenBeret: XD 忘記電腦正在跑一個超大矩陣 等跑完再RUN這個 06/22 21:31
→ zu11tune13: 沒有錢錢嗎$_$ 06/23 21:22
推 GreenBeret: 有的 等下寄出 3Q 06/24 13:17
→ zu11tune13: 收到囉3Q 06/24 23:46
%% ========
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
（為了將中間的矩陣命名）
而使用稀疏矩陣的目的在此可能看不出來
但當在點數多時稀疏矩陣的優勢就很明顯了
以上供您參考～
-----
Sent from JPTT on my iPhone
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 101.15.139.6
※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1529667860.A.921.html