有一个二维 101*100000 的矩阵
是一维热传模型f(x,t)的矩阵
矩阵的Row是长度,Column是时间
在矩阵(1,1)是模型的f(0,0)
在矩阵(2,1)是模型的f(0.002,0)
在矩阵(2,2)是模型的f(0.002,0.05)
我已经可以绘制出整个矩阵的图
但是x轴和t轴都不是我要的长度
我只有办法用矩阵index当轴绘图
没办法用实际的x跟t当轴绘图
请求支援
code:
T = ones(101,1)*100 ;
S = zeros(101,100000) ;
for m=1:100000
T = A * T + B ; %%A,B值都符合矩阵运算
S(:,m) = T ;
end
mesh(S);
另外我想问如何写出一个检查是否收敛的code
没有程式经验的我只能写出
for m=1:100000
S(:,m) = T
if abs(S(50,m)-S(50,m-1))<=0.0001 %%测x在中点是否收敛
break
end
end
但他很明显不work..