作者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值 該怎麼寫
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 61.224.191.213 (臺灣)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/MATLAB/M.1564022889.A.471.html
※ GreenBeret:轉錄至看板 Programming 07/25 10:49
1F：→ book5601: if k % 750 == 0 07/25 12:03
2F：→ book5601: 其實妳可以寫c=f*c 可以省一行 07/25 12:05
※ 編輯: GreenBeret (61.224.191.213 臺灣), 07/26/2019 23:57:41