| 查看: 1302 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
解方程，没有报错，但是一直在运行，感觉计算量不该很大 已有1人参与
|
||
前面主要是在输入一些矩阵，最后希望得到一个fx=。。。的表达式，谢谢各位大神帮帮忙啦 F4F[GMCIHEQ[FORD1U_QJSM.png |
回复此楼» 猜你喜欢
【求调剂】新能源材料本科，一志愿211，初试321
已经有3人回复
-
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
9楼2016-04-06 11:02:31
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
