看板 MATLAB 關於我們 聯絡資訊
%% ======== 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://www.ptt.cc/bbs/MATLAB/M.1564022889.A.471.html GreenBeret:轉錄至看板 Programming 07/25 10:49
book5601: if k % 750 == 0 07/25 12:03
book5601: 其實妳可以寫c=f*c 可以省一行 07/25 12:05
※ 編輯: GreenBeret (61.224.191.213 臺灣), 07/26/2019 23:57:41