当前位置：文档之家 › 复合辛普森求积公式

复合辛普森求积公式

  • 格式：doc
  • 大小：24.71 KB
  • 文档页数：2

复合辛普森求积公式

复合辛普森求积公式，听起来是不是有点让人摸不着头脑？别急，让我慢慢给您讲讲。

咱先来说说这公式到底是干嘛的。简单点说，它就是用来计算一些复杂图形面积或者一些函数积分的好帮手。

比如说，您想知道一个形状不规则的土地面积，直接测量那可太难了。这时候复合辛普森求积公式就派上用场啦。

我记得有一次给学生们讲这个公式的时候，有个小同学瞪着大眼睛问我：“老师，这公式真能算出那么奇怪的形状的面积吗？”我笑着回答：“那当然，只要咱们用对了方法，啥难题都能解决。”

那咱们来具体看看这公式怎么用。它的原理其实就是把一个大的区间分成很多小的区间，然后在每个小区间上用特定的方法去近似计算面积或者积分。

想象一下，把一块大蛋糕切成很多小块，每小块咱们都能大概算出它的大小，加起来不就差不多知道整个蛋糕的大小了嘛。

在实际应用中，可不能马虎。参数要选对，计算要仔细，不然一个小错误可能就导致结果相差十万八千里。

就像有一回，我让学生们自己动手用复合辛普森求积公式计算一个函数的积分。有个粗心的孩子，计算过程中少算了一个区间，结果出来和正确答案差了好多。我就告诉他：“这就好比你做算术题，少加了一个数，能对吗？”

而且，要真正掌握这个公式，得多做练习题。光听我在这儿讲可不行，得自己动手去算，去琢磨。

学习复合辛普森求积公式，就像是一场探险。刚开始可能会觉得有点难，但是只要您坚持，一点点地去理解，去尝试，就会发现其中的乐趣和奥秘。

总之，复合辛普森求积公式虽然有点复杂，但只要咱们有耐心，多练习，就能把它拿下，让它成为我们解决数学问题的有力武器！

合集下载

c语言版本复合梯形法、辛普森法求积公式

c语言版本复合梯形法、辛普森法求积公式

复合梯形法求积公式（f（x）=sinx/x）

#include "stdio.h"

#include "math.h"

void main()

{

int k;

double a,b,n,h,x;

double f,f1,f2,T,F=0.0;

printf("请输入a b n :");

scanf("%lf%lf%lf",&a,&b,&n);

h=(b-a)/n;

for(k=1;k

{

x=a+k*h;

f=sin(x)/x;

printf("\nf=%0.7f",f);

F+=f;

printf("\nF=%0.7f",F);

}

if(a!=0.0)

f1=sin(a)/a;

else(f1=1.0);

printf("\nf1=%0.7f",f1);

f2=sin(b)/b;

printf("\nf2=%0.7f",f2);

T=0.5*h*(f1+2*F+f2);

printf("\nT=%0.7f",T);

}

复合辛普森公式（f（x）=sinx/x）

#include "stdio.h"

#include "math.h"

void main()

{

int k;

double a,b,n,h;

double x,f,f1,f2,F1=0.0,F2=0.0,S=0.0;

printf("请输入a b n :");

scanf("%lf%lf%lf",&a,&b,&n);

h=(b-a)/n;

for(k=0;k

{

x=a+k*h+0.5*h;

f=sin(x)/x;

printf("\nf=%0.7f",f);

F1+=f;

printf("\nF1=%0.7f",F1);

}

for(k=1;k

{

f=0;

x=a+k*h;

f=sin(x)/x;

printf("\nf=%0.7f",f);

(完整版)复合梯形公式与复合辛普森公式对比

(完整版)复合梯形公式与复合辛普森公式对比

1

SHANGHAI JIAO TONG UNIVERSITY

题目名称：复合梯形公式与复合辛普森公式对比

学生姓名:

学生学号:

班 级:

学院(系):

2

目录

1. 概述 .............................................................................................................. 3

2. 问题提出 ...................................................................................................... 4

3. 算法推导 ...................................................................................................... 5

4. 算法框图 ...................................................................................................... 6

4.1复合梯形公式算法流程图 ...................................................................... 6

4.2 复合辛普森公式算法流程图 ................................................................. 7

5. MATLAB源程序 ......................................................................................... 8

复合辛普森求积

复合辛普森求积

第三次实验

实验名称 复合辛普森求积计算积分 实验时间 2012.05.06

姓名 班级 数应二班 学号 成绩

一、实验目的，内容

1.目的：

学习并理解复合辛普森求积计算积分的matlab实现。

2. 内容：

用matlab计算积分104ln9xxdx（ 精确值-0.4444），并求出达到。

二．代码

1.

function Sn=ComSimpson(a,b,n)

%复合辛普森求积

%f表示被积函数，本题中由f.m调用

%a,b分别表示积分上下限

%n表示区间分割次数

%sn表示该方法计算所返回的值

h=(b-a)/n;

for k=0:n

x(k+1)=a+k*h;

x_k(k+1)=x(k+1)+1/2*h;

if (x(k+1)==0)|(x_k(k+1)==0)

x(k+1)=10^(-10);

x_k(k+1)=10^(-10)；%误差精度

end

end

S1=h/6*(f(x(1))+f(x(n+1)));%复合辛普森求积公式第一部分

for i=2:n

F_1(i)=h/3*f(x(i));

end

for j=1:n

F_2(j)=2*h/3*f(x_k(j));

end

S2=sum(F_1)+sum(F_2);%复合辛普森求积公式第二部分

Sn=S1+S2;%算得值

f的表达式由f.m（见附）文件调用。

附：function y=f(x)

y=sqrt(x).*log(x);

三.数值结果:

在命令窗口输入如下指令时，出现如下结果

Sn=ComSimpson(0,1,300)

Sn =-0.4438

Sn=ComSimpson(0,1,700)

Sn =-0.4442

数值计算考题五——复合梯形求积公式与复合辛普森求积公式求积分

数值计算考题五——复合梯形求积公式与复合辛普森求积公式求积分

数值计算考题五

1. 分别用复合梯形求积公式与复合辛普森求积公式求积分I=102xesinx dx的近似值，要求误差不超过=0.510-5.

解：方法一： 复合梯形求积公式

复合梯形求积公式是将积分区间划分为n个很小的区间，然后将各个小区间的面积相加而得到在整个积分区间上的积分，当分成的小区间数n时，求得的面积就等于积分的精确值。

由复合梯形求积公式的余项RnT可得满足精度要求0.510-5时区间ba,被分成的区间数n的最小值为700，所以在编程时循环次数应大于等于这个值，方可满足精度要求。以下是编写的C语言程序：

#include

#include

void main()

{

int n=700,i;

double x,f=0.0,t,h,T=0.0,c=2.0,a=0.0,b=1.0;

h=(b-a)/n;

for(i=0;i

{

x=a+i*h;

f=f+exp(pow(x,c))*sin(x);

}

t=(h/2)*(2*f+sin(1)*exp(1));

printf("T=%f\n",t);

}

输出结果为T=0.778746.

方法二：复合辛普森求积公式：

复合辛普森求积法是将积分区间分割之后，在每个小区间[xi,xi+1]上运用辛普森求积公式。以下是编写的c语言程序：

#include

#include

void main()

{

int n=700,i;

double x1,x2,f1=0.0,f2=0.0,t,h,T=0.0,c=2.0,a=0.0,b=1.0;

h=(b-a)/n;

for(i=0;i

{

x1=a+i*h;

x2=a+(i+0.5)*h;

f1=f1+exp(pow(x1,c))*sin(x1); f2=f2+exp(pow(x2,c))*sin(x2);