| 查看: 1460 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
解方程，没有报错，但是一直在运行，感觉计算量不该很大 已有1人参与
|
||
前面主要是在输入一些矩阵，最后希望得到一个fx=。。。的表达式，谢谢各位大神帮帮忙啦 F4F[GMCIHEQ[FORD1U_QJSM.png |
回复此楼» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有16人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有45人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有11人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有4人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
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
