![]() |
中国面包师贴吧-楼主(阅:2584/回:0)数字计算机傅里叶变换电路2第六部分拉普拉斯变换AD采样电路 根据拉普拉斯变换 +∞ -st F(s)=£[f(t)]= ∫ e f(t)dt, 0 设φ(s)=F`(s) F(s) φ(s)=tgα= s -st φ(s)=tgα=e f(t)dt, F(s) α=arctg s 设w(s)=tgβ -st e f(t) w(t)=tgβ= t -st e f(t) β=arctg t 因为α+β=π/2 -st F(s) e f(t) arctg +arctg =π/2 s t -st F(s) e f(t) arctg =π/2-arctg s t -st F(s) e f(t) =tg[π/2-arctg ] s t -st e f(t) F(s) =tg[π/2 -arctg ]s (1) t 这样就得到一个计算拉普拉斯变换的公式 同理可证, -st e f(t) F(s) arctg =π/2-arctg t s -st e f(t) F(s) =tg[π/2-arctg ] t s -st F(s) e f(t) =tg[π/2-arctg ]t s F(s) -st f(t) =tg[π/2-arctg ]t/e (2) s 这样就得到一个计算拉普拉斯逆变换的公式 下面的电路可参见《数字电路》上册,上海师范大学物理系编,上海人民出版社1975年出版 1 F(p)= -ap p(1+e ) 因为, -st e f(t) F(s) =tg[π/2 -arctg ]s (1) t -st e f(t) 1 tg[π/2-arctg ]s= t -as s(1+e ) 这样就得到一个采集方波信号的函数F(s), 用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,当有方波信号输入时,等式左右两端电路的电压值相等。同时计算机将采集得到的方法信号经过处理后变成函数F(s),存储在寄存器中, 当计算机从端口01发送方波信号,被物体反弹回来以后,经过端口02采集这个方波信号,经过上面的计算如果等式左右两端相等,则证明方波信号就是发射出去的方波信号。发射方波信号的函数如下: F(s) -st f(t) =tg[π/2-arctg ]t/e (2) s 1 -st f(t) =tg[π/2-arctg ]t/e 2 -as s (1+e ) 寄存器A和寄存器B里面的数据相互比较,如果两者相等减法器输出0, 此时,或门输出1,与门输出1,状态寄存器A保存1,证明接收信号就是发射的信号。 下面的电路可参见《数字电路》上册,上海师范大学物理系编,上海人民出版社1975年出版 -ap 1-e F(p)= 2 -ap ap (1+e ) 因为 F(s) -st f(t) =tg[π/2-arctg ]t/e (2) s -as 1-e -st f(t) =tg[π/2-arctg ]t/e 2 -as as (1+e ) 这样就得到一个采集方波信号的函数f(t), 用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,这样就形成了一个产生三角波的电路. 定积分是把函数在某个区间上的图像[a,b]分成n份,用平行于y轴的直线把其分割成无数个矩形,再求当n→∞时所有这些矩形面积的和, 习惯上,我们用等差级数分点,即相邻两端点的间距△x是相等的。但必须指出,△x不相等,积分值仍然相同。我们假设这些“矩形面积和”。 S=f(x )△x +f(x )△x +...+f(x )△x 1 1 2 2 n-1 n-1 那么当n→∞时,△x的最大值趋于0,所以所有的△x趋于0,所以S仍然区域积分值。 利用这个规律,在我们了解牛顿-莱布尼兹公式之前,我们便可以对某些函数进行积分。 例如:证明对于函数 x f(x)=k (k∈Q,k≠-1)有 b a +∞ bk -ak ∫ f(x)dx= 0 k+1 证明:选择等比级数来分点,令公比 a k a-b q= , q=k b k 那么矩形面积和为 a a k 2 S =k (aq -a)+k q (aq -aq) n a 提取k (aq-a),则有 a (k+1) 2(k+1) (n-1)(k+1) S =ak (q-1)[1+q +q +...q ] n 利用等比级数公式,得到 b a q-1 b a bk -ak S = (bk -ak )= n (k+1) N q -1 其中, (k+1) q -1 N= q-1 1/v 设k=u/v(u,v∈Z),令q =s,则 (k+1) (k+1) (k+1) s v-1 s v-1 s v-1 s-1 N= = = v v v s -1 s -1 s -1 s-1 令n增加,则s,q都趋于1,因而N的极限为 u+v u = +1=k+1 v v 根据拉普拉斯变换 +∞ -st F(s)=£[f(t)]= ∫ e f(t)dt, 0 因为, x f(x)=k (k∈Q,k≠-1)有 b a +∞ bk -ak ∫ f(x)dx= 0 k+1 -st 设φ(s)=e f(t) 同理可证: -bs -as +∞ -st be f(t)-ae f(t) ∫ e f(t)dt= 0 -s [e f(t)]+1 所以, -bs -as +∞ -st be f(t)-ae f(t) F(s)=£[f(t)]=∫ e f(t)dt= (3) 0 -s [e f(t)]+1 所以, -bs -as be f(t)-ae f(t) F(s)= -s [e f(t)]+1 -bs -as -s be f(t)-ae f(t) [e f(t)]+1= F(s) -bs -as -s be f(t)-ae f(t) e f(t)- +1=0 F(s) 1 f(t)= -bs -as be -ae -s -e F(s) F(s) f(t)= (4) -bs -as -s be -ae -e F(s) 下面的电路可参见《数字电路》上册,上海师范大学物理系编,上海人民出版社1975年出版 1 F(p)= -ap p(1+e ) -bs -as +∞ -st be f(t)-ae f(t) F(s)=£[f(t)]=∫ e f(t)dt= (3) 0 -s [e f(t)]+1 -bs -as be f(t)-ae f(t) 1 = -s -as [e f(t)]+1 s(1+e ) 这样就得到一个采集方波信号的函数F(s), 用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,当有方波信号输入时,等式左右两端电路的电压值相等。同时计算机将采集得到的方法信号经过处理后变成函数F(s),存储在寄存器中, 当计算机从端口01发送方波信号,被物体反弹回来以后,经过端口02采集这个方波信号,经过上面的计算如果等式左右两端相等,则证明方波信号就是发射出去的方波信号。发射方波信号的函数如下: F(s) f(t)= (4) -bs -as -s be -ae -e F(s) 1 -as s(1+e ) f(t)= -s -bs -as e be -ae - -as s(1+e ) -s -bs -as e be -ae - -as s(1+e ) f(t)= -as s(1+e ) 下面的电路可参见《数字电路》上册,上海师范大学物理系编,上海人民出版社1975年出版 -ap 1-e F(p)= 2 -ap ap (1+e ) F(s) f(t)= (4) -bs -as -s be -ae -e F(s) 2 -ms 1 ms (1+e ) f(t)= -ms -ms -bs -as -s 1-e 1-e be -ae -e 2 -ms ms (1+e ) 这样就得到一个采集方波信号的函数f(t), 用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,这样就形成了一个产生三角波的电路. 第七部分梅林变换 根据梅林变换 +∞ s-1 {Mf}(s)=φ(s)= ∫ x f(x)dx 0 逆变换是: 1 +∞ -s {M φ}(x)=f(x)= ∫ x φ(s)ds 2πi 0 这是在复平面上的垂直线上的线积分。在melin反演定理中给出了这种反演有效的条件。 这个转换以芬兰数学家Hjalmar Mellin命名。 例1,求函数f(t)=1的梅林变换 解:由梅林变换的定义,有 +∞ s-1 φ(s)=M(1)= ∫ x dx 0 1 +∞ = x 2πi 0 1 = (Re(t)>0) s 例2.求函数f(x)=x的梅林变换 解:由梅林变换的定义有 +∞ s-1 φ(s)=M(s)= ∫ x xdx 0 +∞ s = ∫ x dx 0 1 s+1 +∞ = x s+1 0 x = s+1 n 例3.求函数f(x)=x 的梅林变换,其中n是正整数,Re(s)>0 解: n +∞ s-1 n φ(s)=M(s )= ∫ x x dx 0 +∞ s-1+n = ∫ x dx 0 1 s+n +∞ = x s+1 0 n x = s+n ax 例4.求函数f(x)=e 的梅林变换 解: ax +∞ s-1 ax φ(s)=M(e )= ∫ x e dx 0 1 s-1 ax +∞ s-1 +∞ s-2 ax = x e - ∫ x e dx a 0 a 0 s-1 s-2 s-1 s-2 s-3 = M[x ]= =M[x ]=... a a a (s-1)! (s-1)! = M[1]= n n a sa 下面的内容可参见《富利叶变换》,苏联I.N.Sneddon著,何衍睿,张燮译,科学出版社1958年出版。 1.积分变换式 在相当的一段时期中,人们已经认识到,奥里佛*赫维赛德所创造的用来求解物理学与电工学中的瞬时问题的运算微积。在形式上是等于拉普拉斯变换式的系统应用的。这正是大多数运算微积的近代数本所采用的表现方法。 注:例如,参看 H.Jeffreys,"Operational Methods in Mathematical Ohysics"(剑桥,伦敦,1931), J.RCarrson,"Electric Circuit Theory and the Operational Calculus"(MeGraw-Hill,纽约,1926), N.W.MeLachlan,"Complex Variable and the Operator Calculus"(剑桥,伦敦,1942); G.Doetsch,"Theorie und Anwendung der Laplace-Transformation"(Dover,纽约,1944) H.S.Carslaw and J.C.Jaeger,"Operational Methods in Applied Matematics"(牛津,纽约,1941); R.V.Churchill,"Modern Operational Mathematics in Engineering"(MeGraw-Hill,纽约,1944) 例如说,倘若函数f(x)是由一个微分方程与某些边界条件所决定的。那么在某些场合下,可将f(x)的边值问题化为函数。 ∞ -px φ(p)=∫ f(x)e dx (1) 0 -px 的边值问题比较简便,而这个函数是由f(x)乘以e 再对x由0到∞求积分而得的。 用此种方式所定义的函数φ(p)显然是变数p的函数,它叫做函数f(x)的拉普拉斯变换式。 ∞ I (α)=∫ f(x)K(α,x)dx (2) 0 是收敛的,那么由方程(2)便决定变数α的一个函数,这个函数叫做函数f(x)的以K(α,x)为核的积分变换式。此种核的最简单的例子为: -αx K(α,x)=e 由它可以引出拉普拉斯变换式(1)。另一个在这种方式下常常用到的核是 α-1 K(α,x)=x 由它可以引出变换式 ∞ α-1 F(α)=∫ f(x)x dx (3) 0 梅林最先有系统的研究了这种类型的积分变换式,因而由方程(3)所决定的函数F(α)便叫做函数f(x)的梅林变换式。 注:H.Mellin,Aeta Soc.Fennicae,21,1-115(1896);Acta Math.,25,139(1902)。 当核K(α,x)是正弦或余弦,或者是贝塞而微分方程的解案时,便有特殊的变换式产生,以后我们要定义这些变换式,并讨论它们的一些性质。由定义(2)立刻可以看出,倘若f(x),g(x)是两个函数,它们具有以K(α,x)为核的积分变换式,那么f(x),g(x)之和的积分变换式即为 ∞ ∞ ∞ ∫ (f+g)K(α,x)dx=∫ f(x)K(α,x)dx+∫ g(x)K(α,x)dx 0 0 0 也就是说,它等于函数f(x),g(x)的积分变换式之和,其次,倘若c是一个常数,那么 ∞ ∞ ∫ cf(x)K(α,x)dx=c∫ f(x)K(α,x)dx 0 0 这些方程说明了一件事情:将函数f(x)化为它的积分变换式I (f)的算子是一个线性算子。此时我们也可以将方程(2)看作是函数f(x)与I (f)之间的一种变换。并利用巴拿哈空 α 间的性质来建立此种类型的变换的抽象理论,这种方案在数学上是极有趣的。 但若我们的主要目的是在数学物理中的应用时,那么这样的做法并不会有多少成果。因此,这里不易作出这种讨论。我们仅仅考虑积分变换式的某些性质,它们在以后的边值问题的分析中是有用的。 定理10.设F(a)为函数f(x)的富利叶变换式,设 1 ∞ iax F(a)= ∫ f(x)e dx (38) 2π -∞ 那么f(x)便可以由F(a)用下面的关系式表出: 1 ∞ -iax F(a)= ∫ F(a)e dx (38a) 2π -∞ 我们注意一点:函数f(x)与它的富利叶变换式F(a)的关系并不是对称的,而函数与它的富利叶余弦及正弦变换式F (a),F (a)的关系却是对称的。这些结果往往也可以用稍微不 c s 同的形式陈述出来。 例如,我们往往把方程(38)与(38a)写成如此的式样比较便利: 倘若 ∞ iax f(a)=∫ f(x)e dx (39) -∞ 那么 1 ∞ -iax f(x)= ∫ f(x)e dx (39a) 2π -∞ 1 将变换式的定义里面的因子 略去 2π 这种写法显然是常用的,但它却破坏了函数f(x)与它的变换式之间的一部分对称性。必须指出,我们只有在如次的情形下,才有权利来应用这些公式:函数f(x)在适当的区间(-∞,∞)或者(0,∞)内满足狄利克莱条件,而且积分 ∞ ∫ │f(x)│dx -∞ 是收敛的。 6.梅林变换式 6.1梅林变换定理 以前我们已经看到(1),函数f(x)的梅林变换式F(s)是由如次的表达式来定义的: ∞ s-1 F(s)=∫ f(x)x dx (99) 0 倘若在方程(38)中,令 x ξ=e ,s=e+ia, 那么它便化为下面的式样: s-c 1 ∞ -c s-1 F( )= ∫ ξ f(lnξ)ξ dξ i 2π 0 用同样的代换,也可以将(38a)化为 1 ∞ s-c c-s f(lnξ)= ∫ F( )ξ ds i 2π 0 i 因此,如果我们令 -1/2 -c g(ξ)=(2π) ξ f(lnξ ) s-c G(s)= F( ) i 的话,便得到梅林变换的反演公式如下: 定理15.倘若对于某个k>0而言,积分 ∞ k-1 ∫ ξ │g(ξ)│dξ 0 是有界的,并设 ∞ s-1 G(s)=∫ ξ g(ξ)dξ 0 那么, 1 c+i∞ -s g(ξ)= ∫ G(s)ξ ds 2πi c-i∞ 其中c>k, 函数的各阶导数的梅林变换式之间的关系,没有富利叶变换式与拉普拉斯变换式的情形那样简单,例如, r r-1 r-1 ∞ d f d f s-1 ∞ ∞ d f s-2 ∫ x dx=[ x ] -(s-1) ∫ x dx 0 r r-1 0 0 r-1 dx dx dx 因此,倘若我们假设f是如此的函数,使得上式的方括号化为零,那么我们便得到下面的关系式: (r) (r-1) F (s)=-(s-1)F (s-1) r (r) d f 其中F (s)代表导数 f的梅林变换式。 r dx 重复应用这种规则,便可以得到 (r) r Γ(s) F (s)=(-1) F(s-r) (100) Γ(s-r) 注:Γ(s)表示s!,Γ(s-r)表示(s-r)! 这个公式将已知函数的导数的梅林变换式用函数本身的梅林变换式表出。 6.2关于梅林变换的褶合式定理 设F(s)和G(s)分别为函数f(x)和g(x)的梅林变换式,则由定义可知,乘积f(x)g(x)的梅林变换式为 ∞ s-1 ∞ s-1 1 c+i∞ -σ ∫ f(x)g(x)x dx=∫ g(x)x dx ∫ F(σ)x dσ= 0 0 2πi c-i∞ 1 c+i∞ -σ ∞ s-σ-1 = ∫ F(σ)x dσ∫ g(x)x dx 2πi c-i∞ 0 1 c+i∞ = ∫ F(σ)G(s-σ)dσ (101) 2πi c-i∞ 此式有一特殊情形: 1 ∞ c+i∞ -s ∫f(x)g(x)dx= ∫ g(x)dx ∫ F(s)x ds 2πi 0 c-i∞ 1 c+i∞ = ∫ F(s)G(1-s)ds (101a) 2πi c-i∞ 放此,我们也可以得到乘积F(s)G(s)的梅林变换式如下: 1 c+i∞ -s ∫ F(s)G(s)x ds= 2πi c-i∞ 1 c+i∞ -s ∞ s-1 = ∫ F(s)x ds∫ g(u)u du 2πi c-i∞ 0 ∞ du 1 c+i∞ x -s =∫ g(u) ∫ F(s)( ) ds 0 u 2πi c-i∞ u ∞ x du =∫ f( )g(u) (102) 0 u u 1 c+i∞ 1 c+i∞ ∞ s-1 ∫ F(s)G(s)ds= ∫ F(s)ds∫ g(u)u du= 2πi c-i∞ 2πi c-i∞ 0 ∞ du 1 c+i∞ 1 -s =∫ g(u) ∫ F(s)( ) ds 0 u 2πi c-i∞ u ∞ 1 du =∫ f( )g(u) (102a) 0 u u 第八部分梅林变换 根据梅林变换 +∞ s-1 {Mf}(s)=φ(s)= ∫ x f(x)dx 0 逆变换是: 1 c+i∞ -s {M φ}(x)=f(x)= ∫ x φ(s)ds 2πi c-i∞ 这是在复平面上的垂直线上的线积分。在melin反演定理中给出了这种反演有效的条件。这个转换以芬兰数学家Hjalmar Mellin命名。 设φ(s)=F`(s) φ(s) g(s)=tgα= s -st g(s)=tgα= x f(x)dx φ(s) α=arctg s 设w(s)=tgβ s-1 x f(x) w(t)=tgβ= x s-1 x f(x) β=arctg x 因为α+β=π/2 s-1 φ(s) x f(x) arctg +arctg =π/2 s x s-1 φ(s) x f(x) arctg =π/2-arctg s x s-1 φ(s) x f(x) =tg[π/2-arctg ] s x s-1 x f(x) φ(s)=tg[π/2 -arctg ]s (1) x 这样就得到一个计算梅林变换的公式,同理可证, s-1 x f(x) φ(s) arctg =π/2-arctg x s s-1 x f(x) φ(s) =tg[π/2-arctg ] x s s-1 φ(s) x f(x)=tg[π/2-arctg ]x s φ(s) 2-s f(t) =tg[π/2-arctg ]x (2) s 这样就得到一个计算梅林逆变换的公式。 根据梅林变换,得 1 φ(s)= 1 s+ n a 它的原函数f(x的图像如右图所示 s-1 x f(x) φ(s)=tg[π/2 -arctg ]s (1) x s-1 x f(x) 1 tg[π/2 -arctg ]s= x 1 s+ n a 这样就得到一个采集方波信号的函数φ(s),用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,当有方波信号输入时,等式左右两端电路的电压值相等。同时计算机将采集得到的方法信号经过处理后变成函数F(s),存储在寄存器中。当计算机从端口01发送方波信号,被物体反弹回来以后,经过端口02采集这个方波信号,经过上面的计算如果等式左右两端相等,则证明方波信号就是发射出去的方波信号。发射方波信号的函数如下: φ(s) 2-s f(t) =tg[π/2-arctg ]x (2) s 1 2-s f(t) =tg[π/2-arctg ]x 1 1+ n sa 根据梅林变换,得 1 s- n a φ(s)= 2 1 as (s+ ) n a 它的原函数f(x的图像如右图所示 因为, φ(s) 2-s f(t) =tg[π/2-arctg ]x (2) s 1 s- n a 2-s f(x) =tg[π/2-arctg ]x 2 1 as (s+ ) n a 这样就得到一个采集方波信号的函数f(t),用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,这样就形成了一个产生三角波的电路。 定积分是把函数在某个区间上的图像[a,b]分成n份,用平行于y轴的直线把其分割成无数个矩形,再求当n→∞时所有这些矩形面积的和, 习惯上,我们用等差级数分点,即相邻两端点的间距△x是相等的。但必须指出,△x不相等,积分值仍然相同。我们假设这些“矩形面积和”。 S=f(x )△x +f(x )△x +...+f(x )△x 1 1 2 2 n-1 n-1 那么当n→∞时,△x的最大值趋于0,所以所有的△x趋于0,所以S仍然区域积分值。 利用这个规律,在我们了解牛顿-莱布尼兹公式之前,我们便可以对某些函数进行积分。 例如:证明对于函数 k f(x)=x (k∈Q,k≠-1)有 k+1 k+1 +∞ b -a ∫ f(x)dx= 0 k+1 证明:选择等比级数来分点,令公比 n a b q= , =q b a n 2 b=a*q ,x =a, x =a 0 2 (i+1) i 且△x =x -x =aq -aq i i+1 i 那么矩形面积和为 k k k 2 k (n-1)k n (n-1) S =a *(a*q-a)+a *q *(a*q -a*q)+...+a q [aq -aq ] n n k 提取a (aq-a),则有 (k+1) (k+1) 2(k+1) (n-1)(k+1) S =a (q-1)[1+q +q +...+q ] n 利用等比级数公式,得到 (k+1) (k+1) q-1 (k+1) (k+1) b -a S = (b -a )= n (k+1) N q -1 其中, (k+1) q -1 N= q-1 1/v 设k=u/v(u,v∈Z),令q =s,则 u+v (k+1) u+v s v-1 s v-1 s v-1 s-1 N= = = v v v s -1 s -1 s -1 s-1 令n增加,则s,q都趋于1,因而N的极限为 u+v u = +1=k+1 v v 根据梅林变换 +∞ s-1 {Mf}(s)=φ(s)= ∫ x f(t)dt, 0 因为, k f(x)=x (k∈Q,k≠-1)有 k+1 k+1 +∞ b -a ∫ f(x)dx= 0 k+1 s-1 设g(t)=x f(x)dx 同理可证: s s +∞ s-1 b f(x)-a f(x) ∫ x f(x)dx= 0 s-1 [x f(x)]+1 所以, s s +∞ s-1 b f(x)-a f(x) {Mf}(s)=φ(s)=∫ x f(x)dx= (3) 0 s-1 [x f(x)]+1 所以, s s b f(x)-a f(x) {Mf}(s)=φ(s)= s-1 [x f(x)]+1 s s s-1 b f(x)-a f(x) [s f(x)]+1= φ(s) s s s-1 b f(x)-a f(x) x f(x)- +1=0 φ(s) s s s-1 b -a f(x)[x - ]+1=0 φ(s) 1 f(x)= s s b -a s-1 -x φ(s) φ(s) f(x)= (4) s s s-1 b -a -x φ(s) 根据梅林变换,得 1 F(p)= 1 s+ n a 它的原函数f(x的图像如右图所示 因为, s s +∞ s-1 b f(x)-a f(x) {Mf}(s)=φ(s)=∫ x f(x)dx= (3) 0 s-1 [x f(x)]+1 s s b f(x)-a f(x) 1 = s-1 -as [x f(x)]+1 s(1+e ) 这样就得到一个采集方波信号的函数φ(s),用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,当有方波信号输入时,等式左右两端电路的电压值相等。 同时计算机将采集得到的方法信号经过处理后变成函数F(s),存储在寄存器中 当计算机从端口01发送方波信号,被物体反弹回来以后,经过端口02采集这个方波信号,经过上面的计算如果等式左右两端相等,则证明方波信号就是发射出去的方波信号。发射方波信号的函数如下: φ(s) f(x)= (4) s s s-1 b -a -x φ(s) 1 -as s(1+e ) f(t)= s s s-1 1 b -a -x -as s(1+e ) 根据梅林变换,得 1 s- n a φ(s)= 2 1 as (s+ ) n a 它的原函数f(x的图像如右图所示 因为, φ(s) f(x)= (4) s s s-1 b -a -x φ(s) 1 s- n m 2 1 ms (s+ ) n m f(x)= 1 s- n s s s-1 m b -a -x 2 1 ms (s+ ) n m 这样就得到一个采集方波信号的函数f(t), 用加法器,乘法器乘法器,除法器,tg计算电路表示上面的公式,这样就形成了一个产生三角波的电路。 |
| 发帖须知: 1,发帖请遵守《计算机信息网络国际联网安全保护管理办法》、《互联网信息服务管理办法》、 《互联网电子公告服务管理规定》、《维护互联网安全的决定》等法律法规。 2,请对您的言论负责,我们将保留您的上网记录和发帖信息。 3,在此发帖表示认同我们的条款,我们有权利对您的言论进行审核、删除或者采取其他在法律、地方法规等条款规定之内的管理操作。 |