C程序设计：编制求解定积分的通用函数

Posted 2023-05-15

参考技术A #include<stdio.h>
#include<math.h>
float f1(float a)

float z;
z=a+1;
return(z);

float f2(float a)

float z;
z=2*a+3;
return(z);

void integral(float a,float b,float (*fun)(float ))

float i,trace,hight1,hight2,area,sum=0;
trace=(b-a)/1000;
for ( i=a;i<b;i+=trace )

hight1=(*fun)(a);
a+=trace;
hight2=(*fun)(a);
area=(((hight1+hight2)*trace)/2);
sum+=area;

printf("sum=%f\n",sum);


void main()

float a,b,sum1;
printf("*************************\n");
printf("请输入该定积分的上下区间a和b：\n");
scanf("%f %f",&a,&b);
integral(a,b,f1(a));
integral(a,b,f2(a));
追问

但是你定义的函数不对啊

怎么定义lnx和e^x

追答

不会吧？

这个…不用定义吧

好吧，我错了，其实我也刚学，不是很懂

本回答被提问者采纳

人工智能之数学基础篇—微积分

人工智能之数学基础篇—微积分

• 1 微积分的基本思想
• 2 微积分的解释
• 3 定积分
• • 3.1 定积分的定义
  • 3.2 定积分的几何含义
• 4 定积分的性质
• 5 牛顿一莱布尼茨公式
• • 5.1 积分上限的函数及其导数
  • 5.2 牛顿一莱布尼茨公式
  • 5.3 牛顿一莱布尼茨公式的几何解释
• 6 综合实例——Python中常用的定积分求解方法
• 7 Python中数值积分常用函数简介

  微积分作为初等数学和高等数学的分水岭，在现代科学中有着极其重要的作用。17世纪末 ，微积分的概念和技巧不断扩展并被广泛应用于解决天文学、物理学中的各种实际问题。例如，求瞬时速度 、 曲线切线 、 曲线长 、曲线围成的面积 、 由面围成的体积等问题 。

  微积分是高等数学中研究函数的微分、积分及有关概念和应用的数学分支 ， 是数学的一个基础学科，内容包括极限 、 微分学、积分学及其应用 。微分学包括导数运算，是一套关于变化率的理论，在人工智能之数学基础篇—高等数学基础中篇中已经介绍过；积分学包括积分运算，为计算面积、 体积、 水压力等几何物理问题提供一套通用的方法。 本文将主要介绍微积分的一些核心概念及如何利用Python编程工具来解决定积分相关的问题 。

  在介绍微积分理论之前先来看一个小故事，微积分基本定理又称为牛顿—莱布尼茨公式，从名字就可以看出两个人对微积分的创立都做出了很大贡献，但是在到底谁是微积分的创立者这一问题上，二人曾争论不休。莱布尼茨和牛顿曾是好朋友，但因为这件事，莱布尼茨至死没和牛顿和解。那到底是谁创立了微积分呢？后人核实，两人各自独立地创立了微积分，牛顿在研究物理中的运动问题过程中创立了微积分，而莱布尼茨是从几何方面出发也独立发现了微积分，所以又称牛顿—莱布尼茨公式。

1 微积分的基本思想

  微积分以研究函数变化规律为目的，主要运用到的数学工具就是微分和积分。

  微分是对函数局部变化率的线性描述。微分学的基本思想是“无限细分”和“等效替代”，其几何意义如图1所示。设函数 $y=f(x)$，假设函数上有一点 $p$，当点 $p$ 在沿着横坐标 $\Delta x$ 移动时，其在纵坐标上的变化范围为 $\Delta y$，如图1(a)所示。特别的，当 $p$ 的移动范围足够小时（${\lim }_{\Delta x\to 0}dy=0,{\lim }_{\Delta x\to 0}dx=0$，如图1(b)所示），$p$ 点纵坐标的变化值 $\Delta y$ 与该点切线的变化距离 $dy$ 之间的差值 $\left|\Delta y-dy\right|$ 比 $\Delta y$ 要小的多，这时可以用 $p$ 点附近的一个切线段来近似替代原函数。

图1 微分的几何意义

  $p$ 点切线的斜率是导数，而 $f{\left(x\right)}^{{}^{\prime }}=\frac{dy}{dx}$，所以 $dy=f{\left(x\right)}^{{}^{\prime }}dx$。根据无穷小的概念，$dx$，$dy$ 都是微分。

  积分是在确定函数的导数基础上通过一定的数学方法对原函数进行求解的过程。积分的基本思想是通过微分的 “无限求和” 来进行的 。 微分是对函数的求导过程，因此可以将积分看作微分的一个逆向过程，其几何意义如图 2 所示。设函数 $y=f(x)$ 把区间 [$a$，$b$] 均分为 4 份，如图2(a)所示 ，整体等于部分之和，即 $f(b)-f(a)=\sum \Delta y$。继续细分，如图2(b)所示，把区间均分到最细，即间隔为 $dx$，对应的 $\Delta y$ 也变成了 $dy$ ， 所以 $f(b)-f(a)=\sum dy$。根据前面微分的无限细分可知 $dy=f{\left(x\right)}^{{}^{\prime }}dx$，因此 $f(b)-f(a)=\sum f{\left(x\right)}^{{}^{\prime }}dx$。

图2 积分的几何意义

2 微积分的解释

  下面来讨论一个问题，如何求解曲边梯形的面积？

  设曲线方程为 $y=f(x)$，在区间 [$a$，$b$] 上非负、连续 。 由直线 $x=a$ ，$x=b$ , $y=0$ 及曲线 $y=f(x)$ 所围成的平面图形如图 3 所示，称之为曲边梯形，其中曲线弧段称为曲边梯形的曲边。

图3 曲边梯形

  曲边梯形面积不像矩形 、圆等规则图形有特定的公式求解，但可以采用以直代曲的思想来进行求解。

  分别以 4 个小矩形和 9 个小矩形替代，如图4所示 。

图4 曲边梯形的分割

  可以看出 9 个小矩形比 4 个小矩形求面积的精确度高，那么如果无限分割呢？在 $a$ 与 $b$ 之间分割出 $n$ 等份，$n$ 趋近于无穷，是不是会无限接近曲边梯形的面积呢? 曲边梯形的无限分割如图5所示。

图5 曲边梯形的无限分割

  因此，将其极限值
$$A=\underset{n\to \infty }{\lim }\sum _{i=1}^{n}f\left({\xi }_i\right)\Delta x_i,x_{i-1}<{\xi }_i<x_i$$定义为曲边梯形的面积。

  积分的基本思想就是“无限求和”，需要尽可能地将每一个矩形的底边分割成无穷小，这样所求的结果即为曲边梯形的面积。数学家莱布尼茨为了体现求和的感觉，将 Sum 中的第一个字母 s 拉长了，简写成
$$Sum\left[f(x)\Delta x\right]\Rightarrow {\int }_{um}^{}f(x)dx$$就成了常见的积分符号 $\int$。

  下面再来讨论一个问题：求变速直线运动的路程。

  设某物体做直线运动，已知速度 $\nu =\nu (t)$ 是时间间隔 [$t_1,t_2$] 上 $t$ 的连续函数，且 $v(t)⩾0$，计算在这段时间内物体所经过的路程。

  匀速直线运动的公式：路程=速度 x 时间。但是在我们的问题中，速度不是常数，而是随时间变化的变量，因此所求路程 s 不能直接按匀速直线运动的路程公式来计算。物体运动的速度函数 $\nu =\nu (t)$ 是连续变化的，在很短的时间内，速度的变化很小，因此如果把时间间隔分小，在小段时间内，以匀速运动近似代替变速运动，就可以算出各部分路程的近似值，再求和得到整个路程的近似值，最后，通过对时间间隔无限细分的极限过程，求得物体在时间间隔 [$t_1,t_2$] 内的路程结果如下。
s = lim ⁡ n → ∞ ∑ i = 1 n v ( ξ i ) Δ t i ,   t i − 1 < ξ i < t i s=\\lim_{n\\rightarrow \\infty }\\sum_{i=1}^{n}v\\left ( \\xi _{i} \\right )\\Delta t_{i},\\: t_{i-1}< \\xi _{i} < t_{i} s=n→∞lim​i=1∑以上是关于C程序设计：编制求解定积分的通用函数的主要内容，如果未能解决你的问题，请参考以下文章

c语言 求定积分的通用函数

用C语言编写一个求定积分的程序

用C语言求定积分

编写一个求定积分的通用函数，用它分别求以下几个积分：

C语言如何求定积分?

C语言求解定积分（蒙特卡罗方法）