Matlab解热传导方程代码(共3页).doc

Sample MATLAB codes 1.%Newton Cooling Lawclear; close all; clc; h = 1;T(1) = 10; %T(0)error = 1;TOL = 1e-6;k = 0;dt = 1/10;while error TOL, k = k+1; T(k+1) = h*(1-T(k)*dt+T(k); error = abs(T(k+1)-T(k);endt = linspace(0,dt*(k+1),k+1); plot(t,T),hold on, plot(t,1,r-.)xlabel(Time),ylabel(Temperature),title(T_0 = ,num2str(T(1), , T_infty = 1),legend(Cooling Trend,Steady State)2.%Boltzman Cooling Lawclear; close all