24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 390  |  回复: 3
当前主题已经存档。

wuhaoshrek

金虫 (初入文坛)

[交流] 【求助】大家帮我看看这个matlab程序啊?【已解决】

RT~
小弟刚学数值分析,自己编的matlab程序,复合-3点高斯-勒让德求积公式,代码见下面!为什么计算pi总是不对啊?pi=4/(1+x^2),0~1积分。是这样的,这个程序要用很大很大的n才能算稍微准确一点,但是用其他自己编的辛普森公式之类的,只用几段就可以算出来了,应该高斯公式是精度最高的啊?哪位达人抽点时间帮我看一下代码程序?有什么问题啊??谢谢啦!金币相送!!

function I=fguass(fh,a,b,n)
%fh:输入一个函数
%a:积分下限,b:积分上限,n:积分区间区间的分段数
h=(b-a)/n;
x=a:h:b;
x1=0.5*h-sqrt(15)/10*h+x;
x2=0.5*h+x;
x3=0.5*h+sqrt(15)/10*h+x;
fx1=5/9*fh(x1);
fx2=8/9*fh(x2);
fx3=5/9*fh(x3);
fx=h/2*(fx1+fx2+fx3);
I=sum(fx);

[ Last edited by nono2009 on 2009-11-8 at 16:33 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gothicfat

木虫 (著名写手)

鄙视一切名校高傲者!!!

额。。我学的是C,我不知道同学有会的没
鄙视一切名校高傲者!!!
2楼2009-11-08 02:00:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

★ ★ ★ ★ ★ ★
wuhaoshrek(金币+1,VIP+0):非常感谢! 11-8 16:04
wuhaoshrek(金币+4,VIP+0):都给你! 11-8 16:04
nono2009(金币+1,VIP+0):谢谢热心应助! 11-8 16:32
减去求和的最后一项就好了
I=sum(fx)-fx(end);
或者
I=sum(fx(1:end-1));
3楼2009-11-08 09:31:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuhaoshrek

金虫 (初入文坛)

nono2009(金币+0,VIP+0):今后matlab相关交流请到仿真模拟版。呵呵 11-8 16:32
非常感谢楼上的师兄!!谢谢!!
4楼2009-11-08 16:05:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuhaoshrek 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见