| 查看: 1057 | 回复: 8 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
gengbiaolu铜虫 (正式写手)
|
[求助]
求高人指点通过mathmatica编程数值求解耦合方程组: 已有2人参与
|
||
|
求高人指点通过mathmatica编程数值求解耦合方程组: (E-1)x_1=a(x_0+x_2) (E-2)x_2=a(x_1+x_3) (E-3)x_3=a(x_2+x_4) (E-4)x_4=a(x_3+x_5) (E-5)x_5=a(x_4+x_6) (E-6)x_6=a(x_5+x_7) |x_1|^2+|x_2|^2+|x_3|^2+|x_4|^2+|x_5|^2+|x_6|^2=1-|x_0|^2-|x_7|^2 其中a为常数 a，x_0=0.1; x_7=0.2; 求x_1, x_2,...x_6 及相应的 E。 我用Mathmatic解下方程得到六组解 NSolve[{(a - 1)*x1 - 0.7*(x2 + 0.6) == 0, (a - 2)*x2 - 0.7*(x1 + x3) == 0, (a - 3)*x3 - 0.7*(x2 + x4) == 0, (a - 4)*x4 - 0.7*(x3 + 0.6) == 0, x1^2 + x2^2 + x3^2 + x4^2 + 0.6^2 + 0.6^2 == 1}, {x1, x2, x3, x4,x5,x6, a}] 我想设a变化，如a=0:0.1:1，每一个a求得相应的六组解，我如何将所得到的结果全部存下来以便画图呀？ |
回复此楼» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
-
职称评审没过，求安慰
已经有53人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
gengbiaolu
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8.7
- 散金: 90
- 红花: 2
- 帖子: 458
- 在线: 70.2小时
- 虫号: 762835
- 注册: 2009-05-04
- 专业: 原子和分子物理
|
求高人指点通过mathmatica编程数值求解耦合方程组: (E-1)x_1=a(x_0+x_2) (E-2)x_2=a(x_1+x_3) (E-3)x_3=a(x_2+x_4) (E-4)x_4=a(x_3+x_5) (E-5)x_5=a(x_4+x_6) (E-6)x_6=a(x_5+x_7) |x_1|^2+|x_2|^2+|x_3|^2+|x_4|^2+|x_5|^2+|x_6|^2=1-|x_0|^2-|x_7|^2 其中a为常数 a，x_0=0.1; x_7=0.1; 求x_1, x_2,...x_6 及相应的 E。 我用Mathmatic解下方程得到六组解 NSolve[{(E - 1)*x1 -a*(x2 + 0.1) == 0, (E - 2)*x2 - a*(x1 + x3) == 0, (E - 3)*x3 -a*(x2 + x4) == 0, (E - 4)*x4 -a*(x3 + x5) == 0, (E - 5)*x5 - a*(x4 + x6) == 0,(E - 6)*x6 -a*(x5 + 0.1) == 0,x1^2 + x2^2 + x3^2 + x4^2 +x5^2+x6^2+ 0.1^2 + 0.1^2 == 1}, {x1, x2, x3, x4,x5,x6, a}] 我想设a变化，如a=0:0.1:1，每一个a求得相应的六组解，我如何将所得到的结果全部存下来以便画图呀？ 每一个确定的a是可以求出数据的。 |
4楼2015-10-07 11:36:50
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2015-10-07 11:10:59
gengbiaolu
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8.7
- 散金: 90
- 红花: 2
- 帖子: 458
- 在线: 70.2小时
- 虫号: 762835
- 注册: 2009-05-04
- 专业: 原子和分子物理
|
求高人指点通过mathmatica编程数值求解耦合方程组: (E-1)x_1=a(x_0+x_2) (E-2)x_2=a(x_1+x_3) (E-3)x_3=a(x_2+x_4) (E-4)x_4=a(x_3+x_5) (E-5)x_5=a(x_4+x_6) (E-6)x_6=a(x_5+x_7) |x_1|^2+|x_2|^2+|x_3|^2+|x_4|^2+|x_5|^2+|x_6|^2=1-|x_0|^2-|x_7|^2 其中a为常数 a，x_0=0.1; x_7=0.1; 求x_1, x_2,...x_6 及相应的 E。 我用Mathmatic解下方程得到六组解 NSolve[{(a - 1)*x1 - 0.7*(x2 + 0.6) == 0, (a - 2)*x2 - 0.7*(x1 + x3) == 0, (a - 3)*x3 - 0.7*(x2 + x4) == 0, (a - 4)*x4 - 0.7*(x3 + 0.6) == 0, x1^2 + x2^2 + x3^2 + x4^2 +x5^2+x6^2+ 0.1^2 + 0.1^2 == 1}, {x1, x2, x3, x4,x5,x6, a}] 我想设a变化，如a=0:0.1:1，每一个a求得相应的六组解，我如何将所得到的结果全部存下来以便画图呀？ 每一个确定的a是可以求出数据的。 |
3楼2015-10-07 11:31:02
gengbiaolu
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 8.7
- 散金: 90
- 红花: 2
- 帖子: 458
- 在线: 70.2小时
- 虫号: 762835
- 注册: 2009-05-04
- 专业: 原子和分子物理
gengbiaolu: 回帖置顶 2015-10-07 12:00:01
|
求高人指点通过mathmatica编程数值求解耦合方程组: (E-1)x_1=a(x_0+x_2) (E-2)x_2=a(x_1+x_3) (E-3)x_3=a(x_2+x_4) (E-4)x_4=a(x_3+x_5) (E-5)x_5=a(x_4+x_6) (E-6)x_6=a(x_5+x_7) |x_1|^2+|x_2|^2+|x_3|^2+|x_4|^2+|x_5|^2+|x_6|^2=1-|x_0|^2-|x_7|^2 其中a为常数 a，x_0=0.1; x_7=0.1; 求x_1, x_2,...x_6 及相应的 E。 我用Mathmatic解下方程得到六组解 NSolve[{(E - 1)*x1 -a*(x2 + 0.1) == 0, (E - 2)*x2 - a*(x1 + x3) == 0, (E - 3)*x3 -a*(x2 + x4) == 0, (E - 4)*x4 -a*(x3 + x5) == 0, (E - 5)*x5 - a*(x4 + x6) == 0,(E - 6)*x6 -a*(x5 + 0.1) == 0,x1^2 + x2^2 + x3^2 + x4^2 +x5^2+x6^2+ 0.1^2 + 0.1^2 == 1}, {x1, x2, x3, x4,x5,x6, E}] 我想设a变化，如a=0:0.1:1，每一个a求得相应的六组解，我如何将所得到的结果全部存下来以便画图呀？ 每一个确定的a是可以求出数据的。 |
5楼2015-10-07 11:59:12