| 查看: 323 | 回复: 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可能个超越方程，碰到这种怎么办 |
回复此楼» 猜你喜欢
现代”学阀”该如何界定
已经有12人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
申博发邮件
已经有11人回复
-
我的奶奶
已经有3人回复