练习题-18 计算两个积分

2024-06-19 02:20

本文主要是介绍练习题-18 计算两个积分,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

  1. 计算积分 I = ∫ R e − t 4 d t . I=\int_{\mathbb{R}} e^{-t^4} dt. I=Ret4dt.

:令 x = t 4 x=t^4 x=t4. 则 I = 2 ∫ 0 ∞ e − x ⋅ 1 4 ⋅ x − 3 / 4 d x = 1 2 Γ ( 1 4 ) I = 2\int_0^\infty e^{-x} \cdot \frac{1}{4}\cdot x^{-3/4} dx=\frac{1}{2} \Gamma(\frac{1}{4}) I=20ex41x3/4dx=21Γ(41).

:已知 Γ ( 1 4 ) \Gamma(\frac{1}{4}) Γ(41)是超越数.

  1. 计算积分 J = ∫ 0 ∞ 1 1 + t 2 sin ⁡ t t d t . J=\int_0^{\infty} \frac{1}{1+t^2} \frac{\sin t}{t}dt. J=01+t21tsintdt.

J = 1 2 ∫ R 1 1 + t 2 sin ⁡ t t d t J=\frac{1}{2} \int_{\mathbb{R}} \frac{1}{1+t^2} \frac{\sin t}{t}dt J=21R1+t21tsintdt. 用围道积分与留数来计算:

J = I m [ ∫ C R + ∫ C 1 + ∫ C 2 + ∫ C r e i z 2 z ( z 2 + 1 ) d t ] J=\mathrm{Im} \left[\int_{C_R} + \int_{C_1} + \int_{C_2} + \int_{C_r} \frac{e^{iz}}{2z(z^2+1)} dt \right] J=Im[CR+C1+C2+Cr2z(z2+1)eizdt].
其中 R > r > 0 R>r>0 R>r>0, C R C_R CR是上半圆弧 ∣ z ∣ = R |z|=R z=R, C r C_r Cr是上半圆弧 ∣ z ∣ = r |z|=r z=r, C 1 C_1 C1沿着实轴从 − R -R R − r -r r, C 2 C_2 C2沿着实轴从 r r r R R R. 四段曲线构成一个闭合围道 γ \gamma γ,正方向取逆时针方向。

在上半平面(或者在 γ \gamma γ所围成的区域内), 被积函数只有一个极点 z = i z=i z=i. 按照Jordan引理,当 R → ∞ R\to \infty R时,积分趋向于 0 0 0; 在 C 1 , C 2 C_1, C_2 C1,C2上的积分会抵消,当 r → 0 r \to 0 r0时,极点 z = 0 z=0 z=0对积分的“贡献”只有一半。根据柯西积分定理, J J J等于下列“留数贡献”的虚部:
2 π i R e s ( e i z 2 z ( z 2 + 1 ) , z = i ) + π i R e s ( e i z 2 z ( z 2 + 1 ) , z = 0 ) = 2 π i lim ⁡ z → i ( z − i ) e i z 2 z ( z 2 + 1 ) + π i lim ⁡ z → 0 z e i z 2 z ( z 2 + 1 ) = π i 2 ( 1 − 1 e ) . \begin{aligned} &2\pi i\, \mathrm{Res}\left(\frac{e^{iz}}{2z(z^2+1)} , z=i \right)+\pi i \,\mathrm{Res} \left(\frac{e^{iz}}{2z(z^2+1)} , z=0\right) \\ =& 2\pi i \lim_{z \to i} \frac{(z-i)e^{iz}}{2z(z^2+1)} + \pi i \lim_{z\to 0} \frac{ze^{iz}}{2z(z^2+1)}\\ =& \frac{\pi i}{2}(1-\frac{1}{e}). \end{aligned} ==2πiRes(2z(z2+1)eiz,z=i)+πiRes(2z(z2+1)eiz,z=0)2πizilim2z(z2+1)(zi)eiz+πiz0lim2z(z2+1)zeiz2πi(1e1).

所以 J = π 2 ( 1 − 1 e ) . J=\frac{\pi }{2}(1-\frac{1}{e}). J=2π(1e1).

这篇关于练习题-18 计算两个积分的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1073735

相关文章

计算绕原点旋转某角度后的点的坐标

问题: A点(x, y)按顺时针旋转 theta 角度后点的坐标为A1点(x1,y1)  ,求x1 y1坐标用(x,y)和 theta 来表示 方法一: 设 OA 向量和x轴的角度为 alpha , 那么顺时针转过 theta后 ,OA1 向量和x轴的角度为 (alpha - theta) 。 使用圆的参数方程来表示点坐标。A的坐标可以表示为: \[\left\{ {\begin{ar

【云计算 复习】第1节 云计算概述和 GFS + chunk

一、云计算概述 1.云计算的商业模式 (1)软件即服务(SaaS) 有些景区给游客提供烧烤场地,游客需要自己挖坑或者砌烧烤台,然后买肉、串串、烧烤。 (2)平台即服务(PaaS) 有些景区给游客提供烧烤场地,同时搭建好烧烤台,游客只需要自己带食材和调料、串串、烧烤。 (3)基础设施即服务(IaaS) 有些景区给游客提供烧烤场地,同时搭建好烧烤台,还有专门的厨师来烧烤,用户不需要关心前面的所有

剑指offer(C++)--和为S的两个数字

题目 输入一个递增排序的数组和一个数字S,在数组中查找两个数,使得他们的和正好是S,如果有多对数字的和等于S,输出两个数的乘积最小的。 class Solution {public:vector<int> FindNumbersWithSum(vector<int> array,int sum) {vector<int> result;int len = array.size();if(

剑指offer(C++)--两个链表的第一个公共结点

题目 输入两个链表,找出它们的第一个公共结点。 解法一 两个链表一定有交点的话,方法是指向短链表指针先走完,然后指向长链表,指向长链表指针后走完,指向短链表。所以,第二次走过,一定会在交点相遇。 class Solution {public:ListNode* FindFirstCommonNode( ListNode *pHead1, ListNode *pHead2) {ListN

两个基因相关性CPTAC蛋白组数据

目录 蛋白数据下载 ①蛋白数据下载 1,TCGA-选择泛癌数据  2,TCGA-TCPA 3,CPTAC(非TCGA) ②蛋白相关性分析 1,数据整理 2,蛋白相关性分析 PCAS在线分析 蛋白数据下载 CPTAC蛋白组学数据库介绍及数据下载分析 – 王进的个人网站 (jingege.wang) ①蛋白数据下载 可以下载泛癌蛋白数据:UCSC Xena (xena

剑指offer—编程题7(用两个栈实现一个队列)

题目:用两个栈实现一个队列。队列的声明如下,请实现它的两个函数appendTail 和deleteHead,分别完成在队列尾部插入结点和在队列头部删除结点的功能。 代码如下: [java]  view plain copy print ? public class Test07 {       /**       * 用两个栈模拟的队列       *

算法9—两个巨大正整数相加

两个巨大整数相加,可能会造成溢出,或者它的大小已经超出基本数据类型的范围,所以,我们对巨大整数进行相加时,可以把它们转换成字符串,然后通过字符串的处理进行整数相加。 这里有两种做法:第一种,把整数存在一个字符数组里进行处理。代码如下: [java]  view plain copy public static String addThroughArray(String

算法8—不通过比较,找出两个数的最大值

问题: 比如:给定两个值 5和10,不通过比较,直接找出最大值。 分析: 一旦涉及到不用比较找最大值,想都不用想,一般只能通过位运算来实现。  max = a - ((a-b)&((a-b)>>31)) 或者 max = ((a+b)+|a-b|)/2 如果找最小值,我们只需把两个值相加,减去max即可。

算法6— 判断两个链表是否相交

问题: 给出两个单向链表的头指针,比如h1、h2,判断链表是否相交,如果不相交返回NULL;如果相交,返回指向第一个相交节点的指针。时间复杂度控制在O(n)。 分析: 如果两单向链表相交的话,一定是Y型相交,不可能出现X型,弄清楚这点后接下来的工作就是: (1)先找到h1,h2的最后一个节点L1和L2,同时记录节点数量a,b;(这里假设 a > b) (2)判断最后一个节点是否相同

SpringBoot中如何监听两个不同源的RabbitMQ消息队列

spring-boot如何配置监听两个不同的RabbitMQ 由于前段时间在公司开发过程中碰到了一个问题,需要同时监听两个不同的rabbitMq,但是之前没有同时监听两个RabbitMq的情况,因此在同事的帮助下,成功实现了监听多个MQ。下面我给大家一步一步讲解下,也为自己做个笔记; 详细步骤: 1. application.properties 文件配置: u.rabbitmq.ad