一​维​非​稳​态​导​热​方​程-PID温度反演

%% 已知参数 lamda = 10; % 导热系数 cp = 440; % 热容 rou = 7800; % 密度 qw = 500000; % 热流 a = lamda/rou/cp; c = qw/lamda; xspan = [0 0.012]; tspan = [0 10]; ngrid = [1000 20]; n = ngrid(1); m = ngrid(2); x = linspace(xspan(1), xspan(2), m); t = linspace(tspan(1), tspan(2), n); T0_real = 5*x; %% 调用函数计算T(x,tao) T = HeatTrans(a,c,T0_real,xspan,tspan,ngrid); Tref = T; N = zeros(n,m); Treal = Tref + N; %% 试凑法初步确定PID参数 % 这里采用的试凑法的方法是迭代20步看哪组参数效果更好