In[]:=
Clear[y0]
In[]:=
Solve0==
10125
14818ϕ
+150y0(-6308.297007407409`-148180.`ϕ)+150946244.5511111114`+
13000.850814965874`
ϕ
+1.11135`*^7ϕ-
0.022395039313571143`
3/2
(6959.`+(337500.`-2250.`y0)ϕ)
ϕ
,y0,​​Assumptions->ϕ==0.00145
Solve
:Solve 无法求解具有不精确系数的系统. 答案是通过求解相应的精确系统并且将结果数值化处理得到的.
Out[]=
{{y0-13.9237},{y013.8775}}
​​​​
In[]:=
Solve0==
10125
14818ϕ
+150y0(-6308.297007407409`-148180.`ϕ)+150946244.5511111114`+
13000.850814965874`
ϕ
+1.11135`*^7ϕ-
0.022395039313571143`
3/2
(6959.`+(337500.`-2250.`y0)ϕ)
ϕ
,y0
Solve
:可能存在使某些解或所有解都无效的参数值.
Out[]=
y0-
5.66966×
-41
10
(-1.35977×
40
10
-7.21673×
41
10
ϕ+2.25969×
43
10
2
ϕ
)
ϕ
+(3.57166×
-41
10
-6.1863×
-41
10
)-5.42084×
83
10
-
1.86357×
80
10
2
ϕ
-
1.7613×
82
10
ϕ
-2.71686×
85
10
ϕ-5.10621×
86
10
2
ϕ

1/3
3.60415×
125
10
+
5.08803×
120
10
3
ϕ
+
7.10064×
122
10
2
ϕ
+
4.68908×
124
10
ϕ
-6.12471×
127
10
ϕ-1.84178×
129
10
2
ϕ
-2.30769×
130
10
3
ϕ
+
1.
-3.68844×
14
10
-3.62032×
28
10
ϕ
6.99024×
51
10
-5.24268×
54
10
ϕ+1.31067×
57
10
2
ϕ
3/2
(7.67757×
53
10
+5.40501×
55
10
ϕ+1.31067×
57
10
2
ϕ
)
3
ϕ
-(2.25001×
-41
10
+3.89713×
-41
10
)
1/3
3.60415×
125
10
+
5.08803×
120
10
3
ϕ
+
7.10064×
122
10
2
ϕ
+
4.68908×
124
10
ϕ
-6.12471×
127
10
ϕ-1.84178×
129
10
2
ϕ
-2.30769×
130
10
3
ϕ
+
1.
-3.68844×
14
10
-3.62032×
28
10
ϕ
6.99024×
51
10
-5.24268×
54
10
ϕ+1.31067×
57
10
2
ϕ
3/2
(7.67757×
53
10
+5.40501×
55
10
ϕ+1.31067×
57
10
2
ϕ
)
3
ϕ
,y0-
5.66966×
-41
10
(-1.35977×
40
10
-7.21673×
41
10
ϕ+2.25969×
43
10
2
ϕ
)
ϕ
+(3.57166×
-41
10
+6.1863×
-41
10
)-5.42084×
83
10
-
1.86357×
80
10
2
ϕ
-
1.7613×
82
10
ϕ
-2.71686×
85
10
ϕ-5.10621×
86
10
2
ϕ

1/3
3.60415×
125
10
+
5.08803×
120
10
3
ϕ
+
7.10064×
122
10
2
ϕ
+
4.68908×
124
10
ϕ
-6.12471×
127
10
ϕ-1.84178×
129
10
2
ϕ
-2.30769×
130
10
3
ϕ
+
1.
-3.68844×
14
10
-3.62032×
28
10
ϕ
6.99024×
51
10
-5.24268×
54
10
ϕ+1.31067×
57
10
2
ϕ
3/2
(7.67757×
53
10
+5.40501×
55
10
ϕ+1.31067×
57
10
2
ϕ
)
3
ϕ
-(2.25001×
-41
10
-3.89713×
-41
10
)
1/3
3.60415×
125
10
+
5.08803×
120
10
3
ϕ
+
7.10064×
122
10
2
ϕ
+
4.68908×
124
10
ϕ
-6.12471×
127
10
ϕ-1.84178×
129
10
2
ϕ
-2.30769×
130
10
3
ϕ
+
1.
-3.68844×
14
10
-3.62032×
28
10
ϕ
6.99024×
51
10
-5.24268×
54
10
ϕ+1.31067×
57
10
2
ϕ
3/2
(7.67757×
53
10
+5.40501×
55
10
ϕ+1.31067×
57
10
2
ϕ
)
3
ϕ
,y0-
5.66966×
-41
10
(-1.35977×
40
10
-7.21673×
41
10
ϕ+2.25969×
43
10
2
ϕ
)
ϕ
-7.14333×
-41
10
-5.42084×
83
10
-
1.86357×
80
10
2
ϕ
-
1.7613×
82
10
ϕ
-2.71686×
85
10
ϕ-5.10621×
86
10
2
ϕ

1/3
3.60415×
125
10
+
5.08803×
120
10
3
ϕ
+
7.10064×
122
10
2
ϕ
+
4.68908×
124
10
ϕ
-6.12471×
127
10
ϕ-1.84178×
129
10
2
ϕ
-2.30769×
130
10
3
ϕ
+
1.
-3.68844×
14
10
-3.62032×
28
10
ϕ
6.99024×
51
10
-5.24268×
54
10
ϕ+1.31067×
57
10
2
ϕ
3/2
(7.67757×
53
10
+5.40501×
55
10
ϕ+1.31067×
57
10
2
ϕ
)
3
ϕ
+4.50002×
-41
10
1/3
3.60415×
125
10
+
5.08803×
120
10
3
ϕ
+
7.10064×
122
10
2
ϕ
+
4.68908×
124
10
ϕ
-6.12471×
127
10
ϕ-1.84178×
129
10
2
ϕ
-2.30769×
130
10
3
ϕ
+
1.
-3.68844×
14
10
-3.62032×
28
10
ϕ
6.99024×
51
10
-5.24268×
54
10
ϕ+1.31067×
57
10
2
ϕ
3/2
(7.67757×
53
10
+5.40501×
55
10
ϕ+1.31067×
57
10
2
ϕ
)
3
ϕ

◼
  • 由于其他解带有虚数，所以取第三个解
    • ◼
  • Since the other solutions carry imaginary numbers, take the third solution
    • In[]:=
    y0[ϕ_]=-
    5.669663743463788`*^-41(-1.359771270037994`*^40-7.216733621951227`*^41ϕ+2.2596928476699755`*^43
    2
    ϕ
    )
    ϕ
    -7.143328696215792`*^-41-5.420835445021492`*^83-
    1.8635733069744022`*^80
    2
    ϕ
    -
    1.7612965310011027`*^82
    ϕ
    -2.7168576814593794`*^85ϕ-5.106211765810845`*^86
    2
    ϕ
    
    1/3
    3.60415106260003`*^125+
    5.088030633716089`*^120
    3
    ϕ
    +
    7.100637217932908`*^122
    2
    ϕ
    +
    4.689075160885056`*^124
    ϕ
    -6.124711541755471`*^127ϕ-1.8417791612792977`*^129
    2
    ϕ
    -2.3076940411782083`*^130
    3
    ϕ
    +
    1
    3
    ϕ
    1.`
    -3.68843711342438`*^14-3.620317657903336`*^28ϕ
    6.990239970258592`*^51-5.242679977663929`*^54ϕ+1.3106699944126697`*^57
    2
    ϕ
    3/2
    (7.677572311801298`*^53+5.405008883625345`*^55ϕ+1.3106699944126697`*^57
    2
    ϕ
    )
    +4.500015095340188`*^-41
    1/3
    3.60415106260003`*^125+
    5.088030633716089`*^120
    3
    ϕ
    +
    7.100637217932908`*^122
    2
    ϕ
    +
    4.689075160885056`*^124
    ϕ
    -6.124711541755471`*^127ϕ-1.8417791612792977`*^129
    2
    ϕ
    -2.3076940411782083`*^130
    3
    ϕ
    +
    1
    3
    ϕ
    1.`
    -3.68843711342438`*^14-3.620317657903336`*^28ϕ
    6.990239970258592`*^51-5.242679977663929`*^54ϕ+1.3106699944126697`*^57
    2
    ϕ
    3/2
    (7.677572311801298`*^53+5.405008883625345`*^55ϕ+1.3106699944126697`*^57
    2
    ϕ
    )
    Out[]=
    -
    5.66966×
    -41
    10
    (-1.35977×
    40
    10
    -7.21673×
    41
    10
    ϕ+2.25969×
    43
    10
    2
    ϕ
    )
    ϕ
    -7.14333×
    -41
    10
    -5.42084×
    83
    10
    -
    1.86357×
    80
    10
    2
    ϕ
    -
    1.7613×
    82
    10
    ϕ
    -2.71686×
    85
    10
    ϕ-5.10621×
    86
    10
    2
    ϕ
    
    1/3
    3.60415×
    125
    10
    +
    5.08803×
    120
    10
    3
    ϕ
    +
    7.10064×
    122
    10
    2
    ϕ
    +
    4.68908×
    124
    10
    ϕ
    -6.12471×
    127
    10
    ϕ-1.84178×
    129
    10
    2
    ϕ
    -2.30769×
    130
    10
    3
    ϕ
    +
    1.
    -3.68844×
    14
    10
    -3.62032×
    28
    10
    ϕ
    6.99024×
    51
    10
    -5.24268×
    54
    10
    ϕ+1.31067×
    57
    10
    2
    ϕ
    3/2
    (7.67757×
    53
    10
    +5.40501×
    55
    10
    ϕ+1.31067×
    57
    10
    2
    ϕ
    )
    3
    ϕ
    +4.50002×
    -41
    10
    1/3
    3.60415×
    125
    10
    +
    5.08803×
    120
    10
    3
    ϕ
    +
    7.10064×
    122
    10
    2
    ϕ
    +
    4.68908×
    124
    10
    ϕ
    -6.12471×
    127
    10
    ϕ-1.84178×
    129
    10
    2
    ϕ
    -2.30769×
    130
    10
    3
    ϕ
    +
    1.
    -3.68844×
    14
    10
    -3.62032×
    28
    10
    ϕ
    6.99024×
    51
    10
    -5.24268×
    54
    10
    ϕ+1.31067×
    57
    10
    2
    ϕ
    3/2
    (7.67757×
    53
    10
    +5.40501×
    55
    10
    ϕ+1.31067×
    57
    10
    2
    ϕ
    )
    3
    ϕ
    In[]:=
    y0[0.00145]
    Out[]=
    1712.28+1.77636×
    -14
    10
    
    ◼
  • 这边代入相同的参数ϕ，却无法得到一样的解
    • ◼
  • This side substitutes the same parameters ϕ, but can not get the same solution