| 查看: 478 | 回复: 2 | ||
fgwei2000铜虫 (正式写手)
|
[求助]
请大家帮忙看看我这个Mathematica程序有什么问题，为什么给不出一个完整的波形？
|
|
请大家帮忙看看我这个程序有什么问题，为什么给不出一个完整的波形？ 我这个用的是夹逼原理，就是通过调节参数EE来最终达到一个完整的正弦或者余弦波形，为什么我调节EE，波形都不发生变化呢？ 是不是这种类型的方程，我的程序本身有问题。 第一步：这一步表示先定义一个方程； Equation[\[Sigma]0_, \[Delta]0_, R_, a_, M_, fv_, k_, EE_] := G[r] (-((-1 + k) k)/ r^2 - (M + EE - \[Delta]0/(1 + Exp[(r - R)/a])) (M - EE + \[Sigma]0/(1 + Exp[(r - R)/a])) - 2 k (fv/( 4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] + D[\[Delta]0/(1 + Exp[(r - R)/a]), r]))/r - D[fv/(4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] + D[\[Delta]0/(1 + Exp[(r - R)/a]), r]), r] - (fv/( 4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] + D[\[Delta]0/(1 + Exp[(r - R)/a]), r]))^2 + (k/r + fv/(4 M) (D[\[Sigma]0/(1 + Exp[(r - R)/a]), r] + D[\[Delta]0/(1 + Exp[(r - R)/a]), r])) D[\[Sigma]0/(1 + Exp[(r - R)/a]), r]/( M - EE + \[Sigma]0/(1 + Exp[(r - R)/a]))) - D[\[Sigma]0/(1 + Exp[(r - R)/a]), r]/( M - EE + \[Sigma]0/(1 + Exp[(r - R)/a])) D[G[r], r] + D[G[r], {r, 2}]; 第二步：用NDSolve数值求解该方程。 sol = NDSolve[{Equation[-66.596, 706.4, 7.030, 0.584, 939, 0, 2, -10] == 0, G[0.001] == 0, G'[0.001] == 1}, G[r], {r, 0.001, 2}] 第三步：Plot图象 g1 = Plot[(G[r])/1 /. sol, {r, 0.001, 2}, PlotStyle -> {Thin, RGBColor[1, 0, 0]}, FrameLabel -> {r, "\!\(\*FormBox[ RowBox[{\" \", RowBox[{\"G\", \"(\", \"r\", \"  \"}]}],TraditionalForm]\)"}, Frame -> True, PlotRange -> Automatic] 输出的图象不随参数EE改变，让人失望。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有11人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 52194.9
- 散金: 1007
- 红花: 45
- 沙发: 92
- 帖子: 6975
- 在线: 1034小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
2楼2011-06-10 10:56:37
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 52194.9
- 散金: 1007
- 红花: 45
- 沙发: 92
- 帖子: 6975
- 在线: 1034小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
3楼2011-06-10 11:00:54
