【CQF Math Class 数学笔记】
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)1dx=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)1dx=∫(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)1dx=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) 1exp(−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) 1exp(−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)231∗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) 1exp(−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) 1exp(−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 数学笔记】相关推荐
- Kill Math: 让数学不只是符号
Kill Math: 让数学不只是符号 Author: Bret Victor / April 11, 2011 译者:神奇的战士 Bret Victor 是苹果公司的前 IPAD 的交互设计师, 大 ...
- latex 上标_1700页数学笔记火了!全程敲代码,速度飞快易搜索,硬核小哥教你上手LaTeX+Vim...
晓查 乾明 发自 凹非寺 转自: 量子位 | 公众号 QbitAI 又出现一位"神仙"本科生! 数学课上,全程键盘手打1700页笔记. 速度紧追老师板书,公式.图形一个不落. 效 ...
- 1700 页数学笔记火了!全程敲代码,速度飞快易搜索,硬核小哥教你上手 LaTeX+Vim...
点击"小詹学Python",选择"星标"公众号 重磅干货,第一时间送达 本文转载自量子位,禁二次转载 又出现一位"神仙"本科生! 数学课上, ...
- LaTeX+Vim 数学笔记
又出现一位"神仙"本科生! 数学课上,全程键盘手打1700页笔记. 速度紧追老师板书,公式.图形一个不落. 效果?请看下图: 不仅排版媲美教科书,而且还能够批注,检索关键词-- 笔 ...
- 读书笔记:《思考的乐趣:Matrix67数学笔记》第4章 统计数据的陷阱
<思考的乐趣:Matrix67数学笔记>第4章讲了几个统计学上的陷阱,由于现在流行的大数据与统计学很有渊源,所以认真读了这一章,在<大数据时代>中指出只考虑相关性就够了,而不考 ...
- 数学笔记3——导数3(隐函数的导数)
数学笔记3--导数3(隐函数的导数) 幂函数的扩展形式 f(x) = xn的导数:f'(x) = nxn-1,n是整数,该公式对f(x) = xm/n, m,n 是整数同样适用. 推导过程: 什么是隐 ...
- 程序员的数学笔记3--迭代法
第三节课程,介绍的是迭代法. 前两节笔记的文章: 程序员的数学笔记1–进制转换 程序员的数学笔记2–余数 03 迭代法 什么是迭代法 迭代法,简单来说,其实就是不断地用旧的变量值,递推计算新的变量值. ...
- 程序员的数学笔记2--余数
上一节程序员的数学笔记1–进制转换是介绍了进制,特别是十进制和二进制之间的转换,移位操作和逻辑操作. 今天介绍的是余数,看完本节笔记,你会发现生活中有很多东西都有余数的影子. 余数 余数的特性 整数是 ...
- 清华美女学霸数学笔记曝光, 精美程度无与伦比
数学学习是从基础开始的,所以在自学数学基础知识有了一定的基础和能力后,我们可以对自己提出新的更高的要求:写一点数学读书笔记,目的是培养自己的创造性思维.这样,可以逼着自己主动去发现问题.思考问题.解决 ...
最新文章
- 旋转矩阵、欧拉角、轴角
- flask使用debug模式时,存在错误时,会占用设备内存直至服务重启才释放;debug模式会开启一个守护进程(daemon process)...
- access 一亿条数据_循环运算数据溢出
- 云存储领域的那些安全难题与解决途径
- 【小白学PyTorch】16.TF2读取图片的方法
- autopep8规范你的python代码
- How to change the status of Prepayment invoice
- glDrawElements参数在新旧版本传最后一个参数的不同
- Java基础学习总结(72)——提升 java 代码的运行效率
- 浅议.NET、ASP.NET和C#的关系
- 【超清视频】SCCM2012精讲系列课程03:SCCM2012的基础环境初始化(一)
- POJ 1625 Censored!(AC自动机 + DP + 大数 + 拓展ASCII处理)题解
- kotlin写的音视频app(已开源)
- 科普:什么是人工智能
- 怎么解决计算机键盘驱动,键盘驱动不兼容,详细教您键盘驱动不兼容怎么解决...
- linux 下numpy安装包,Linux下安装numpy
- android开发教程!看懂这些帮你轻松解决就业问题!知乎上转疯了!
- pxe自动装机利用tfp,http,nfs服务实现。
- Java实现的一个简易网络画板
- 啤酒每罐2.3元,饮料每罐1.9元