[5000P] 这个程式 我要每隔750次 储存一次 怎么写

楼主: GreenBeret (绿扁帽)   2019-07-25 10:48:07
%% ========
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值 该怎么写
作者: book5601 (苍澜野兽)   2019-07-25 12:03:00
if k % 750 == 0其实妳可以写c=f*c 可以省一行

Links booklink

Contact Us: admin [ a t ] ucptt.com