使用NDSolve出错,国赛2018a第一问求解
mathematica吧
全部回复
仅看楼主
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
level 5
你是不是露了一个边界条件呀,你这个偏微分方程求解后有三个不同的未知常量,也就是说需要三个不同的边界条件来求解。
2020年09月01日 16点09分 2
吧务
level 15
反正帖子挖上来你们也不看是吧:https://tieba.baidu.com/p/3979151769
2020年09月05日 06点09分 3
1