24小时热门版块排行榜    

查看: 1632  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

坠落天石

木虫 (小有名气)

[求助] 用mathematics解方程组

NSolve[{-d^2 - (4 H \[Pi]^2 Sqrt[36 + R^4/H^2] r0^3)/(d^2 m R^2) + (
    4 H \[Pi]^2 r0^3 Sqrt[R^4/H^2 + 4 r0^2])/(
    d^2 m R^2) + (-2 R^2 + 4 Sqrt[R^4 + 4 H^2 r0^2])/(
    4 H k \[Pi] r0^2 + m Sqrt[R^4 + 4 H^2 r0^2]) +
    1/R^2 8 H \[Pi]^2 ((2 Sqrt[2]
           H k \[Pi] r0^6 (d^2 m R^4 - Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]) ArcTan[(
           6 Sqrt[2] d^2 H m)/Sqrt[
           d^4 m^2 R^4 -
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6] Sqrt[
          d^4 m^2 R^4 -
           d^2 m Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]) - (2 Sqrt[2]
           H k \[Pi] r0^6 (d^2 m R^4 + Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]) ArcTan[(
           6 Sqrt[2] d^2 H m)/Sqrt[
           d^4 m^2 R^4 +
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6] Sqrt[
          d^4 m^2 R^4 +
           d^2 m Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]) + (32 Sqrt[2]
           H^4 k^2 \[Pi]^2 Sqrt[(36 H^2 + R^4)/H^2]
           r0^9 ArcTanh[(Sqrt[2] d^2 m Sqrt[36 H^2 + R^4])/Sqrt[
           
           d^4 m^2 R^4 -
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          36 H^2 + R^4] Sqrt[d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]
           Sqrt[d^4 m^2 R^4 -
           d^2 m Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]) - (32 Sqrt[2]
           H^4 k^2 \[Pi]^2 Sqrt[(36 H^2 + R^4)/H^2]
           r0^9 ArcTanh[(Sqrt[2] d^2 m Sqrt[36 H^2 + R^4])/Sqrt[
           d^4 m^2 R^4 +
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          36 H^2 + R^4] Sqrt[d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]
           Sqrt[d^4 m^2 R^4 +
           d^2 m Sqrt[d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]])) -
    1/R^2 8 H \[Pi]^2 ((2 Sqrt[2]
           H k \[Pi] r0^6 (d^2 m R^4 - Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]) ArcTan[(
           2 Sqrt[2] d^2 H m r0)/Sqrt[
           d^4 m^2 R^4 -
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6] Sqrt[
          d^4 m^2 R^4 -
           d^2 m Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]) - (2 Sqrt[2]
           H k \[Pi] r0^6 (d^2 m R^4 + Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]) ArcTan[(
           2 Sqrt[2] d^2 H m r0)/Sqrt[
           d^4 m^2 R^4 +
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6] Sqrt[
          d^4 m^2 R^4 +
           d^2 m Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]) + (32 Sqrt[2]
           H^4 k^2 \[Pi]^2 r0^9 Sqrt[(R^4 + 4 H^2 r0^2)/H^2]
           ArcTanh[(Sqrt[2] d^2 m Sqrt[R^4 + 4 H^2 r0^2])/Sqrt[
           d^4 m^2 R^4 -
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          R^4 + 4 H^2 r0^2] Sqrt[
          d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6] Sqrt[
          d^4 m^2 R^4 -
           d^2 m Sqrt[
            d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]) - (32 Sqrt[2]
           H^4 k^2 \[Pi]^2 r0^9 Sqrt[(R^4 + 4 H^2 r0^2)/H^2]
           ArcTanh[(Sqrt[2] d^2 m Sqrt[R^4 + 4 H^2 r0^2])/Sqrt[
           d^4 m^2 R^4 +
            d^2 m Sqrt[
             d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]]])/(d^2 m Sqrt[
          R^4 + 4 H^2 r0^2] Sqrt[
          d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6] Sqrt[
          d^4 m^2 R^4 +
           d^2 m Sqrt[d^4 m^2 R^8 + 256 H^4 k^2 \[Pi]^2 r0^6]])) +
    1/(d^6 m^3 R^8)
      8 H^2 k \[Pi]^2 r0^3 (d^2 m r0 R^2 (d^2 m r0 R^2 +
          2 \[Pi] (-4 H k + m R^2) r0^3) +
       2 \[Pi]^2 (-4 H k + m R^2)^2 r0^6 Log[
         d^2 m r0 R^2 + \[Pi] (4 H k - m R^2) r0^3]) == 0,
  r0^2 - (144500 R^2)/(240000 a \[Pi] + 5 \[Pi] R^2) == 0}, {a, r0}]
解上面两个方程组,我想得出a与r0之间的关系,但是运行不了结果。我也尝试求数值解,但是
x = {1, 2, 3};
z+x==0;Print[z]
结果不是z={-1,-2,-3}而是
z
z
z
不知道咋回事,求哪位高手帮帮忙!
回复此楼

» 收录本帖的淘帖专辑推荐

matlab

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)

Clear["Global`*"];
eq[x_] := {x + y == 0};
var = {y};
r1[x_?NumericQ] := (y /. NSolve[eq[x], var])[[1]]
xx = {0.1, 0.2, 0.3, 0.4, 0.5};
Map[r1[#] &, xx]
Plot[r1[x], {x, 0, 1}]
--------------------
用到你代码的时候 ,把 方程和 变量换下 前三行
多个解的时候 最好打印出来看看 换下[[1]]
或者用FindRoot 速度快些 不搜索全部解
18楼2013-10-03 17:54:01
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 坠落天石 的主题更新
信息提示
请填处理意见