| 查看: 1247 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
解方程，没有报错，但是一直在运行，感觉计算量不该很大 已有1人参与
|
||
前面主要是在输入一些矩阵，最后希望得到一个fx=。。。的表达式，谢谢各位大神帮帮忙啦 F4F[GMCIHEQ[FORD1U_QJSM.png |
回复此楼» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有5人回复
-
青年基金C终止
已经有3人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
限项规定
已经有7人回复
-
救命帖
已经有8人回复
-
招博士
已经有4人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
baghnac
金虫 (小有名气)
- 应助: 14 (小学生)
- 金币: 1706.4
- 红花: 6
- 帖子: 177
- 在线: 44小时
- 虫号: 1365578
- 注册: 2011-08-11
- 性别: GG
- 专业: 宇宙学
4楼2016-04-01 17:55:37
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
5楼2016-04-02 12:30:18