| 查看: 3425 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
matlab 积分计算 算不出结果，求帮助！！！
|
||||
|
q =1.0e+008 *[2.3916; -3.6135; 1.8065; -0.6109]; syms y ; F2=int(exp(q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3)),y,0.9,1.1); F3=log(F2); F4=vpa(F3) 可以计算出F4 = 149272173.68701522388239715123227，但是当q =1.0e+009 *[2.3916; -3.6135; 1.8065; -0.6109] 时，运算不出结果，其实结果大概为149272173*10，怀疑是不是因为F2数值太大，matlab处理不了？q的取值为1.0e+008数量级时F2=2*10^64828081，q的取值为1.0e+009数量级时F2就计算不出来了，提示 numeric::int(1/exp(152725000*y^4 - (5051340117333333*y^3)/8388608 + 890400000*y^2 - (4908174540800001*y)/8388608 - 1348575000), y = 9/10..11/10) 但其实取对数以后得到的F3并不大，请各位大侠帮忙解决一下~~~非常感激~~~ |
回复此楼» 猜你喜欢
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有5人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有3人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
11楼2012-03-23 22:50:43
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢应助 2012-03-23 02:57:02
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢应助 2012-03-23 02:57:02
|
q =1.0e+009 *[2.3916; -3.6135; 1.8065; -0.6109]; syms y ; F2=int(exp((q(1)*y+1/2*q(2)*(y.^2-1)+1/3*q(3)*(y.^3-3*y)+1/4*q(4)*(y.^4-6*y.^2+3))/1e10),y,0.9,1.1); F3=log(F2); F4=vpa(F3)+1e10 Warning: Explicit integral could not be found. F4 = 9999999998.5398333654497435804148 |
2楼2012-03-22 16:42:16
ldoop 
铁杆木虫 (著名写手)
- 应助: 7 (幼儿园)
- 金币: 3485.8
- 散金: 3151
- 红花: 16
- 沙发: 1
- 帖子: 2225
- 在线: 285.5小时
- 虫号: 533053
- 注册: 2008-03-25
- 专业: 群落生态学
3楼2012-03-22 17:16:39
4楼2012-03-22 18:47:28
