sh函数

定义

双曲函数（HyperbolicFunction）包括下列六种函数：

sinh/双曲正弦：sinh(x)=[e^x-e^(-x)]/2

cosh/双曲余弦：cosh(x)=[e^x+e^(-x)]/2

tanh/双曲正切：tanh(x)=sinh(x)/cosh(x)=[e^x-e^(-x)]/[e^x

+e^(-x)]

coth/双曲余切：coth(x)=cosh(x)/sinh(x)=[e^x+e^(-x)]/[e^(x)

-e^(-x)]

ch/双曲正割：ch(x)=1/cosh(x)=2/[e^x+e^(-x)]

csch/双曲余割：csch(x)=1/sinh(x)=2/[e^x-e^(-x)]

cosh^2(t)-sinh^2(t)=1

和性质t>0对于所有的t。

参数t不是圆角而是双曲角，它表示在x轴和连接原点和双曲线上的点

(cosht,sinht)的直线之间的面积的两倍。

函数coshx是关于y轴对称的偶函数。

函数sinhx是奇函数，就是说-sinhx=sinh(-x)且sinh0=0。[3]

实变双曲函数

y=sh(x)，定义域：R，值域：R，奇函数，函数图像为过原点并且穿越Ⅰ、

Ⅲ象限的严格单调递增曲线，当x->+∞时是（1/2）e^x的等价无穷大，函数图

像关于原点对称。

y=ch(x)，定义域：R，值域：[1，+∞），偶函数，函数图像是悬链线，最

低点是（0，1），在Ⅰ象限部分是严格单调递增曲线，当x->+∞时是（1/2）e^x

的等价无穷大，函数图像关于y轴对称。

y=th(x)，定义域：R，值域：（-1，1），奇函数，函数图像为过原点并且

穿越Ⅰ、Ⅲ象限的严格单调递增曲线，其图像被限制在两渐近线y=1和y=-1之

间，lim[x->+∞，tanh(x)=1],lim[x->-∞，tanh(x)=-1]。

y=cth(x)，定义域：{x|x≠0}，值域：{x||x|>1}，奇函数，函数图像分为

两支，分别在Ⅰ、Ⅲ象限，函数在（-∞，0）和（0，+∞）分别单调递减，垂直

渐近线为y轴，两水平渐近线为y=1和y=-1，lim[x->+∞，coth(x)=1]，

lim[x->-∞，coth(x)=-1]。

y=sch(x)，定义域：R，值域：（0，1]，偶函数，最高点是（0，1），函数

在（0，+∞）严格单调递减，x轴是其渐近线，lim[x->∞，ch(x)]=0。

y=xh(x)，定义域：{x|x≠0}，值域：{x|x≠0}，奇函数，函数图像分为两

支，分别在Ⅰ、Ⅲ象限，函数在（-∞，0）和（0，+∞）分别单调递减，垂直渐

近线为y轴，两水平渐近线为x轴，lim[x->∞，csch(x)]=0。

双曲函数名称的变更：sh也叫sinh，ch也叫cosh，th也叫tanh，cth也叫

coth，sch也叫ch，xh也叫csch。

双曲正弦：sh(z)=[e^z-e^(-z)]/2

双曲余弦：ch(z)=[e^z+e^(-z)]/2

解析性：shz，chz是全平面的解析函数。

周期性：shz，chz是周期函数，周期为2πi，这是完全不同于实变函数中

的性质。

反双曲函数

反双曲函数是双曲函数的反函数.，它们的定义为：

arcsh(x)=ln[x+sqrt(x^2+1）]

arcch(x)=ln[x+sqrt(x^2-1）]

arcth(x)=ln[sqrt(1-x^2)/(1-x)]=ln[(1+x)/(1-x)]/2

arccth(x)=ln[sqrt(x^2-1)/(x-1)]=ln[(x+1)/(x-1)]/2

arcsch(x)=±ln[1+sqrt(1-x^2)/x]

arcxh(x)=ln[1-sqrt(1+x^2)/x]，如果x<0

ln[1+sqrt(1+x^2)/x]，如果x>0

其中，

sqrt为squareroot的缩写，即平方根

3三角函数编辑

双曲函数与三角函数有如下的关系：

*sinhx=-i*sin(i*x)

*coshx=cos(i*x)

*tanhx=-i*tan(i*x)

*cothx=i*cot(i*x)

*chx=c(i*x)

*cschx=i*csc(i*x)

i为虚数单位，即i*i=-1

ch^2(x)-sh^2(x)=1

cth^2(x)-xh^2(x)=1

th^2(x)+sch^2(x)=1

加法公式

sinh(x+y)=sinh(x)*cosh(y)+cosh(x)*sinh(y)

cosh(x+y)=cosh(x)*cosh(y)+sinh(x)*sinh(y)

tanh(x+y)=[tanh(x)+tanh(y)]/[1+tanh(x)*tanh(y)]

coth(x+y)=（1+coth(x)*coth(y))/(coth(x)+coth(y))

减法公式

sinh(x-y)=sinh(x)*cosh(y)-cosh(x)*sinh(y)

cosh(x-y)=cosh(x)*cosh(y)-sinh(x)*sinh(y)

tanh(x-y)=[tanh(x)-tanh(y)]/[1-tanh(x)*tanh(y)]

coth(x-y)=（1-coth(x)*coth(y))/(coth(x)-coth(y))

二倍角公式

sinh（2x)=2*sinh(x)*cosh(x)

cosh（2x)=cosh^2(x)+sinh^2(x)=2*cosh^2(x)-1=2*sinh^2(x)

+1

tanh（2x)=2tanh(x)/（1+tanh^2(x))

coth（2x)=（1+coth^2(x))/2coth(x)

半角公式

cosh^2(x/2）=(cosh(x)+1）/2

sinh^2(x/2）=(cosh(x)-1）/2

tanh(x/2）=(cosh(x)-1）/sinh(x)=sinh(x)/(cosh(x)+1）

coth(x/2）=sinh(x)/(coth(x)-1）=(coth(x)+1）/sinh(x)

三倍角公式

sin（3*x)=3*sin(x)+4*sin^3(x)

sinh（3*x)=3*sinh(x)+4*sinh^3(x)

德莫佛公式

（cosh(x）±sinh(x))^n=cosh(nx)±sinh(nx)

双曲函数的恒等式都在圆三角函数有相应的公式。Osborn'srule指出：将

圆三角函数恒等式中，圆函数转成相应的双曲函数，有两个sinh的积时（包括

coth^2(x),tanh^2(x),csch^2(x),sinh(x)*sinh(y））则转换正负号，则可得

到相应的双曲函数恒等式。如

5导数编辑

（sinh(x))'=cosh(x)

（cosh(x))'=sinh(x)

（tanh(x))'=ch^2(x)

（coth(x))'=-csch^2(x)

（ch(x))'=-ch(x)tanh(x)

（csch(x))'=-csch(x)coth(x)

（arcsinh(x))'=1/sqrt(x^2+1）

（arccosh(x))'=1/sqrt(x^2-1）(x>1）

（arctanh(x))'=1/（1-x^2）(|x|<1）

（arccoth(x))'=1/（1-x^2）(|x|>1）

6不定积分编辑

∫sinh(x)dx=cosh(x)+c

∫cosh(x)dx=sinh(x)+c

∫ch^2(x)dx=tanh(x)+c

∫csch^2(x)dx=-coth(x)+c

∫ch(x)tanh(x)dx=-ch(x)+c

∫csch(x)coth(x)dx=-csch(x)+c

∫tanh(x)dx=ln(cosh(x))+c

∫coth(x)dx=ln|sinh(x)|+c

∫ch(x)dx=arctan(sinh(x))+c=2arctan(e^x)+c1=2arctan(tanh(x/2））

+c2

∫csch(x)dx=ln|coth(x)-csch(x)+c=ln|tanh(x/2）|+c

∫[1/sqrt(x^2+1）]dx=arcsinh(x)+c=ln(x+sqrt(x^2+1））+c

∫[1/sqrt(x^2-1）]dx=sgn(x)arccosh|x|+c=ln|x+sqrt(x^2-1）|+c

（sgn是符号函数.sgn(x)=x/|x|,x≠0;sgn(x)=0,x=0)

7级数表示编辑

sinh(z)=z+z^3/3!+z^5/5!+z^7/7!+...+z^（2k-1）/（2k-1）！+...(z∈C)

cosh(z)=1+z^2/2!+z^4/4!+z^6/6!+...+z^（2k)/（2k)!+...(z∈C)

arcsinh(z)=z-（1/6）z^3+（3/40）z^5-（5/112）z^7+...+(-1）^k[（2k-1）！！

/（2k)!!][z^（2k+1）/（2k+1）]+...(|z|<1）

arctanh(z)=z+z^3/3+z^5/5+z^7/7+...+z^（2k-1）/（2k-1）+...(|z|<1）

近似：

z→+∞，sinh(z)≈exp(z)/2；cosh(z)≈exp(z)/2；tanh(z)→1；

z→-∞，sinh(z)≈-exp(-z)/2；cosh(z)≈exp(-z)/2；tanh(z)→-

1；

z→0，sinh(z)≈z；cosh(z)≈1+z^2/2；

8实际应用编辑

双曲函数并非单纯是数学家头脑中的抽象，在物理学众多领域可找到丰富的

实际应用实例。

阻力落体

在空气中由静止开始下落的小石块既受重力的作用又受到阻力的作用。设小

石块的质量为m，速度为v，重力加速度为g，所受空气阻力假定与v2正比，阻

尼系数为μ。设初始时刻小石块静止。求其小石块运动速度与时间的关系。

解：

小石块遵循的运动方程为

mdv/dt=mg―μv2⑻

这是Riccati方程，它可以精确求解。

依标准变换方式，设

v=（m/μ）（z′/z）⑼

代入⑻式，再作化简，有

z''―（gμ/m)z=0⑽

⑽式的通解是

z=C1exp（√gμ/mt）+C2exp（－√gμ/mt）⑾

其中，C1和C2是任意常数。

由于小石块在初始时刻是静止的，初始条件为

v（0）=0⑿

这等价于

z′（0）=0⒀

因此，容易定出

C2=C1⒁

将⒁式代入⑾式，再将⑾式代入⑼式，就可得

满足初始条件的解

v=√mg/μtanh（√μg/mt)⒂

我们可以作一下定性的分析。小石块初始时刻静止。因此，随着时间增加，

开始时小石块速度较小，小石块所受的阻力影响较小，此时，小石块与不受阻力

的自由落体运动情况相类似，小石块加速度几乎是常数。反映在图1中，起始段

t和v的关系是直线。当小石块速度很大时，重力相对于阻力来说可以忽略，阻

力快速增加到很大的数值，导致小石块的速度几乎不再增加。此时，小石块加速

度接近零，v几乎不随时间而变化。从图1中可以看到，一段时间后，v相不多

是一平行于t轴的直线。

导线电容

真空中两条圆柱形无穷长平行直导线，横截面的半径分别为R1和R2，中心

线相距为d(d>R1+R2）。试求它们间单位长度的电容。

解：

设这两条导线都带电，单位长度的电荷量分别是为λ和―λ。

我们可以用电像法精确求解。电像法的思路是：

由于在静电平衡情况时，导线是等势体，因而我们可设想用偶极线来取代这

两条圆柱形带电导线，适当地选择偶极线的位置，使它们所产生的两个等势面恰

好与原来两导线的表面重合。这样就满足了边界条件。这里采用的偶极线是两条

无穷长的均匀带电平行直线，它们单位长度的电荷量也分别为λ和―λ。这偶

极线便是原来两带电导线的电像。于是就可以计算电势，从而求出电容来。为此

先求偶极线的等势面。

以偶极线所在的平面为z-x平面，取笛卡儿坐标系，使偶极线对称地处在z

轴的两侧，它们到z轴的距离都是a。如图2所示。这偶极线所产生的电势便为

φ=φ1+φ2

=（λ/2πε0）In（r1′/r1）+（―λ/2πε0）In（r2′/r2）

=（λ/2πε0）In[（r2/r1）（r1′/r2′）]⒃

y

P

r2r1

R2―λ+λR1x

O

aa

a2a1

图2：带电导线与其镜像

式中r1′和r2′分别是偶极线λ和―λ到某个电势参考点的距离。为方便

起见，我们取z轴上的电势为零，这样，r1′=r2′=a，于是，⒃式便化为

φ=（λ/2πε0）In（r2/r1）⒄

由于对称性，平行于z轴的任何一条直线都是偶极线的等势线。所以，我们

只须考虑z-y平面内任意一点P（z，y）的电势即可。于是

φ=（λ/4πε0）In{[（x2+a2）+y2]/[（x2―a2）+y2]}⒅

故偶极线的等势面方程便为

[（x2+a2）+y2]/[（x2―a2）+y2]=k2⒆

式中

k2=e4πε0φ/λ⒇

令

c=[（k2+1）/（k2―1）]a（21）

则⒆式可化为

（x―c）2+y2=[4k2/（k2―1）2]a2（22）

这表明，偶极线的等势面都是轴线平行于z轴的圆柱面，它们的轴线都在z

轴上z=c处，其横截面的半径为

R=∣2k/（k2―1）∣a（23）

这个结果启示，我们可以找到偶极线的两个等势面，使它们分别与原来两导

线的表面重合。这只要下列等式成立就可以了：

a1=∣c1∣=[（k12+1）/（k12―1）]a（24）

R1=∣2k1/（k12―1）∣a（25）

a2=∣c2∣=[（k22+1）/（k22―1）]a（26）

R2=∣2k2/（k22―1）∣a（27）

d=a1+a2（28）

由（24）至（27）式得

a12―R12=a2=a22―R22（29）

原来两导线表面的方程是

R1：（x―a1）2+y2=R12（30）

R2：（x+a2）2+y2=R22（31）

利用（29）式，可以把（30）和（31）式分别化为

x2+y2+a2=2a1x（32）

x2+y2+a2=―2a2x（33）

利用（32）和（33）两式，由⒅式得出，半径为R1和R2的两导线的电势分

别为

φ1=（λ/4πε0）In[（a1+a）/（a1―a）]（34）

φ2=―（λ/4πε0）In[（a2+a）/（a2―a）]（35）

于是两导线的电势差便为

U=φ1+φ2=（λ/2πε0）In[（a1+a）（a2―a）/R1R2]（36）

用已知的量消去未知数，可以得出

U=（λ/2πε0）In[（d2―R12―R2）/2R1R2+√[（d2―R12―R2）/2R1R2]2―1]

（37）

最后得出原来两导线为l一段的电容为

C=Q/U=2πε0l/In[（d2―R12―R22）/2R1R2+√[（d2―R12―R22）/

2R1R2]2―1]（38）

单位长度的电容为

c=2πε0/In[（d2―R12―R22）/2R1R2+√[（d2―R12―R22）/2R1R2]

2―1]（39）

利用反两曲余弦关系式

archx=In[（x+√x2―1）]（40）

对本题的精确解表示作简洁表示

c=2πε0/arch[（d2―R12―R22）/2R1R2]（41）

最后一式可以在一般手册上查到。

粒子运动

一电荷量为q、静质量为m0的粒子从原点出发，在一均匀电场E中运动，E=Eez

沿z轴方向，粒子的初速度沿y轴方向，试证明此粒子的轨迹为

x=（W0/qE）[cosh（qEy/p0c）―1]（42）

式中p0是粒子出发时动量的值，W0是它出发时的能量。

解：

带有电荷量q的粒子在电磁场E和B中的相对论性的运动方程为

dp/dt=q（E+v×B）（43）

式中v是粒子的速度，p是粒子的动量

p=mv=mv0/√1－v2/c2（44）

本题运动方程的分量表示式为

dpx=qE

dpy=0

dpz=0（45）

解之，有

px=qEt+C1

py=C2

pz=C3（46）

代入t=0时初始条件

px（0）=0

py（0）=p0

pz（0）=0（47）

定出积分常数后，可知

px=qEt

py=p0

pz=0（48）

粒子的能量为

W=mc2

=√p2c2+m02c4

=√（px2+py2+pz2）c2+m02c4

=√q2E2c2t2+W02（49）

因dx/dt=qEt/m=qEc2t/√q2E2c2t2+W02（50）

积分得

x=∫[qEc2t/√q2E2c2t2+W02]dt

=[√q2E2c2t2+W02－W02]/qE（51）

又由（48）式得

dy/dt=p0/m=p0c2/√q2E2c2t2+W02（52）

积分得

y=∫[p0c2/√q2E2c2t2+W02]dt

=（p0c/qE）arsh（qEct/W0）（53）

或（qEct/W0）=sinh（qEy/p0c）（54）

在（51）式和（54）式中消去t，有

x=（W0/qE）[√1+sinh2（qEy/p0c）－1]（55）

利用恒等变换公式

cosh2x―sinh2x=1（56）

（55）式可以写成

x=（W0/qE）[cosh2（qEy/p0c）－1]（57）

（57）式是一种悬链线。

图3：匀强电场中粒子的悬链线运动轨迹

讨论：

因双曲余弦泰勒级数展开式是

cosh（x）=1+x2/2!+x4/4!+x6/6!+……（58）

当v/c→0时，保留前2项，得

x=（qE/2mv02）y2（59）

（59）式是抛物线轨迹。《普通物理学》教材用经典牛顿力学求解，普遍会

给有这个结果。这表示，非相对论确是相对论在v/c→0时的极限。或者说，（59）

式成立的条件是v/c<<1，这也是牛顿力学的适用范围。

非线性方程

如著名的KdV（Korteweg-deVries）方程的形式为

ux+uux+βuxxx=0（60）[3]

它是非线性的频散方程，其中β是频散系数。用双曲函数展开法求其某些

特殊精确解。

解：

考虑其行波解

u（x，t）=φ（ξ）（61）

其中，

ξ=kx－ωt+ξ0（62）

KdV方程成为

－ωφξ+kφφξ+k3βφξξξ=0（63）

记

f=1/（coshξ+r），g=sinhξ/（coshξ+r）（64）

尝试

φ=a0+a1f+a2g（65）

注意存在关系式

df/dξ=－fg

dg/dξ=1－g2－rg

g2=1－2rf+（r2－1）f2（66）

将（65）式代入（63）式，并在（66）式的帮助下使所得方程中各项只含有

f和g的幂次项，且g的幂次项不大于1。合并f和g的同次幂项并取其系数为

零，就得到方程（63）对应的非线性代数方程组

－6βk3b1（r2－1）2=0,

－6βk3a1（r2－1）=0,

－2kb1（r2－1）（－6βk2r+a1）=0,

－k（－6βk2ra1+a12－b12+b12r2）=0,

b1（4βk3+ka0－ka0r2+3ka1r－7βk3r2+cr2－c）=0,

ωa1+kb12r－βk3a1－ka0a1=0,

－b1（ka1+ωr－βk3r－ka0r）=0（67）

用计算机代数系统Maple对此超定方程组进行运算，可求得k≠0，ω≠0时

的一个非平凡精确解

φ=（ω－βk3）/k+6βk2/（coshξ+1）=0（68）

其中，k、ω、ξ0为任意常数。

（68）式是孤波解，图4绘出了其函数图像形状（作图时取了β=1/6k2，

ω=βk3）。

图4：KdV方程的孤波解

从以上的讨论中可知，无论是在经典或近代的物理学内容中，还是在正在发

展中的物理学内容中，双曲函数起着不可或缺的重要作用。

悬链线

形如y=acosh(x/a)(a为常数）的函数的图象又叫悬链线，可以由柔软的绳

子得到，有点象抛物线，但其实两者差距很大.据说莱布尼兹（Leibniz）于1690

年最先解出悬链线方程，惠更斯（Huygens）和伯努利兄弟（Jacob

Bernoulli,JohannBernoulli）随其后.惠更斯在1691年把悬链线命名为

catenary.悬链线与抛物线有这样的关系：悬链线是直线上滚动的抛物线的焦

点的运动轨迹.悬链线的顶点的渐开线是曳物线（tractrix).这条曳物线的渐进

线称为悬链线的准线，悬链线绕准线旋转形成的曲面叫做悬链面.

数学证明

设最低点A处受水平向左的拉力H，右悬挂点处表示为C点，在AC弧线区段

任意取一段设为B点，则B受一个斜向上的拉力T，设T和水平方向夹角为θ，

绳子的质量为m，受力分析有：Tsinθ=mg；Tcosθ=H，

tanθ=dy/dx=mg/H，mg=ρs，，其中s是右段AB绳子的长度，ρ是绳子

线重量密度，代入得微分方程dy/dx=ρs/H；利用弧长公式ds=√（1+dy^2/dx^2）

*dx；所以s=∫√（1+dy^2/dx^2）*dx;所以把s带入微分方程得dy/dx=ρ∫√

（1+dy^2/dx^2）*dx/H;.....⑴对于⑴设p=dy/dx微分处理得

p'=ρ/H*√（1+p^2）......⑵p'=dp/dx;对⑵分离常量求积分∫dp/√

（1+p^2）=∫ρ/H*dx得ln[p+√（1+p^2）]=ρx/H+C，即asinhp（反双曲正

弦）=ρx/H+C当x=0时，dy/dx=p=0；带入得C=0；整理得asinhp=ρx/H另

祥解：（ln[p+√（1+p^2）]=ρx/H）；p=sh（ρx/H)（1+p^2=e^（2ρx/H)-2pe^

（ρx/H)+p^2）；（p=[e^（ρx/H)-e^(-ρx/H)]/2=dy/dx）；y=ch（ρx/H)*

H/ρ（y=H/（2ρ）*[e^（ρx/H)+e^(-ρx/H)]）；令a=H/ρ：y=a*cosh

(x/a)（y=a[e^(x/a)+e^(-x/a)]/⑵=a*cosh(x/a））。