| 查看: 2527 | 回复: 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并不大，请各位大侠帮忙解决一下~~~非常感激~~~ |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有10人回复
-
2025年遐想
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求matlab 解方程的问题
已经有5人回复
- 关于Matlab算符号积分
已经有6人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- 求助大侠帮忙，有关matlab模拟数据
已经有29人回复
- 求助使用MATLAB计算绿色叶片的覆盖度的问题
已经有24人回复
- matlab 计算过程中工作空间变量保存及读取 求解决内存不足的方法
已经有5人回复
- 【求助】matlab 计算一个二重积分，高手们帮帮忙吧
已经有6人回复
- 求一定积分的解析解或数值解(数值解则要求给出matlab程序)
已经有5人回复
- 怎样用MATLAB计算出图中白色圆形的面积
已经有6人回复
- 求助matlab积分的问题
已经有4人回复
- 【求助】催化反应动力学matlab计算各基元反应的速率常数时，该如何避免较小量被忽略？
已经有3人回复
- 【求助】matlab中保存计算数据的问题
已经有8人回复
4楼2012-03-22 18:47:28
5楼2012-03-22 18:50:14
|
感谢版主解答！可是又发现问题了，如果把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
9楼2012-03-22 22:11:09
11楼2012-03-23 22:50:43