※ [本文轉錄自 MATLAB 看板 #1TEHXfHn ]
作者: GreenBeret (綠扁帽) 看板: MATLAB
標題: [討論] 這個程式 我要每隔750次 儲存一次 怎麼寫
時間: Thu Jul 25 10:48:07 2019
%% ========
clear;clc;
%% ========
% 假設係數(done)
% 假設初始條件
C=[一個3800*1的矩陣];
% 其他假設
t=750;% 計算次數
%% ========
% 建立F矩陣
gg=[0.4664 0.1072 0.4264]; hh=eye(1700,3800);
aa=conv2(hh,gg,'same');
aa=[aa];
gg=[0.4668 0.1072 0.426]; hh=eye(601,3800);
bb=conv2(hh,gg,'same');
bb=[bb];
bb(1,:)=[];
bb=circshift(bb,[0,1699]);
gg=[0.40328 0.2336 0.36312]; hh=eye(1501,3800); %要刪除第一行
(4400-1087-448=2867)
cc=conv2(hh,gg,'same');
cc=[cc];
cc(1,:)=[];
cc=circshift(cc,[0,2299]);
F=sparse([aa;bb;cc]);
%% ========
% 開始計算
for k=1:t
D=F*C% 計算D矩陣
C=D;% 用D矩陣取代C矩陣
end
請問如果我要每隔計算750次儲存一次D值 該怎麼寫