跪求,急!用matlab 在分点数同样多的条件下,用梯形和辛普森公式求sinx在...

说明,用下面程序时,必须要先确定m的值

这两个程序都有误差估计设n=2m+1(这是因为,辛普森公式要求的)梯形公式:h=π/2/2m=π/4m;a=0;b=pi/2;x=a:h:b;y=sin(x);z1=(y(1)+y(n))h/2;

z2=sum(y(2:n-1))h;

z=z1+z2,syms

tf=sin(t);

intf=int(f,t,a,b)(计算精确值),

Fs=double(intf),Juewucha8000=abs(z8000-Fs)(计算误差)辛普森公式a=0;b=π/2;;

h=(b-a)/(2m);

x=a:h:b;

y=sin(x);z1=y(1)+y(2m+1);

z2=2sum(y(2:2:2m));

z3=4sum(y(3:2:2m));z=(z1+z2+z3)h/3,

syms

t,f=exp((-t^2)/2)/(sqrt(2pi));intf=int(f,t,a,b),

Fs=double(intf);

Juewucha=abs(z-Fs)

function T = myInt3(f,a,b,e) 

h = b-a;

T1 = h(f(a)+f(b))/2;

T2 = T1/2 + h/2f(a+h/2);

n = 2;

while abs(T2-T1) >= e

    T1 = T2;

    S = 0;    

    for j=0:n-1

        x = a + (2j+1)(b-a)/(2n);

        S = S+ f(x);

    end

    T2 = T1/2 +S(b-a)/(2n);

    n = n2;   

end

T = T2;

测试如下:

对函数3x^2+4sin(x)在区间[-2,2]上积分

>> f=@(x)(3x^2+4sin(x)');

>> myInt3(f,-2,2,1e-4)

ans =

 160000

>> clear >> syms a b c d e; >> y=dsolve('aD2y+bDy+cy=0','y(0)=d','Dy(0)=e') y = (2ae + bd + d(b^2 - 4ac)^(1/2))/(2exp((t(b - (b^2 - 4ac)^(1/2)))/(2a))(b^2 - 4ac)^(1/2)) - (2ae + bd - d(b^2 - 4ac)^(1/2))/(7795

以上就是关于跪求,急!用matlab 在分点数同样多的条件下,用梯形和辛普森公式求sinx在...全部的内容，包括:跪求,急!用matlab 在分点数同样多的条件下,用梯形和辛普森公式求sinx在...、基于matlab的变步长复化梯形积分公式、在matlab中用隐式梯形积分法对二阶微分方程（例如电路）进行编程等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！