| 查看: 2812 | 回复: 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并不大，请各位大侠帮忙解决一下~~~非常感激~~~ |
回复此楼» 猜你喜欢
急招9月入学博士，要有4级、最晚7月硕士毕业。精密电机驱控课题；学位材料
已经有6人回复
-
有没有快的中文核心比较快录用的，纳米材料光催化
已经有4人回复
-
本人42，博士刚毕业，现在找不到工作，怎么办？:(
已经有21人回复
-
河北省自然基金
已经有6人回复
-
博士申请
已经有5人回复
-
有人投过CCC中国控制会议吗？
已经有3人回复
-
3,4-二羟基苯乙酮如何纯化？
已经有5人回复
-
国基评审
已经有10人回复
-
2026-博士申请
已经有4人回复
-
考研调剂
已经有3人回复
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