CQF Math Class Notes

Bayes’ Probability

P ( B ∣ A ) = P ( A ∣ B ) ∗ P ( B ) P ( A ) P(B \vert A) = {P(A \vert B) * P(B) \over P(A)} P(B∣A)=P(A)P(A∣B)∗P(B)​
or
P ( B ∣ A ) ∗ P ( A ) = P ( A ∩ B ) = P ( B ∩ A ) = P ( A ∣ B ) ∗ P ( B ) P(B \vert A) * P(A) = P(A \cap B) = P(B \cap A) = P(A \vert B) * P(B) P(B∣A)∗P(A)=P(A∩B)=P(B∩A)=P(A∣B)∗P(B)
我们的假设条件是A事件基于B的发生的概率与B事件基于A发生的概率是一样的。
举个例子:从一副牌连续翻出两张K的概率是多少 (大小鬼被剔除)?
K ∩ K = 4 / 52 ∗ 3 / 51 = 1 / 221 ≈ 0.5 % K \cap K = 4/52 * 3/51 = 1/221 \approx 0.5 \% K∩K=4/52∗3/51=1/221≈0.5%
我们从上面的例子可以看出,两个条件互换所得出的结果是一样的。从而得出上面贝叶斯概率公式的结论。


The Crazy Intergral and its Rules and Formulas

我们如何来证明这道题?可以说我们基本上把能用的微积分知识和定理都撸了一遍,这是一道非常好的一个练习题:
∫ 1 ( x 2 + 1 ) d x = t a n − 1 x + C \int {1 \over(x^2+1)} dx = tan^{-1}x + C ∫(x2+1)1​dx=tan−1x+C
= ∫ 1 ( x + i ) ( x − i ) d x = ∫ ( 1 / 2 i ) ( x + i ) d x + ∫ ( − i / 2 i ) ( x − i ) d x = \int {1 \over {(x + i)(x - i)}} dx = \int {(1/2i) \over (x+i)}dx + \int {(-i/2i) \over (x - i)}dx =∫(x+i)(x−i)1​dx=∫(x+i)(1/2i)​dx+∫(x−i)(−i/2i)​dx
= 1 2 i ( l n ( x − i ) − l n ( x + i ) ) = {1 \over 2i} \bigl(ln(x-i)-ln(x+i)\bigr) =2i1​(ln(x−i)−ln(x+i))
这个答案看着好像已经结束了,当时这是不可以的也不应该!因为我们原问题并不包含虚数,所以答案里包含虚数是不合理的,所以要继续简化!
下面我们要用到欧拉恒等式(Euler’s Equation) e i θ = cos ⁡ θ + i sin ⁡ θ e^{i\theta} = \cos \theta + i \sin \theta eiθ=cosθ+isinθ来进行运算:
a + b i = r e i θ a + bi = re^{i \theta} a+bi=reiθ
r = ( a 2 + b 2 ) r = \sqrt{(a^2+b^2)} r=(a2+b2) ​
θ = t a n − 1 b a \theta = tan^{-1}{b \over a} θ=tan−1ab​
然后我们继续进行替换:
= 1 2 i ( l n ( x − i ) − l n ( x + i ) ) = 1 2 i ( l n ( ( x 2 + ( − 1 ) 2 ) ∗ e i t a n − 1 − 1 x ) − l n ( ( x 2 + ( 1 ) 2 ) ∗ e i t a n − 1 1 x ) ) = {1 \over 2i} \bigl(ln(x-i)-ln(x+i)\bigr) = {1 \over 2i} \Bigl(ln(\sqrt{(x^2+(-1)^2)}*e^{itan^{-1}{-1 \over x}})-ln(\sqrt{(x^2+(1)^2)}*e^{itan^{-1}{1 \over x}})\Bigr) =2i1​(ln(x−i)−ln(x+i))=2i1​(ln((x2+(−1)2) ​∗eitan−1x−1​)−ln((x2+(1)2) ​∗eitan−1x1​))
= 1 2 i ( l n ( x 2 + 1 ) + l n ( e i t a n − 1 − 1 x ) − l n ( x 2 + 1 ) − l n ( e i t a n − 1 1 x ) ) = {1 \over 2i} \Bigl(ln{\sqrt{(x^2+1)}} + ln({e^{itan^{-1}{-1 \over x}}}) - ln{\sqrt{(x^2+1)}} - ln({e^{itan^{-1}{1 \over x}})}\Bigr) =2i1​(ln(x2+1) ​+ln(eitan−1x−1​)−ln(x2+1) ​−ln(eitan−1x1​))
l n ( x 2 + 1 ) 抵消掉 ln{\sqrt{(x^2+1)}}抵消掉 ln(x2+1) ​抵消掉, l n ( e ) 抵消掉 ln{(e)}抵消掉 ln(e)抵消掉,
t a n − 1 tan^{-1} tan−1是一个奇函数(odd function),所以我们可以将(-1/x)的提出负号出来,即:
= 1 2 i ( − 2 i ∗ t a n − 1 1 x ) = − t a n − 1 1 x + C 1 = {1 \over 2i}{\bigl(-2i * tan^{-1}{1 \over x}}\bigr) = -tan^{-1}{1 \over x} + C1 =2i1​(−2i∗tan−1x1​)=−tan−1x1​+C1
其实做到这里我们可以算是做完了,我个人觉得老师未必会扣你的分数。但是呢,其实还是可以有一丢丢的优化空间,因为:
t a n − 1 x + c o t − 1 x = π 2 tan^{-1}x + cot^{-1}x = {\pi \over 2} tan−1x+cot−1x=2π​,且 cot ⁡ − 1 x = tan ⁡ − 1 1 x \cot^{-1}x = \tan^{-1}{1 \over x} cot−1x=tan−1x1​,左右移动项得出:
= tan ⁡ − 1 x − π 2 + C 1 = \tan^{-1}x - {\pi \over 2} + C1 =tan−1x−2π​+C1
π 2 + C 1 {\pi \over 2} + C1 2π​+C1他们其实都可以看作是常数项,所以可以合并,最终最优解,即:
∫ 1 ( x 2 + 1 ) d x = t a n − 1 x + C \int {1 \over(x^2+1)} dx = tan^{-1}x + C ∫(x2+1)1​dx=tan−1x+C


Trignometry Formula

cos ⁡ ( α + β ) = cos ⁡ α ∗ cos ⁡ β − sin ⁡ α ∗ sin ⁡ β \cos(\alpha+\beta) = \cos\alpha*\cos\beta - \sin\alpha*\sin\beta cos(α+β)=cosα∗cosβ−sinα∗sinβ
cos ⁡ ( α − β ) = cos ⁡ α ∗ cos ⁡ β + sin ⁡ α ∗ sin ⁡ β \cos(\alpha-\beta) = \cos\alpha*\cos\beta + \sin\alpha*\sin\beta cos(α−β)=cosα∗cosβ+sinα∗sinβ
sin ⁡ ( α ± β ) = sin ⁡ α ∗ cos ⁡ β ± cos ⁡ α ∗ sin ⁡ β \sin(\alpha\pm\beta) = \sin\alpha*\cos\beta \pm \cos\alpha*\sin\beta sin(α±β)=sinα∗cosβ±cosα∗sinβ

tan ⁡ ( α + β ) = tan ⁡ α + tan ⁡ β 1 − tan ⁡ α ∗ tan ⁡ β \tan(\alpha+\beta) = {{\tan\alpha+\tan\beta} \over {1-\tan\alpha*\tan\beta}} tan(α+β)=1−tanα∗tanβtanα+tanβ​

tan ⁡ ( α − β ) = tan ⁡ α − tan ⁡ β 1 − tan ⁡ α ∗ tan ⁡ β \tan(\alpha-\beta) = {{\tan\alpha-\tan\beta} \over {1-\tan\alpha*\tan\beta}} tan(α−β)=1−tanα∗tanβtanα−tanβ​


Forward Kolmogorov Equation

The equation is given by

δ p δ t ′ = C 2 δ 2 p δ y ′ 2 {\delta p \over \delta t'} = C^2{\delta^2p \over \delta y'^2} δt′δp​=C2δy′2δ2p​

for the transition density function p ( y , t ; y ′ , t ′ ) ; c 2 ∈ R p(y,t;y',t');c^2\in \mathbb{R} p(y,t;y′,t′);c2∈R, The states (y,t) are past and are fixed while (y’,t’) refers t ofuture ones and are variables. By simple substitution show that

p ( y , t ; y ′ , t ′ ) = 1 2 c π ( t ′ − t ) e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) p(y,t;y',t') = {1 \over 2c\sqrt{\pi(t'-t)}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) p(y,t;y′,t′)=2cπ(t′−t) ​1​exp(−4c2(t′−t)(y′−y)2​) ---- (1.2) satisfies the KFE

show that 1.2.satisfies
∫ − ∞ ∞ p ( y , t ; y ′ , t ′ ) d y ′ = 1 \int_{-\infty}^\infty p(y,t;y',t')dy' = 1 ∫−∞∞​p(y,t;y′,t′)dy′=1

let’s solve this problem.

LHS:

δ p δ t ′ = 1 2 c π ∗ e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) ∗ − 1 2 ( t ′ − t ) 3 2 + 1 2 c π ( t ′ − t ) e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) ∗ ( y ′ − y ) 2 4 c 2 ( t ′ − t ) 2 {\delta p \over \delta t'} = {1 \over 2c\sqrt{\pi}} * exp(-{(y'-y)^2 \over 4c^2(t'-t)}) * -{{1\over2}(t'-t)^{3\over2}} + {1 \over 2c\sqrt{\pi(t'-t)}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) * {(y'-y)^2 \over 4c^2(t'-t)^2} δt′δp​=2cπ ​1​∗exp(−4c2(t′−t)(y′−y)2​)∗−21​(t′−t)23​+2cπ(t′−t) ​1​exp(−4c2(t′−t)(y′−y)2​)∗4c2(t′−t)2(y′−y)2​ ---- product rule

= ( y ′ − y ) 2 − 2 c 2 ( t ′ − t ) 8 c 3 π ( t ′ − t ) 5 2 e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) = {(y'-y)^2-2c^2(t'-t) \over 8c^3\sqrt{\pi}(t'-t)^{5\over2}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) =8c3π ​(t′−t)25​(y′−y)2−2c2(t′−t)​exp(−4c2(t′−t)(y′−y)2​)

RHS:

δ p δ y ′ = 1 2 c π ( t ′ − t ) ∗ e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) ∗ − 2 ( y ′ − y ) 4 c 2 ( t ′ − t ) {\delta p \over \delta y'} = {1 \over 2c\sqrt{\pi(t'-t)}} * exp(-{(y'-y)^2 \over 4c^2(t'-t)}) * {-2(y'-y) \over 4c^2(t'-t)} δy′δp​=2cπ(t′−t) ​1​∗exp(−4c2(t′−t)(y′−y)2​)∗4c2(t′−t)−2(y′−y)​

= − ( y ′ − y ) ( 4 c 3 π ( t ′ − t ) 3 2 ∗ e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) =-{(y'-y)\over(4c^3{\sqrt\pi}(t'-t)^{3\over2}} * exp(-{(y'-y)^2 \over 4c^2(t'-t)}) =−(4c3π ​(t′−t)23​(y′−y)​∗exp(−4c2(t′−t)(y′−y)2​)

δ 2 p δ y ′ 2 = − 1 ( 4 c 3 π ( t ′ − t ) 3 2 ∗ e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) − ( y ′ − y ) ( 4 c 3 π ( t ′ − t ) 3 2 ∗ e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) ∗ − 2 ( y ′ − y ) 4 c 2 ( t ′ − t ) {\delta^2 p \over \delta y'^2}=-{1\over(4c^3{\sqrt\pi}(t'-t)^{3\over2}} * exp(-{(y'-y)^2 \over 4c^2(t'-t)}) - {(y'-y)\over(4c^3{\sqrt\pi}(t'-t)^{3\over2}} * exp(-{(y'-y)^2 \over 4c^2(t'-t)}) * {-2(y'-y) \over 4c^2(t'-t)} δy′2δ2p​=−(4c3π ​(t′−t)23​1​∗exp(−4c2(t′−t)(y′−y)2​)−(4c3π ​(t′−t)23​(y′−y)​∗exp(−4c2(t′−t)(y′−y)2​)∗4c2(t′−t)−2(y′−y)​

= ( y ′ − y ) 2 − 2 c 2 ( t ′ − t ) ( 8 c 5 π ( t ′ − t ) 5 2 ∗ e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) ={(y'-y)^2-{2c^2(t'-t)}\over(8c^5{\sqrt\pi}(t'-t)^{5\over2}} * exp(-{(y'-y)^2 \over 4c^2(t'-t)}) =(8c5π ​(t′−t)25​(y′−y)2−2c2(t′−t)​∗exp(−4c2(t′−t)(y′−y)2​)

Hence: L H S = δ p δ t ′ = ( y ′ − y ) 2 − 2 c 2 ( t ′ − t ) 8 c 3 π ( t ′ − t ) 5 2 e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) = C 2 δ 2 p δ y ′ 2 = R H S LHS = {\delta p \over \delta t'} = {(y'-y)^2-2c^2(t'-t) \over 8c^3\sqrt{\pi}(t'-t)^{5\over2}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) = C^2{\delta^2p \over \delta y'^2} = RHS LHS=δt′δp​=8c3π ​(t′−t)25​(y′−y)2−2c2(t′−t)​exp(−4c2(t′−t)(y′−y)2​)=C2δy′2δ2p​=RHS

we know that ( y ′ − y ) 2 − 2 c 2 ( t ′ − t ) 8 c 3 π ( t ′ − t ) 5 2 e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) {(y'-y)^2-2c^2(t'-t) \over 8c^3\sqrt{\pi}(t'-t)^{5\over2}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) 8c3π ​(t′−t)25​(y′−y)2−2c2(t′−t)​exp(−4c2(t′−t)(y′−y)2​) is a common solution to FKE, then let’s see if its PDF is equal to 1 as all PDF should sum up to 1.

∫ − ∞ ∞ p ( y , t ; y ′ , t ′ ) d y ′ = ∫ − ∞ ∞ 1 2 c π ( t ′ − t ) e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) d y ′ \int_{-\infty}^\infty p(y,t;y',t')dy' = \int_{-\infty}^\infty {1 \over 2c\sqrt{\pi(t'-t)}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) dy' ∫−∞∞​p(y,t;y′,t′)dy′=∫−∞∞​2cπ(t′−t) ​1​exp(−4c2(t′−t)(y′−y)2​)dy′

from observation, it is better to use chain rule in this situation, and it’ll make things a lot easier.

let u = ( y ′ − y ) 2 c ( t ′ − t ) u ={(y'-y) \over 2c{\sqrt(t'-t)}} u=2c( ​t′−t)(y′−y)​ then d u d y ′ = 1 2 c ( t ′ − t ) {du \over dy'} = {1 \over {2c {\sqrt{(t'-t)}}}} dy′du​=2c(t′−t) ​1​, thus d y ′ = 2 c ( t ′ − t ) d u dy' = {2c {\sqrt{(t'-t)}}}du dy′=2c(t′−t) ​du

∫ − ∞ ∞ 1 2 c π ( t ′ − t ) e x p ( − ( y ′ − y ) 2 4 c 2 ( t ′ − t ) ) d y ′ ⇒ 1 2 c π ( t ′ − t ) ∗ ∫ − ∞ ∞ e − u 2 d u ∗ 2 c ( t ′ − t ) = 1 π ∗ ∫ − ∞ ∞ e − u 2 d u = 1 \int_{-\infty}^\infty {1 \over 2c\sqrt{\pi(t'-t)}}exp(-{(y'-y)^2 \over 4c^2(t'-t)}) dy' \Rightarrow {1 \over 2c\sqrt{\pi(t'-t)}} * {\int_{-\infty}^\infty e^{-u^2}}du * {2c {\sqrt{(t'-t)}}} = {1 \over \sqrt\pi} * {\int_{-\infty}^\infty e^{-u^2}}du = 1 ∫−∞∞​2cπ(t′−t) ​1​exp(−4c2(t′−t)(y′−y)2​)dy′⇒2cπ(t′−t) ​1​∗∫−∞∞​e−u2du∗2c(t′−t) ​=π ​1​∗∫−∞∞​e−u2du=1

why? because ∫ − ∞ ∞ e − u 2 d u = π {\int_{-\infty}^\infty e^{-u^2}}du = \sqrt\pi ∫−∞∞​e−u2du=π ​ is known by definition. and we are done with this quesiton.

— 待更 —

【CQF Math Class 数学笔记】相关推荐

  1. Kill Math: 让数学不只是符号

    Kill Math: 让数学不只是符号 Author: Bret Victor / April 11, 2011 译者:神奇的战士 Bret Victor 是苹果公司的前 IPAD 的交互设计师, 大 ...

  2. latex 上标_1700页数学笔记火了!全程敲代码,速度飞快易搜索,硬核小哥教你上手LaTeX+Vim...

    晓查 乾明 发自 凹非寺  转自: 量子位 | 公众号 QbitAI 又出现一位"神仙"本科生! 数学课上,全程键盘手打1700页笔记. 速度紧追老师板书,公式.图形一个不落. 效 ...

  3. 1700 页数学笔记火了!全程敲代码,速度飞快易搜索,硬核小哥教你上手 LaTeX+Vim...

    点击"小詹学Python",选择"星标"公众号 重磅干货,第一时间送达 本文转载自量子位,禁二次转载 又出现一位"神仙"本科生! 数学课上, ...

  4. LaTeX+Vim 数学笔记

    又出现一位"神仙"本科生! 数学课上,全程键盘手打1700页笔记. 速度紧追老师板书,公式.图形一个不落. 效果?请看下图: 不仅排版媲美教科书,而且还能够批注,检索关键词-- 笔 ...

  5. 读书笔记:《思考的乐趣:Matrix67数学笔记》第4章 统计数据的陷阱

    <思考的乐趣:Matrix67数学笔记>第4章讲了几个统计学上的陷阱,由于现在流行的大数据与统计学很有渊源,所以认真读了这一章,在<大数据时代>中指出只考虑相关性就够了,而不考 ...

  6. 数学笔记3——导数3(隐函数的导数)

    数学笔记3--导数3(隐函数的导数) 幂函数的扩展形式 f(x) = xn的导数:f'(x) = nxn-1,n是整数,该公式对f(x) = xm/n, m,n 是整数同样适用. 推导过程: 什么是隐 ...

  7. 程序员的数学笔记3--迭代法

    第三节课程,介绍的是迭代法. 前两节笔记的文章: 程序员的数学笔记1–进制转换 程序员的数学笔记2–余数 03 迭代法 什么是迭代法 迭代法,简单来说,其实就是不断地用旧的变量值,递推计算新的变量值. ...

  8. 程序员的数学笔记2--余数

    上一节程序员的数学笔记1–进制转换是介绍了进制,特别是十进制和二进制之间的转换,移位操作和逻辑操作. 今天介绍的是余数,看完本节笔记,你会发现生活中有很多东西都有余数的影子. 余数 余数的特性 整数是 ...

  9. 清华美女学霸数学笔记曝光, 精美程度无与伦比

    数学学习是从基础开始的,所以在自学数学基础知识有了一定的基础和能力后,我们可以对自己提出新的更高的要求:写一点数学读书笔记,目的是培养自己的创造性思维.这样,可以逼着自己主动去发现问题.思考问题.解决 ...

最新文章

  1. 旋转矩阵、欧拉角、轴角
  2. flask使用debug模式时,存在错误时,会占用设备内存直至服务重启才释放;debug模式会开启一个守护进程(daemon process)...
  3. access 一亿条数据_循环运算数据溢出
  4. 云存储领域的那些安全难题与解决途径
  5. 【小白学PyTorch】16.TF2读取图片的方法
  6. autopep8规范你的python代码
  7. How to change the status of Prepayment invoice
  8. glDrawElements参数在新旧版本传最后一个参数的不同
  9. Java基础学习总结(72)——提升 java 代码的运行效率
  10. 浅议.NET、ASP.NET和C#的关系
  11. 【超清视频】SCCM2012精讲系列课程03:SCCM2012的基础环境初始化(一)
  12. POJ 1625 Censored!(AC自动机 + DP + 大数 + 拓展ASCII处理)题解
  13. kotlin写的音视频app(已开源)
  14. 科普:什么是人工智能
  15. 怎么解决计算机键盘驱动,键盘驱动不兼容,详细教您键盘驱动不兼容怎么解决...
  16. linux 下numpy安装包,Linux下安装numpy
  17. android开发教程!看懂这些帮你轻松解决就业问题!知乎上转疯了!
  18. pxe自动装机利用tfp,http,nfs服务实现。
  19. Java实现的一个简易网络画板
  20. 啤酒每罐2.3元,饮料每罐1.9元

热门文章

  1. 每周分享第 31 期
  2. 对于多目标跟踪论模型CTracker的解析
  3. Android开发 Studio4.0 APP logo 适配
  4. 什么叫死区时间_关于pwm死区时间的介绍
  5. Billy Belceb病毒编写教程---Win32篇
  6. DICOM医学影像协议
  7. About kettle carte
  8. resultful 风格
  9. 文华学院大学计算机基础考试数据库ip,文华学院大学计算机基础模拟试卷
  10. sql中on和where的区别在哪