24小时热门版块排行榜    

查看: 1997  |  回复: 7
本帖产生 1 个 数学EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

keepflying

至尊木虫 (著名写手)

[求助] 求带积分项(积分项除了积分变量还有一个符号常量)的非线性方程

问题描述

求变量T的值。除了T之外,其他参数均为已知。
c1=3.741832E-16;
c2=1.438786E-2;
sigma=5.6688e-008;            %斯特潘常数σ
Q=1000;                              %Q的值还有5000,10000,50000   
lambda1=2.5*1e-6;              %积分下限 λ1

积分上限如果不方便取无穷大(inf)的话,取100*1e-6。若不能获取解析解的话,能求出数值解也行。初学Matlab,水平太菜,请各位不吝赐教。
搞了两天还没算出来,快疯掉了。
积分项,纠结....
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

热分析,热设计,热测试,热变形仿真分析及测试
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

keepflying

至尊木虫 (著名写手)

keepflying: 回帖置顶 2012-04-13 20:37:08
引用回帖:
2楼: Originally posted by xiuyouxu at 2012-04-13 17:23:46:
建议采用牛顿迭代法,不过积分部分的计算还是比较麻烦,没办法,只能用数值积分方法求解。

下载了个maple,费了牛劲把积分项搞出来了。也不知道积分的结果对否:
4.189007356*10^(-10)*T*(-3.971367012*10^12-7.19400*10^5*T^2-2.070145440*10^9*T-125.*T^3+125.*T^3*(exp(1/T))^(28776/5))/(exp(1/T))^(28776/5)

可有谁有空帮忙验证下?
热分析,热设计,热测试,热变形仿真分析及测试
3楼2012-04-13 20:33:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

xiuyouxu

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
建议采用牛顿迭代法,不过积分部分的计算还是比较麻烦,没办法,只能用数值积分方法求解。
忘记自己,忘记一切烦恼(欢迎访问我的网站兆字节:http://www.mathbeta.com/)
2楼2012-04-13 17:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylberyl

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
1stopt貌似可解此类方程。
4楼2012-04-14 04:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

1stOpt的确很容易求解:
CODE:
LoopConstant Q=[1000,5000,10000,50000];
Constant c1=3.741832E-16,c2=1.438786E-2, sigma=5.6688e-008,
         lambda1=2.5*1e-6;
Function Q=sigma*T^4-int(c1*x^(-5)/(exp(c2/(x*T))-1),x=lambda1,1E8);

Q        T
1000        -364.440988765583
5000        544.966388377671
10000        -648.077906496062
50000        969.102507615287
5楼2012-10-16 16:27:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 15/750 2026-11-09 05:37 by l48hFd8Nwae2
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 05:33 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 11/550 2026-11-09 05:15 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 21:29 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-05 10/500 2026-11-06 12:32 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
信息提示
请填处理意见