C语言使用矩形法求定积分的通用函数

本文为大家分享了C语言使用矩形法求定积分的通用函数，供大家参考，具体内容如下

要求：

写一个用矩形法求定积分的通用函数，分别求：sin(x),cos(x),e^x 。

分析：

矩形法，学过高等数学就知道化曲为直的思想。将定积分化为多个函数连续的和。基本思想是将区间[a，b]化成n等分，当n越大的时候结果越准确。图形化成一小块一小块的矩形。底边长都为(b-a)/n.高为每个等分点的函数值。然后将每个矩形的面积相加即为所求。

如：

y=x;

可以通过矩形的方法来无限逼近定积分的求解，如下：

因为被分成n等分，就可以认为每一等分是一个矩形，那么每一矩形的面积为： 每一个矩形面积为：***Sn=f(x)(b-a)/n 总面积为：****S=S1+S2+…+Sn

#include <stdio.h>
#include <stdlib.h>
#include <math.h>//引入sin x,cos x,e^x的库
//使用指向函数的指针变量来复用一个通用函数
int main()
{
  int i,j,k,n,m;
  float res;
  float (*fun)(float);//定义指向函数的指针变量
  float integral (float a,float b,float (*fun)(float),int n);//n为将积分区间（b-a)分成n等分，当n的值越大的时候结果越精确，数学上的定义是取无穷大
  float fsin(float);//sin x函数的声明
  float fcos(float);//cos x函数的声明
  float fexp(float);//e^x 函数的声明
  printf ("请输入积分的下限\n");
  scanf ("%d",&m);
  printf ("请输入积分的上限\n");
  scanf ("%d",&n);
  printf("请输入你要计算的函数的具体函数\n");
  printf("1.sin(x) 2.cos(x) 3.e^x\n");
  scanf ("%d",&i);
  switch(i)
  {
  case 1:
    fun=fsin;//函数地址（入口）交给指针变量，灵活性强
    break;

  case 2:
    fun=fcos;
    break;
  case 3:
    fun=fexp;
  }
  res=(fun)(2.00);
  //printf("xxxx=%f\n",res);
  printf ("计算的结果为\n");
  res=integral(m,n,fun,200000);
  printf("res=%f",res);
  return 0;
}
float fsin(float x)
{
  // printf("fsinx=%f\n",x);
  return sin(x);
}
float fcos(float x)
{
  //printf("fcosx=%f\n",x);
  return cos(x);
}
float fexp(float x)
{
  //getchar();
  //printf("fexp=%f\n",x);
  return exp(x);
}
float integral(float a,float b,float (*fun)(float),int n)
{
  //矩形法计算,定积分转换为连续求和的形式
  int i=0;
  float x=a,s=0;
  float h=(b-a)/n;

  for (i=1;i<=n;i++)
  {
    x=x+h;
    s=s+((*fun)(x)*h);

  }
  return s;
}

结果：

以上就是本文的全部内容，希望对大家的学习有所帮助，也希望大家多多支持我们。

赞 (0)

相关推荐

  • C语言求解定积分的方法

    本文实例为大家分享了C语言求解定积分的具体方法,供大家参考,具体内容如下 题目要求: 求下面函数的定积分: 思路: 求一个函数的定积分,其实就是求它的面积,如对函数 求积分. 也就是要求出绿色部分的面积,如下: 我们可以通过矩形的方法来无限逼近定积分的求解,如下: 因为被分成n等分,就可以认为每一等分是一个矩形,那么每一矩形的面积为: 每一个矩形面积为:***Sn=f(x)(b-a)/n 总面积为:****S=S1+S2+-+Sn 通过这样的一个思路就可以完成,定积分的求解. 这样这三个定积分的

  • C语言实现黎曼和求定积分

    本文实例为大家分享了C语言程序实现黎曼和求定积分,供大家参考,具体内容如下 通过黎曼和解定积分既是把在xy平面中函数曲线与x轴区间区域划分成多个矩形并求它们的面积之和,矩形数量越多,得出的面积越精确. #include <stdio.h> #include <string.h> #include <stdlib.h> #include <math.h> int main(){ float function1(float); //函数f(x)1 float f

  • C语言实现求定积分的方法

    本文实例讲述了C语言实现求定积分的方法.分享给大家供大家参考.具体实现方法如下: 复制代码 代码如下: #include <cmath>  #include <cstdio> #define ACC 1000 float solve(float (*p)(float),float up,float down,int acc); float fun_exp(float x); float fun_qua(float x); void main(){ char selection; f

  • C语言:利用指针编写程序,用梯形法计算给定的定积分实例

    题目要求 利用指针编写程序,用梯形法计算下列公式中的定积分: 参考代码 首先说明一下指针的用处:因为所传递的参数均为数字,并不需要使用指针提高效率,故这里使用指针指向函数. 请注意calc()函数中的这一语句: double(*pfunction)() = &function; 即实现了我所描述的过程. 代码如下: #include <stdio.h> #include <math.h> double function(double x){ return (double)(

  • C语言使用矩形法求定积分的通用函数

    本文为大家分享了C语言使用矩形法求定积分的通用函数,供大家参考,具体内容如下 要求: 写一个用矩形法求定积分的通用函数,分别求:sin(x),cos(x),e^x . 分析: 矩形法,学过高等数学就知道化曲为直的思想.将定积分化为多个函数连续的和.基本思想是将区间[a,b]化成n等分,当n越大的时候结果越准确.图形化成一小块一小块的矩形.底边长都为(b-a)/n.高为每个等分点的函数值.然后将每个矩形的面积相加即为所求. 如: y=x; 可以通过矩形的方法来无限逼近定积分的求解,如下: 因为被分

  • 如何用矩形法(梯形法)求定积分

    如何用矩形法(梯形法)求定积分

    分析:高中的时候,我们学习过,可以通过矩形法或者矩形法来求定积分. 思路就是将积分区间划分成n等份,然后将这n等份近似看成矩形(或梯形),然后对所有的矩形(或梯形)的面积进行求和. 简单的例子:求函数X^2在的定积分 矩形法: 复制代码 代码如下: #include<iostream>#include<math.h>using namespace std;int main(){ float fun(float x); float a,b; cout<<"请输入

  • Go语言实现选择法排序实例

    本文实例讲述了Go语言实现选择法排序的方法.分享给大家供大家参考.具体实现方法如下: 复制代码 代码如下: package main import "fmt" func select_sort(a []int) {  len := len(a)  for i:=0; i < len-1; i++ {   k := i   j:= i + 1     for ; j < len; j++ {    if a[j] < a[k] { k = j }   }   if k

  • java递归法求字符串逆序

    本文实例讲述了java递归法求字符串逆序的方法.分享给大家供大家参考.具体实现方法如下: public static String reverseString(String x) { if(x==null || x.length()<2) return x; return reverseString(x.substring(1,x.length()))+ x.charAt(0); } 希望本文所述对大家的java程序设计有所帮助.

  • JS 用6N±1法求素数 实例教程

    用6N±1法求素数 任何一个自然数,总可以表示成为如下的形式之一: 6N,6N+1,6N+2,6N+3,6N+4,6N+5 (N=0,1,2,-) 显然,当N≥1时,6N,6N+2,6N+3,6N+4都不是素数,只有形如6N+1和6N+5的自然数有可能是素数.所以,除了2和3之外,所有的素数都可以表示成6N±1的形式(N为自然数). 根据上述分析,我们可以构造另一面筛子,只对形如6 N±1的自然数进行筛选,这样就可以大大减少筛选的次数,从而进一步提高程序的运行效率和速度. 以下代码需要自然数大于

  • 浅谈C/C++ 语言中的表达式求值

    经常可以在一些讨论组里看到下面的提问:"谁知道下面C语句给n赋什么值?" m = 1; n = m+++m++; 最近有位不相识的朋友发email给我,问为什么在某个C++系统里,下面表达式打印出两个4,而不是4和5: a = 4; cout << a++ << a; C++ 不是规定 << 操作左结合吗?是C++ 书上写错了,还是这个系统的实现有问题? 注:运行a = 4; cout << a++ << a; 如在Visua

  • python、Matlab求定积分的实现

    python求定积分 计算 from sympy import * x = symbols('x') print(integrate(sin(2*x)/(1+x**2), (x, 0, 3))) sympy库中integrate函数 integrate(f, (x, lower_bound, upper_bound)) # f-函数,x-变量,lower_bound-下限,upper_bound-上限 但是发现求不出来,如果是sin(2*x)就可以,为什么? matlab求定积分 syms x

  • python 求定积分和不定积分示例

    求f(x) = sin(x)/x 的不定积分和负无穷到正无穷的定积分 sin(x)/x 的不定积分是信号函数sig ,负无穷到正无穷的定积分为pi import math import numpy as np import matplotlib.pyplot as plt from sympy import * #用于求导积分等科学计算 def draw_plot_set():#设置画图格式 ax = plt.gca() #改变坐标轴位置 ax.spines['right'].set_color

随机推荐