应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » 【求助】20金币求助解方程
北京石油化工学院2026年研究生招生接收调剂公告
81/1返回列表
查看: 516  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

mltgood

铜虫 (正式写手)

[交流] 【求助】20金币求助解方程

本人不善数学，无奈建模出来这样一个方程，详见附件，现在想求出m(t)，不知如何解，请高手赐教。能提供有效帮助的，本人送20金币，谢谢。

[ Last edited by luoqiquan on 2009-12-19 at 17:47 ]
回复此楼

» 猜你喜欢
1楼 2009-04-13 12:06:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)
★ ★ ★ ★ ★
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 ]
赞一下(10人) 回复此楼
2楼2009-04-13 12:35:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mltgood

铜虫 (正式写手)
★ ★
fegg7502(金币+2,VIP+0):thank you very much! 12-21 01:11
这么复杂啊，不知道下面的方程有解吗?求n（L,t)
赞一下(10人) 回复此楼
3楼2009-04-13 16:57:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
mltgood(金币+20,VIP+0):虽然还有很多问题要请教，但是还是非常感谢您的帮助。另看您不在线，刚发了封email给你，上面有这个问题更深一步的信息。 4-13 17:37
这样左边积出来是无穷大，没法算，其实上面的结果就挺好，复杂点可以用软件慢慢算啊，如果数据少就发帖求助吧，若数据多又不会MATLAB什么的就EXCEL吧，呵呵

[ Last edited by fspdlh on 2009-4-13 at 17:21 ]
赞一下(4人) 回复此楼
4楼2009-04-13 17:14:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mltgood

铜虫 (正式写手)
应该不是无穷大，虽然上限是无穷，但是L是有一定大小的，所以当L超过一定值后n（L,t)就是零了的。
赞一下 回复此楼
5楼2009-04-13 17:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)
★ ★
mingdong(金币+2,VIP+0):感谢您的帮助！ 4-14 15:15
没收到lz的邮件啊，我又想了想，我算的这个还不大行呢，因为要对L进行积分，积分式中就不能有含L的不确定的关系式，所以lz那个n(L)必须展开，我想zl还是详细点把那些关系式都写上去吧，好让大家帮您想啊
赞一下(10人) 回复此楼
6楼2009-04-13 20:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mltgood

铜虫 (正式写手)
★ ★ ★
fegg7502(金币+3,VIP+0):thank you very much! 12-21 01:11
呵呵，好的。其实这个就是求t时刻晶体粒度分布。由于我的实验只能获得最终的粒度分布n(L,tn)，而要根据晶体粒度分布进行动力学分析，必须知道粒度随时间的分布情况。我根据其它关系，得出了以上关系（首先第二个关系式）。粒度分布是有一定范围的，并且随着时间的不同这个范围也是不一样的，所以只能将积分上限设为无穷大，只不过超过晶体的最大粒度，n（L,t)为零，最终左边积分不会为无穷。
赞一下(10人) 回复此楼
7楼2009-04-13 21:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luoqiquan

木虫 (著名写手)
斑内整理，帖子挪个位，谢谢！
赞一下 回复此楼
8楼2009-12-19 17:48:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mltgood 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312 化工或制药调剂 +7 小小墨123 2026-04-02 8/400 2026-04-02 22:15 by ZXlzxl0425
[考研] 材料调剂 +8 懒羊羊轻置玉臀 2026-04-02 8/400 2026-04-02 22:03 by liu823948201
[考研] 一志愿北京工业大学，324分求调剂 +7 零八# 2026-03-28 7/350 2026-04-02 21:09 by 1104338198
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +8 双马尾痞老板2 2026-04-02 9/450 2026-04-02 14:45 by 5896
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +4 哇呼哼呼哼 2026-04-02 4/200 2026-04-02 13:10 by yulian1987
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 377求调剂 +3 RASKIN 2026-04-02 3/150 2026-04-02 09:45 by zzchen2000
[考研] 266求调剂 +10 阳阳哇塞 2026-04-01 10/500 2026-04-02 09:07 by vgtyfty
[考研] 求调剂 +5 零八# 2026-03-27 5/250 2026-04-01 14:40 by yulian1987
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 化学0703 调剂 306分 一志愿211 +12 26要上岸 2026-03-28 12/600 2026-04-01 11:10 by chemdavid
[考研] 08工科，295，接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[硕博家园] 博一被送出联培感觉不适应怎么办 +3 全村的狗 2026-03-31 3/150 2026-04-01 10:44 by 328838485
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 266求调剂 +3 哇呼哼呼哼 2026-03-29 3/150 2026-03-31 10:06 by cal0306
[考研] 抱歉 +3 田洪有 2026-03-30 3/150 2026-03-30 19:11 by 迷糊CCPs
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
信息提示
关闭
请填处理意见
关闭