| 查看: 452 | 回复: 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改变，让人失望。 |
回复此楼» 猜你喜欢
博士申请都是内定的吗？
已经有8人回复
-
博士读完未来一定会好吗
已经有32人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有8人回复
-
读博
已经有5人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 50412.5
- 散金: 1007
- 红花: 45
- 沙发: 90
- 帖子: 6745
- 在线: 1006.3小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
2楼2011-06-10 10:56:37
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 50412.5
- 散金: 1007
- 红花: 45
- 沙发: 90
- 帖子: 6745
- 在线: 1006.3小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
3楼2011-06-10 11:00:54