| 查看: 451 | 回复: 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改变，让人失望。 |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有24人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 50405.3
- 散金: 1007
- 红花: 45
- 沙发: 90
- 帖子: 6743
- 在线: 1006小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
2楼2011-06-10 10:56:37
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 50405.3
- 散金: 1007
- 红花: 45
- 沙发: 90
- 帖子: 6743
- 在线: 1006小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
3楼2011-06-10 11:00:54