有一個二維 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..