matlab课下作业(2)--2009.12.25 10:10

 matlab课下作业--2009.12.25 10:10

1.梯形法求pi

 

2.复合辛普森发法求pi

二次差值插值函数x(2k)至x(2k+1)对应的sk(x)的积分为
(h/3)*(f(2k)+4*f(2k+1)+f(2k+2));
推导之后，整个区间上的s(k)的和为

(h/3)(f(0)+f(2m)+4*(f1+f3+f5...)+2*(f2+f4+f6...));
这里的h为(b-a)/2m;

 

function sinpsen(f,a,b,m) % function sinpsen % f=4./(1+x.^2); % a=0;b=1;m=2000; s1=0;s2=0; i=0:m-1; i1=a+(2*i+1)*(b-a)/2/m; s1=4./(1+i1.^2); i=1:m-1; i2=a+(2*i)*(b-a)/2/m; s2=4./(1+i2.^2); s=((b-a)/6/m)*(4/(1+a^2)+4/(1+b^2)+4*sum(s1)+2*sum(s2)) % s =3.1416;
3.蒙特卡洛方法求pi

%概率的方法
N=10000000; x=rand(N,1); y=rand(N,1); count=0; for i=1:N if (x(i)^2+y(i)^2<=1) %随机投点落在圆内产生的概率为count/n count=count+1; end end PI=4*count/N PI =3.1413
                                                                             2009.12.25 10:10