| 查看: 309 | 回复: 0 | ||
[求助]
符号运算，想求一个解析解，用solve还是reduce或者别的函数
|
|
\[Alpha] = Sqrt[rn*(rn + Sqrt[fx/d])]; \[Beta] = Sqrt[rn*(Sqrt[fx/d] - rn)]; md1 = {{\[Beta]^2 + u*rn^2, 0}, {0, \[Alpha]^3 + (u - 2)*rn^2*\[Alpha]}}; md3 = {{\[Beta]^2 + u*rn^2, 0}, {0, \[Alpha]^3 + (u - 2)*rn^2*\[Alpha]}}; ma1 = {{\[Alpha]^2 - u*rn^2, 0}, {0, \[Beta]^3 - \[Beta]*(u - 2)*rn^2}}; ma3 = {{\[Alpha]^2 - u*rn^2, 0}, {0, \[Beta]^3 - \[Beta]*(u - 2)*rn^2}}; md = ArrayFlatten[{{md1, 0}, {0, md3}}]; ma = ArrayFlatten[{{ma1, 0}, {0, ma3}}]; mdni = Inverse[md]; ms = mdni.ma; mp = {{0, -Coth[\[Alpha]*b], 0, 1/Sinh[\[Alpha]*b]}, {-Cot[\[Beta]*b], 0, -1/Sin[\[Beta]*b], 0}, {0, 1/Sinh[\[Alpha]*b], 0, -Coth[\[Alpha]*b]}, {-1/Sin[\[Beta]*b], 0, -Cot[\[Beta]*b], 0}}; mr = ms.mp; mf = IdentityMatrix[4] - mr; df = FullSimplify[Det[mf]]; Solve[df == 0, fx] 过程都没有问题，最后希望是一个fx=。。。的表达式，其中n是正整数，b,d,rn,u都是实数 用solve一直解不出来，是不是因为df==0可能个超越方程，碰到这种怎么办 |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有13人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
20