- 在线时间
- 2 小时
- 最后登录
- 2014-2-28
- 注册时间
- 2014-2-24
- 听众数
- 1
- 收听数
- 0
- 能力
- 0 分
- 体力
- 5 点
- 威望
- 0 点
- 阅读权限
- 10
- 积分
- 2
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 1
- 主题
- 1
- 精华
- 0
- 分享
- 0
- 好友
- 3
升级   40% 该用户从未签到 - 自我介绍
- 希望和大家多多交流
 |
如题,我的代码大概是:1 D/ M$ ^0 o1 M/ H
pde={XXXX};(*偏微分方程*)! @4 X9 N* |4 g' G2 _, |' b2 Q4 W
ic={XXX};(*初始条件*)
5 r0 ?! n- j* p4 K2 V! lbc={XXX};(*边界条件*)* ]2 H. D5 s3 E+ g; b! ]* n8 u
s=NDSolve[{pde, ic, bc}, {u}, {t, 0, 10}, {x, 0, 100}, {y, 0, 100}];- U* K) Q8 ]) }9 H. n2 l2 s
5 |( {' A9 j, G. O) e1 Q) u4 f' n' y也就是要解一个关于x、y的二维偏微分。; R! t2 D) m4 G' ]' C6 K+ n3 X
: p' ~0 {+ @! E3 A$ ]" J6 E3 G
我自己写了一个程序可以计算出该偏微分的数值解,但是需要用NDSolve来验证我的程序是不是正确,这就需要知道NDSolve得出的数值解在 {x, 0, 100}, {y, 0, 100} 每个网格点上的具体的值是多少。但是我不知道如何去获得具体的数值,比如怎样才能在mathematica中打印一个 t=10时刻的100*100的二维矩阵,矩阵的内容就是对应的u(10,x,y)的值。希望大家帮帮我!在这里小妹先谢谢啦!!!~~~0 c/ j6 U# u! ]9 q3 a* y% Z1 N
我才开始用mathematica,我觉得数值解就保存在s中,但是不知道那个s到底是个什么数据结构,希望大神能够详细讲解。, V/ Y$ k& e4 Y" A9 i
' q. F5 E) l8 K }. Q ?
9 c& ^) y, r7 c, j3 k/ m3 W |
zan
|