Introductory university-level calculus, linear algebra, abstract algebra, probability, statistics, and stochastic processes.
Riemann integral
Get link
Facebook
X
Pinterest
Email
Other Apps
-
You may have learned that the integral is the inverse operation of differentiation. Here, we define the integral as the calculation of area. This approach has an advantage that can be easily extended to higher dimensions.
Definition (Partition of an interval)
The partition of the closed interval \([a,b]\) is a finite sequence \(\Delta = \{x_n\}\) of the form
Let \(f(x)\) be a bounded function on \([a,b]\). Let us define the following quantities:
\[
\begin{eqnarray}
M_i &=& \sup\{f(x)\mid x_{i} \leq x \leq x_{i+1}\},\\
m_i &=& \inf\{f(x)\mid x_{i} \leq x \leq x_{i+1}\}.
\end{eqnarray}
\]
Remark. If \(f(x)\) is continuous on \([a,b]\), \(f(x)\) has maximum and minimum values on each \([x_i, x_{i+1}]\) (the Extreme Value Theorem). Thus, \(M_i\) and \(m_i\) are the maximum and minimum values. If \(f(x)\) is not continuous, it may not have maximum or minimum values. Nevertheless, since we are assuming \(f(x)\) is bounded, supremum and infimum do exist. Thus \(M_i\) and \(m_i\) are well-defined. □
are called the upper Riemann sum and lower Riemann sum, respectively, with respect to the partition \(\Delta\).
Remark. Upper or lower Riemann sums are also known as upper and lower Darboux sums. □
Each \(m_i\cdot(x_{i+1}-x_{i})\) (or \(M_i\cdot(x_{i+1}-x_{i})\)) represents the area of a rectangle with the ``width'' of \(x_{i+1} - x_{i}\) and the ``height'' of \(m_i\) (or \(M_i\)). Note that \(m_i\) (or \(M_i\)) can be negative. So this area is a signed area.
The Riemann sums approximate the signed area enclosed by the graph of \(y=f(x)\), the \(x\)-axis, \(x = a\) and \(x = b\). Clearly,
\[s_{\Delta} \leq S_{\Delta}.\]
Definition (Refinement of a partition)
Let \(\Delta = \{x_0, x_1, \cdots, x_n\}\) and \(\Delta' = \{x_0',x_1',\cdots, x_m'\}\) be partitions of \([a,b]\). \(\Delta'\) is said to be a refinement of \(\Delta\) if each \(a_i\in\Delta\) is equal to some \(a_j'\in \Delta'\).
Remark. If \(x_i = x_j'\) and \(x_{i+1} = x_{j+k}'\), then the sub-interval \([x_{i}, x_{i+1}]\) is partitioned into smaller sub-intervals
where \(s_{\Delta}\) and \(s_{\Delta'}\) are lower Riemann sums of \(f(x)\) with respect to \(\Delta\) and \(\Delta'\), respectively, and \(S_{\Delta}\) and \(S_{\Delta'}\) are the corresponding upper Riemann sums.
Proof. Let \(I_i = [x_{i}, x_{i+1}]\) be a sub-interval of \(\Delta\), and \(I_j' = [x_{j}', x_{j+1}']\) be a sub-interval of \(\Delta'\) such that \(I_j'\subset I_i\).
Let \(m_i = \inf\{f(x)\mid x \in I_i\}\) and \(m_j' = \inf\{f(x)\mid x \in I_j' \}\).
Then, \[m_i \leq m_j'\] because \(m_i\) is at least as low as \(m_j'\) (since \(I_i\) covers a wider range than \(I_j'\)).
If \(x_i = x_j'\) and \(x_{i+1} = x_{j+k}'\), then
Therefore the sequence of the lower Riemann sums is monotone increasing with respect to refinements of partitions, and is bounded above by upper Riemann sums (e.g., \(S_{\Delta}\), etc). By the continuity axiom of real numbers, \(\sup_{\Delta} s_{\Delta}\) exists (\(\sup_{\Delta}\) means supremum over all possible partitions \(\Delta\)). Similarly, \(\inf_{\Delta}S_{\Delta}\) exists. Clearly,
Based on the above observation, we define the definite integral as follows.
Definition (Definite integral)
Let \(f(x)\) be a bounded function on the closed interval \([a,b]\). Let \(s_{\Delta}\) and \(S_{\Delta}\) be the lower and upper Riemann sums of \(f(x)\) with respect to the partition \(\Delta\). The function \(f(x)\) is said to be Riemann-integrable, or simply integrable, if \(\sup_{\Delta}s_{\Delta} = \inf_{\Delta}S_{\Delta}\) and we write
This function is not Riemann-integrable. Take any sub-interval \([a, b]\) of any partition \(\Delta\). Since rational numbers are dense, there is always at least one rational number in \([a,b]\). Hence \(S_{\Delta} = 1\). Similarly, irrational numbers are dense so that \(s_{\Delta} = 0\). Therefore
Defining the birth process Consider a colony of bacteria that never dies. We study the following process known as the birth process , also known as the Yule process . The colony starts with \(n_0\) cells at time \(t = 0\). Assume that the probability that any individual cell divides in the time interval \((t, t + \delta t)\) is proportional to \(\delta t\) for small \(\delta t\). Further assume that each cell division is independent of others. Let \(\lambda\) be the birth rate. The probability of a cell division for a population of \(n\) cells during \(\delta t\) is \(\lambda n \delta t\). We assume that the probability that two or more births take place in the time interval \(\delta t\) is \(o(\delta t)\). That is, it can be ignored. Consequently, the probability that no cell divides during \(\delta t\) is \(1 - \lambda n \delta t - o(\delta t)\). Note that this process is an example of the Markov chain with states \({n_0}, {n_0 + 1}, {n_0 + 2}...
Generational growth Consider the following scenario (see the figure below): A single individual (cell, organism, etc.) produces \(j (= 0, 1, 2, \cdots)\) descendants with probability \(p_j\), independently of other individuals. The probability of this reproduction, \(\{p_j\}\), is known. That individual produces no further descendants after the first (if any) reproduction. These descendants each produce further descendants at the next subsequent time with the same probabilities. This process carries on, creating successive generations. Figure 1. An example of the branching process. Let \(X_n\) be the random variable representing the population size (number of individuals) of generation \(n\). In the above figure, we have \(X_0 = 1\), \(X_1=4\), \(X_2 = 7\), \(X_3=12\), \(X_4 = 9.\) We shall assume \(X_0 = 1\) as the initial condition. Ideally, our goal would be to find how the population size grows through generations, that is, to find the probability \(\Pr(X_n = k)\) for e...
The birth-death process Combining birth and death processes with birth and death rates \(\lambda\) and \(\mu\), respectively, we expect to have the following differential-difference equations for the birth-death process : \[\begin{eqnarray}\frac{{d}p_0(t)}{{d}t} &=& \mu p_1(t),\\\frac{{d}p_n(t)}{{d}t} &=& \lambda(n-1)p_{n-1}(t) - (\lambda + \mu)np_n(t) + \mu(n+1)p_{n+1}(t),~~(n \geq 1).\end{eqnarray}\] You should derive the above equations based on the following assumptions: Given a population with \(n\) individuals, the probability that an individual is born in the population during a short period \(\delta t\) is \(\lambda n \delta t + o(\delta t)\). Given a population with \(n\) individuals, the probability that an individual dies in the population is \(\mu n \delta t + o(\delta t)\). The probability that multiple individuals are born or die during \(\delta t\) is negligible. (The probability of one birth and one death during \(\delta t\) is also negligible.) Consequ...
Comments
Post a Comment