一个计算问题请大家帮助！！！急！！！_excel吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：69,697贴子：
一个计算问题请大家帮助！！！急！！！收藏
p=b+(a1-r)*b1/(1+r)+(a2-r)*b2/(1+r)^2+...+(a10-r)*b10/(1+r)^10 已知除了r以外的所有值怎么倒退计算出r?是那种批量计算的。请大家帮帮忙吧谢谢！！！
用 "单变量求解" 。请你自行搜索。
可是用单变量得一个一个输入值吧？不能批量计算吧
你真逗。五次方程已麻烦得要命，你这十次方程就乖乖的逐个求解罢。以下摘录自维基百科：寻找五次方程的解一直是个重要的数学问题。一次方程和二次方程很早就找到了公式解，经过数学家们的努力，后来三次方程及四次方程也有了解答，但是之后很长的一段时间里没有人知道五次方程是否存在公式解。相形之下，解五次方程显得格外的困难。
后来，保罗·鲁菲尼（Paolo Ruffini）和尼尔斯·阿贝尔证明了一般的五次方程，不存在统一的根式解（即由方程的系数通过有限次的四则运算及根号组合而成的公式解）。认为一般的五次方程没有公式解存在的看法其实是不正确的。事实上，利用一些超越函数，如Θ函数或戴德金η函数即可构造出五次方程的公式解。另外，若只需求得数值解，可以利用数值方法（如牛顿法）得到相当理想的解答。
证明一般五次以上的方程式无根式解的人是埃瓦裏斯特·伽罗瓦，他巧妙地利用群论处理了上述的问题。
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或急求帮助！！计算时一直显示正在运行，没有任何反应怎么解决！…_mathematica吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：5,800贴子：
急求帮助！！计算时一直显示正在运行，没有任何反应怎么解决！…收藏
角度包含在里面了
上完整代码吧
WorkingPrecision设置为3的代码及其运行结果的图形：Subscript[F, H] = Monitor[Table[NIntegrate[Subscript[f, HD], {z, 0, d + \[Eta]},
WorkingPrecision -& 3] + Subscript[F, HI], {\[Theta], -\[Pi], \[Pi], 0.1}],ProgressIndicator[\[Theta], {-Pi, Pi}]]ListPlot[%]我的电脑只用了5分钟的样子。我不了解模型所以不知道这误差是多少，但是目测这计算质量还行？
嗯，又试了一下，你的这个积分，直接使用符号计算（也就是Integrate）其实是积不出来的吧？症结在此。因为你使用的是Integrate，所以软件默认先尝试得到该积分的符号解，这一尝试相当耗时。（我在不赋值的情况下用了数分钟，把不同的theta值赋进去估计是导致更糟糕的情况，结果你才老算不完。）如果一开始就指定使用数值算法（NIntegrate）就会好很多。就算不指定精度，只要使用NIntegrate的话也还是相当快的——写帖的这点时间里已经算完了，结果图形和之前的几乎一样所以不贴了。
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或求高人帮助计算预产期！_百度宝宝知道扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
安全计算帮助文件
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口&& 查看话题
matlab 积分计算 算不出结果，求帮助！！！
q =1.0e+008&&*;
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 = ，但是当q =1.0e+009&&* 时，运算不出结果，其实结果大概为，怀疑是不是因为F2数值太大，matlab处理不了？q的取值为1.0e+008数量级时F2=2*10^，q的取值为1.0e+009数量级时F2就计算不出来了，提示 numeric::int(1/exp(*y^4 - (3333*y^3)/8388608 + *y^2 - (0001*y)/8388608 - ), y = 9/10..11/10) 但其实取对数以后得到的F3并不大，请各位大侠帮忙解决一下~~~非常感激~~~
q =1.0e+009&&*;
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.&&
. 应该是值溢出了，建议稍微处理下表达式，比如分子都是 10^9乘以一个数，然后相加后除以10^10，那么直接10就可以了 : Originally posted by dbb627 at
q =1.0e+009&&*;
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 ... 你好，非常感谢你的回复。但是你的变形是不等价的，因此得出的结果是有问题的，请再帮我考虑一下，谢谢！ : Originally posted by ldoop at
应该是值溢出了，建议稍微处理下表达式，比如分子都是 10^9乘以一个数，然后相加后除以10^10，那么直接10就可以了 你好，非常感谢你的回复。你看到的可能是一楼作过变形后的表达式，但是那个变形和我的表达式是不等价的，原函数的形式请看主贴，貌似不太好处理...请再帮我考虑一下，谢谢！ : Originally posted by tianaylk at
你好，非常感谢你的回复。但是你的变形是不等价的，因此得出的结果是有问题的，请再帮我考虑一下，谢谢！ 确实是有问题
下面变形应该等价
q =1.0e+008&&*;
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.&&
>> q =1.0e+009&&*;
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.&&
log(numeric::int(1/exp(*y^4 - (3333*y^3)/8388608 + *y^2 - (0001*y)/8388608 - ), y = 9/10..11/10)) +
>> q =1.0e+009&&*;
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.&&
. : Originally posted by dbb627 at
确实是有问题
下面变形应该等价
q =1.0e+008&&*;
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
Warnin ... 感谢版主解答！可是又发现问题了，如果把q的数值改变，例如把q(4)从-0.6109改为-5.6109，就又算不出结果了
q =1.0e+009&&*;
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
log(numeric::int(1/exp(*y^4 - (3333*y^3)/8388608 - *y^2 - (0001*y)/8388608 + ), y = 9/10..11/10)) +
q =1.0e+009&&*;
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
during evaluation of 'mlfloat'
Error in ==> sym.vpa at 42
& && & r = mupadfeval('mlfloat',s);
也就是说随着q的变化，指数中应该减去10的多少次方不是固定的，应该怎么办呢？还请版主再帮忙解答一下，谢谢！ : Originally posted by tianaylk at
感谢版主解答！可是又发现问题了，如果把q的数值改变，例如把q(4)从-0.6109改为-5.6109，就又算不出结果了
q =1.0e+009&&*;
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 ... log(exp(100))
>> log(exp(1000))
这个是 你计算的数字实在很大，matlab没法表示，而你想通过减去一个数把它控制100左右，比较困难，我暂时还想不出较好的方法，也许你考虑将你的参数归一化，然后积分，或者考虑换个软件 : Originally posted by dbb627 at
log(exp(100))
&& log(exp(1000))
这个是 你计算的数字实在很大，matlab没法表示，而你想通过减去一个数把它控制100左右，比较困难，我暂时还想不出较好的方法 ... 是啊，计算的确有点大，但其实取对数以后并不大，不知道该如何解决这个问题，已经采用符号积分，但还是只能计算到q为1e8数量级，希望版主有空再帮忙考虑一下，非常感谢~ 问题还没完全解决，还请高人指点！}