物体が力
![{\displaystyle f(x)=-kx}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b99171922918ac66fd4e6a446d5eb7d5a18997a4)
を受けるときの運動を考える。このとき運動方程式は、
![{\displaystyle m{\ddot {x}}(t)+kx(t)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3ac404f8fdd967a088bdc779b2271896ef5ccc68)
で表される。変数k を、k = m ω2 の関係を使ってωに置き換えると、
![{\displaystyle {\ddot {x}}(t)+\omega ^{2}x(t)=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9cd3752d8569ef6534eecd60ecaa0bdf20af2158)
が得られる。この式は、定数係数線形2階微分方程式であるので簡単に解くことができ、解は、
![{\displaystyle x(t)=A\cos \omega t+B\sin \omega t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e7cb7390640ccf2af0ae88b5416e879a82a022a9)
となる。ここでA , B は任意定数であり、これらを決めるためには2つの初期条件が必要だが、ここでは特に、t = 0 で、
![{\displaystyle x(0)=x_{0},\quad v(0)\equiv {\dot {x}}(t)=v_{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b42ec01780f7dabd5e783292c20668d3fdddd220)
とすると
![{\displaystyle A=x_{0},\quad B={\frac {v_{0}}{\omega }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cfedbd5ea5acd780c9e483e762416a2869c562fb)
で表される。よって、この式の解は
![{\displaystyle {\begin{aligned}x(t)&=x_{0}\cos \omega t+{\frac {v_{0}}{\omega }}\sin \omega t\\&={\sqrt {x_{0}^{2}+\left({\frac {v_{0}}{\omega }}\right)^{2}}}\cos(\omega t-\theta ),\quad \tan \theta ={\frac {v_{0}}{x_{0}\omega }}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e83956234d85976f6c8f9ee0d42b904337db2bc9)
となる。この運動を単振動、ωを固有角振動数と呼ぶ。
次に力が
![{\displaystyle f(x,{\dot {x}})=-kx-m\gamma {\dot {x}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ab419597ab4a21ea9dd99a4c691426f425bce22)
で与えられる場合を考える。第2項は、速度に比例した力である。この様な力は空気抵抗などに見られる。この場合の運動方程式は
![{\displaystyle {\begin{aligned}m{\ddot {x}}(t)+m\gamma {\dot {x}}(t)+kx(t)&=0\\\therefore \quad {\ddot {x}}(t)+\gamma {\dot {x}}(t)+\omega ^{2}x(t)&=0\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/93a07ee7814820c6ff765ebe8b438cbf4e10a847)
となる。
これも定数係数の2階微分方程式なので解けるのだが、ここでは具体的に計算する。まず、
![{\displaystyle x(t)=e^{at}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bfe4299d064f002554b12ee643f1ae3acfab4895)
と仮定する。これを運動方程式に代入すると
![{\displaystyle a^{2}+\gamma a+\omega ^{2}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2717fbbb111740058b0476b25ebd830cab84e137)
が得られ、a は
![{\displaystyle a_{\pm }={\frac {1}{2}}(-\gamma \pm {\sqrt {\gamma ^{2}-4\omega ^{2}}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f9a1fb4d075fbc2fcbcd5bbf9ef0771e2b6421a6)
で与えられることが分かる。ここで、抵抗力の係数であるγが小さい数であるとするなら、根号の中身は負になり、この解は複素数によって与えられる。実部をとると一般解は
![{\displaystyle {\begin{aligned}x(t)&=Ae^{ia_{+}t}+Be^{ia_{-}t}\\&=e^{-\gamma t/2}\{\alpha \sin(t{\sqrt {\omega ^{2}-\gamma ^{2}/4}})+\beta \cos(t{\sqrt {\omega ^{2}-\gamma ^{2}/4}})\}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd6dea58764f417fde4ee6b74ae64a7ea7764acc)
となる(A , B , α, βは任意係数)。この解で、係数 exp(-γt /2) は粒子の運動が抵抗力を受けて時間的に減衰していく様子を示している。一方
![{\displaystyle \alpha \sin(t{\sqrt {\omega ^{2}-\gamma ^{2}/4}})+\beta \cos(t{\sqrt {\omega ^{2}-\gamma ^{2}/4}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/af14975f1acaf16d07e16fe92ad02fb9945e38ea)
の項は、この物体がωに近い角振動数で振動していることを示している。この振動を減衰振動と呼ぶ。
物体が単振動の力に加えて、周期的な外力を受けている場合を考える。このとき、
![{\displaystyle f(x,t)=-kx+ml\sin \omega _{0}t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7a0589a2ba6c777502af9063b291e2b53dee93a8)
となる(外力の大きさを表すパラメータとして、ここでは、後の計算を簡単にするために、パラメータをml という形に置いた)。
このとき運動方程式は、
![{\displaystyle {\begin{aligned}m{\ddot {x}}(t)+m\omega ^{2}x(t)&=ml\sin \omega _{0}t\\\therefore \quad {\ddot {x}}(t)+\omega ^{2}x(t)&=l\sin \omega _{0}t\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7c18696abb6dabcd95e8ca34a5d223db7eb37ec3)
となる。この式は、定数係数2階常微分方程式に加えて、右辺に関数項が加わった形をしている。このとき特別な解の形を予測して、この方程式の特解を求めることが出来る。ここでは、解を
![{\displaystyle x(t)=C\sin \omega _{0}t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f2daffec5f194895b625aa38c78d4e052b462bb)
と仮定する。これを運動方程式に代入すると
![{\displaystyle {\begin{aligned}&(-C\omega _{0}^{2}+\omega ^{2}C)\sin \omega _{0}t=l\sin \omega _{0}t\\&\therefore \quad C={\frac {l}{\omega ^{2}-\omega _{0}^{2}}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3deadc0707c01a0d5718ec7e70c5d02890597c95)
が得られる。よって、
![{\displaystyle x(t)={\frac {l}{\omega ^{2}-\omega _{0}^{2})}}\sin \omega _{0}t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d6919005729f46dea115425ef545ce184faa9f47)
は、この方程式の特解となる。左辺の線形方程式に対する一般解を加えると、この方程式全体に対する解は、
![{\displaystyle x(t)={\frac {l}{\omega ^{2}-\omega _{0}^{2}}}\sin \omega _{0}t+A\sin \omega t+B\cos \omega t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55daf8402d33faffae5d8f90cf08faa07dbd8a07)
と表される。
この式は、2つの部分に分かれている。まず後半の
![{\displaystyle A\sin \omega t+B\cos \omega t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6b3b3b1942045c15510dd9f9049ec70bd816b12d)
の部分は、物体が外力を受けていても、単振動と同じ角振動数の運動を繰りかえすことを表しており、次に前半の
![{\displaystyle {\frac {l}{-\omega _{0}^{2}+\omega ^{2}}}\sin \omega _{0}t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9bd0dd70b48266364bfac3181600206c11ab3f59)
の部分は、外力と同じ角振動数
の周期的な運動が引き起こされることを表している。
この式を見ると、ω = ω0 のとき、x (t ) の値が無限大になるように見える。実際、系の固有角振動数ωと外力の角振動数ω0 が近いとき、物体にきわめて大きな振動が引き起こされることが知られている。この現象を共鳴と呼ぶ。
共鳴の場合にx が無限大になるという結果がでたのは、もちろん物理的におきることではなく、上の解法(特に特解として仮定した形)が不十分だったことを意味する。正しい解を求めるには解として仮定する形をより広げなければならないが、それとは別に物理的によりわかりやすい方法もある。アイデアは「上の解法で得た解はωとω0 が少しでも違えばうまくいく。一方、パラメータを連続的に変えると現象も(多分)連続的に変わるはずだから、ω = ω0 の場合というのは、ω0→ω の極限として得られるはずである。だから上の解でその極限を取れば共鳴の場合の正しい解も得られる。」ということである。但し、上の解で単純にその極限をとってももちろんうまくいかない。係数A , B が意味を失うからである。そこでA , B の代りに常にはっきりした物理的意味を持つ量、すなわち初期条件(t = 0 でのx , dx /dt の値)を使って上の解を書き直し、その上でω0→ω の極限を取る。t = 0 でx (0) = X , dx (0)/dt = V とすると
![{\displaystyle A={\frac {1}{\omega }}\left(V+{\frac {l\omega _{0}}{\omega _{0}^{2}-\omega ^{2}}}\right),\quad B=X}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe79b2fafd79be615f964de5f4b77a8093c3ae0e)
になるので解は
![{\displaystyle {\begin{aligned}x(t)&={\frac {l}{-\omega _{0}^{2}+\omega ^{2}}}\sin \omega _{0}t+{\frac {1}{\omega }}\left(V+{\frac {l\omega _{0}}{\omega _{0}^{2}-\omega ^{2}}}\right)\sin \omega t+X\cos \omega t\\&={\frac {l}{\omega (\omega _{0}+\omega )}}{\frac {\omega \sin \omega _{0}t-\omega _{0}\sin \omega t}{-\omega _{0}+\omega }}+{\frac {V}{\omega }}\sin \omega t+X\cos \omega t\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ec6ae7acf56783ad24ae855aad060c06fbdf9aba)
となる。こうした上でω0→ω の極限を取ると、極限は確かに存在して次のようになる:
![{\displaystyle x(t)={\frac {l}{2\omega ^{2}}}(\sin \omega t-\omega t\cos \omega t)+{\frac {V}{\omega }}\sin \omega t+X\cos \omega t}](https://wikimedia.org/api/rest_v1/media/math/render/svg/106f3468aa6d4df62a56ea3b94f59306d4b19799)
特に重要なのは第2項で、cosにt が掛かっている点である。振動の振幅が時間に比例して増大していくのである。この結果から共鳴の本質が分かる。つまり共鳴では外力から与えられるエネルギーが蓄積されていくというのが本質で、そのため振幅が時間とともに果てしも無く大きくなっていくのである(一休さんの有名なエピソードで、お寺の鐘を指一本で大きく動かすというのがある。正確な周期で指で押すことを続けると共鳴により少しづつエネルギーが蓄積されて最後には大きな振動になる)。
なお、上の解法では現象のパラメータ変化に対する連続性を前提としている。この前提はたいてい成り立つが、どんな場合でもとまではいかない。今の場合にどうかを確認するには得られた解を元の微分方程式(でω0 = ωとしたもの)に代入すればよい。確かに解になっていることが分かる。また、ω0 がωに近づくにつれ振動のグラフがどう変わって行くかを見るのも興味深い。