| 查看: 284 | 回复: 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可能个超越方程,碰到这种怎么办 |
» 猜你喜欢
国自然申请面上模板最新2026版出了吗?
已经有12人回复
售中科院一区文章,我:8 O 5 5 1 O 5 4,含JCR-TOP
已经有4人回复
Materials Today Chemistry审稿周期
已经有5人回复
溴的反应液脱色
已经有7人回复
推荐一本书
已经有12人回复
基金申报
已经有4人回复
计算机、0854电子信息(085401-058412)调剂
已经有4人回复
纳米粒子粒径的测量
已经有7人回复
常年博士招收(双一流,工科)
已经有4人回复
参与限项
已经有5人回复












回复此楼