24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tianaylk

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 141.5
帖子: 16
在线: 9小时
虫号: 865837
注册: 2009-10-08
专业: 电力系统

[求助] 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并不大，请各位大侠帮忙解决一下~~~非常感激~~~

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

求matlab 解方程的问题已经有5人回复
关于Matlab算符号积分已经有6人回复
求回归方程参数估计的最优化算法matlab代码已经有9人回复
求助大侠帮忙，有关matlab模拟数据已经有29人回复
求助使用MATLAB计算绿色叶片的覆盖度的问题已经有24人回复
matlab 计算过程中工作空间变量保存及读取求解决内存不足的方法已经有5人回复
【求助】matlab 计算一个二重积分，高手们帮帮忙吧已经有6人回复
求一定积分的解析解或数值解(数值解则要求给出matlab程序) 已经有5人回复
怎样用MATLAB计算出图中白色圆形的面积已经有6人回复
求助matlab积分的问题已经有4人回复
【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？已经有3人回复
【求助】matlab中保存计算数据的问题已经有8人回复

1楼2012-03-22 15:57:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tianaylk

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 141.5
帖子: 16
在线: 9小时
虫号: 865837
注册: 2009-10-08
专业: 电力系统

引用回帖:

2楼: Originally posted by dbb627 at 2012-03-22 16:42:16:
q =1.0e+009 *;
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 b ...

你好，非常感谢你的回复。但是你的变形是不等价的，因此得出的结果是有问题的，请再帮我考虑一下，谢谢！

赞一下

回复此楼

4楼2012-03-22 18:47:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 11 个回答

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

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

赞一下

回复此楼

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.