关于偏微分方程组mathematica或matlab编程求解,有没有简单点的办法,方程见图谢啦

题目详情

关于偏微分方程组mathematica或matlab编程求解,有没有简单点的办法,方程见图

谢啦

▼优质解答

答案和解析

只会Mathematica.数值解的话可以用NDSolve.你给的条件还不足,需要补充边界条件,以及具体参数.一个小例子：
eqn = D[u[t, r], t] == D[u[t, r], r, r] + D[u[t, r], r]/r;
int = u[0, r] == 1;
bon1 = u[t, 0.1] == 1;
bon2 = u[t, 1] == 1 + t;
sol = NDSolve[{eqn, int, bon1, bon2}, u, {t, 0, 1}, {r, 0.1, 1}];
Plot3D[u[t, r] /. sol, {t, 0, 1}, {r, 0.1, 1}]
注意这里把0点给挖掉了,因为那里是奇点,而Mathematica对偏微分数值计算的此类地方要求比较严,不去掉的话运算会出错.这个就是个最简单的圆形分布的,外边界的温度随时间t线性增加的温度随时间变化图了.非齐次项的引入是完全类似的.读懂了我的代码的话,你就该知道怎么加了.

看了关于偏微分方程组mathem...的网友还看了以下：

高等数学-常微分方程解法e^-lncosx〔∫secx*e^lncosxdx+C〕=x+C/cos 　2020-06-02 …

分部积分法求不定积分用分部积分法求,∫e^xsinxdx∫e^-xcosxd× 　2020-06-12 …

∫e^2xsinxdx=?网上有2种答案.1.(1/5)(e^2x)(2sinx-cosx)+c2 　2020-08-03 …

∫xe^x/(√e^x-2)dx这个我要t=(√e^x-2)最后得到∫（t^2+2）ln(t^2+ 　2020-08-03 …

求救成本会计形成性测评求大神多选题11.在下列方法中，属于成本计算的基本方法有（）A品种法B分批法C 　2020-11-02 …

编译原理题目5、设有文法（E）：E→E+T|TT→T*F|FF→(E)|i1)该文法含有左递归吗?若　2020-11-23 …

植树节到了，某市举行大型植树活动，共有1430人参加植树，要把人数分成相等的若干队，且每队人数在10 　2020-11-26 …

七张卡片正面分别标有1,2,3,4,5,6,7七个数字,后面分别标有a,b,c,d,e,f、g七个字　2020-12-26 …

部分积分法的时候是取求导数比较简单的项是u吗?e^x乘以cosxdx的时候教科书是取u=e^xdv= 　2020-12-31 …