mathematica求解一阶偏微分方程

题目如图用mathematica写的代码如下,pde = Pi^2*D[u[x, t], t] == D[u[x, t], x, x];bcs = {u[x, 0] == Sin (Pi*x), u[0, t] == 0};acs = {Pi*E^(-t) + Derivative[1, 0][u][1, t] == 0};NDSolveValue[{pde, bcs, acs}, u[x, t], {x, 0, 1}, {t, 0, 1}]总是会报错:NDSolveValue::ndnum: 在 t == 0.` 处碰到一个导数的非数值量.跪求各位大佬帮忙解答,改不了重新写一下这个题的解法也行,跪谢跪谢!!!

第1个回答  2018-03-23

mathematica求解一阶偏微分方程
mathematica求解一阶偏微分方程 10 题目如图用mathematica写的代码如下,pde=Pi^2*D[u[x,t],t]==D[u[x,t],x,x];bcs={u[x,0]==Sin(Pi*x),u[0,t]==0};acs={Pi*E^(-t)+Derivative[1,0][u][1,t]==0};NDSolveValue[{pde,bcs,... 题目如图用mathematica写的代码如下,pde = Pi^2*...

mathematica解一个微分方程,急!
首先,最明显的一点,你语法错了。把你的单独的T全部改成T[h, t]。然后,在Mathematica里最好不要用大写字母开头的自定义变量,这可能和Mathematica的内置函数撞车,当然,在你的代码里这并没引发问题。语法的问题就这些。然后是方程本身的问题……你为什么给了三个边界啊?一般而言在数值求解偏微分方...

急求!!Mathematica中解偏微分方程的命令
In[1]=DSolve[{y'[x] + y[x] == a Sin[x], y[0] == 0}, y[x], x]Out[1]={{y[x] -> -(1\/2) a \\\\[ExponentialE]^-x (-1 + \\\\[ExponentialE]^ x Cos[x] - \\\\[ExponentialE]^x Sin[x])}} In[1]=linearequation = y'[t] - 3*t*y[t] ==...

mathematica基础——解方程和微分方程(组)
Mathematica提供了多种解微分方程的方法,包括使用DSolve、DSolve[{eqn1,eqn2,…},{u1,u2,…},…]来求解微分方程组等。在使用DSolve时,确保正确添加变量参数,并注意方程组和偏微分方程的特定格式。此外,了解一些基本的函数,如Root函数,对于理解解的表示和性质至关重要。解决微分方程后,您可以通过...

关于偏微分方程组mathematica或matlab编程求解,有没有简单点的办法,方...
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对偏微分数值计算...

关于mathematica求解偏微分方程组的边界条件
用Solve 或者DSolve[{方程,条件},函数,自变量]里面可以规定条件比如 再比如这个微分方程

用mathematica求解具有边界条件的连续性方程(偏微分方程),要可执行的...
不过,对于你这个问题,还存在更准确的处理方法,因为它的解析解是可求的,只要用上LaplaceTransform:eq = With[{n = n[t, x]}, D[n, t] == D[n, x, x]];bc = {n[t, 0] == 0, D[n[t, x], x] == 0 \/. x -> 1};ic = n[0, x] == n0 DiracDelta[x - 1];...

Mathematica能不能运用拉普拉斯变换解偏微分方程??如果能,如何运用啊...
确切地说,应该是它的反变换不能用初等函数(以及Mathematica内置的某些特殊函数)的有限次复合来表达。这很正常啊,就和无理数发明之前我们没法表达根号2,虚数发明之前我们没法表达x^2+1==0的根,不定义误差函数 。如果你依旧指的是解析式的话,那么,很可能。当然了,这种事最好还是试试,不过我刚...

这种偏微分方程定解问题怎么用Mathematica的NDSolve求解?
跟具体方程形式有很大关系。一种较简单的情况是1和2满足的方程可以通过适当的方法(引入分段函数等)统一成一个方程,这方面的实例还是比较多的,比如Stackexchange的帖子《How to solve different PDE defined in different regions coupled through boundary condition》(编号115594,链接不贴了,怕被吞)。

mathematica 调用前面得出的方程
需要 手动 替换。NDSolve[{x''[t] == -r*d^2*Exp[y[t]\/-10000]*x'[t]^2\/m, x[0] == 0,x'[0] == vx0} \/. s1[[1]], x[t], {t, 0, 35}]

相似回答