| 查看: 493 | 回复: 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改变，让人失望。 |
回复此楼» 猜你喜欢
285求调剂
已经有10人回复
-
353求调剂
已经有7人回复
-
求调剂
已经有7人回复
-
材料化工306分找合适调剂
已经有11人回复
-
调剂
已经有8人回复
-
工科求调剂
已经有13人回复
-
复试调剂
已经有5人回复
-
一志愿北京化工085600 310分求调剂
已经有11人回复
-
一志愿北京2，材料与化工308求调剂
已经有10人回复
-
0817化学工程与技术求调剂，一志愿中海洋319
已经有9人回复
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 53463.7
- 散金: 1007
- 红花: 48
- 沙发: 94
- 帖子: 7195
- 在线: 1059.1小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
3楼2011-06-10 11:00:54
ygl820
至尊木虫 (知名作家)
- 应助: 133 (高中生)
- 金币: 53463.7
- 散金: 1007
- 红花: 48
- 沙发: 94
- 帖子: 7195
- 在线: 1059.1小时
- 虫号: 81941
- 注册: 2005-07-22
- 性别: GG
- 专业: 有机合成
2楼2011-06-10 10:56:37
10