| 查看: 1287 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
解方程，没有报错，但是一直在运行，感觉计算量不该很大 已有1人参与
|
||
前面主要是在输入一些矩阵，最后希望得到一个fx=。。。的表达式，谢谢各位大神帮帮忙啦 F4F[GMCIHEQ[FORD1U_QJSM.png |
回复此楼» 猜你喜欢
11408 321分求调剂
已经有4人回复
-
301求调剂
已经有12人回复
-
085801 总分275 本科新能源 求调剂
已经有9人回复
-
085600 295分求调剂
已经有17人回复
-
289求调剂
已经有23人回复
-
291求调剂
已经有23人回复
-
304求调剂
已经有15人回复
-
285求调剂
已经有11人回复
-
材料调剂
已经有11人回复
-
292分，材料与化工，申请调剂
已经有19人回复
4楼2016-04-01 17:55:37
zdzqy
主管区长
- 应助: 4 (幼儿园)
- 金币: 176.3
- 散金: 199
- 红花: 1
- 帖子: 245
- 在线: 136.9小时
- 虫号: 1163354
- 注册: 2010-12-06
- 性别: MM
- 专业: 固体力学
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
