数学建模社区-数学中国
标题:
求助mathematica大神,如何打印NDSolve的结果
[打印本页]
作者:
longnight_121
时间:
2014-2-24 21:03
标题:
求助mathematica大神,如何打印NDSolve的结果
如题,我的代码大概是:
3 t# T$ p; V$ Q2 \# h! `% ]
pde={XXXX};(*偏微分方程*)
$ G9 v7 b F9 H) v7 h2 y! v
ic={XXX};(*初始条件*)
! J6 v3 [; @- u/ v+ i
bc={XXX};(*边界条件*)
; S0 _* M5 ^9 z/ J! s
s=NDSolve[{pde, ic, bc}, {u}, {t, 0, 10}, {x, 0, 100}, {y, 0, 100}];
* _% S; |6 t, r, {+ }% q& ~3 M
6 D7 q9 N4 Y0 n+ i
也就是要解一个关于x、y的二维偏微分。
8 _% v; g1 F, N0 W) s9 N
; T! v- N* P- n' p0 Z1 Z
我自己写了一个程序可以计算出该偏微分的数值解,但是需要用NDSolve来验证我的程序是不是正确,这就需要知道NDSolve得出的数值解在 {x, 0, 100}, {y, 0, 100} 每个网格点上的具体的值是多少。但是我不知道如何去获得具体的数值,比如怎样才能在mathematica中打印一个 t=10时刻的100*100的二维矩阵,矩阵的内容就是对应的u(10,x,y)的值。希望大家帮帮我!在这里小妹先谢谢啦!!!~~~
" F( D; k) [; c+ f( c: h- k
我才开始用mathematica,我觉得数值解就保存在s中,但是不知道那个s到底是个什么数据结构,希望大神能够详细讲解。
$ G6 s, Q. L: {2 @2 s! j7 d
F2 u: u' h% B1 ?1 l
0 r6 w" T( R6 Y; C
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5