理论教育 拉普拉斯变换与微通道间流体流动

拉普拉斯变换与微通道间流体流动

时间:2023-11-07 理论教育 版权反馈
【摘要】:定义8.2 设函数 f 在区间 [0,∞) 上有定义,如果含复参数 p 的无穷积分对 p 的某一取值范围内是收敛的,则称为函数 f 的拉普拉斯变换,f 称为原函数,F 称为象函数. 称由象函数求原函数的积分为拉普拉斯逆变换,即有 L1[F]=f.微分定理8.2 设 f (n=1,2,…

拉普拉斯变换与微通道间流体流动

定义8.2 设函数 f(t) 在区间 [0,∞) 上有定义,如果含复参数 p 的无穷积分

对 p 的某一取值范围内是收敛的,则称

为函数 f(t) 的拉普拉斯变换,f(t) 称为原函数,F(p) 称为象函数. 称由象函数求原函数的积分

为拉普拉斯逆变换,即有 L−1[F(p)]=f(t).

微分定理8.2 设 f(n)(t) (n=1,2,…) 是分段连续函数,则

例: 非齐次波动方程

解:对方程和边界条件作关于t 的拉普拉斯变换,得

方程 (8.16) 的通解为

由 (8.16) 知,A=B=0,所以

对 (8.18) 作逆拉普拉斯变换,得

Matlab 解算:在程序中用 U(x,p) 代替 ũ(x,p). 对 ksin(πx⁄a) 进行拉普拉斯变换:

>>syms t a x positive,syms k p

>> F=laplace(k*sin(pi*x/a),t,p) %回车,得(www.daowen.com)

F =k*sin(pi*x/a)/p

利用拉普拉斯变换的性质和上面的Matlab结果,对原方程和边界条件作 Laplace变换后,可得 (8.15) – (8.16). 用Matlab 求解问题 (8.15) – (8.16):

>> syms x p c p a k

>> U=dsolve('D2U–p^2*U/c^2= –k*sin(pi*x/a)/p','U(0)=0,U(a)=0','x');

>> disp('U(x,p)= '),pretty(U) %回车,得

U(x,p)=

对 U(x,p)=ũ(x,p) 作拉普拉斯逆变换:

>> syms t positive

>> u=ilaplace(U,p,t)

u =k*c^2*a^2*sin(pi*x/a) * (1/pi^2/c^2-1/pi^2/c^2*cosh(pi*(-a^2*c^2)^(1/2)/a^2*t))

>> pretty(u)

于是原问题的解为

免责声明:以上内容源自网络,版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。

我要反馈