level 1
Kim_Xerxes
楼主
代码如下:
D1 = 0.082/(1377*300);
D2 = 0.37/(2100*862);
NDSolve[{D[T1[x, t], t] == D1*D[T1[x, t], x, x], T1[x, 0] == 37,
T1[0, t] == 75}, T1, {x, 0, 0.06}, {t, 0, 5400}]
g1 = Plot3D[Evaluate[T1[x, t] /. %], {x, 0, 0.01}, {t, 0, 5400},
PlotRange -> All]
NDSolve[{D[T2[x, t], t] == D2*D[T2[x, t], x, x], T2[x, 0] == 37,
Evaluate[T2[x, t] /. {x -> 0.06}] ==
Evaluate[T1[x, t] /. {x -> 0.06}],
0.37*Evaluate[D[T2[x, t], x] /. {x -> 0.06}] ==
0.082*Evaluate[D[T1[x, t], x] /. {x -> 0.06}]}, T2, {x, 0.06,
0.66}, {t, 0, 5400}]
g2 = Plot3D[Evaluate[T2[x, t] /. %], {x, 0.06, 0.66}, {t, 0, 5400},
PlotRange -> All]
Show[g1, g2]
错误:
NDSolve::ibcinc: Warning: boundary and initial conditions are inconsistent.
NDSolve::ndnum: Encountered non-numerical value for a derivative at t == 0.`.
如果大佬有国赛2018a第一问求解的方法麻烦讲解一下,萌新人傻了
2020年08月29日 08点08分
1
D1 = 0.082/(1377*300);
D2 = 0.37/(2100*862);
NDSolve[{D[T1[x, t], t] == D1*D[T1[x, t], x, x], T1[x, 0] == 37,
T1[0, t] == 75}, T1, {x, 0, 0.06}, {t, 0, 5400}]
g1 = Plot3D[Evaluate[T1[x, t] /. %], {x, 0, 0.01}, {t, 0, 5400},
PlotRange -> All]
NDSolve[{D[T2[x, t], t] == D2*D[T2[x, t], x, x], T2[x, 0] == 37,
Evaluate[T2[x, t] /. {x -> 0.06}] ==
Evaluate[T1[x, t] /. {x -> 0.06}],
0.37*Evaluate[D[T2[x, t], x] /. {x -> 0.06}] ==
0.082*Evaluate[D[T1[x, t], x] /. {x -> 0.06}]}, T2, {x, 0.06,
0.66}, {t, 0, 5400}]
g2 = Plot3D[Evaluate[T2[x, t] /. %], {x, 0.06, 0.66}, {t, 0, 5400},
PlotRange -> All]
Show[g1, g2]
错误:
NDSolve::ibcinc: Warning: boundary and initial conditions are inconsistent.
NDSolve::ndnum: Encountered non-numerical value for a derivative at t == 0.`.
如果大佬有国赛2018a第一问求解的方法麻烦讲解一下,萌新人傻了