偏微分一維熱傳導問題-資料下載頁
2025-03-27 00:21本頁面
【正文】 for k = 1:M2 j = 2:M1。 B(j1,1) = U(k,j)。 x = A\B。 for j = 2:M1 U(k+1,j) = x(j1)。 %k+1時刻的不同位置的溫度分布 end。end。 %作圖x = 0:h:1。y = 0:t:2。[xx,yy]=meshgrid(x,y)。figure(1)。surf(xx,yy,U)。shading flattitle(39。一維熱傳導方程數(shù)值解溫度分布圖39。)。xlabel(39。位置x39。)。ylabel(39。時間t39。)。zlabel(39。溫度T39。)。 figure(2)s = 0。for i= 1:8000 s = s+(200*(1(1)^i))/(i*pi)*sin(i*pi*xx).*exp(i^2*pi^2*yy)。end。surf(xx,yy,s)。title(39。一維熱傳導方程解析解溫度分布圖39。)。xlabel(39。位置x39。)。ylabel(39。時間t39。)。zlabel(39。溫度T39。)。 figure(3)x = 0:h:1。y = 0:t:2。[xx,yy] = meshgrid(x,y)。dd = Us。surf(xx,yy,dd)。title(39。一維熱傳導方程誤差溫度分布圖39。)。xlabel(39。位置x39。)。ylabel(39。時間t39。)。zlabel(39。誤差（數(shù)值解減解析解）39。)。 figure(4)z = zeros(M2+1,1)。if mod((M1+1),2)~=0 i = 1:M2+1。 z(i,1) = U(i,M1/2)。else i = 1:M2+1。 z(i,1) = U(i,(M1+1)/2)。end。plot(z)。title(39。溫度隨時間增加的趨勢圖39。)。xlabel(39。時間t39。)。ylabel(39。溫度T39。)。17
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1