数学建模社区-数学中国
标题:
求助mathematica大神，如何打印NDSolve的结果
[打印本页]
作者:
longnight_121
时间:
2014-2-24 21:03
标题:
求助mathematica大神，如何打印NDSolve的结果
如题，我的代码大概是：
/ G$ y! F" v8 Q. D" Z5 X- o
pde={XXXX};(*偏微分方程*)
% \9 h( z* P) A
ic={XXX};(*初始条件*)
" y. a" w7 G& X- T& S1 t
bc={XXX};(*边界条件*)
9 z4 y+ u4 ]! F# ^" m9 Z, p
s=NDSolve[{pde, ic, bc}, {u}, {t, 0, 10}, {x, 0, 100}, {y, 0, 100}];
% E, A! L# q L+ c" j4 a- p# b
! D% o& H+ P8 i' @
也就是要解一个关于x、y的二维偏微分。
3 [4 @- j- ~0 l' C# h2 l- g" ~7 i
( {$ F2 d& N& c6 F2 W9 c
我自己写了一个程序可以计算出该偏微分的数值解，但是需要用NDSolve来验证我的程序是不是正确，这就需要知道NDSolve得出的数值解在 {x, 0, 100}, {y, 0, 100} 每个网格点上的具体的值是多少。但是我不知道如何去获得具体的数值，比如怎样才能在mathematica中打印一个 t=10时刻的100*100的二维矩阵，矩阵的内容就是对应的u(10,x,y)的值。希望大家帮帮我！在这里小妹先谢谢啦！！！~~~
. o$ I4 S: x2 v$ o* r& S6 `* W) M
我才开始用mathematica，我觉得数值解就保存在s中，但是不知道那个s到底是个什么数据结构，希望大神能够详细讲解。
# f3 F6 u, [, g" V$ [6 T Z# I
0 L M- ?+ ~4 `7 b+ Z1 X
8 u/ B) H* W! H4 g
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5