| 查看: 464 | 回复: 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改变，让人失望。 |
回复此楼» 猜你喜欢
限项规定
已经有6人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
救命帖
已经有6人回复
-
存款400万可以在学校里躺平吗
已经有37人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
最失望的一年
已经有18人回复
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 51110.5
- 散金: 1007
- 红花: 45
- 沙发: 92
- 帖子: 6870
- 在线: 1023.6小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
2楼2011-06-10 10:56:37
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 51110.5
- 散金: 1007
- 红花: 45
- 沙发: 92
- 帖子: 6870
- 在线: 1023.6小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
3楼2011-06-10 11:00:54