| 查看: 1465 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
解方程，没有报错，但是一直在运行，感觉计算量不该很大 已有1人参与
|
|||
前面主要是在输入一些矩阵，最后希望得到一个fx=。。。的表达式，谢谢各位大神帮帮忙啦 F4F[GMCIHEQ[FORD1U_QJSM.png |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
 |
10楼2016-04-06 15:58:30
2楼2016-03-31 18:16:09
|
这是代码 a=20 b=1 u=0.33 rn=n*Pi/a \[Alpha]=Sqrt[rn*(rn+Sqrt[fx/d])] \[Beta]=Sqrt[rn*(Sqrt[fx/d]-rn)] md1={{u*rn^2-\[Alpha]^2,0},{\[Alpha]^3,\[Beta]^3-(2-u)*rn^2*\[Beta]}} md3={{u*rn^2-\[Alpha]^2,0},{\[Alpha]^3,\[Beta]^3-(2-u)*rn^2*\[Beta]}} md2={{-1,0,-1,0},{\[Alpha],-\[Beta],-\[Alpha],\[Beta]},{+\[Alpha]^2-u*rn^2,0,-\[Alpha]^2+u*rn^2,0},{rn^2-\[Alpha]^3,-rn^2*\[Beta]-\[Beta]^3,-\[Alpha]^3+rn^2*\[Alpha],-rn^2*\[Beta]-\[Beta]^3}} ma1={{\[Alpha]^2-u*rn^2,-\[Beta]^2-u*rn^2},{\[Alpha]^3-(2-u)*\[Alpha]*rn^2,0}} ma3={{\[Alpha]^2-u*rn^2,-\[Beta]^2-u*rn^2},{\[Alpha]^3-(2-u)*\[Alpha]*rn^2,0}} ma2={{1,1,1,1},{\[Alpha],0,-\[Alpha],0},{-\[Alpha]^2+u*rn^2,\[Beta]^2+u*rn^2,\[Alpha]^2-u*rn^2,-\[Beta]^2-u*rn^2},{rn^2*\[Alpha]-\[Alpha]^3,0,rn^2*\[Alpha]-\[Alpha]^3,0}} md=ArrayFlatten[{{md1,0,0},{0,md2,0},{0,0,md3}}] ma=ArrayFlatten[{{ma1,0,0},{0,ma2,0},{0,0,ma3}}] mdni=Inverse[md] ms=mdni.ma mp1={{0,0,-Exp[-\[Alpha]*b],0},{0,Cot[\[Beta]*b],0,-1/Sin[\[Beta]*b]},{-Exp[-\[Alpha]*b],0,0,0},{0,-1/Sin[\[Beta]*b],0,Cot[\[Beta]*b]}} mp=ArrayFlatten[{{mp1,0},{0,mp1}}] mr=ms.mp mf=IdentityMatrix[8]-mr Solve[Det[mf]==0,fx] |
3楼2016-04-01 10:40:22
baghnac
金虫 (小有名气)
- 应助: 14 (小学生)
- 金币: 1706.4
- 红花: 6
- 帖子: 177
- 在线: 44小时
- 虫号: 1365578
- 注册: 2011-08-11
- 性别: GG
- 专业: 宇宙学
4楼2016-04-01 17:55:37
