%% ======== clear;clc; %% ======== % 假设系数 a=1; b=0.8; c=1; % 假设初始条件 C=[1:7]'; % 其他假设 t=10;% 计算次数 %% ======== % 建立F矩阵 F=sparse(7,7);% 建立稀疏矩阵(节省储存空间) F(1,1:2)=[b,c];% 头 for i=2:6 F(i,[i-1:i+1])=[a,b,c];% 中 end F(7,6:7)=[a,b];% 尾 %% ======== % 开始计算 for k=1:t D=F*C% 计算D矩阵 C=D;% 用D矩阵取代C矩阵 D=nan(7,1);% 将D矩阵清除(这行可有可无) end S=full(F)% 列印F 以上是我之前写数值的写法 其中我将D与C矩阵的关系写为D=FC （为了将中间的矩阵命名） 而使用稀疏矩阵的目的在此可能看不出来 但当在点数多时稀疏矩阵的优势就很明显了 以上供您参考～ ----- Sent from JPTT on my iPhone --

※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 101.15.139.6 ※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1529667860.A.921.html