作者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/cn.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