24小时热门版块排行榜

返回列表

silencezy

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 714.3
帖子: 43
在线: 8小时
虫号: 4547153
注册: 2016-03-28
专业: 固体力学

[求助] 符号运算，想求一个解析解，用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可能个超越方程，碰到这种怎么办

回复此楼

» 猜你喜欢

1楼 2016-05-18 14:24:38

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

相关版块跳转我要订阅楼主 silencezy 的主题更新

返回列表