| 查看: 1529 | 回复: 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 ] |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 亲们，求助一下，这个非线性方程组怎么解？
已经有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
20