14 Apr

费曼积分法(8):求高斯积分

自从了解了费曼积分法之后,我就一直想着用费曼积分法来求高斯积分$\int_0^{\infty} e^{-x^2}dx=\frac{\sqrt{\pi}}{2}$这个神奇的积分,但一直无果。在《数学桥》里边,作者是通过将其转变为二重积分来解决的,简洁而巧妙。但是为了显示费曼积分法的威力,我一直想找到高斯积分的其他求法。上星期在《数学物理方法》中看到作者用拉普拉斯变换求出了该积分,眼睛不禁为之一亮,不过这属于积分变换内容,属于“积分符号内取积分”的技巧,在此不作讨论。今天在网上查找资料时,在“赵洁”的一篇论文《含参变量积分》中,看到了一种属于费曼积分法范畴内的方法,特与大家分享。

从“事后分析”来看,高斯积分的结果涉及到了$\sqrt{\pi}$这个量,一般来说我们常见的公式出现$\pi$的不少,可是几乎没有出现$\sqrt{\pi}$的,所以一般来说我们都将它平方。我们引入
$$f(x)=(\int_0^x e^{-t^2}dt)^2$$

点击阅读全文...

9 Apr

2^29365247-1也不是素数!

第二个数字也测试完了,$2^{29365247}-1$也不是一个素数。继续努力!

[Comm thread Apr 9 12:22] Sending result to server: UID: bojone/bojone, M29365247 is not prime. Res64: AFA532C54A91F89B. We4: D8A82DF8,3429750,00000000, AID: C
[Comm thread Apr 9 12:22]
[Comm thread Apr 9 12:22] PrimeNet success code with additional info:
[Comm thread Apr 9 12:22] LL test successfully completes double-check of M29365247
[Comm thread Apr 9 12:22] CPU credit is 29.2021 GHz-days.
[Comm thread Apr 9 12:22] Done communicating with server.
8 Apr

2^29363731-1不是素数!

2^29363731-1

2^29363731-1

很小的时候就开始对素数感兴趣了,后来是在一本《未解之谜》上看到了梅森素数、完全数、孪生素数等等东西,觉得甚是好玩。在初中买了计算机之后,就关注到了Prime 95这个梅森素数的分布式计算程序,以前也尝试过运行它,不过由于那时候计算机配置较低,一般都是运行到20%左右就没有坚持下去了。

上大学入手了一台四核的笔记本,就在去年10月份左右再次运行了这个程序,由于是四核,一次性可以同时测试四个数字。经过半年的运行,今天终于测试完了第一个数字:$2^{29363731}-1$。正如预料中的,这不是一个素数。不管怎样,它是我第一个完成的测试,也算是自己的一个独立的成果啦,呵呵,自娱自乐一番。

点击阅读全文...

27 Mar

费曼积分法(7):欧拉数学的综合

在本系列的第五篇文章中,BoJone导出了一些看似不合理的公式,而且并没有说明它的应用和来源。其实,这些都是我在研究以下积分的时候总结出来的:

$$\int_{-\infty}^{+\infty} \frac{\cos x}{a^2+x^2}dx$$

点击阅读全文...

24 Mar

费曼积分法(6):教科书上的两道练习题

我们的《数学分析》教程上有两道比较有趣的定积分,经测试可以用费曼积分法的思路解决。

$$\begin{aligned}\int_0^1 \frac{\ln(1+x)}{1+x^2}dx \\ \int_0^{\pi} \frac{x \sin x}{1+\cos^2 x}dx\end{aligned}$$

No.1

点击阅读全文...

24 Mar

费曼积分法(5):欧拉数学的传承

在大学第二学期,我们的《数学分析》终于龟速地爬行到了定积分这一章节。对于一些比较复杂的定积分,我总想用自己的方法来解决它,这就重新燃起了我对“费曼积分法——积分符号内取微分”的热情。尤其是我用费曼积分法解决了几道比较有趣复杂的定积分问题时,成就感高涨,遂在此总结,与大家共勉。

这和欧拉数学有什么关系呢?之前已经提到过,欧拉数学是用一种不严谨却极具创造性的方式,给予我们对数学的介乎感性和理性的直观理解。我觉得费曼积分法也属于这个范畴内,它着眼于用一种特殊的视角解决问题,而暂时忽略掉数学严密性。在读费曼的故事中,我感觉到这种思想是贯穿他一生的研究之中的。

本文继续对费曼积分法的研究,得出一些不是很严谨的结论,为以后的应用奠下基础。

一、不成立的函数

首先我们重新考虑$\int_0^{\infty} \frac{\sin x}{x}dx$。这一次我们将它引入复数范畴内,考虑:
$$\int_0^{\infty}\frac{\cos x+i \sin x}{x}dx=\int_0^{\infty}\frac{e^{ix}}{x}dx$$

点击阅读全文...

7 Mar

轻微的扰动——摄动法简介(3)

微分方程领域大放光彩

虽然微分方程在各个计算领域都能一展才华,不过它最辉煌的光芒无疑绽放于微分方程领域,包括常微分方程和偏微分方程。海王星——“笔尖上发现的行星”——就是摄动法的著名成果,类似的还有冥王星的发现。天体力学家用一颗假设的行星的引力摄动来解释已知行星的异常运动,并由此反推未知行星的轨道。我们已不止一次提到过,一般的三体问题是混沌的,没有精确的解析解。这就要求我们考虑一些近似的方法,这样的方法发展起来就成为了摄动理论。

跟解代数方程一样,摄动法解带有小参数或者大参数的微分方程的基本思想,就是将微分方程的解表达为小参数或大参数的幂级数。当然,这是最直接的,也相当好理解,不过所求得的级数解有可能存在一些性态不好的情况,比如有时原解应该是一个周期运动,但是级数解却出现了诸如$t \sin t$的“长期项”,这是相当不利的,因此也发展出各种技巧来消除这些项。可见,摄动理论是一门应用广泛、集众家所大成的实用理论。下面我们将通过一些实际的例子来阐述这个技巧。

点击阅读全文...

21 Feb

[问题解答]有多少位数字?

解决完上一题《有多少个5?》后,子瑞表示看到一道类似的题目,当然,这道题比上一道难一些:

一个数,各个数字加起来等于900,乘以2后各个数字加起来还是等于900,已知这个数字只有3、4、5、6组成,请问满足条件的最大数与最小数的积有多少位数?

要解答这个问题,我们只需要知道最大数和最小数分别有多少位即可。因为最大数必然是6...3的形式,而最小数只能是3...6的形式,它们的位数之和就是所求的位数。

怎样比较两个数的大小呢?显然,在不同位数的数时,位数多的数要大,同样位数才从高到低逐位比较。因此,我们应当考虑位数的最大与最小。

点击阅读全文...