| 查看: 3262 | 回复: 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并不大，请各位大侠帮忙解决一下~~~非常感激~~~ |
回复此楼» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有10人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab 解方程的问题
已经有5人回复
- 关于Matlab算符号积分
已经有6人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 求助大侠帮忙，有关matlab模拟数据
已经有29人回复
- 求助使用MATLAB计算绿色叶片的覆盖度的问题
已经有24人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 【求助】matlab 计算一个二重积分，高手们帮帮忙吧
已经有6人回复
- 求一定积分的解析解或数值解(数值解则要求给出matlab程序)
已经有5人回复
- 怎样用MATLAB计算出图中白色圆形的面积
已经有6人回复
- 求助matlab积分的问题
已经有4人回复
- 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
已经有3人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
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
5楼2012-03-22 18:50:14
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★
jjdg: 金币+1, 感谢应助 2012-03-23 02:57:21
jjdg: 金币+1, 感谢应助 2012-03-23 02:57:21
|
确实是有问题 下面变形应该等价 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))-1e8),y,0.9,1.1); F3=log(F2); F4=vpa(F3)+1e8 Warning: Explicit integral could not be found. F4 = 149272173.68701522388239715123227 >> 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))-1e8),y,0.9,1.1); F3=log(F2); F4=vpa(F3)+1e8 Warning: Explicit integral could not be found. F4 = log(numeric::int(1/exp(152725000*y^4 - (5051340117333333*y^3)/8388608 + 890400000*y^2 - (4908174540800001*y)/8388608 - 1248575000), y = 9/10..11/10)) + 100000000 >> 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))-1e9),y,0.9,1.1); F3=log(F2); F4=vpa(F3)+1e9 Warning: Explicit integral could not be found. F4 = 1492721786.6977684879699111309597 |
6楼2012-03-22 19:34:19
|
感谢版主解答！可是又发现问题了，如果把q的数值改变，例如把q(4)从-0.6109改为-5.6109，就又算不出结果了 q =1.0e+009 *[2.3916; -3.6135; 1.8065; -5.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))-1e9),y,0.9,1.1); F3=log(F2); F4=vpa(F3)+1e9 Warning: Explicit integral could not be found. > In sym.int at 64 F4 = log(numeric::int(1/exp(1402725000*y^4 - (5051340117333333*y^3)/8388608 - 6609600000*y^2 - (4908174540800001*y)/8388608 + 3401425000), y = 9/10..11/10)) + 1000000000 q =1.0e+009 *[2.3916; -3.6135; 1.8065; -5.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. > In sym.int at 64 ??? Error using ==> mupadfeval at 28 Error: Singularity [ln]; during evaluation of 'mlfloat' Error in ==> sym.vpa at 42 r = mupadfeval('mlfloat',s); 也就是说随着q的变化，指数中应该减去10的多少次方不是固定的，应该怎么办呢？还请版主再帮忙解答一下，谢谢！ |
7楼2012-03-22 20:28:48
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 程序强帖: 2
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
8楼2012-03-22 21:33:10
9楼2012-03-22 22:11:09
10楼2012-03-23 19:30:00
