这个微分方程怎样通过mathematica 全微分求出来，已经求了好几天了，求解答

来源：蜘蛛抓取(WebSpider) 时间：2016-12-31 02:20 标签： mathematica 微分方程
君，已阅读到文档的结尾了呢~~
如何利用mathematica求微分变换,mathematica求积分,利用微分求近似值,傅里叶变换微分性质,微分方程拉氏变换,微分变换,求微分,微分怎么求,全微分怎么求,matlab求微分方程
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
如何利用mathematica求微分变换
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口 上传我的文档
 下载
 收藏
免责声明：1，本账号发布文档来源于互联网和个人收集，仅用于技术分享交流用，版权为原作者所有。 2，文档内容纯属来自网络意见，与本账号立场无关。涉及政治言论一律相应删除，请大家监督。 3. 本账号发布文档均来源于个人收集和互联网，仅用于分享、学习、交流等使用.如在使用下载、浏览过程中出现问题，请及时站内留言或留下信箱等联系方式. 4. 本人不对文档合法性承担任何法律责任。 5. 如有文档存在侵权行为请告知，经核实将予以删除. 6，如有侵犯原您的版权,请提出指正，我们将立即删除相关资料有其它问题也欢迎与本人联系
 下载此文档
正在努力加载中...
Mathematica 求解微分方程的解、通解和初值问题
下载积分：1500
内容提示：Mathematica 求解微分方程的解、通解和初值问题
文档格式：PDF|
浏览次数：95|
上传日期： 10:37:36|
文档星级：
该用户还上传了这些文档
Mathematica 求解微分方程的解、通解和初值问题
官方公共微信(輕虹微瀾)
(魔法佑斗)
第三方登录：小木虫 --- 500万硕博科研人员喜爱的学术科研平台
&&查看话题
如何使用mathematica求解非线性微分方程组
如何采用matnematica求解非线性微分方程组，其中的微分方程为抛物线型的！
主要用到什么函数，或者有什么例子没有。
我看了NDsolve的用法了，在mathmatica帮助里面给出了他的几种用法，
但是跟我的要求有点出入
我不知道NDsolve能不能用在下面的情况呢？
还请大家指点。
就是这个吧？
恩，是很像，关键是我的边界条件并不是周期性的啊！
试试看吧，既然mathematica支持这样的格式，没准儿其他边界条件也可以解出来呢。如果不能解就会给出提示的。
这个格式也是帮助文件里面的，东西很多，我也没有全看，你也可以再找找
我的mathematica里面help中关于NDsolve中没有你上面列出的用法，会不会是版本问题，如果帮助里没有这个功能，是不是意味着我的程序中采用你给的格式也会出错。因为我现在程序运行不成，错误很奇怪，我读不懂。
rpn里面的p和n是什么？p？n？
首先应保证方程的正确性，输出来看看，没有问题了再用NDSolve求解
再有就是看看边界条件是否正确
目前我copy你的code，按照我的理解改了下，提示说
NDSolve::overdet: There are fewer dependent variables, {n,p}, than equations, so the system is overdetermined.
大约是边界不合理
rpn和gpn就是上面表达式的代号，类似变量名之类的。
我把rpn和gpn直接代到方程组里面了，运行倒是没有什么错误了，关键是图上面空空，应该是怎么回事儿呢？
up = 400;& && && && && && && && && && && && && &&&
un = 8000;& && && && && && && && && && && && && && &
Kb = 1.38*10^(-23);& && && && && && && &
h = 6.63*10^(-34);& && && && && && && &
e = 1.6*10^(-19);& && && && && && && && &
r = 0.35;& && && && && && && && && && && && && && &
y = 0.6;& && && && && && && && && && && && && && && && &
p0 = 1*10^16;& && && && && && && && && && && && &
n0 = 1*10^16;& && && && && && && && && && && && &
tn = 1*10^(-9);& && && && && && && && && && && && &
tp = 1*10^(-9);& && && && && && && && && && && && &
p1 = p0;& && && && && && && && && && && && && && && && && &&&
n1 = n0;& && && && && && && && && && && && && && && && && && &
T = 300;& && && && && && && && && && && && && && && && && && &
a = 1.65*10^4;& && && && && && && && && && && && && &&&
v = (3/532)*10^17;& && && && && && && && &&&
X = 10;& && && && && && && && && && && && && && && && && && && &&&
Pt = 100;& && && && && && && && && && && && && && && && && && && &
c1 = up*Kb*T/e;
c2 = un*Kb*T/e;
(*gpn=a*y*(1-r)*Pt*Exp/(h*v);*)
(*rpn=(p*n-p0*n0)/(tn (p+p1)+tp (n+n1));*)
NDsolve, t] ==
& & c1*D, x, x] +
& &&&a*y*(1 - r)*Pt*
& && &Exp/(h*v) - (p*n - p0*n0)/(tn (p + p1) +
& && &&&tp (n + n1)),& && &&&
& && && && && &
& &D, t] ==
& & c2*D, x, x] +
& &&&a*y*(1 - r)*Pt*
& && &Exp/(h*v) - (p*n - p0*n0)/(tn (p + p1) + tp (n + n1)),
& && && && && & p == p0, n == n0,
& && && && && & ND, x, 0] == 0, ND, x, 10] == 0,
& && && && && & ND, x, 0] == 0, ND, x, 10] == 0},
& && && && && & {p, n}, {t, 0, 1000}, {x, 0, 10}];
谢谢你，我也发现是边界条件的问题了。
浏览器进程
打开微信扫一扫
随时随地聊科研

我要回帖

更多关于 mathematica 微分方程 的文章

 

随机推荐