| 查看: 1531 | 回复: 3 | |||
ghost118木虫 (小有名气)
|
[求助]
matlab龙贝格求积 已有1人参与
|
|
小弟按照书上的算法程序输入了matlab，被积函数是ln(1+x)/x,积分区间是0到1，下面是算法程序： a=0;b=1;eps=0.5*10^-6; M=1;h=b-a;err=1;J=0;R=zeros(4,4); R(1,1)=h*(1+f(b))/2; while(err>eps) J=J+1; h=h/2; x=a+h:2*h:b-h; R(J+1,1)=R(J,1)/2+h*sum(f(x)); for K=1:min(3,J) R(J+1,K+1)=R(J+1,K)+(R(J+1,K)-R(J,K))/(4^K-1); end if(J>3) err=abs(R(J+1,4)-R(J,4)); end end quad=R(J,4); 运行出来的结果和实际相差很大，而且用这个程序做书上的例题，可以做出来，是不是被积函数是反常积分的缘故呢，希望能有大神从中指点，不胜感激 [ 发自手机版 http://muchong.com/3g ] |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 亲们，求助一下，这个非线性方程组怎么解？
已经有11人回复
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2014-11-12 18:49:39
ghost118
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1773.7
- 散金: 148
- 红花: 1
- 帖子: 286
- 在线: 136.4小时
- 虫号: 1968273
- 注册: 2012-09-02
- 性别: GG
- 专业: 钢铁冶金
3楼2014-11-13 00:16:51
ghost118
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1773.7
- 散金: 148
- 红花: 1
- 帖子: 286
- 在线: 136.4小时
- 虫号: 1968273
- 注册: 2012-09-02
- 性别: GG
- 专业: 钢铁冶金
4楼2014-11-13 00:17:33
