※ [本文转录自 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值 该怎么写