| 查看: 560 | 回复: 7 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】20金币求助解方程
|
|||
|
本人不善数学，无奈建模出来这样一个方程，详见附件，现在想求出m(t)，不知如何解，请高手赐教。能提供有效帮助的，本人送20金币，谢谢。 [ Last edited by luoqiquan on 2009-12-19 at 17:47 ] |
回复此楼» 猜你喜欢
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
物理学I论文润色/翻译怎么收费?
已经有198人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
-
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有0人回复
★ ★ ★ ★ ★
mingdong(金币+5,VIP+0):谢谢！ 4-13 15:31
mingdong(金币+5,VIP+0):谢谢！ 4-13 15:31
|
设nL=n(L),mt=m(t),Ft=F(t),则问题转化为解积分方程，比较复杂，建议代入数据再算。 附件1为原题，附件2为结果，有一个实数解和两个复数解，实数解应该是lz想要的，即： -1/6*(-3^(1/3)*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(2/3)+Ln^4*nL^2*3^(2/3)-3*nL*Ln^2*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3))/nL/Ln/((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3) MATLAB计算过程如下： >> Ln=sym('Ln') Ln = Ln >> left=int('nL*(x-mt)^3',0,Ln) left = 1/4*nL*Ln^4-nL*mt*Ln^3+3/2*nL*mt^2*Ln^2-nL*mt^3*Ln >> mt=solve('1/4*nL*Ln^4-nL*mt*Ln^3+3/2*nL*mt^2*Ln^2-nL*mt^3*Ln=Ft','mt') mt = 1/6/nL/Ln*((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)-1/2*Ln^3*nL/((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/2*Ln -1/12/nL/Ln*((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/4*Ln^3*nL/((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/2*Ln+1/2*i*3^(1/2)*(1/6/nL/Ln*((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/2*Ln^3*nL/((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)) -1/12/nL/Ln*((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/4*Ln^3*nL/((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/2*Ln-1/2*i*3^(1/2)*(1/6/nL/Ln*((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+1/2*Ln^3*nL/((-108*Ft+3*3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)) >> mt=simplify(mt) mt = -1/6*(-3^(1/3)*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(2/3)+Ln^4*nL^2*3^(2/3)-3*nL*Ln^2*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3))/nL/Ln/((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3) 1/12*(-3^(1/3)*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(2/3)+Ln^4*nL^2*3^(2/3)+6*nL*Ln^2*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+3^(5/6)*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(2/3)*i+3*i*3^(1/6)*Ln^4*nL^2)/nL/Ln/((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3) -1/12*(3^(1/3)*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(2/3)-Ln^4*nL^2*3^(2/3)-6*nL*Ln^2*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3)+3^(5/6)*((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(2/3)*i+3*i*3^(1/6)*Ln^4*nL^2)/nL/Ln/((-36*Ft+3^(1/2)*(Ln^8*nL^2+432*Ft^2)^(1/2))*nL^2*Ln^2)^(1/3) [ Last edited by fspdlh on 2009-4-13 at 12:49 ] |
2楼2009-04-13 12:35:56
3楼2009-04-13 16:57:09
4楼2009-04-13 17:14:06
5楼2009-04-13 17:36:02
6楼2009-04-13 20:20:29
7楼2009-04-13 21:25:54
8楼2009-12-19 17:48:53
